用matlab解矩阵一般方法 用matlab解矩阵
\u7528matlab\u89e3\u77e9\u9635\u65b9\u7a0b1\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
d =
0.8101 0.8484 0.8484 0.6870
0.4940 0.1044 + 0.3755i 0.1044 - 0.3755i -0.6611
0.3042 -0.3138 + 0.1397i -0.3138 - 0.1397i 0.2986
0.0841 -0.0182 - 0.1004i -0.0182 + 0.1004i -0.0416
v =
4.0724 0 0 0
0 -0.0137 + 0.5507i 0 0
0 0 -0.0137 - 0.5507i 0
0 0 0 -0.0450
矩阵分析是解决很多问题的好方法,但是很多时候矩阵的运算比较繁琐,特别是高阶矩阵运算。这时候如果用matlab来计算就方便快捷得多。下面我将介绍一些基本的矩阵运算方法。如加,减,乘,除,转置,求逆。
约定:
a=[1,3,5;2,4,6;7,9,8] b=[9,6,4;3,4,5;2,3,4]
工具/原料
matlab
加和减:
加减法的命令很简单,直接用加或者减号就可以了。如:
c=a+b
d=a-b
乘法:
一般乘法:c=a*b,要求a的列数等于b的行数。
如果a,b是一般的向量,如a=[1,2,3] b=[3,4,5]
点积: dot(a,b),
叉积: cross(a,b)
卷积: conv(a,b)
除法:一般在解线性方程组时会用到。
x=a\b 如果ax=b,则 x=a\b是矩阵方程的解。
x=b/a 如果xa=b, 则x=b/a是矩阵方程的解。
转置:
转置时,矩阵的第一行变成第一列,第二行变成第二列,。。。
x=a.'
求逆:
要求矩阵为方阵。这在矩阵运算中很常用。
x=inv(a)
方法/步骤
绛旓細闄ゆ硶锛氫竴鑸湪瑙g嚎鎬ф柟绋嬬粍鏃朵細鐢ㄥ埌銆倄=a\b 濡傛灉ax=b锛屽垯 x=a\b鏄煩闃垫柟绋嬬殑瑙c倄=b/a 濡傛灉xa=b, 鍒檟=b/a鏄煩闃垫柟绋嬬殑瑙c傝浆缃細杞疆鏃讹紝鐭╅樀鐨勭涓琛屽彉鎴愮涓鍒楋紝绗簩琛屽彉鎴愮浜屽垪锛屻傘傘倄=a.'姹傞嗭細瑕佹眰鐭╅樀涓烘柟闃点傝繖鍦ㄧ煩闃佃繍绠椾腑寰堝父鐢ㄣ倄=inv(a)
绛旓細4銆佹暟鎹浆鎹銆傝嫢鈥滅洊0鈥濈嚎鐨勬暟鐩瓑浜庣煩闃电殑缁存暟鍒欒烦杩囨姝ワ紝鑻モ滅洊0鈥濈嚎鐨勬暟鐩皬浜庣煩闃电殑缁存暟鍒欒繘琛屾暟鎹浆鎹傛湰棰樺睘浜庡悗鑰咃紝鍒欑洿鎺ユ眰鏈浼樿В銆傚n缁寸煩闃碉紝鎵惧嚭涓嶅悓琛屻佷笉鍚屽垪鐨刵涓0锛屽姣忎釜0鐨勪綅缃唬琛ㄤ竴瀵归厤缃叧绯伙紝鍏蜂綋姝ラ濡備笅銆傦紙1锛夊厛鎵惧彧鍚湁涓涓0鐨勮锛堟垨鍒楋級锛屽皢璇ヨ锛堟垨鍒楋級涓殑0...
绛旓細1銆侀鍏堬紝鎵撳紑MATLAB鍛戒护琛岀獥鍙o紝杈撳叆a = [1 2 3; 4 5 6; 7 8 9; 6 4 7]鍒涘缓涓涓猘鐭╅樀锛屽涓嬪浘鎵绀猴紝鐒跺悗杩涘叆涓嬩竴姝ャ2銆佸叾娆★紝瀹屾垚涓婅堪姝ラ鍚庯紝杈撳叆a锛1锛2,1锛3锛変互鎻愬彇鐭╅樀鐨勭1,2琛岀1鍒3鍒楋紝濡備笅鍥炬墍绀猴紝鐒跺悗杩涘叆涓嬩竴姝ャ3銆佹帴鐫锛屽畬鎴愪笂杩版楠ゅ悗锛屾寜Enter閿紝鍙互鐪嬪埌鐭...
绛旓細瑙e喅杩欎釜闂鐨鏂规硶濡備笅锛1銆佺涓姝ユ墦寮matlab鍛戒护琛岀獥鍙o紝杈撳叆a锛濓蓟123锛456锛789锛647锛斤紝鍒涘缓涓涓猘鐭╅樀銆2銆佺浜屾杈撳叆a锛1:2锛1:3锛夛紝杩涜鎻愬彇a鐭╅樀鐨勭1锛2琛岀1鍒3鍒椼3銆佺涓夋鎸夊洖杞﹂敭涔嬪悗锛屽彲浠ョ湅鍒板凡缁忔彁鍙栦簡a鐭╅樀鐨勭1锛2琛岀1鍒3鍒椼4銆佺鍥涙杈撳叆a锛2:3锛1:2锛夛紝杩涜鎻...
绛旓細1銆侀鍏堟墦寮鐢佃剳涓婄殑鈥matlab鈥濊蒋浠讹紝鍦鍛戒护琛岀獥鍙h緭鍏ヤ竴涓悜閲弉锛屽悜閲弉鏈4涓厓绱犮傛帴鐫浣跨敤diag鍑芥暟鐢熸垚瀵硅鐭╅樀銆2銆乨iag(n,k)鍙互鎶婂悜閲忔斁鍦ㄧk鏉″瑙掔嚎涓婏紝k涓烘鍊硷紝琛ㄧず鍙充笂銆俴涓鸿礋鍊硷紝鍒欒〃绀哄乏涓嬨傚湪鍛戒护琛岃緭鍏iag(n,2)锛屽嵆鍙妸鍚戦噺鏀惧湪鍙充笂鐨勭浜屾潯瀵硅绾夸笂銆3銆佸湪鍛戒护琛岃緭鍏iag(n,-...
绛旓細涓銆鐭╅樀鐨勮〃绀鏂规硶 1銆佺煩闃靛厓绱犲繀椤鍦鈥漑]鈥濆唴锛2銆佺煩闃电殑鍚岃鍏冪礌涔嬮棿鐢ㄧ┖鏍硷紙鎴栤,鈥濓級闅斿紑锛3銆佺煩闃电殑琛屼笌琛屼箣闂寸敤鈥;鈥濓紙鎴栧洖杞︾锛夐殧寮锛4銆佺煩闃电殑鍏冪礌鍙互鏄暟鍊笺佸彉閲忋佽〃杈惧紡鎴栧嚱鏁帮紱5銆佺煩闃电殑灏哄涓嶅繀棰勫厛瀹氫箟銆備簩銆佺煩闃电殑杩愮畻 1銆佺畻鏈繍绠 MATLAB鐨鍩烘湰绠楁湳杩愮畻鏈夛細锛(鍔)銆侊紞(鍑)銆...
绛旓細>> 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銆佹墦寮matlab杞欢锛2銆佽緭鍏ヤ竴涓嬩唬鐮侊紝娓呭睆鍜屽伐浣滅┖闂达紱3銆佽緭鍏ヨ嚜宸辩殑鐭╅樀锛屽鏈緥涓缃竴涓煩闃礎=[1 2 ;3 4 ];4銆佽緭鍏姹傝В鐭╅樀閫嗙殑绋嬪簭锛屽鍥炬墍绀猴紱B=inv(A)5銆侀獙璇佽嚜宸辨眰瑙g殑閫嗭紝涓や釜鐭╅樀鐨勪箻绉负鍗曚綅闃碉紝姹傞嗘纭紱
绛旓細鐪嬭捣鏉瑙f硶鏄鐨勩備綘鍙互鐢ㄥ緱鍒扮殑X杩涜楠岃瘉锛氱湅鐪婣*X*B=C鏄惁鎴愮珛銆傚鏋滈獙璇佹纭紝閭h偗瀹氬氨涓嶆槸鐭╅樀璁$畻鐨勯棶棰樹簡銆
绛旓細瑙o細E(Y)=1*(0.12+0.03+0.15)+3*(0.05+0.25+0.20)+5*(0.15+0.02+0.03)锛汦(X)=1*(0.12+0.05+0.15)+2*(0.03+0.25+0.02)+3*(0.15+0.20+0.03)锛汦(XY)=1*1*0.12+1*2*0.03+1*3*0.15 +3*1*0.05+3*2*0.25+3*3*0.20 +5*1*0.15+5*2*...