第4章-Matlab的数值计算.ppt

6.条件数cond(A，1)计算A的1-范数下的条件数cond(A)或cond(A，2)计算A的2-范数下的条件数cond(A，inf)计算A的-范数下的条件数*7.行阶梯形R=rref(A)求矩阵A的减少行阶梯形式8.有理基z=nullz的列向量为方程组的正交规范基z=null(A,’r’)z的列向量是方程AX=0的有理基*9.子空间夹角subspace(A,B)返回两个列指定的A,B子空间的夹角例4-36求魔方矩阵的子空间夹角。H=magic(8);A=H(:,2:4);B=H(:,5:8);theta=subspace(A,B)*4.3.3矩阵的特征值分析eig(A)求包含矩阵A的特征值的向量。[X,D]=eig(A)产生一个矩阵A的特征值在对角线上的对角矩阵D和矩阵X，它们的列是相应的特征向量，满足AX=XD。*4.3.3矩阵的分解1.LU分解MATLAB的语法为：[L,U]=lu(A)%A=LU求上三角矩阵U和交换下三角矩阵L，使LU=A。[L,U,P]=lu(A)%PA=LU求上三角矩阵U、有单位对角线的下三角矩阵L和交换矩阵P，满足LU=PA。*高斯消去法（LU矩阵分解法）将任意一个方阵分解成一个上三角矩阵与一个下三角矩阵的积。 A=LU其中L为主对角线元素都是1的下三角矩阵（或其行变换形式），U为一个上三角矩阵。*2.Cholesky因式分解chol(A)求正定矩阵A的Cholesky因子，是一个上三角矩阵。如果A不是一个正定矩阵，则给出一个错误信息。[G,err]=chol(A)求矩阵A的Cholesky因子G。如果A不是一个正定矩阵，则不给出错误信息，而是将err设为非零值。*3.正交分解[Q,R]=qr(A)求得m×m的矩阵Q和上三角矩阵R，Q的列形成了一个正交基，Q和R满足A=Q×R。[Q,R,P]=qr(A)求矩阵Q、上三角矩阵R和交换矩阵P。Q的列形成一个正交基，R的对角线元素按大小降序排列，它们满足A×P=Q×R。*4.SVD分解计算奇异值分解的函数svd的格式为：svd(A)返回一个包含矩阵A的奇异值的向量。[U,S,V]=svd(A)返回奇异值为对角线的矩阵S和A同型的矩阵，m×m和n×n的酉矩阵U和V，使A=U×S×V′和U′×A×V=S成立。S中奇异值是非负数且按降序排列。*第4章Matlab的数值计算(1)——Matlab具有出色的数值计算能力，占据世界上数值计算软件的主导地位。——Matlab数值计算是使用Matlab的基础，是Matlab强大计算功能的体现。*主要内容：①数据分析；②矩阵分析；③多项式运算；④数值插值和拟合；⑤数值积分和微分。*4.1概述数值计算的特点：1）数值计算是定义在数值数组的基础上，在计算之前必须定义变量并赋值。2）数值计算存在误差。3）数值计算的计算速度较快。4）数值计算除了算术运算、关系运算和逻辑运算三种，还提供了许多运算函数。*4.2数据分析（1）了解数据分析函数：max,min,mean,sum,prod等（2）协方差和相关系数函数：cov(x),cov(A),corrcoef(x)等（3）频谱分析傅立叶变换函数：y=fft(x)，fftn，ifft，ifftn等*例4-10计算x的协方差、y的协方差、

x与y的互协方差。x=[12345];y=[246810];cx=cov(x)cy=cov(y)cxy=cov(x,y)**例4-11求行向量（0，0.1，0.2…，3.14）的离散傅立叶变换，并通过傅立叶反变换验证结果是正确的。4.3矩阵操作4.3.1矩阵的结构变换1.转置对于实矩阵用（’）符号或（.’）求转置结果是一样的；然而对于含复数的矩阵，则（’）将同时对复数进行共轭处理，而（.’）则只是将其排列形式进行转置。*例4-12分析下列语句的执行结果。a=[123];b=a′,c=a.′%行向量的转置 a=[1;2;3];b=a′%列向量的转置a=[