第2章MATLAB数值计算21.pptVIP

一些测试函数 B 多项式求值 求多项式f(x)在某点或某些点的函数值的函数是polyval(f,x)。 若x为一数值，则求多项式在该点的值；若x为向量或矩阵，则对向量或矩阵中的每个元素求其多项式的值。 例 已知一个多项式， (1)计算f(x)=0 的全部根。 (2)由方程f(x)=0根构造一个多项式g(x)，并与f(x)进行对比。 (3)计算f(5)、f(7.8)、f(9.6)、f(12.3)的值。 命令如下： P=[3,0,4,-5,-7.2,5]; %输入多项式系数；缺项补0; X=roots(P) %求方程f(x)=0的根 G=poly(X) %求多项式g(x) X0=[5,7.8,9.6,12.3]; f=polyval(P,X0) %求多项式f(x)在给定点的值 C 多项式的四则运算 ? 由于多项式是利用向量来表示，多项式的四则运算可以转化为向量的运算。 多项式的加减为对应项系数的加减，因此可以通过向量的加减来实现。但是在向量的加减中两个向量需要有相同的长度，因此在进行多项式加减时，需要将短的向量前面补0。 多项式的乘法实际上是多项式系数向量之间的卷积，用Conv函数。 多项式的除法为乘法的拟运算，可以通过反卷积函数deconv来实现。 命令如下： P1=[3,0,4,-5,-6,7]; P2=[0,0,3,-2,-3,0]; %输入多项式系数；缺项补0;补短; f1+f2=P1+P2 %多项式求和 f1+f2=P1+P2 %多项式求差 f1*f2=conv(P1,P2) %多项式乘 f1/f2=deconv(P1,P2) %多项式f1(x)除以f2(x) 例 已知两个多项式， 试求： D 多项式的导函数 对多项式求导数的函数是： p=polyder(P) 求多项式P的导函数 p=polyder(P,Q) 求P*Q的导函数 ? 补充：2.6 数据统计与分析 1 求矩阵最大和最小元素 求向量的最大最小元素 ①y=max(X) 返回向量X的最大元素存入y。 ②[y,I]=max(X) 返回向量X的最大元素存入y，最大元素的序号存入I。 求矩阵的最大和最小元素 ①max(A) 返回一个行向量，向量的第i个元素是A矩阵的第i列上的最大元素。 ②[Y,U]=max(A) 返回两个行向量，Y向量记录A的每列的最大元素，U向量记录每列最大元素的行号。 两个向量或矩阵对应元素的比较 ①U=max(A,B) A,B是两个同型的向量或矩阵。结果U是与A,B同型的向量或矩阵，U的每个元素等于A,B对应元素的较大者。 ②U=max(A,n) n是一个标量。结果U是与A同型的向量或矩阵，U的每个元素等于A对应元素和n中的较大者。 min函数的用法和max完全相同。 例 求矩阵A的每行及每列的最大和最小元素，并求整个矩阵的最大和最小元。 程序如下 A=[13,-56,78;25,63,-235;78,25,563;1,0,-1]; max(A,[],2) %求每行最大元素 min(A,[],2) %求每行最小元素 max(A) %求每列最大元素 min(A) %求每列最小元素 max(max(A)) %求整个矩阵的最大元素 min(min(A)) %求整个矩阵的最小元素 2 求矩阵的平均值和中值 求矩阵或向量元素的平均值的函数是 mean，求中值的函数是 median。 mean（X）返回向量 X 的算术平均值。 mean（A）返回一个行向量，其第 i个元素是矩阵 A 的第 i列的算术平均值。 mean（A，dim） 当 dim 为 1 时，该函数等同于 mean（A），当 dim 为 2 时，返回一个列向量，其第 i个元素是矩阵 A 的第 i行的算术平均值。 median（X）返回向量 X 的中值。 median（A）返回一个行向量，其第 i个元素是矩阵 A 的第 i列的中值。 median（A，dim）当 dim 为 1 时，该函数等同于 median（A），当 dim 为 2 时，返回一个列向量，其第 i个元素是矩阵 A 的第 i行的中值。 中间值就是将数组按大小顺序排列，取中间的那个数作为中值，如果数组有偶数个，那么就取中间那两个数的平均值 3 矩阵元素求和与求积 矩阵和向量求和与求积的函数是 sum 和 prod sum（X）返回向量 X 各元素的和。 sum（A）返回一个行向量，其第 i个元素是矩阵 A