實驗1 矩陣的建立和基本運算

實驗 矩陣的建立和基本運FileExit，MATLABExit52->>b=2-2-1.1輸入“>>a=5>>b=[1,3,5,7,9,11]%1.222，c=(1.3,2.5,7.6,2,-3)，d=(23,20,17,14,11,8,5,2)，>>A=[2,3,5;1,3,5;6,9,4]%235135694問題：A(2,3)A(7)又如何？體會以上輸入的結(jié)果，注意，>>d=numel(A)%A的元素數(shù)，5.x>>[n,m]=size(A)%A的行(n)、列(m)4>>[i,j]=find(A>3)；%A3%”后面是注釋語句，被忽略而不執(zhí)行；對一個數(shù)組可用n=length(A),A若是矩陣，nA的行、列數(shù)的最大值。>>A(2,：)；%A2>>A（[1，3]，：；%>>A（2：3，1：2）%A2，31，2 >>A([1，3]，：)=A（[3，1]，：； %A13問題1.4：A2，3>>A(2,：)=4；%A24>>A（find(A==3)）=–3；%A3的所有元素換為>>A(2，：)=[%A2356359379>>reshape(A,2,6)%A2×6ans323355576899>>A(4,5)=3;%A的維數(shù)，A4ⅹ5 ans=2343366553251232436365352505 Ak>>tril(A,k);Ak>>triu(A,k);Ak注意:“:”表示“全部”>>flipud(A);%A>>fliplr(A);%A>>rot90(A);%A問題1.5:rot90(A,2)rot90(A,-2)>>A= n>>A=ones(n,m);nⅹm1>>A=zeros(n,m);nⅹm0>>A=rand(n,m);nⅹm0～１問題1.6：產(chǎn)生一個在區(qū)間[10,20]4>>randn(m,n);mⅹn>> 1～n>>randperm(6)Ans= >>logspace(a,b,n);%在(10a,10bn>>diag(a,b,n);a,b,c,d,… nHilbertH(i,j)=1/(i+j->>magic(n);n>>4/2;%42, %42, %43>>sqrt(4);%4>>exp(3);%e3次方,>>log(4);%4的自然對數(shù),log10(4)10為底,log2(4)21.1 %A>>det(A);%A的行列式,A %A %A %A>>[X,D]=eig(A);%AX %A的跡,A %>>A+B;%A,B必須是同維矩陣,3+A>>A- %A,B必須是同維矩陣,3-A %A.*B %(A./B進行比較>>A\B;%(A.\B進行比較 %A^2A*A(A.^2進行比較注意:矩陣的加減乘除按相關(guān)規(guī)則運算,否則給出警告信息;“.*”,“./”,“.\”,“.^”稱為點運算(或稱數(shù)組運算,又稱元素群運算),點運算是前后矩陣1.7:A的本征矢量和本征值,2^A(A必須是方陣)2.^A的區(qū)別.1.3.>>savedate %a,b,cdate.mat>>load %date.mat %將列出工作空間所有變量的變量名﹑大小﹑字節(jié)數(shù)﹑>>helpsqrtsqrt三﹑MATLAB2～32～3AA21A1，3，4A13A3×414×24四、1.2>>A1=[1,2,3;4,5,6;7,8,9];A2=A1‵;A3=A1-123456789︰,︰,2 7123456789︰,︰,2 7258369︰,︰,3 --20-420說明:AB的張量積定義為 M a aB mn其中ABBA均為mp×nq矩陣,但一般ABBA 【例1.3】A

B 6

AB>>A=[12;34];B=[123;456;78>>C=kron(A,B)C=123246456878936948 %求A的2-范數(shù),和norm(A)相同. %求行范數(shù)（無窮大范數(shù)A最大列之和 %AFrobenius范數(shù)

22 ALULUL和一個U的乘積,A=LU.格式:[L,U]=lu(X)%U為上三