如何用matlab解方程组

 我来答
棺欢目烫cl
2017-06-13 · TA获得超过340个赞
知道小有建树答主
回答量:321
采纳率:0%
帮助的人:163万
展开全部
matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法:
(1)x=inv(A)*b — 采用求逆运算解方程组;

(2)x=A\B — 采用左除运算解方程组
PS:使用左除的运算效率要比求逆矩阵的效率高很多~
例:
x1+2x2=8
2x1+3x2=13
>>A=[1,2;2,3];b=[8;13];
>>x=inv(A)*b
x =
2.00
3.00
>>x=A\B
x =
2.00
3.00;
即二元一次方程组的解x1和x2分别是2和3。

对于同学问到的用matlab解多次的方程组,有符号解法,方法是:先解出符号解,然后用vpa(F,n)求出n位有效数字的数值解.具体步骤如下:
第一步:定义变量syms x y z ;
第二步:求解[x,y,z,]=solve('eqn1','eqn2',,'eqnN','var1','var2','varN');
第三步:求出n位有效数字的数值解x=vpa(x,n);y=vpa(y,n);z=vpa(z,n);。
如:解二(多)元二(高)次方程组:
x^2+3*y+1=0
y^2+4*x+1=0
解法如下:
>>syms x y;
>>[x,y]=solve('x^2+3*y+1=0','y^2+4*x+1=0');
>>x=vpa(x,4);
>>y=vpa(y,4);
结果是:
x =
1.635+3.029*i
1.635-3.029*i
-.283
-2.987
y =
1.834-3.301*i
1.834+3.301*i
-.3600
-3.307。
二元二次方程组,共4个实数根;

解答如下:
基本方法是:solve(s1,s2,…,sn,v1,v2,…,vn),即求表达式s1,s2,…,sn组成的方程组,求解变量分别v1,v2,…,vn。
具体例子如下:
x^2 + x*y + y = 3
x^2 - 4*x + 3 = 0
解法:
>> [x,y] = solve('x^2 + x*y + y = 3','x^2 - 4*x + 3 = 0')
运行结果为
x =
1 3
y =
1 -3/2

即x等于1和3;y等于1和-1.5


>>[x,y] = solve('x^2 + x*y + y = 3','x^2 - 4*x + 3= 0','x','y')
x =
1 3
y =
1 -3/2
结果一样,二元二方程都是4个实根。

通过这三个例子可以看出,用matlab解各类方程组都是可以的,方法也有多种,只是用到解方程组的函数,注意正确书写参数就可以了,非常方便。

2、变参数非线性方程组的求解
对于求解非线性方程组一般用fsolve命令就可以了,但是对于方程组中某一系数是变化的,该怎么求呢?

%定义方程组如下,其中k为变量
function F = myfun(x,k)
H=0.32;
Pc0=0.23;W=0.18;
F=[Pc0+H*(1+1.5*(x(1)/W-1)-0.5*(x(1)/W-1)^3)-x(2);
x(1)-k*sqrt(x(2))];

%求解过程
H=0.32;
Pc0=0.23;W=0.18;
x0 = [2*W; Pc0+2*H]; % 取初值
options = optimset('Display','off');
k=0:0.01:1; % 变量取值范围[0 1]
for i=1:1:length(k)
kk=k(i);
x = fsolve(@(x) myfun(x,kk), x0, options);%求解非线性方程组
x1(i)=x(1);
x2(i)=x(2);
end
plot(k,x1,'-b',k,x2,'-r');
xlabel('k')
legend('x1','x2')
lhmhz
高粉答主

2019-10-25 · 专注matlab等在各领域中的应用。
lhmhz
采纳数:7263 获赞数:17004

向TA提问 私信TA
展开全部

用matlab求解方程组,可以用solve函数求得其解析值。

这里用a—表示α,b—表示β,c—表示θ,d—表示η

实现求解此方程组的代码:

syms D Pr a Pd b c s w k d

[Pr,s]=solve((1-c)*D-Pr+a*Pd-(1-b)*s-Pr+w==0,(1-b)*(Pr-w)-k*d*s==0,'Pr,s')

运行结果

已赞过 已踩过<
你对这个回答的评价是?
评论 收起
兄弟连教育
2018-07-28 · 百度知道合伙人官方认证企业
兄弟连教育
兄弟连教育成立于2006年,11年来专注IT职业教育,是国内专业的IT技术培训学校。2016年成功挂牌新三板(股票代码:839467)市值过亿。开设专注程序员培训专注php、Java、UI、云计算、Python、HTML5、
向TA提问
展开全部
S=solve('2*x*y=1,x+2=y+z,x+y-z=4','x,y,z'); %前面的参数是方程组列表,后面是未知变量列表
S.x %输出未知数x的值
S.y %输出未知数y的值
S.z %输出未知数z的值
f=@(x)2*x; %定义一个匿名函数y=2x,其中@(x)表示x是匿名函数的自变量
fplot(f,[-10,10]) %画图函数,第一个参数是函数名,第二参数要画函数的区间[-10,10]
本回答被网友采纳
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
乐就乐520
2016-08-11 · TA获得超过8971个赞
知道大有可为答主
回答量:7977
采纳率:91%
帮助的人:1543万
展开全部
S=solve('2*x*y=1,x+2=y+z,x+y-z=4','x,y,z'); %前面的参数是方程组列表,后面是未知变量列表
S.x %输出未知数x的值
S.y %输出未知数y的值
S.z %输出未知数z的值
f=@(x)2*x; %定义一个匿名函数y=2x,其中@(x)表示x是匿名函数的自变量
fplot(f,[-10,10]) %画图函数,第一个参数是函数名,第二参数要画函数的区间[-10,10]
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
ranymondyaobai
2017-11-08 · TA获得超过1.5万个赞
知道大有可为答主
回答量:3022
采纳率:86%
帮助的人:366万
展开全部
用solve函数
先举一例,解方程"x^2+100*x+99=0"
在matlab ”Command Window"中输入如下命令:
x=solve('x^2+100*x+99=0','x')
首先来求一个二元一次方程组
9x+8y=10 式1
13x+14y=12 式2
[x,y]=solve('9*x+8*y=10','13*x+14*y=12','x','y')
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
收起 更多回答(7)
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式