Matlab7.x程序设计第4章.pptVIP

第4章 MATLAB 数值计算(2) ;4.4 多项式函数;4.4.2多项式的算术运算;4;2 乘法运算 ;3 除法运算 [Q,r]=deconv(P1,P2) 对多项式P1和P2作除法运算。其中Q返回多项式P1除以P2的商式，r返回P1除以P2的余式。 注意deconv是conv的逆函数，即有P1=conv(P2,Q)+r。 ;例4-45 计算 Q=deconv([1 8 0 0 -10],[2 -1 3]) Q = 0.5000 4.2500 1.3750 [Q,r]=deconv([1 8 0 0 -10],[2 -1 3]) 执行结果如下： Q = 0.5000 4.2500 1.3750 r = 0 0 0 -11.3750 -14.1250 由执行结果可知商是: 余式是: ;4.4.3 导函数;例4-46 求的导数;a=[3 -2 1];b=[4 5 6];polyder(a,b) 执行结果为： ans = 48 21 24 -7 结果是 ;例4-48 求 的导数。 ;例4-49 求有理分式 的导数。;4.4.4 多项式求根;例4-50 求多项式的 根;例4-51 求方程 的根。 ;例 4-52已知 (1) 计算 的全部根。(2) 由方程 的根构造一个多项式并与 进行对比。 ;执行结果为： X = -0.3046 + 1.6217i -0.3046 - 1.6217i -1.0066 1.0190 0.5967 G = 1.0000 0.0000 1.3333 -1.6667 -2.4000 1.6667 注意：构造的多项式的首项系数为1。 ;4.4.5多项式估值;例4-53 已知 分别计算 和 时 的值。;2 矩阵多项式求值;例4-54 当x取 时求 的值。;polyval(p,a) 执行结果： ans = -4 -10 8 8 200 4 106 44 332 ;例4-55 当x=8时求(x-1)(x-2) (x-3)(x-4) 的值。 ;4.4.6 部分分式函数;例4-56 求 的部分分式。 ;执行结果为： res = -6.6667 5.0000 1.6667 poles = -4.0000 -3.0000 -1.0000 k = [ ] ;4.4.7多项式积分;例4-57 求 。 ;4.5插值和拟合;method是插值方法，可用的方法有： ‘linear’ 默认方法，线性插值。 ‘nearest’ 最邻近插值。 ‘spline’ 三次样条插值。 ‘cubic’ 三次插值，要求x的值等距离。 所有插值方法均要求x是单调的。 注意：X1的取值范围不能超出X的给定范围，否则，会给出“NaN”错误。 ;例4-58 求正弦、余弦函数在区间[0 1]内间 隔为0.25的各点的值。;2 二维数值插值;4.5.2 数据拟合;例4-59已知 在[1，3]区间10个采样点的函数值，求 的4次拟???多项式p(x)。 ;例4-60 已知数表如下 x … 9 10 11 12 … f(x)… 2.6093 1.7586 1.3979 1.9483 … ?求4次拟合多项式f(x)，并计算f(10.5)的近似值。 ;4.6数值微分与积分;例4-61 求矩阵 的差分。;4.6.2数值积分 1数值积分基本原理 高等数学中求解定积分的数值方法有简单的梯形法、辛普生(Simpson)法、牛顿－柯特斯(Newton-Cotes)法等，基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1]，i=1,2,…,n，其中x1=a，xn+1=b。这样求定积分问题就分解为求和问题。 ;2 一重积分;[I,n]=quad8(fname,a,b,tol,trace) 基于牛顿－柯特斯法求定积分。其中参数的含义和quad函数相似，只是tol的缺省值取10-6。 ;例4-62 求 ;例4