基于MATLAB的数值分析.pptxVIP

第三章线性代数3.1常用矩阵函数norm：矩阵或向量范数Forvectors...norm(V,P)=sum(abs(V).^P)^(1/P).norm(V)=norm(V,2).norm(V,inf)=max(abs(V)).norm(V,-inf)=min(abs(V)).例x=[12345];x=[120300405];x=[100002000300405];[norm(x,1),norm(x,2),norm(x,3),norm(x,inf)]

01:不同范数意义下的单位圆03描绘norm(x,1)=1;norm(x,2)=1;norm(x,inf)=102运行以下Matlab程序，文件名为：normpolt.m04的图形

Formatrices...norm(X)isthelargestsingularvalueofX,max(svd(X)).norm(X,2)isthesameasnorm(X).norm(X,1)isthe1-normofX,thelargestcolumnsum,=max(sum(abs(X))).norm(X,inf)istheinfinitynormofX,thelargestrowsum,=max(sum(abs(X))).norm(X,fro)istheFrobeniusnorm,sqrt(sum(diag(X*X))).norm(X,P)isavailableformatrixXonlyifPis1,2,inforfro.例x=[120300;400506;708009];[norm(x,1),norm(x,2),norm(x,inf),norm(x,fro)]

dot(x,y)向量的内积det(A)方阵的行列式；rank(A)矩阵的秩；trace(A)矩阵的迹；rref(A)初等变换化矩阵A为阶梯矩阵inv(A)矩阵的逆；即A-1pinv(A)矩阵的广义逆A+null(A)零空间的基阵roth(A)值空间的基阵orth(A)将A标准正交化cond(A,flag)矩阵的条件数,flag=2,1,inf,fro;

例：分别求x=[1378-2],y=[393-39]的长度与它们的夹角。x=[1378-2];y=[393-39];1xx=norm(x,2);2yy=norm(y,2);3theta=acos(dot(x,y)/(xx*yy));s=[xx,yy,theta]4例：给定一组线性无关的向量，将其标准正交化a=magic(5);b=orth(a)5

d=eig(A)：方阵的特征值；[V,D]=eig(A):A*V=V*D[V,J]=jordan(A):A*V=V*Jc=condeig(A)：向量c中包含矩阵A关于各特征值的条件数[V,D,c]=condeig(A)：例:A=[100;120;123],d=eig(A),[V,D]=eig(A),C=condeig(A),[V,D,C]=condeig(A),

例：观察7阶随机矩阵特征值的分布a=rands(7,7)%产生7阶随机矩阵e=eig(a)title(特征值的分布);plot(real(e),imag(e),o)xlabel(实轴);ylabel(虚轴);注：本例验证了如下定理：实方阵的特征值或为实数或呈共轭对出现。

例：观察正交矩阵的特征值分布a=rands(7,7);b=orth(a);%构造一个正交矩阵theta=0:0.01:2*pi;e=eig(b);plot(real(e),imag(e),r*,cos(theta),sin(theta));axisequaltitle(正交矩阵特征值的分布);xlabel(实轴);ylabel(虚轴);注：本例验证了正交矩阵的特征值分布在复平面的单位圆上。