matlab中矩阵方程 matlab中矩阵方程组求解。
matlab\u5982\u4f55\u89e3\u77e9\u9635\u65b9\u7a0b\uff1f1. \u8bbeAx = b\uff0c\u6c42x\uff0c\uff08x\u548cb\u90fd\u662f\u5411\u91cf\uff09\u5219\uff1ax = A\b
\u5c31\u8fd9\u4e48\u7b80\u5355\u3002\u60f3\u4e0d\u51fa\u66f4\u590d\u6742\u7684\u65b9\u6cd5\u4e86\u3002
2. \u4f60\u4e0a\u9762\u7684\u770b\u8d77\u6765\u662f\u4e2a\u9f50\u6b21\u65b9\u7a0b\uff0c\u4f46\u662f\u56e0\u4e3ab1 b2 b3\u5df2\u77e5\uff0c\u53ef\u4ee5\u79fb\u5230\u7b49\u53f7\u53e6\u4e00\u8fb9\uff1b\u518d\u628aV3=V4\u6574\u7406\u5230\u77e9\u9635\u65b9\u7a0b\u91cc\u6216\u8005\u5e72\u8106\u4ece\u65b9\u7a0b\u4e2d\u53bb\u6389V4\uff0c\u5c31\u53ef\u4ee5\u7528\u4e0a\u9762\u65b9\u6cd5\u89e3\u4e86\u3002
\u8bf7\u53c2\u8003 \u6700\u5c0f\u4e8c\u4e58\u6cd5\uff0c\u5c06\u4f60\u7684\u7b49\u5f0f\u91cd\u65b0\u6392\u5217\uff0c\u7136\u540e\u6700\u5c0f\u4e8c\u4e58\u6c42\u89e3\u3002\u5982\u679c\u4e0d\u662f\u7ebf\u6027\u7b49\u5f0f\uff0c\u90a3\u4e48\u5148\u9700\u8981\u7ebf\u6027\u5316
若f=0（齐次方程），则用克莱姆法则求便可求其解（求行列式用命令det())；
若f不为0（非齐次），则用高斯消元法（即初等变换法）可求其解。
Matlab中提供命令lu对矩阵进行LU分解，如果是稀疏矩阵，则可使用命令lunic对矩阵进行LU分解。
你的4元一次方程若表示成：AX=f
则A=[a,b,c,d],X=[x,y,z,u]',在Matlab中输入
>>A=[.......]; %%把a,b,c,d都输入
>>f=[.......]; %%输入你的f
>>X=A\f; %%列出线性方程
>>X %%显式方程的解
绛旓細棰樹富鐨闂鍙互鐢╢or鍙屽惊鐜鍙ュ拰solve锛堟垨vpasolve锛夊嚱鏁版眰鍑篎I(i)銆丟I(i)鐨勫笺傚疄鐜版柟娉曪細k=1;for i=1:4 for j=1:4 syms FI0銆丟I0 [FI0,GI0]=vpasolve(FI0*dGI(i,j)-GI0*dSI(i,j)==0,FI0*dSI(i,j)-GI0*dGI(i,j)==0)FI(k)=FI0,GI(k)=GI0 k=k+1;end end ...
绛旓細鐢╯=solve('A*x=B')鏄緱涓嶅埌s鍊鐨銆傚簲璇ヤ娇鐢ㄤ笅鍒楀懡浠 s=A\B %鐭╅樀闄ゆ硶锛屽緱鍒版渶灏忎簩涔樿繎浼艰В銆
绛旓細绋嬪簭锛欰=[1 2 3;1 3 1; 0 1 2];B=[1 0; 0 1;-1 0];X=(A+2)\B%鍖栫畝鐭╅樀鏂圭▼銆侫X=B瀵瑰簲X=A\B銆俋A=B瀵瑰簲X=B/A缁撴灉锛歑 = 6.0000 -0.2500 -3.0000 0.5000 -1.0000 -0.2500
绛旓細鍙互锛岃鐢ㄧ鍙疯繍绠楃殑鍔熻兘
绛旓細鍏堝皢xP=0涓よ竟鍙栬浆缃紝寰桺'x'=0,姹傚嚭x'鍐嶈浆缃洖鏉ユ眰鍑簒.褰鐭╅樀鏂圭▼P'x'=0,rank(P)=r<n鏃,鍗砅鏄寮傜殑锛岃繖鏃舵柟绋嬫湁鏃犵┓澶氫釜瑙,鐢∕ATLAB鍙眰鍑哄畠鐨勪竴涓熀鏈В绯伙紝 鍩烘湰瑙ｇ郴鐨勭嚎鎬х粍鍚堝氨鏄畠鐨勯氳В锛堜竴鑸В锛夎屾眰鍩烘湰瑙ｇ郴鐢matlab 涓鐨勫懡浠 x=null(P', r )鍗冲彲.鍏朵腑:r=rank(P)....
绛旓細鍙傝冧唬鐮侊細鍙傛暟瀹氫箟m1 = 1; m2 = 2; m3 = 3;k1 = 4; k2 = 5; k3 = 6;c1 = 7; c2 = 8; c3 = 9;Qy = 10;寰垎鏂圭▼K = [k1 -k1 0; -k1 k1+k2 -k2; 0 -k2 k2+k3];C = [c1 -c1 0; -c1 c1+c2 -c2; 0 -c2 c2+c3]dy = @(t,y)[y(4:6); (-K*y(...
绛旓細涓ゅ紡鐩稿姞鍙緱锛孉=(C+D)/2 鎯冲噺涓猴細B=(C-D)/2
绛旓細璇鐭╅樀鏂圭▼锛屽彲浠olve鍑芥暟鏉ユ眰瑙ｃ備緥濡傦細a+b+c=3,d+e+f=5,F1=10;F2=20;鍒 X=[i1;i2];A=[3+i1;5+i2];B=[10;20];[i1,i2]=solve(X-A.*B)姹傝В缁撴灉 i1=-10/3;i2=-100/3
绛旓細浣鐨闂鎻愬緱濂藉儚鏈夌偣闂锛屼袱涓彉閲忎笁涓鏂圭▼锛岀敤solve鍑芥暟鐩存帴瑙ｅ彧鑳芥槸鏃犺В锛岃岀敤fsolve鍑芥暟姹傝В锛屽垯寰楀埌瑙ｆ槸浠や汉璐硅В鐨勶紙涓嶆亽绛夛級銆傚浣犳妸闂杞崲鎴愪紭鍖栬В锛屽嵆宸茬煡P鍊硷紝姹侭銆丩銆丯鍊笺傚埌鍙互鐢╨sqnonlin锛堬級鍑芥暟姹傚緱銆傚綋P=[30;20;40]鏃讹紝B=142.57掳锛孡=47.86掳锛孨=[1.9012;5.5885;-4...
绛旓細clear a0=[2 3;6 4];a1=[-10 1;-20 3];a2=[1 3;5 2];syms p1 p2 p3 p4 p;p=[p1 p2;p3 p4]浜屾鐭╅樀鏂圭▼鏄細p^2*a2+p*a1+a0=0;eq=p^2*a2+p*a1+a0;[p1,p2,p3,p4]=solve(eq(1,1),eq(1,2),eq(2,1),eq(2,2),p1,p2,p3,p4);p1=double(p1);p2=double(...
