第6章MATLAB数值计算.ppt

第6章 MATLAB数值计算 6.1.1 数据统计与分析 1. 求矩阵最大元素和最小元素 （1）求向量的最大值和最小值 y=max(X)：返回向量X的最大值存入y，如果X中包含复数元素，则按模取最大值。 [y,I]=max(X)：返回向量X的最大值存入y，最大值的序号存入I，如果X中包含复数元素，则按模取最大值。 求向量X的最小值的函数是min(X)，用法和max(X)完全相同。 例 求向量x的最大值。 命令如下： x=[-43,72,9,16,23,47]; y=max(x) %求向量x中的最大值 [y,l]=max(x) %求向量x中的最大值及其该元素的位置 （2）求矩阵的最大值和最小值 max(A)：返回一个行向量，向量的第i个元素是矩阵A的第i列上的最大值。 [Y,U]=max(A)：返回行向量Y和U，Y向量记录A的每列的最大值，U向量记录每列最大值的行号。 （3）两个向量或矩阵对应元素的比较 函数max和min还能对两个同型的向量或矩阵进行比较，调用格式为： U=max(A,B)：A,B是两个同型的向量或矩阵，结果U是与A,B同型的向量或矩阵，U的每个元素等于A,B对应元素的较大者。 U=max(A,n)：n是一个标量，结果U是与A同型的向量或矩阵，U的每个元素等于A对应元素和n中的较大者。 例 求两个2×3矩阵x, y所有同一位置上的较大元素构成的新矩阵p。 2. 求矩阵的平均值和中值 mean(X) 返回向量X的算术平均值 median(X) 返回向量X的中值 mean(A) 返回一个行向量，其第i个元素是A的第i列的算术平 均值 median(A)返回一个行向量，其第i个元素是A的第i列的中值 mean(A,dim) 当dim=1时等同于mean(A); 当dim=2时返回一个列向量，其第i个元素是A的第i行的算术平均值 median(A,dim)当dim=1时等同于median(A); 当dim=2时返回一个列向量，其第i个元素是A的第i行的中值 例6.11 用一个3次多项式在区间[0，2π]内逼近函数。 命令如下： X=linspace(0,2*pi,50); Y=sin(X); P=polyfit(X,Y,3) %得到3次多项式的系数和误差 Y1 = polyval(P,X); plot(X,Y,:o,X,Y1,-*); legend(原始曲线,拟合曲线); 1. 多项式的四则运算 （1）多项式的加减运算 （2）多项式乘法运算 因为多项式是用其系数构成的离散序列表示的，因此多项式的乘法和除法就对应于卷积和去卷积的操作。 2. 多项式的导函数 3. 多项式求值 4. 多项式求根 6.2 数值微积分 6.2.1 数值微分 1. 数值差分与差商 2. 数值微分的实现 在MATLAB中，没有直接提供求数值导数的函数，只有计算向前差分的函数diff，其调用格式为： DX=diff(X)：计算向量X的向前差分 DX=diff(X,n)：计算X的n阶向前差分。 6.2 数值微积分 (1)被积函数是一个解析式 6.4.1 直接解法 1．利用左除运算符的直接解法 对于线性方程组Ax=b，可以利用左除运算符“\”求解： x=A\b 6.4.3 求线性方程组的通解 线性方程组的求解分为两类：一类是求方程组的惟一解即特解，另一类是求方程组的无穷解即通解。这里对线性方程组 Ax=b的求解理论作一个归纳：p172 常微分方程： clear;clc; p = [1 2 1]; x=-1:0.01:1; dp = polyder(p); dpx = polyval(dp,x); y = x.^2 + 2*x +1; dx = diff(y)/0.01; x1 = x(1:end-1); gx = 2*x+2; plot(x,dpx,r,x1,dx,g,x,gx,b); figure(2) subplot(1,3,1); plot(x,dpx,r); subplot(1,3,2); plot(x1,dx,g); subplot(1,3,3); plot(x,gx,b); 例：f(x) = x2+2x+1 　　　用不同方法求其倒数，并画出图象． 例：求以下的定积分。其精确值为0.7468204… %（1）使用字符串表示被处理函数 fun=exp(-x.*x); %（2）调用积分指令求积分值 E1 = quad (fun,0,1) E2 = quadl (fun,0,1) 区别：求数值积分是采取不