矩陣運(yùn)算與Matlab命令

矩陣運(yùn)算與Matlab命令
1.1 知識(shí)要點(diǎn)與背景：日常矩陣及其運(yùn)算
【????? A=[4 2 3;1 3 2;1 3 3;3 2 2],???????? % 表1-1、表1-2的數(shù)據(jù)分別寫(xiě)成矩陣形式
? B=[35 20 60 45;10 15 50 40;20 12 45 20]???????? 】
【?????? C=A*B????????? %矩陣乘法，求各訂單所對(duì)應(yīng)的原材料和勞動(dòng)力 。???? 】??
【?? whos?????????? % 查看Matlab工作空間中變量及其規(guī)模? 】????
1.2實(shí)驗(yàn)與觀(guān)察：矩陣和Matlab語(yǔ)言
1.2.1 向量的生成和運(yùn)算
? 【????? x=linspace(0,4*pi,100);?? %將[0,4π]區(qū)間100等分，產(chǎn)生了一個(gè)100維向量
?? y=sin(x);?????????????? %計(jì)算函數(shù)值，產(chǎn)生了一個(gè)與x同維的100維函數(shù)向量y
?? y1=sin(x).^2;???????? %計(jì)算函數(shù)向量，注意元素群運(yùn)算
??? y2=exp(-x).*sin(x);???????????????
?? %以x為橫坐標(biāo)，y為縱坐標(biāo)畫(huà)函數(shù)的圖用不同的線(xiàn)型將函數(shù)曲線(xiàn)繪制在一個(gè)圖上
??? plot(x,y,'-',x,y1,'-',x,y2,'.-')??????????? 】
1. 向量的創(chuàng)建
???????????? ◆直接輸入向量。
【x1=[1 2 4],x2=[1,2,1],x3=x1' 】
????? ◆冒號(hào)創(chuàng)建向量 。
【???? x1=3.4:6.7
?x2=3.4:2:6.7
?x3=2.6:-0.8:0?? 】
?????? ◆生成線(xiàn)性等分向量。
【???? x=linspace(0,1,5)?? 】
2. 向量的運(yùn)算
?【?? y=sin(x)??? 】
???????? 【? y1=sin(x).^2;?? y2=exp(-x).*sin(x);?? 】????
1.2.2.矩陣創(chuàng)建和運(yùn)算
1.創(chuàng)建矩陣
（1）數(shù)值矩陣的創(chuàng)建
??????? ◆直接輸入法創(chuàng)建簡(jiǎn)單矩陣。
【?? A=[1 2 3 4; 5 6 7 8; 9 10 11 12] 】
【?? B=[-1.3,sqrt(3);(1+2)*4/5,sin(5);exp(2),6] 】
（2）符號(hào)矩陣的創(chuàng)建
?????? ◆
【?? syms a11 a12 a13 a14 a21 a22 a23 a24 a31 a32 a33 a34 …
??????? b11 b12 b13 b14 b21 b22 b23 b24 b31 b32 b33 b34
A1=[a11 a12 a13 a14 ;a21 a22 a23 a24; a31 a32 a33 a34],
B1=[b11 b12 b13 b14 ;b21 b22 b23 b24; b31 b32 b33 b34]????????? 】?
2.矩陣的運(yùn)算
????????
【?? C=A1+B1,D=A1-B1? 】
???????
【???? syms c
? cA=c*A1? 】
???????
【?? C=A1*B1?? 】
{??? ??? Error using ==> sym/mtimes,? Inner matrix dimensions must agree.??? }
???????
【???? A2=A1(:,1:3), B1? 】
【? G=A2*B1??????? 】
【?? g11=A2(1,:)*B1(:,1)? 】
??????? 【??? A,?? A_trans=A'?? 】
???????
【?? H=[1 2 3 ; 2 1 0 ; 1 2 3 ], K=[1 2 3 ; 2 1 0 ; 2 3 1]
h_det=det(H), k_det=det(K),H_inv=inv(H),K_inv=K^-1?????????????? 】
????
【????? A=[3 0 1; 1 1 0;0 1 4];
?? B=inv(A-2*eye(3))*A, B=(A-2*eye(3))\A??? 】
3.分塊矩陣：矩陣的裁剪、分割、修改與抽取
（1）
【??? A=[1 0 1 1 2;0 1 -1 2 3;3 0 5 1 0;2 3 1 2 1], vr=[1,3];vc=[1,3];
A1=A(vr,vc)?????? %取出A的1、3行和1、3列的交叉處元素構(gòu)成新矩陣A1 】
???? ◆將上面的矩陣A分為四塊，并把它們賦值到矩陣B中，觀(guān)察運(yùn)行后的結(jié)果。
【?? A11=A(1:2,1:2),A12=A(1:2,3:5),A21=A(3:4,1:2),A22=A(3:4,3:5)
B=[A11 A12;A21 A22]?????????? 】
（2）矩陣的修改和提取
??????? ◆? 【?? A=[1 0 1 1 2;0 1 -1 2 3;3 0 5 1 0;2 3 1 2 1]
? A(1,:)=[0 0 0 0 0];?????? A????????????????? 】
?????? ◆ 觀(guān)察：
【??? B(:,[2,4])=[ ]???? %刪除矩陣B的第2、4列?? 】
（3）矩陣元素的抽取
4.生成特殊矩陣
?????? 。
??????? ◆
【?? y1=rand(1,5), y2=rand(1,5),
rand('seed',3), x1=rand(1,5),? rand('seed',3), x2=rand(1,5)? 】
5. 常用矩陣函數(shù)
6. 數(shù)據(jù)的簡(jiǎn)單分析
??????? ◆
【????? rand('seed',1);A=rand(3,6),
?Asort=sort(A), Amax=max(A), Asum=sum(A)? 】
1.2.3? Matlab工作環(huán)境和編程
2.Matlab的基本設(shè)計(jì)
1.3應(yīng)用、思考與練習(xí)
1.3.1? 關(guān)系矩陣
1.3.2 投入產(chǎn)出
1.3.3 循環(huán)比賽的名次
? 【??? A=[0 1 1 0; 0 0 1 1; 0 0 0 1; 1 0 0 0],
???? e=ones(4,1); c=A*e;??
????? s=c'????????????????????????????? 】

★ 畫(huà)矩陣結(jié)構(gòu)圖的gplot指令。
?????? ◆（3）????????
【????? clf, A=[0 1 1 0;0 0 1 1;0 0 0 1;1 0 0 0];? xy=[0 1;0 0;-1 –0.5;1 –0.5];
?? graphy_plot(A,xy,1,0.5),? % gplot(A,xy)?? 】
1.3.4? 參考程序
graphy_plot.m
?【? function y=graphy_plot(A,xy,l,p)
?????? %畫(huà)矩陣的有向結(jié)構(gòu)圖。 A為鄰接矩陣，xy為頂點(diǎn)坐標(biāo)，l控制參數(shù)，l=0，畫(huà)無(wú)向圖；????
?????? %l~=0,畫(huà)有向圖。p為控制箭頭大小的參數(shù)。
a=-max(abs(xy(:,1)))*1.1;b=max(abs(xy(:,1)))*1.1;
c=-max(abs(xy(:,2)))*1.1;d=max(abs(xy(:,2)))*1.1;
if l=0
?? gplot(A,xy),axis([a b c d]),hold on,
elseif l~=0
?? U=[];V=[];X=[];Y=[];
?? n=length(A(:,1)) ;
?? for i=1:n
????? k=find(A(i,:)~=0);
????? m=length(k);
????? if(m~=0)
???????? for j=1:m
??????????? u(1)=(xy(k(j),1)-xy(i,1)); v(1)=(xy(k(j),2)-xy(i,2));
??????????? u(2)=eps;???? v(2)=eps;???? U=[u;U];V=[v;V];
??????????? X=[[xy(i,1) xy(k(j),1)];X]; Y=[[xy(i,2) xy(k(j),2)];Y];
????????? end
????? text(xy(i,1),xy(i,2),['\bullet\leftarrow\fontsize{16}\it{V}',…
????? um2str(i)]);??? hold on,
?????? end
????????? end
?? gplot(A,xy),axis([a b c d]),hold on,
?? h=quiver(X,Y,U,V,p);set(h,'color','red');hold on,
?? plot(xy(:,1),xy(:,2),'k.','markersize',12),hold on,
end? , hold off?????????????? 】

200多個(gè)MATLAB經(jīng)典教程和MATLAB論文請(qǐng)查看：matlab教程