06- 值计算功能.pptVIP

第6章 数值计算功能 6.3 微分和积分 x=0:0.1:10; y=sin(x); z=cumsum(y)*0.1; plot(x,y,r-,x,z,k*) 6.4 MATLAB 7与线性代数 矩阵的特征值函数及其功能 矩阵的分解运算 三角分解(lu) 正交分解(qr) 特征值分解(eig) Chollesky分解(chol) 奇异值分解(svd) 矩阵的翻转 A=randn(3) A = -0.4326 0.2877 1.1892 -1.6656 -1.1465 -0.0376 0.1253 1.1909 0.3273 B=fliplr(A) B = 1.1892 0.2877 -0.4326 -0.0376 -1.1465 -1.6656 0. 7. MATLAB 7 与概 率 统 计 均匀分布的随机数据的产生 （unifrnd） R=unifrnd(A,B)返回区间[A,B]上的连续型均匀分布 R=unifrnd(A,B,M,N)返回一个M*N的矩阵 * MATLAB 7 简明教程 清华大学出版社 数值微分 有限差分 1．求元素之差函数diff( ) 格式一：A=diff(x) 功能：计算x中相邻元素之间的差值或近似导数。如果x为向量，则返回一个比x少一个元素的向量，其元素值为[x(2)-x(1)，x(3)-x(2)，…，x(n)-x(n-1)]；如果x为矩阵，则返回一个列间差值的矩阵：[x(2:n,:) - x(1:n-1,:)]。 格式二：A=diff(x, n) 功能：使用diff函数递归 n 次，计算第n阶差值。例如，diff(x,2) = diff(diff(x))。 例6-21 使用diff函数求数值微分。 解：在命令窗口中输入如下命令，并按Enter键确认。 h = .002; x = 0:h:pi; diff(sin(x.^2))/h; %is an approximation to 2*cos(x.^2).*x X=[6,9,3,4,0;5,4,1,2,5;6,7,7,8,0;7,8,9,10,0] X = 6 9 3 4 0 5 4 1 2 5 6 7 7 8 0 7 8 9 10 0 A=diff(X) A = -1 -5 -2 -2 5 1 3 6 6 -5 1 1 2 2 0 B=diff(X,2) B = 2 8 8 8 -10 0 -2 -4 -4 5 2．求数值梯度函数gradient( ) 两变量函数F(x, y)的梯度定义为 对N个变量函数 F(x,y,z, … ) 其梯度为 梯度可看作指向F增加方向的向量集。 格式一：Fx=gradient(F) 功能：F为一向量，返回F的一维数值梯度，Fx与 一致，表示x方向的差分。 格式二：[Fx, Fy]=gradient(F) 功能：F为一矩阵，返回二维数值梯度的x和y分量。 Fx与 表示x(列)方向的差分 Fy与 一致，表示y(行)方向的差分。每个方向点间距离设为1。 函数的数值积分 数值积分实际是计算函数曲线下的面积，即定积分 基本思想是将整个积分区域划分为若干子区间。 矩形求积函数cumsum cumsum (x)对向量x，返回一个向量，该向量的第N个元素是x 的前N个元素的和 x1=[1 2 3 4 5 6 7 8] x1 = 1 2 3 4 5 6 7 cumsum(x1) ans = 1 3 6 10 15 21 28 36 x2=[1 2 3 ;4 5 6;7 8 9] x2 = 1 2 3 4 5 6 7 8 9 cumsum(x2) ans = 1 2 3 5