matlab求矩阵方程组
重新整理一下你的方程。你的第一个方程是
（pi0,pi1）=(pi0,pi1)*[A00,A01;A10,R*A2]
取个转置，就是
（pi0,pi1）(T)
=
[A00,A10;A01,R*A2]*(pi0,pi1)(T)
相当于求解线性方程组
Ax
=
0
A
=
[A00
-
1,A10;A01,R*A2
-
1]
x
=
(pi0,pi1)(T)
你还有一个约束条件，对吧：
1=pi0*e+pi1*inv(I-R)*e;
这相当于线性方程组的第三行。
所以相当于：
//
代码如下
A
=
[A00
-
1,A10;
A01,R*A2
-
1;
e,*inv(I-R)*e];
b
=
[0;
0;
1];
x
=
A\b;
//
Matlab自动求解线性方程组的语句，没记错的话，对于这种超定方程组，应该也能自动搞定，算出最小二乘解
//
代码结束
最后，求出来的x是个2*1的向量，两个分量分别对应pi0和pi1
你试试吧
先将xp=0两边取转置，得p'x'=0,求出x'再转置回来求出x.
当矩阵方程p'x'=0,rank(p)=r
评论
0
0
0
加载更多
绛旓細杩欎釜鐭╅樀绛変簬 [ R 0 R -X][ 0 R X R][ R -X R 0][ X R 0 R]鐨勯嗛樀涔樹互[a b c d]'matlab鐨璇濓紝灏辨眰涓婇潰4x4鐭╅樀鐨勯嗗氨琛岋紝鍗崇煩闃电殑-1娆℃柟锛岀劧鍚庡啀涔樹互[a b c d]'锛屽緱鍒颁竴涓4x1鐭╅樀灏辨槸[x1 x2 x3 x4]'銆傚鏋滅郴鏁版病鏈夐嗛樀锛岃鏄庤鏂圭▼涓嶅畬澶囷紝鏃犺В ...
绛旓細MATLAB锛圡atrix Laboratory鍦ㄧ嚎鎬т唬鏁颁腑锛孧ATLAB鎻愪緵浜嗚澶氬己澶х殑宸ュ叿鍜屽嚱鏁帮紝浣垮緱澶嶆潅鐨璁＄畻鍜屽垎鏋愬彉寰楃畝鍗曘備互涓嬫槸涓浜汳ATLAB鍦ㄧ嚎鎬т唬鏁颁腑鐨勫簲鐢細鐭╅樀杩愮畻锛歁ATLAB涓虹煩闃佃繍绠楁彁渚涗簡涓撻棬鐨勫嚱鏁帮紝濡傜煩闃典箻娉曪紙*锛夈佺煩闃靛姞娉曪紙+锛夈佺煩闃佃浆缃紙'锛夈佺煩闃电殑閫嗭紙inv锛夌嚎鎬鏂圭▼缁勬眰瑙锛歁ATLAB浣跨敤宸﹂櫎杩愮畻绗︼紙\锛夊拰...
绛旓細1. 璁続x = b锛屾眰x锛岋紙x鍜宐閮芥槸鍚戦噺锛夊垯锛歺 = A\b 灏辫繖涔堢畝鍗曘傛兂涓嶅嚭鏇村鏉傜殑鏂规硶浜嗐2. 浣犱笂闈㈢殑鐪嬭捣鏉ユ槸涓綈娆℃柟绋嬶紝浣嗘槸鍥犱负b1 b2 b3宸茬煡锛屽彲浠ョЩ鍒扮瓑鍙峰彟涓杈癸紱鍐嶆妸V3=V4鏁寸悊鍒鐭╅樀鏂圭▼閲屾垨鑰呭共鑴嗕粠鏂圭▼涓幓鎺塚4锛屽氨鍙互鐢ㄤ笂闈㈡柟娉曡В浜嗐
绛旓細鏍规嵁杈撳叆鐨鐭╅樀绫诲瀷鍜屽ぇ灏忥紝MATLAB浼氳嚜鍔ㄩ夋嫨鏈鍚堥傜殑绠楁硶鏉姹傝В鏂圭▼缁锛屼粠鑰屾彁楂璁＄畻鏁堢巼銆2銆佺█鐤忕煩闃垫敮鎸侊細瀵逛簬绋鐤忕煩闃碉紙鐭╅樀涓ぇ閮ㄥ垎鍏冪礌涓洪浂锛夌殑绾挎鏂圭▼缁勬眰瑙锛孧ATLAB閲囩敤浜嗕笓闂ㄤ紭鍖栫殑绠楁硶锛屽埄鐢ㄧ█鐤忕煩闃电殑鐗规畩缁撴瀯锛岄伩鍏嶄簡涓嶅繀瑕佺殑璁＄畻鍜屽瓨鍌ㄥ紑閿锛屼粠鑰屽姞蹇眰瑙ｉ熷害銆3銆佸绾跨▼璁＄畻锛歁ATLAB鍙互鍦ㄥ鏍...
绛旓細鍏堝皢xP=0涓よ竟鍙栬浆缃紝寰桺'x'=0,姹傚嚭x'鍐嶈浆缃洖鏉ユ眰鍑簒.褰鐭╅樀鏂圭▼P'x'=0,rank(P)=r<n鏃,鍗砅鏄寮傜殑锛岃繖鏃舵柟绋嬫湁鏃犵┓澶氫釜瑙,鐢∕ATLAB鍙眰鍑哄畠鐨勪竴涓熀鏈В绯伙紝 鍩烘湰瑙ｇ郴鐨勭嚎鎬х粍鍚堝氨鏄畠鐨勯氳В锛堜竴鑸В锛夎屾眰鍩烘湰瑙ｇ郴鐢╩atlab 涓鐨勫懡浠 x=null(P', r )鍗冲彲.鍏朵腑:r=rank(P)....
绛旓細浣跨敤matlab姹傝В绾挎鏂圭▼缁锛屽彲浠ヨ繖鏍疯В銆傞鍏堬紝鍐欏嚭绾挎ф柟绋嬬粍鐨勭郴鏁般傚嵆 A=[1 -9 -10;-9 1 -5;8 7 1];鍏朵簩锛屽啓鍑虹嚎鎬ф柟绋嬬粍鐨勫父鏁伴」绯绘暟銆傚嵆 B=[1;0;4];鐒跺悗锛岀敤鐭╅樀宸﹂櫎娉曪紝姹傚嚭X=[x1,x2,x3] 鐨勮В銆傚嵆 X=A\B 杩愯缁撴灉 ...
绛旓細b)闄ゆ硶锛氫竴鑸湪瑙ｇ嚎鎬鏂圭▼缁鏃朵細鐢ㄥ埌銆倄=a\b 濡傛灉ax=b锛屽垯 x=a\b鏄鐭╅樀鏂圭▼鐨勮В銆倄=b/a 濡傛灉xa=b, 鍒檟=b/a鏄煩闃垫柟绋嬬殑瑙ｃ傝浆缃細杞疆鏃讹紝鐭╅樀鐨勭涓琛屽彉鎴愮涓鍒楋紝绗簩琛屽彉鎴愮浜屽垪锛屻傘傘倄=a.'姹傞嗭細瑕佹眰鐭╅樀涓烘柟闃点傝繖鍦ㄧ煩闃佃繍绠椾腑寰堝父鐢ㄣ倄=inv(a)
绛旓細浣犳妸y褰撴垚鑷彉閲忥紝鐢诲浘锛屽彲浠ユ柟渚夸粠y=-1:.1:1銆傚彂鐜皒鍙栧硷紝鏈澶ф槸0.2602銆傛墍浠ュ彲浠ュ緢鏄庢樉锛屽彂鐜板彧鑳藉彇0.1803 鍜 0.2067 鎵浠ヤ綘鐨勭▼搴忓彧鐢ㄥ啓涓涓惊鐜姞鍒ゆ柇鍊硷紝鍗冲皬浜0.2602鍗冲彲銆傝皟鐢╯olve鍑芥暟 瑙ｅ嚭0.1882847702315369927482666845352鍜-0.22384699275570656092042225557401 鍏跺疄浣犺繖涓〃杈惧紡杩樻槸鏈夐棶棰樼殑锛...
绛旓細A=[ 2 3 1 3 2 5 1 3 2];b1=[5 7 6]';b2=[8 9 10]';閫鐭╅樀 rats(inv(A)*b1)ans = -1/9 13/9 8/9 rats(inv(A)*b2)ans = -5/6 17/6 7/6 solve:res=solve('2*x+3*y+z=5','3*x+2*y+5*z=7','x+3*y+2*z=6');[res.x res.y res.z]ans = [...
绛旓細缁欎綘涓涓緥瀛愶紝鏉ヨ鏄庡浣鐢∕atlab姹绾挎鏂圭▼缁鐨勯氳В銆>> a=[1 -1 1 -1;-1 1 1 -1;2 -2 -1 1]; %绾挎ф柟绋嬬粍鐨勭郴鏁鐭╅樀 >> b=[1;1;-1]; % 甯稿垪鍚戦噺 >> [rank(a) rank([a,b])]ans = 2 2 %绉╃浉绛変笖灏忎簬4锛岃鏄庢湁鏃犵┓澶氳В >> rref([a,b]) %绠鍖栬闃舵...
