matlab求解矩阵方程 用matlab解矩阵方程
matlab\u6c42\u89e3\u77e9\u9635\u65b9\u7a0b\u7ec4\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).
1\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
相当于求解KA=Y
K的维数是6*9
然后转置 A^T K^T=Y^T
得到K^T=Y^T\A^T
所以K=(Y^T\A^T)^T
你可以把数据给我，我给你算
K=Y*(X^-1)*(B^-1)*(V^-1)
或者直接用除法
直接用除法：
H = V*B*X;
K = Y/H;
绛旓細1. 璁続x = b锛屾眰x锛岋紙x鍜宐閮芥槸鍚戦噺锛夊垯锛歺 = A\b 灏辫繖涔堢畝鍗曘傛兂涓嶅嚭鏇村鏉傜殑鏂规硶浜嗐2. 浣犱笂闈㈢殑鐪嬭捣鏉ユ槸涓綈娆℃柟绋嬶紝浣嗘槸鍥犱负b1 b2 b3宸茬煡锛屽彲浠ョЩ鍒扮瓑鍙峰彟涓杈癸紱鍐嶆妸V3=V4鏁寸悊鍒鐭╅樀鏂圭▼閲屾垨鑰呭共鑴嗕粠鏂圭▼涓幓鎺塚4锛屽氨鍙互鐢ㄤ笂闈㈡柟娉曡В浜嗐
绛旓細鏈変袱绉嶆柟娉,绗竴绉嶆柟娉曟槸鎵嬬畻鐨勬柟娉(鑰冭瘯鐨勬椂鍊欑敤鐨),鐢卞師寮忓瓙,鍙煡X=BA^-1灏咮鍐欎笂闈,A鍐欏湪涓嬭竟,鐒跺悗閫氳繃鍒楀彉鎹㈡妸A鍙樻垚鍗曚綅鐭╅樀E,鍙樻崲鏃禕涔熻窡鐫杩涜鍒楀彉鎹,褰揂鍙樻垚E鏃,B鐨勫尯鍩熷氨鏄墍姹傜殑X,杩欑鏂规硶瑕佹眰瀛︾敓蹇呴』鎺屾彙,鑰冭瘯鏃剁敤鐨勫氨鏄繖涓柟娉.鏂规硶浜:鍊熺敤璁＄畻杞欢鐨勬柟娉,鍦MATLAB閲,瀹氫箟A=(2...
绛旓細瑙傚療鍥1锛岄暱鏂瑰舰A锛氶暱锛14-2y锛屽锛3x锛孲A闃村奖锛(14-2y)3x锛42x-6xy锛1锛夋眰鍥1涓12涓皬闀挎柟褰y闈㈢НSxy=xy锛屾牴鎹鎰忥紝12Sxy锛婼A闃村奖锛婼B闃村奖锛14x10 12xy锛42x-6xy锛婼B闃村奖锛140 SB闃村奖锛140-6xy-42x锛2锛夋牴鎹鎰忥紝SA闃村奖锛峉B闃村奖锛112銆锛3锛夛紙1锛夛紙2锛変唬鍏ワ紙3锛夊緱锛42x-6xy...
绛旓細鐢╯=solve('A*x=B')鏄緱涓嶅埌s鍊肩殑銆傚簲璇ヤ娇鐢ㄤ笅鍒楀懡浠 s=A\B %鐭╅樀闄ゆ硶锛屽緱鍒版渶灏忎簩涔樿繎浼艰В銆
绛旓細浣犵殑闂鎻愬緱濂藉儚鏈夌偣闂锛屼袱涓彉閲忎笁涓鏂圭▼锛岀敤solve鍑芥暟鐩存帴瑙ｅ彧鑳芥槸鏃犺В锛岃岀敤fsolve鍑芥暟姹傝В锛屽垯寰楀埌瑙ｆ槸浠や汉璐硅В鐨勶紙涓嶆亽绛夛級銆傚浣犳妸闂杞崲鎴愪紭鍖栬В锛屽嵆宸茬煡P鍊硷紝姹侭銆丩銆丯鍊笺傚埌鍙互鐢╨sqnonlin锛堬級鍑芥暟姹傚緱銆傚綋P=[30;20;40]鏃讹紝B=142.57掳锛孡=47.86掳锛孨=[1.9012;5.5885;-4...
绛旓細鍦matlab涓紝濡傛灉A鏄彲閫鐭╅樀 AX=B鐨勮В鏄疉宸﹂櫎B,鍗 X=A\B XA=B鐨勮В鏄疉鍙抽櫎B, 鍗砐=B/A銆傚叿浣撳埌杩欓噷锛欰,B閮芥槸鍙嗙煩闃碉紝X=(A\C)/B
绛旓細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鍊
绛旓細e1=[1,0,0];e2=[0,1,0];e3=[0,0,1];E1=[sqrt(3)/2,1/2,0];E2=[-1/2,sqrt(3)/2,0];E3=[0,0,1];A=[E1',E2',E3',zeros(3,6);zeros(3,3),E1',E2',E3',zeros(3,3);zeros(3,6),E1',E2',E3';];b=[e1,e2,e3]';C=A\b ...
绛旓細鏍规嵁杈撳叆鐨鐭╅樀绫诲瀷鍜屽ぇ灏忥紝MATLAB浼氳嚜鍔ㄩ夋嫨鏈鍚堥傜殑绠楁硶鏉姹傝В鏂圭▼缁勶紝浠庤屾彁楂樿绠楁晥鐜囥2銆佺█鐤忕煩闃垫敮鎸侊細瀵逛簬绋鐤忕煩闃碉紙鐭╅樀涓ぇ閮ㄥ垎鍏冪礌涓洪浂锛夌殑绾挎ф柟绋嬬粍姹傝В锛孧ATLAB閲囩敤浜嗕笓闂ㄤ紭鍖栫殑绠楁硶锛屽埄鐢ㄧ█鐤忕煩闃电殑鐗规畩缁撴瀯锛岄伩鍏嶄簡涓嶅繀瑕佺殑璁＄畻鍜屽瓨鍌ㄥ紑閿锛屼粠鑰屽姞蹇眰瑙ｉ熷害銆3銆佸绾跨▼璁＄畻锛歁ATLAB鍙互鍦ㄥ鏍...
绛旓細warning:explicit solution could not be found.杩欏彞璇濈殑鎰忔濇槸涓嶈兘鎵惧埌绮剧‘瑙ｏ紝涔熷氨鏄matlab閲岄潰solve鍑芥暟鐨勫唴閮ㄧ畻娉曞浜庝綘杩欎釜鏂圭▼缁勪笉閫傜敤锛屾墍浠ヤ綘杩欎釜鏂圭▼缁勪笉鑳界敤solve姹傝В 鑰屼笖锛岃瀵熶綘鐨勪笁涓柟绋嬶紝鏄剧劧鏄棤瑙ｇ殑锛岀涓寮忛櫎浠ョ浜屽紡鍙互寰楀嚭exp(k0)=533.34/534.48 绗簩寮忛櫎浠ョ涓夊紡鍙互寰楀嚭 exp...
