matlab解方程组实例 matlab解n元一次方程组

在求解线性方程组时,会遇到以下几种情形:定解方程组、不定方程组、超定方程组、奇异方程组.作为示例,首先以定解线性方程组为例:在分析如上方程组时,需要知道,方程中有3个未知数,而方程也有3个,所以可以求出(x,y,z)值,转化为矩阵即为: AX = B,其中A为系数矩阵,B为右边值向量.而X即为未知数构成的向量,转化后即为:>> A = [2,3,1; 4,2,3; 7,1,-1]; 如上为系数矩阵;>> B = [4;17;1]; 如上为右边值矩阵;利用矩阵除法:>> X = A\B 求得结果如下图红色箭头所示:x = 1.0000 ; y = -1.0000 ; z = 5.0000;

A=[16 2 3 13;5 11 10 8;9 7 6 12;4 14 15 1];b=[1;3;4;7];det(A) %因为对应的行列式的值为0,有无穷多个解.X=null(A,'r') %得到的结要为1列,说明有一个自由变量x0=pinv(A)*b%通解为X*x4+x0

用solve命令 eg syms x y>> [x y]=solve('(x+3)*(y-2)=(x-3)*(y+10)','(x-1)*(y+3)=(x+2)*(y+12)')

syms x y v t; [x y v t ]=solve(v*t-sqrt((x-1700).^2+(y-700).^2),v*t-sqrt((x-2900).^2+(y-900).^2),v*(t+1.14)-sqrt((x-1400).^2+(y-2200).^2),v*(t-1.37)-sqrt((x-2500).^2+(y-1900).^2)); A=[x y v t]; eval(A) 结果 ans = 1.0e+003 * 0.0040 0.2181 2.1976 1.4146

你这个方程组是不是有点问题,四元一次方程组,怎么还多了q,m,l,a这些变量?如果单纯把上述变量当做待设参数 syms fa fb fc ma q m l a; [fa,fb,fc,ma]=solve('fa-q*2*l-fc=0','-ma-m+fb*l-q*2*l*2*l-fc*3*l=0',. '-ma-fa*l-m-q*2*l*l-fc*a*l=0','-ma-fa*3*l-m-fb*2*l+q*2*l*l=0','fa,fb,fc,ma') 这样就可以了.

K=[5.4000 2.6000 -3.3000 -0.5000 0.7000 0.7000 0 0 -2.6000 5.4000 0.5000 0.6000 -0.7000 2.0000 0 0 -3.3000 -0.5000 4.7000 1.9000 0.3500 0.3500 -0.3500 -0.3500 0..

syms a b Rsd RspRal=72.8;Ralp=51;Rald=21.8;Rbl=50.8;Rblp=2.3;Rbld=48.5;f1=Ral*(1+cos(a))-2*(sqrt(Ralp*Rsp)+sqrt(Rald*Rsd));f2=Rbl*(1+cos(b))-2*(sqrt(Rblp*Rsp)+sqrt(.

用matlab 中的反向斜线运算符(backward slash) 分析: 方程组可化为 2*x-y=-3; 3*x-y=7; ax=b (*) a=[2,-1;3,-1]; b=[-3,7]; x=a\b %可以看成将(*)式左边都除以系数矩阵a >>a=[2,-1;3,-1]; >>b=[-3,7]; >>x=a\b x = 10.0000 % x = 10.0000 23.0000 % y = 23.0000 使用这个方法时,要注意方程组未知数的顺序是有序的;

>> A = [1 1.5 2 9 7;0 3.6 0.5 -4 -4;7 10 -3 22 33;3 7 8.5 21 6; 3 8 0 90 -20];>> b = [3 -4 20 5 16]';>> x=inv(A)*bx = 3.5056 -0.8979 -0.2745 0.1438 0.0137

这种形式是不被matlab接受的.你必须将矩阵方程展开,然后列写为方程组的形式,然后用solve函数求解.