数学实验与建模精品课件：第5章　MATLAB数值计算.pptVIP

第5章　MATLAB数值计算 在科学和工程应用中，往往要进行大量的数学计算。这些运算一般来说难以用手工精确和快捷地进行，而要借助计算机编制相应的程序做近似计算并不断更新和扩充。MATLAB的数值分析功能十分强大中，本章主要讲述MAYLAB在函数、插值和曲线似合分析、微积分和线性方程系统方面的应用。 5.1 特殊矩阵 5.1.1对角阵与三角阵 1. 矩阵的对角元素 (1)提取矩阵的对角线元素 设A为m×n矩阵，diag(A)函数用于提取矩阵A主对角线元素产生一个具有min(m,n)个元素的列向量。 diag(A)函数还有更进一步的形式diag(A,k)，其功能是提取第k条对角线的元素。 例5.1 先建立5×5矩阵A，然后将A的第1行元素乘以1，第2行乘以2，…，第5行乘以5。 2. 矩阵的三角阵 (1)下三角矩阵 求矩阵A的下三角阵的MATLAB函数是tril(A) tril(A)函数也有更进一步的一种形式tril(A,k)，其功能是求矩阵A的第k条对角线以下的元素。 (2)上三角矩阵 在MATLAB中，提取矩阵A的上三角矩阵的函数是triu(A)和triu(A,k)，其用法与提取下三角矩阵的函数tril(A)和tril(A,k)完全相同。 5.1.2 特殊矩阵的生成 1. 魔方矩阵 　魔方矩阵是n*n元素所构成的方阵，其每个元素由不同的1~n2的整数所组成，它的每行、每列以及对角线元素之和均相等，并等于n(1+n2)/2.函数格式为 　　　　　　　　magic(n) 3. 希尔伯特矩阵 Hilbert矩阵的每个元素的值，由行数i和列数j决定，等于1/(i+j-1)，生成希尔伯特矩阵的函数是hilb(n)。MATLAB中，有一个专门求希尔伯特矩阵的逆的函数invhilb(n)，其功能是求n阶的希尔伯特矩阵的逆矩阵。 5.2 矩阵分析 5.2.1 矩阵结构变换 1. 矩阵的转置 转置运算符是单撇号()。 2. 矩阵的旋转 矩阵的旋转利用函数rot90(A,k)，功能是将矩阵A旋转90o的k倍，当k为1时可省略。 3. 矩阵的左右翻转 对矩阵A实施左右翻转的函数是fliplr(A)。 4. 矩阵的上下翻转 对矩阵A实施上下翻转的函数是flipud(A)。 5.2.2 矩阵的逆与伪逆 1. 矩阵的逆 求一个矩阵的逆非常容易。求方阵A的逆可调用函数inv(A)。 例5.4 用求逆矩阵的方法解线性方程组。 命令如下： 一般情况下，用左除比求矩阵的逆的方法更有效，即x=A\b。 2. 矩阵的伪逆 MATLAB中，求一个矩阵伪逆的函数是pinv(A)。 例5.5 求A的伪逆，并将结果送B。 命令如下： A=[3,1,1,1;1,3,1,1;1,1,3,1]; B=pinv(A) 例5.6 求矩阵A的伪逆。 在MATLAB命令窗口，输入命令： A=[0,0,0;0,1,0;0,0,1]; pinv(A) 5.2.3 方阵的行列式 求方阵A所对应的行列式的值的函数是det(A)。 例5.7用克莱姆(Cramer)方法求解线性方程组。 程序如下： D=[2,2,-1,1;4,3,-1,2;8,5,-3,4;3,3,-2,2]; %定义系数矩阵 b=[4;6;12;6]; %定义常数项向量 D1=[b,D(:,2:4)]; %用方程组的右端向量置换D的第1列 D2=[D(:,1:1),b,D(:,3:4)]; %用方程组的右端向量置换D的第2列 D3=[D(:,1:2),b,D(:,4:4)]; %用方程组的右端向量置换D的第3列 D4=[D(:,1:3),b]; %用方程组的右端向量置换D的第4列 DD=det(D); x1=det(D1)/DD; x2=det(D2)/DD; x3=det(D3)/DD; x4=det(D4)/DD; [x1,x2,x3,x4] 5.2.4 矩阵的秩 MATLAB中，求矩阵秩的函数是rank(A)。 例如，求例5.7中方程组系数矩阵D的秩，命令是： 说明D是一个满秩矩阵。 5.2.5 向量和矩阵的范数(不讲) 1. 计算向量3种常用范数的函数 (1)norm(V)或norm(V,2) 计算向量V的2—范数=sum(abs(V).^2)^(1/2) (2)norm(V,1)=sum(abs(V)) 计算向量V的1—范数 (3)norm(V,inf) 计算向量