第5节数值计算功能.ppt

第5节 数值计算功能 教学重点 矩阵的生成和基本的数值运算 稀疏型矩阵的生成和基本操作 数值微分的求解方法 教学过程 矩阵的生成和基本的数值运算 特殊矩阵的生成 稀疏型矩阵的生成和基本操作 稀疏型矩阵和满矩阵的相互转换 微分和积分的定义 数值微分的求解方法 对函数进行数值积分的方法 使用MATLAB 7进行线性代数的运算 使用MATLAB 7进行概率统计方面的运算 1. 矩阵的生成 矩阵的生成有多种方式，通常使用的有4种方法： 在命令窗口中直接输入矩阵 通过语句和函数产生矩阵 在M文件中建立矩阵 从外部的数据文件中导入矩阵 matrix=[1 ,1, 1, 1;2, 2, 2, 2;3, 3, 3, 3;4, 4, 4, 4] matrix = 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 2. 矩阵的基本数值运算 矩阵与常数的四则运算 矩阵之间的四则运算 A=[2 1 -1;2 1 0;1 -1 1] A = 2 1 -1 2 1 0 1 -1 1 B=[1 -1 3;4 3 2]； X=B/A X = -2.0000 2.0000 1.0000 -2.6667 5.0000 -0.6667 3. 特殊矩阵的生成 零矩阵和全1矩阵的生成 （zeros 、ones） 对角矩阵的生成 （diag） 随机矩阵的生成 （rand、randn ） 范德蒙德矩阵的生成 （vander） 魔术矩阵的生成 （magic） Hilbert矩阵和反Hilbert矩阵的生成 （hilb、invhilb ） 特殊矩阵的生成举例 rand(5) ans = 0.9501 0.7621 0.6154 0.4057 0.0579 0.2311 0.4565 0.7919 0.9355 0.3529 0.6068 0.0185 0.9218 0.9169 0.8132 0.4860 0.8214 0.7382 0.4103 0.0099 0.8913 0.4447 0.1763 0.8936 0.1389 4. 稀疏型矩阵 稀疏矩阵的生成 稀疏矩阵与满矩阵的相互转换 稀疏矩阵的操作 （1） 稀疏矩阵的生成 在MATLAB 7中，生成稀疏矩阵用特殊的函数来进行，这些函数有speye、spones、spdiags、sparse、find、full、spalloc、sprand和sprandn 等。 Speye函数应用举例 A=eye(5) A = 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 speye(size(A)) ans = (1,1) 1 (2,2) 1 (3,3) 1 (4,4) 1 (5,5) 1 （2）稀疏矩阵与满矩阵的相互转换 用来将稀疏矩阵和满矩阵相互转换的函数有sparse、full和find等3个函数 。 S(10,50)=82; S(32,14)=82; S(251,396)=25; I=find(S) ％生成S中非零元素的位置 I = 3295 12309 99396 （3）稀疏矩阵的操作 对稀疏矩阵进行操作，主要由nnz、nonzeros、nzmax、sponse、spalloc、isspase、spyfun和spy等函数来实现 。 v=[6 2 7 8 ]; S=diag(v,1) R=sparse(S); N=issparse(S) N = 0 Y=issparse(R) Y = 1 5. 微分和积分 数值微分 使用diff函数求数值微分 使用gradient函数求近似梯度 jacobian函数求多元函数的导数 函数的数值积分 矩形求积 trapz函数(梯