哪位高手会用matlab解矩阵方程AXB+CXD=E,万分感激其中A,B,C,D...
先把方程变型成如下形式：aX+Xb=-c
这个是Sylvester方程,其中的a,b,c你自己转换吧,因为有矩阵逆的符号,我不好表示.（注意c前面有个负号）
然后你用MATLAB自带的工具箱中的lyap()函数求解.命令输入如下：
X=lyap(a,b,c)
绛旓細1.寤虹珛浣犵殑鐭╅樀锛2.姹傜壒寰佸煎拰鐗瑰緛鍚戦噺锛屾渶澶х壒寰佸 [u v]=eig(A);vmax=max(v)3.姹侰R锛孋I锛孯I CI=(umax-6)/(6-1);RI鏈夌幇鎴愮殑琛紝鍘绘煡琛 CR=CI/RI;4.鏉冮噸璁＄畻 n=find(v==max(v))w=u(:,n)/sum(u(:,n))
绛旓細姣忎釜鐭╅樀鍒嗗埆璁＄畻锛>> A=[1 1/9 1/3 1/5;9 1 7 8;3 1/7 1 5;5 1/8 1/5 1]; [V,D]=eig(A); [maxD,idx]=max(diag(D)); maxD,V(:,idx)maxD = 4.6389 ans = -0.0608 -0.9527 -0.2674 -0.1309 >> A=[1 4 3;1/4 1 1;1/3 1 1];>> [V...
绛旓細0.3805 0.2400 0.1202 0.1202 0.0848 0.0543 5.7620 姝鐭╅樀鐨勪竴鑷存у彲浠ユ帴鍙!CI= -0.0476 CR= -0.0378 渚濇锛涙潈閲嶃佹渶澶х壒寰佹牴銆丆I锛孋R锛 鍛靛懙 浣犲仛鐨 灞傛鍒嗘瀽娉曞惂锛 鏈杩戞瘯涓氳璁″垰鎺ヨЕ杩囷紝鍛靛懙
绛旓細鐢卞叕寮忔帹涓涓嬪彲浠ュ緱鍒帮紝鐢ㄤ互涓鐭╅樀宸︿箻A鍙互寰楀埌B:
绛旓細A = 1.0000 5.0000 3.0000 0.2000 1.0000 0.3333 0.3333 3.0000 1.0000 鏈澶х殑鐗瑰緛鍊间负锛歵bmax = 3.0385 鏈澶х殑鐗瑰緛鍊煎搴旂殑鐗瑰緛鍚戦噺涓猴紙鏍囧噯鍖栧悗鐨勶級锛歵bvector = 0.6370 0.1047 0.2583 >> s A1 = 1.0000 5.0000 0.2000 1.0000 鏈澶х殑鐗瑰緛鍊间负...
绛旓細>> syms x1 x2 x3 >> A =[1,2,3;4,5,6;3,2,1];B=[x1;4;2];C=[4;x2;x3];[x1,x2,x3]=solve(A*B-C,'x1','x2','x3')x1 = -10 x2 = -8 x3 = -20
绛旓細1 a=[1,2, 0;3,0,5;7 ,8 ,9]锛沘(2,2)=4;b=a;c=a(2:3,:);2 a=[1, 2, 3;4, 5, 6;7, 8, 9];b=[2 4 6;1 3 5;7 9 10];x=a.*2 y=a.*b z=a*b x = 2 4 6 8 10 12 14 16 18 y = 2 8 18 4 ...
绛旓細褰撶劧鎴戜滑涔熷彲浠ョ敤ezplot浜х敓鍥惧儚 浣嗘敞鎰廵zplot鎵杈撳叆鐨勫弬鏁版槸绗﹀彿锛屼笉鏄鐭╅樀鎴栧悜閲忋俿yms x;杩欓噷瀹氫箟浜嗙鍙峰彉閲弜 ezplot(x^2)鍥犱负x鏄崟涓殑绗﹀彿锛屼笉鏄悜閲忔垨鐭╅樀锛屼篃灏辨槸璇翠粬灏辨槸涓涓崟鐙殑鍏冪礌銆傛墍浠ヤ笉鐢ㄥ姞鐐癸紝鐩存帴骞虫柟灏卞彲浠ヤ簡銆傚綋鐒跺鏋滃啓鎴恱.^2涔熷彲浠ワ紝灏卞拰3^2鍜3.^2涓鏍凤紝鍥犱负鍗曚釜鍏冪礌瀵瑰簲...
绛旓細鍋鐭╅樀 (A,B)锛屽瀹冭繘琛屽垵绛夎鍙樻崲, 灏嗗乏杈瑰寲鎴愬崟浣嶇煩闃, 鍒欏彸杈瑰氨鏄疿锛屽嵆 (E, A^(-1)B)銆傜粰涓よ竟宸︿箻A鐨勯嗛樀锛屽緱鍒扮殑灏辨槸X銆傚彲浠鐢∕ATLAB寰堟柟渚跨殑绠楀嚭鏉ャ倄=(A-1)*B(-1鏄笂鏍) 娉ㄦ剰锛氫竴瀹氭槸宸︿箻銆傝浆鎹㈡垚 AX=B 鐨勫舰寮.XA=B 涓よ竟鍙栬浆缃緱 A^duTX^T = B^T 瀵(A^T,B^T)鐢...
绛旓細璇ュ嚱鏁颁骇鐢熶竴涓暟缁勫厓绱犲间粙浜0鍒版帴杩戜絾涓嶇瓑浜1涔嬮棿鐨勬诞鐐规暟銆傛渶鍚庯紝`rand`瀹氫箟浜鐭╅樀鐨勪笁缁村昂瀵革紝鍗崇煩闃垫湁3灞傦紝姣忓眰鏈4琛屽拰5鍒椼備綘鍙互鏍规嵁闇瑕佽皟鏁磋繖浜涙暟瀛楁潵鍒涘缓涓嶅悓灏哄鐨勪笁缁寸煩闃点傚垱寤哄畬鎴愬悗锛屼綘鍙互浣跨敤MATLAB鐨勫悇绉嶅嚱鏁板拰鍛戒护鏉ユ搷浣滆繖涓笁缁寸煩闃碉紝濡傛煡鐪嬪叾澶у皬銆佸鍏惰繘琛屽垏鐗囨搷浣滄垨杩涜鍏朵粬鏁板...
