第4章Matlab例题.docVIP

例 （1）zeros(3) （2）zeros(2,3) （3）A=[1 2 3;4 5 6]; %产生一个2×3阶矩阵A zeros(size(A)) %一个与矩阵A同样大小的零矩阵 例a=10; b=30; x=a+(b-a)*rand(4) y=0.6+sqrt(0.1)*randn(4) 例 M=100+magic(5) 例format rat %以有理形式输出 H=hilb(4) H=invhilb(4) format short %恢复默认输出格式 例 pascal(5) 例4.6 A=[1:5;2:6;3:7;4:8;5:9] D=diag(1:5); D*A %用D左乘A，对A的每行乘以一个指定常数 例4. p=[3,-7,0,5,2,-18]; A=compan(p); %构造伴随矩阵x1=eig(A) %求A的特征值 x2=roots(p) %直接多项式p的零点 A=[1 -1 1;5 -4 3;2 1 1]; B=inv(A) A*B B*A 例4.9 A=[1,-2,3;3,-1,5;2,1,5]; b=[1;2;3]; x=inv(A)*b 例4.10 解法1： A=[1,-1,1;5,-4,3;2,1,1]; b1=[2;-3;1]; b2=[3;4;-5]; x=A\b1 y=A\b2 解法2： A=[1,-1,1;5,-4,3;2,1,1]; b=[2,3;-3,4;1,-5]; xy=A\b xy = 例4.A=[1,-1,1;5,-4,3;2,1,1]; [L,U]=lu(A) LU=L*U 利用第2种格式对矩阵A进行LU分解[L,U,P]=lu(A) LU=L*U %这种分解其乘积不为A inv(P)*L*U %考虑矩阵P后其乘积等于A 例4. clear A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; b=[13,-9,6,0]; tic x2=A\b; %用左除运算求解 toc tic x1=inv(A)*b; %用求逆运算求解 toc tic [L,U]=lu(A); %LU分解 x3=U\(L\b); %用LU分解求解 toc x1=x1 x2=x2 x3=x3 例4. A=[2,1,1,4;1,2,-1,2;1,-1,3,3]; [Q,R]=qr(A) QR=Q*R [Q,R,E]=qr(A) Q*R/E %验证A=Q*R*inv(E) 例4. A=[2,1,1;1,2,-1;1,-1,3]; R=chol(A) R*R 利用第2种格式对矩阵A进行Cholesky分解： [R,p]=chol() 例4.15 X=[2,0,0,0,0;0,0,0,0,0;0,0,0,5,0;0,1,0,0,-1;0,0,0,0,-5]; A=sparse(X) 例4.16 A=[2,2,1;3,1,-1;4,3,3;5,3,8;6,6,12]; B=spconvert(A) 例17 B=[1,1,2,1,0;2,4,6,6,1;0,3,1,4,2]; %取A对角线上元素构成B d=[-1,0,1]; %产生带状位置向量 A=spdiags(B,d,5,5); %产生稀疏存储的系数矩阵 b=[0,3,2,1,5]; %方程右边参数向量 x=(inv(A)*b) %求解 2 3