matlab求矩阵方程的数值解 用matlab解矩阵方程
MATLAB\u4e2d\u89e3\u65b9\u7a0b\u65f6\u6c42\u6570\u503c\u89e3\u600e\u4e48\u6c42\u4f60\u7684solve\u51fd\u6570\u91cc\u9762\u5df2\u7ecf\u662f\u5b57\u7b26\uff0c\u6ca1\u6709\u5fc5\u8981\u518d\u7528eval
\u5982\u679c\u4f60\u8981\u7cbe\u786e\u89e3\uff0c\u5c31\u76f4\u63a5\u628ax0=eval(a)\u5220\u6389\uff0c\u7ed3\u679c\u5c31\u662fa,\uff0c\u4f46\u662f\u8868\u8fbe\u5f0f\u5f88\u590d\u6742\uff0c\u4f60\u53ef\u4ee5\u52a0\u4e2adouble(a)\u5c31\u53ef\u4ee5\u5f97\u5230\u6570\u503c\u89e3\u4e86
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
你的问题提得好像有点问题，两个变量三个方程，用solve函数直接解只能是无解，而用fsolve函数求解，则得到解是令人费解的（不恒等）。
如你把问题转换成优化解，即已知P值，求B、L、N值。到可以用lsqnonlin（）函数求得。
当P=[30;20;40]时，B=142.57°，L=47.86°，N=[1.9012;5.5885;-4.8164]。
绛旓細宸茬煡6<x<30锛15<y<36,姹備笉瀹鏂圭▼ 2x+5y=126鐨勬暣鏁拌В銆MATLAB绋嬪簭濡備笅锛歺=5:29;y=14:35;[x,y]=meshgrid(x,y);z=2*x+5*y;k=find(z==126);x(k),y(k)杈撳嚭涓猴細ans= 8 13 18 23 28 ans= 22 20 18 16 14 鍗虫柟绋嬫湁5缁勮В锛(8,22),(13,20),(18,18),(23,16),(28,...
绛旓細鏂规硶涓銆佸埄鐢ㄤ笁寰幆璇彞鍜屽垽鏂鍙ワ紝姹傝В鍏舵墍鏈夋暣鏁拌В 鏂规硶浜屻佹暣鏁扮嚎鎬ц鍒掑垎鏀畾鐣屾硶锛屾眰瑙ｅ叾鎵鏈夋暣鏁拌В 涓嬮潰鎸夌涓绉嶆柟娉曟眰瑙ｏ紝鍏惰繃绋嬩负 m=50;n=50;l=50;disp('婊¤冻鏉′欢锛2x+3y+4z=100鐨勮В')for i=1:m for j=1:n for k=1:l x=i;y=j;z=k;xyz=100-(2*x+3*y+4*z);if x...
绛旓細浣犵殑闂鎻愬緱濂藉儚鏈夌偣闂锛屼袱涓彉閲忎笁涓鏂圭▼锛岀敤solve鍑芥暟鐩存帴瑙ｅ彧鑳芥槸鏃犺В锛岃岀敤fsolve鍑芥暟姹傝В锛屽垯寰楀埌瑙ｆ槸浠や汉璐硅В鐨勶紙涓嶆亽绛夛級銆傚浣犳妸闂杞崲鎴愪紭鍖栬В锛屽嵆宸茬煡P鍊硷紝姹侭銆丩銆丯鍊笺傚埌鍙互鐢╨sqnonlin锛堬級鍑芥暟姹傚緱銆傚綋P=[30;20;40]鏃讹紝B=142.57掳锛孡=47.86掳锛孨=[1.9012;5.5885;-4...
绛旓細3銆亁=a\b濡傛灉ax=b锛屽垯 x=a\b鏄煩闃垫柟绋嬬殑瑙ｃ倄=b/a濡傛灉xa=b, 鍒檟=b/a鏄煩闃垫柟绋嬬殑瑙銆4銆佽浆缃椂锛岀煩闃电殑绗竴琛屽彉鎴愮涓鍒楋紝绗簩琛屽彉鎴愮浜屽垪锛屻傘傘倄=a銆5銆佹眰閫嗭細瑕佹眰鐭╅樀涓烘柟闃点傝繖鍦ㄧ煩闃佃繍绠椾腑寰堝父鐢ㄣ倄=inv(a)銆傝繖鍑犵鏂瑰紡閮藉彲浠ヨВ鐭╅樀鏂圭▼銆
绛旓細matlab鏂圭▼姹傝В鏈変袱绉嶆柟娉曪紝绗﹀彿姹傝В鍜屾暟鍊兼眰瑙ｃ1.solve 閫氬父鍦ㄤ笉纭畾鏂圭▼鏄惁鏈夌鍙疯В鐨勬椂鍊欙紝鎺ㄨ崘鍏堜娇鐢╯olve杩涜灏濊瘯锛屽洜涓簊olve鐩告瘮浜庢暟鍊兼眰瑙ｆ潵璇达紝瀹冧笉闇瑕佹彁渚涘垵鍊硷紝骞朵笖涓鑸儏鍐典笅鑳藉寰楀埌鏂圭▼鐨勬墍鏈夎В銆傚浜庝竴浜涚畝鍗曠殑瓒呰秺鏂圭▼锛宻olve杩樿兘澶熻嚜鍔ㄨ皟鐢ㄦ暟鍊艰绠楃郴缁熺粰鍑轰竴涓暟鍊艰В銆俿olve鐨勮皟鐢ㄥ舰寮忥細S =...
绛旓細棰樹富鐨勯棶棰樺彲浠ョ敤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 ...
绛旓細璇鐭╅樀鏂圭▼锛屽彲浠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
绛旓細棣栧厛锛岃鎴戜滑鑱氱劍浜庤繛缁郴缁熺殑绋冲畾瀹堟姢鑰呪斺擫yapunov鏂圭▼銆傚畠婧愪簬涓涓叧閿殑浜屾鍨嬶紙鑳介噺锛夊嚱鏁帮紝鍗砎 = x'Px锛屽叾涓璓涓烘瀹鐭╅樀銆傚V鐨勫鏁帮紝鍗˙V = x'(A'PA + Q)x锛岃嫢鎭掑皬浜庨浂锛岄偅涔堢郴缁熶究澶勪簬绋冲畾鐘舵併傜ǔ瀹氱殑鏍囧織鏄瓨鍦ㄥ绉版瀹氱殑鐭╅樀Q锛屼娇寰桝'PA - Q鏈夎В锛岃繖灏辨槸钁楀悕鐨勮繛缁璍yapuno...
绛旓細.1:1銆傚彂鐜皒鍙栧硷紝鏈澶ф槸0.2602銆傛墍浠ュ彲浠ュ緢鏄庢樉锛屽彂鐜板彧鑳藉彇0.1803 鍜 0.2067 鎵浠ヤ綘鐨勭▼搴忓彧鐢ㄥ啓涓涓惊鐜姞鍒ゆ柇鍊硷紝鍗冲皬浜0.2602鍗冲彲銆傝皟鐢╯olve鍑芥暟 瑙ｅ嚭0.1882847702315369927482666845352鍜-0.22384699275570656092042225557401 鍏跺疄浣犺繖涓〃杈惧紡杩樻槸鏈夐棶棰樼殑锛屾瘮濡倅=0.閭ｄ箞浣犵殑鏋侀檺鏄棤绌峰ぇ銆
绛旓細鐢╯=solve('A*x=B')鏄緱涓嶅埌s鍊肩殑銆傚簲璇ヤ娇鐢ㄤ笅鍒楀懡浠 s=A\B %鐭╅樀闄ゆ硶锛屽緱鍒版渶灏忎簩涔樿繎浼艰В銆
