matlab求解矩阵方程，兄弟们帮帮小弟吧
matlab \u89e3\u77e9\u9635\u65b9\u7a0b\u8fd9\u4e2a\u95ee\u9898\u4ffa\u4e5f\u611f\u5174\u8da3\uff0c\u63d0\u4f9b\u4e00\u79cd\u65b9\u6cd5\uff0c\u4ec5\u4f9b\u53c2\u8003\u3002
clc;clear
P=sym('[p11 p12 p13;p21 p22 p23;p31 p32 p33]')
A=rand(3,3)
At=A'
Q=diag(diag(A))
f=At*P+P*A+Q
\u7ed3\u679c\uff1a
P =
[ p11, p12, p13]
[ p21, p22, p23]
[ p31, p32, p33]
A =
0.3311 0.5676 0.1526
0.6586 0.9805 0.8330
0.8636 0.7918 0.1919
At =
0.3311 0.6586 0.8636
0.5676 0.9805 0.7918
0.1526 0.8330 0.1919
Q =
0.3311 0 0
0 0.9805 0
0 0 0.1919
f =
[ .6622*p11+.6586*p21+.8636*p31+.6586*p12+.8636*p13+.3311, 1.312*p12+.6586*p22+.8636*p32+.5676*p11+.7918*p13, .5230*p13+.6586*p23+.8636*p33+.1526*p11+.8330*p12]
[ .5676*p11+1.312*p21+.7918*p31+.6586*p22+.8636*p23, .5676*p12+1.961*p22+.7918*p32+.5676*p21+.7918*p23+.9805, .5676*p13+1.172*p23+.7918*p33+.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]
\u7136\u540e\u7528solve\u89e39\u4e2a\u65b9\u7a0b\u7ec4\uff0c\u6c42\u5f97P\uff0c\u4e0d\u8fc7\uff0c\u4ffa\u89c9\u7740\uff0c\u8fd9\u65b9\u6cd5\u592a\u7b28\u3002
X=b\u4e58\u4ee5A\u7684\u9006\uff0cA\u7684\u9006\u53ef\u4ee5\u7528inv(A)\u89e3\u51fa
秩 rank(A)行列式 det(A)
迹 trace(A)
逆 inv(A)
行最简形矩阵 rref(A)
特征值和特征向量 [d,v] = eig(A)
解矩阵方程 AX=B A\B
哥帮不了你
绛旓細鍦matlab涓紝濡傛灉A鏄彲閫嗙煩闃 AX=B鐨勮В鏄疉宸﹂櫎B,鍗 X=A\B XA=B鐨勮В鏄疉鍙抽櫎B, 鍗砐=B/A銆傚叿浣撳埌杩欓噷锛欰,B閮芥槸鍙鐭╅樀锛X=(A\C)/B
绛旓細浣犳妸y褰撴垚鑷彉閲忥紝鐢诲浘锛屽彲浠ユ柟渚夸粠y=-1:.1:1銆傚彂鐜皒鍙栧硷紝鏈澶ф槸0.2602銆傛墍浠ュ彲浠ュ緢鏄庢樉锛屽彂鐜板彧鑳藉彇0.1803 鍜 0.2067 鎵浠ヤ綘鐨勭▼搴忓彧鐢ㄥ啓涓涓惊鐜姞鍒ゆ柇鍊硷紝鍗冲皬浜0.2602鍗冲彲銆傝皟鐢╯olve鍑芥暟 瑙ｅ嚭0.1882847702315369927482666845352鍜-0.22384699275570656092042225557401 鍏跺疄浣犺繖涓〃杈惧紡杩樻槸鏈夐棶棰樼殑锛...
绛旓細瑙x=y涓轰粈涔堣鐢ㄧ鍙峰瀷鍙橀噺锛涚洿鎺=锛圓^-1锛墆涓嶅氨濂戒簡 X = 0.4979 0.1445 0.0629 -0.0813
绛旓細end 杩愯缁撴灉 matlab瑙ｆ柟绋嬬粍lnx琛ㄧず鎴恖og(x)鑰宭gx琛ㄧず鎴恖og10(x)1-exp(((log(y))/x^0.5)/(x-1))1銆佽В鏂圭▼鏈杩戞湁澶氫汉闂浣鐢╩atlab瑙ｆ柟绋缁勭殑闂锛屽叾瀹炲湪matlab涓В鏂圭▼缁勮繕鏄緢鏂逛究鐨勶紝渚嬪锛屽浜庝唬鏁版柟绋嬬粍Ax=b(A涓虹郴鏁鐭╅樀锛闈炲寮)鐨姹傝В锛MATLAB涓湁涓ょ鏂规硶锛(1)x=inv(A)*b...
绛旓細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(...
绛旓細warning:explicit solution could not be found.杩欏彞璇濈殑鎰忔濇槸涓嶈兘鎵惧埌绮剧‘瑙ｏ紝涔熷氨鏄matlab閲岄潰solve鍑芥暟鐨勫唴閮ㄧ畻娉曞浜庝綘杩欎釜鏂圭▼缁勪笉閫傜敤锛屾墍浠ヤ綘杩欎釜鏂圭▼缁勪笉鑳界敤solve姹傝В 鑰屼笖锛岃瀵熶綘鐨勪笁涓鏂圭▼锛鏄剧劧鏄棤瑙ｇ殑锛岀涓寮忛櫎浠ョ浜屽紡鍙互寰楀嚭exp(k0)=533.34/534.48 绗簩寮忛櫎浠ョ涓夊紡鍙互寰楀嚭 exp...
绛旓細clear;clc;e=zeros(3);k=9;%k璁惧灏戝氨鏄疕锛坘锛夛紝N锛坘锛墂hile k~=0 k=k-1;ek=[0,0,0;1,0,0;0,1,0]*e*[0,0,0;1,0,0;0,1,0]'+eye(3);H=ek*[0;0;1]*[0,0,1]*ek'*[0;0;1]+1;e=(eye(3)-H*[0,0,1])*ek;N=[0,0,0;1,0,0;0,1,0]-H*[0,...
绛旓細浣跨敤`solve`鍑芥暟姹傝В鏂圭▼缁勶紝渚嬪`sol = solve`姹傝В寰楀埌鐨勮В灏嗗瓨鍌ㄥ湪缁撴瀯浣揱sol`涓2. 浣跨敤绾挎т唬鏁板嚱鏁拌В鏁板兼柟绋嬬粍 褰撴柟绋嬬粍鍖呭惈鍏蜂綋鏁板兼椂锛屽彲浠ヤ娇鐢MATLAB鐨绾挎т唬鏁板嚱鏁版潵姹傝В銆備緥濡傦紝浣跨敤`\`鎿嶄綔绗︽垨`linsolve`鍑芥暟鏉ユ眰瑙ｇ嚎鎬ф柟绋嬬粍銆傛楠ゅ涓嬶細鍒涘缓绯绘暟鐭╅樀鍜屽父鏁板悜閲忥紝渚嬪浣跨敤`A = [a b; ...
绛旓細X0=zeros(5,1);X=fsolve(@(X) hanshu(X),X0);
绛旓細姹傝В鐭╅樀鏂圭▼锛鍙互鍒╃敤鐭╅樀鐨勬ц川鏉ユ眰锛屽嵆 syms alpha1 alpha2 x y A =[1 5*cos(x) 0 -1;5 1 -1 0;3 6*cos(x) -1 0;6 3 0 -1];b=[0;0;1.0472;0];X=[alpha1;alpha2;x;y];S=vpasolve(A*X==b);alpha1=S.alpha1 alpha2=S.alpha2 x=S.x y=S.y 杩愯缁撴灉 ...
