怎么用matlab解这种矩阵方程？ 用matlab解矩阵方程组的代码怎么写？
\u7528matlab\u89e3\u77e9\u9635\u65b9\u7a0b1\u3001\u52a0\u51cf\u6cd5\u7684\u547d\u4ee4\u5f88\u7b80\u5355\uff0c\u76f4\u63a5\u7528\u52a0\u6216\u8005\u51cf\u53f7\u5c31\u53ef\u4ee5\u4e86\u3002\u5982\uff1ac=a+bd=a-b\u3002
2\u3001\u4e00\u822c\u4e58\u6cd5\uff1ac=a*b,\u8981\u6c42a\u7684\u5217\u6570\u7b49\u4e8eb\u7684\u884c\u6570\u3002\u5982\u679ca,b\u662f\u4e00\u822c\u7684\u5411\u91cf\uff0c\u5982a=[1,2,3] b=[3,4,5]\u70b9\u79ef\uff1adot(a,b), \u53c9\u79ef\uff1across\uff08a,b)\u5377\u79ef\uff1aconv(a,b)\u3002
3\u3001x=a\b\u5982\u679cax=b\uff0c\u5219 x=a\b\u662f\u77e9\u9635\u65b9\u7a0b\u7684\u89e3\u3002x=b/a\u5982\u679cxa=b, \u5219x=b/a\u662f\u77e9\u9635\u65b9\u7a0b\u7684\u89e3\u3002
4\u3001\u8f6c\u7f6e\u65f6\uff0c\u77e9\u9635\u7684\u7b2c\u4e00\u884c\u53d8\u6210\u7b2c\u4e00\u5217\uff0c\u7b2c\u4e8c\u884c\u53d8\u6210\u7b2c\u4e8c\u5217\uff0c\u3002\u3002\u3002x=a\u3002
5\u3001\u6c42\u9006\uff1a\u8981\u6c42\u77e9\u9635\u4e3a\u65b9\u9635\u3002\u8fd9\u5728\u77e9\u9635\u8fd0\u7b97\u4e2d\u5f88\u5e38\u7528\u3002x=inv(a)\u3002\u8fd9\u51e0\u79cd\u65b9\u5f0f\u90fd\u53ef\u4ee5\u89e3\u77e9\u9635\u65b9\u7a0b\u3002
\u5148\u5c06xP=0\u4e24\u8fb9\u53d6\u8f6c\u7f6e\uff0c\u5f97P'x'=0,\u6c42\u51fax'\u518d\u8f6c\u7f6e\u56de\u6765\u6c42\u51fax.
\u5f53\u77e9\u9635\u65b9\u7a0bP'x'=0,rank(P)=r<n\u65f6,\u5373P\u662f\u5947\u5f02\u7684\uff0c\u8fd9\u65f6\u65b9\u7a0b\u6709\u65e0\u7a77\u591a\u4e2a\u89e3,\u7528MATLAB\u53ef\u6c42\u51fa\u5b83\u7684\u4e00\u4e2a\u57fa\u672c\u89e3\u7cfb\uff0c
\u57fa\u672c\u89e3\u7cfb\u7684\u7ebf\u6027\u7ec4\u5408\u5c31\u662f\u5b83\u7684\u901a\u89e3\uff08\u4e00\u822c\u89e3\uff09
\u800c\u6c42\u57fa\u672c\u89e3\u7cfb\u7528matlab
\u4e2d\u7684\u547d\u4ee4
x=null(P',
r
)\u5373\u53ef.\u5176\u4e2d:r=rank(P).
w1+2*w2+w3 = 0
w1+3*w2+2*w3 = 0
2*w1+w2+4*w3 = 0
之后用fsolve函数
[w1, w2, w3] = fsolve(w1+2*w2+w3 == 0, w1+3*w2+2*w3 == 0, 2*w1+w2+4*w3 == 0, w1, w2, w3);
绛旓細1銆佸姞鍑忔硶鐨勫懡浠ゅ緢绠鍗曪紝鐩存帴鐢ㄥ姞鎴栬呭噺鍙峰氨鍙互浜銆傚锛歝=a+bd=a-b銆2銆佷竴鑸箻娉曪細c=a*b,瑕佹眰a鐨勫垪鏁扮瓑浜巄鐨勮鏁般傚鏋渁,b鏄竴鑸殑鍚戦噺锛屽a=[1,2,3] b=[3,4,5]鐐圭Н锛歞ot(a,b), 鍙夌Н锛歝ross锛坅,b)鍗风Н锛歝onv(a,b)銆3銆亁=a\b濡傛灉ax=b锛屽垯 x=a\b鏄煩闃垫柟绋嬬殑瑙ｃ...
绛旓細涔嬪悗鐢╢solve鍑芥暟 [w1, w2, w3] = fsolve(w1+2*w2+w3 == 0, w1+3*w2+2*w3 == 0, 2*w1+w2+4*w3 == 0, w1, w2, w3);
绛旓細鐪嬭捣鏉瑙ｆ硶鏄鐨勩備綘鍙互鐢ㄥ緱鍒扮殑X杩涜楠岃瘉锛氱湅鐪婣*X*B=C鏄惁鎴愮珛銆傚鏋滈獙璇佹纭紝閭ｈ偗瀹氬氨涓嶆槸鐭╅樀璁＄畻鐨勯棶棰樹簡銆
绛旓細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(...
绛旓細鍏堝皢xP=0涓よ竟鍙栬浆缃紝寰桺'x'=0,姹傚嚭x'鍐嶈浆缃洖鏉ユ眰鍑簒.褰鐭╅樀鏂圭▼P'x'=0,rank(P)=r<n鏃,鍗砅鏄寮傜殑锛岃繖鏃舵柟绋嬫湁鏃犵┓澶氫釜瑙,鐢∕ATLAB鍙眰鍑哄畠鐨勪竴涓熀鏈В绯伙紝 鍩烘湰瑙ｇ郴鐨勭嚎鎬х粍鍚堝氨鏄畠鐨勯氳В锛堜竴鑸В锛夎屾眰鍩烘湰瑙ｇ郴鐢╩atlab 涓殑鍛戒护 x=null(P', r )鍗冲彲.鍏朵腑:r=rank(P)....
绛旓細棰樹富鐨勯棶棰樺彲浠ョ敤for鍙屽惊鐜鍙ュ拰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 ...
绛旓細1526*p21+.8330*p22][ .1526*p11+.8330*p21+.5230*p31+.6586*p32+.8636*p33, .1526*p12+.8330*p22+1.172*p32+.5676*p31+.7918*p33, .1526*p13+.8330*p23+.3837*p33+.1526*p31+.8330*p32+.1919]鐒跺悗鐢╯olve瑙9涓柟绋嬬粍锛屾眰寰桺锛屼笉杩囷紝淇鸿鐫锛岃繖鏂规硶澶銆
绛旓細matlab涓殑sqrtm(Y)鏄壘杩欐牱鐨鐭╅樀X,瀹冩弧瓒砐*X=Y锛屽苟涓擷鐨勭壒寰佸奸兘鍏锋湁闈炶礋瀹為儴锛岃屽鏋淴鐨勭壒寰佸间腑鏈変换鎰忎竴涓湁璐熺殑瀹炴暟鏍癸紝缁撴灉灏嗕細寰堝鏉傦紝鎵浠atlab鍙墽琛屽苟涓旀樉绀虹壒寰佸煎潎鏈夐潪璐熷疄閮ㄧ殑X銆傚鏋淵鏄竴涓寮傜煩闃碉紝瀹冨彲鑳芥病鏈夋柟鏍广傚洜鑰屽鏋MATLAB鎺㈡煡鍒颁竴涓寮傜煩闃碉紝灏变細鏄剧ずwarning....
绛旓細>> syms x1 x2 x3 >> A =[1,2,3;4,5,6;3,2,1];B=[x1;4;2];C=[4;x2;x3];[x1,x2,x3]=solve(A*B-C,'x1','x2','x3')x1 = -10 x2 = -8 x3 = -20
绛旓細a=[1 0;0 2];b=[1;3];syms x1 x2;y=a*[x1;x2]+b;s=solve(y(1),y(2),'x1','x2');s.x1 %w2鍊 s.x2 %w2鍊
