[急]matlab求解矩阵方程！！我看不懂matlab的英文不会编程！！！大神快来救我！！！ MATLAB7.0求解！！！！大神快来 跪求啊没时间了
\u6c42\u89e3\u7b2c\u4e03\u9898\uff01Matlab\u5927\u795e\u5feb\u6765\uff01%%
clear;
clc;
a = [0.5:0.5:4.5, 5:8]; %% a=5 ,error
figure
ezplot('x^2/0.5+y^2/(25-0.5^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/1+y^2/(25-1^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/1.5^2+y^2/(25-1.5^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/2^2+y^2/(25-2.5^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/2.5^2+y^2/(25-2.5^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/3^2+y^2/(25-3^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/3.5^2+y^2/(25-3.5^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/4.0^2+y^2/(25-4.0^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/4.5^2+y^2/(25-4.5^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/6^2+y^2/(25-6^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/7^2+y^2/(25-7^2) =1',[-10,10,-10,10])
hold on
ezplot('x^2/8^2+y^2/(25-8^2) =1',[-10,10,-10,10])
\u95ee\u9898\u51fa\u5728p\u4e0a\uff0c\u4e24\u79cd\u6539\u6cd5\uff1a
\u65b9\u6cd51:
m\u6587\u4ef6\uff1a\u7b2c1\u884c\u52a0\u81ea\u53d8\u91cfp function [f,k]=conv_m(f1,k1,f2,k2,p)
m\u6587\u4ef6\uff1a\u53bb\u6389\u7b2c10\u884cp=input('p=');
\u8f93\u5165\u4e0d\u53d8\uff1ap=0.5;k1=0:p:2;f1=0.5*k1;k2=k1;f2=f1;[f,k]=conv_m(f1,f2,k1,k2,p)
\u65b9\u6cd52:
m\u6587\u4ef6\uff1a\u7b2c1\u884c\u4e0e\u7b2c10\u884c\u4e0d\u53d8
\u8f93\u5165\u6539\u4e3a\uff1ak1=0:p:2;f1=0.5*k1;k2=k1;f2=f1;[f,k]=conv_m(f1,f2,k1,k2)
\u8fd0\u884c\u540e\u518d\u8f93\u5165p\u503c
使用solve即可解决解方程组的问题。
[u2 u3 i1]=solve('(1/2+1/4)*6-1/2*u2-1/4*u3=i1','-1/2*6+(1/2+1+1/3)*u2-1/3*u3=0','-1/4*6-1/3*u2+(1/3+1/4)*u3=1/3*i1')
运行结果如图：
程序见附件。
绛旓細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鏄煩闃垫柟绋嬬殑瑙ｃ...
绛旓細棣栧厛锛岃鎴戜滑鑱氱劍浜庤繛缁郴缁熺殑绋冲畾瀹堟姢鑰呪斺擫yapunov鏂圭▼銆傚畠婧愪簬涓涓叧閿殑浜屾鍨嬶紙鑳介噺锛夊嚱鏁帮紝鍗砎 = x'Px锛屽叾涓璓涓烘瀹鐭╅樀銆傚V鐨勫鏁帮紝鍗˙V = x'(A'PA + Q)x锛岃嫢鎭掑皬浜庨浂锛岄偅涔堢郴缁熶究澶勪簬绋冲畾鐘舵併傜ǔ瀹氱殑鏍囧織鏄瓨鍦ㄥ绉版瀹氱殑鐭╅樀Q锛屼娇寰桝'PA - Q鏈夎В锛岃繖灏辨槸钁楀悕鐨勮繛缁璍yapuno...
绛旓細鍋囪杩欓噷闇瑕姹傝В鐨鏂圭▼缁勫涓嬶細Matlab璁＄畻婧愮爜锛堜唬鏁板舰寮忥級 res=solve('5*a+5327*b=271.4','5327*a+7277699*b=369320','a','b'); res.a % 杈撳嚭璁＄畻缁撴灉a res.b % 杈撳嚭璁＄畻缁撴灉b Matlab璁＄畻婧愮爜锛鐭╅樀褰㈠紡锛 A=[5,5327;5327,7277699]; B=[271.4;369320]; X=inv(A)*B % inv鏄煩...
绛旓細棰樹富鐨勯棶棰樺彲浠ョ敤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 ...
绛旓細1銆佹墦寮matlab锛岄鍏堝畾涔夊彉閲弜锛歴yms x;2銆乵atlab涓璼olve鍑芥暟鐨勬牸寮忔槸solve(f(x), x)锛姹傝В鐨勬槸f(x) = 0鐨勮В銆傜涓涓緥瀛愶紝姹傝В鏈甯歌鐨勪竴鍏冧簩娆鏂圭▼x^2-3*x+1=0:solve(x^2-3*x+1,x)锛岃В鍑虹殑缁撴灉鐢ㄧ簿纭殑鏍瑰紡琛ㄧず銆3銆matlab瑙鍑虹殑鏍逛笉浠呭寘鍚疄鏍癸紝涔熷寘鍚鏍癸紝渚嬪姹傝В涓夋鏂圭▼x^...
绛旓細浣犵殑绗竴涓鏂圭▼鏄 锛坧i0,pi1锛=(pi0,pi1)*[A00,A01;A10,R*A2]鍙栦釜杞疆锛屽氨鏄 锛坧i0,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...
绛旓細璇鐭╅樀鏂圭▼锛屽彲浠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
绛旓細渚嬶細姹傝В鏂圭▼缁剎+y=1, x-11y=5 eq=@(x)[x(1)+x(2)-1;x(1)-11*x(2)-5][sol,fval]=fsolve(eq,[1,1])杩欓噷瀵逛簬鏂圭▼鐨勭殑杈撳叆闇瑕侀噰鐢鐭╅樀鐨勫舰寮忥紝鍏朵腑x(1)浠ｈ〃x锛寈(2)浠ｈ〃y銆傛湁鏃跺欏彉閲忚緝澶氭椂鍙兘浼氬鏄撴贩娣嗭紝杩欓噷鎻愪緵鍙︿竴绉嶆柟娉曪紝閲囩敤绗﹀彿鍙橀噺褰㈠紡鍐嶅埄鐢matlabFunction杞寲涓哄嚱鏁...
绛旓細.1:1銆傚彂鐜皒鍙栧硷紝鏈澶ф槸0.2602銆傛墍浠ュ彲浠ュ緢鏄庢樉锛屽彂鐜板彧鑳藉彇0.1803 鍜 0.2067 鎵浠ヤ綘鐨勭▼搴忓彧鐢ㄥ啓涓涓惊鐜姞鍒ゆ柇鍊硷紝鍗冲皬浜0.2602鍗冲彲銆傝皟鐢╯olve鍑芥暟 瑙ｅ嚭0.1882847702315369927482666845352鍜-0.22384699275570656092042225557401 鍏跺疄浣犺繖涓〃杈惧紡杩樻槸鏈夐棶棰樼殑锛屾瘮濡倅=0.閭ｄ箞浣犵殑鏋侀檺鏄棤绌峰ぇ銆
绛旓細鐪嬭捣鏉ヨВ娉曟槸瀵圭殑銆備綘鍙互鐢ㄥ緱鍒扮殑X杩涜楠岃瘉锛氱湅鐪婣*X*B=C鏄惁鎴愮珛銆傚鏋滈獙璇佹纭紝閭ｈ偗瀹氬氨涓嶆槸鐭╅樀璁＄畻鐨勯棶棰樹簡銆
