[理学]第四章 多项式与插值.ppt

二、 多项式计算 1. 多项式求根 求多项式 p(x)的根的函数是roots(P)，这里，P是 p(x)的系数向量，该函数返回方程 p(x)=0 的全部根 (含重根，复根)。 2. 多项式求值 求多项式 p(x)在某点或某些点的函数值的函数是 polyval(P, x)。若x为一数值，则求多项式在该点的 值；若x为向量或矩阵，则对向量或矩阵中的每个元 素求其多项式的值。 例1 已知一个多项式 (1)计算f (x)=0 的全部根。 (2)由方程f (x)=0的根构造一个多项式 g(x)，并与 f (x) 进行对比。 (3)计算f (5)、f (7.8)、f (9.6)、f (12.3)的值。 3. 多项式的四则运算 (1)多项式的加减法 例2 设有两个多项式，计算： (1)求f(x)+g(x)、f(x)-g(x)。 (2)求f(x)·g(x)、f(x)/g(x)。 4. 多项式的微分与积分 （1）对多项式求导数的函数是： p=polyder(P) 求多项式P的导函数 p=polyder(P,Q) 求P*Q的导函数 [p,q]=polyder(P,Q) 求P/Q的导函数，导函数的 分子存入p，分母存入q。 （2）对多项式的积分函数： d=poly_itg(c) d是多项式c积分后的系数，但 不包括积分常数 例3 求有理分式的导数。 5 函数的最大值与最小值 MATLAB中用于求最小值的函数是： fmin(f,a,b) 求单变量函数f(x)在区间(a,b)上的最小值点。 fmins(F,X0) 求多变量函数F(x)在估计值X0附近的最小值点。 MATLAB没有专门提供求函数最大值点的函数，但只要注意到-f(x)在区间(a,b)上的最小值点就是f(x)在(a,b)的最大值点，所以fmin(-f,a,b)返回函数f(x)在区间(a,b)上的最大值。 例 求函数f(x)在区间(-10，1)和(1，10)上的最小值点。 首先建立函数文件fx.m： function f=f(x) f=x-1/x+5; return 再在MATLAB命令窗口，输入命令： fmin(fx,-10,-1) %求函数在区间(-10，-1)内的最小值点 fmin(f,1,10) %求函数在区间(1，10)内的最小值点。注意函数名f不用加 例 设有函数f(x,y,z)，求函数f在(0.5,0.5,0.5)附近的最小值。 建立函数文件fxyz.m： function f=f(u) x=u(1);y=u(2);z=u(3); f=x+y.^2./x/4+z.^2./y+2./z; return 在MALAB命令窗口，输入命令： U=fmins(fxyz,[0.5,0.5,0.5]) %求函数的最小值点 fxyz(U) %求函数的最小值 可见，误差峰值出现在端点附近的区间里，这是由于 的局部峰值在端点附近。 二维插值是基于与一维插值同样的基本思想。然而，正如名字所隐含的，二维插值是对两变量的函数z=f(x, y) 进行插值。为了说明这个附加的维数，考虑一个问题。设人们对海水的温度分布估计感兴趣，给定的温度值取自海水表面均匀分布的格栅。 采集了下列的数据： width=1:5; % index for width of plate (i.e.,the x-dimension) depth=1:3; % index for depth of plate (i,e,,the y-dimension) temps=[82 81 80 82 84; 79 63 61 65 81; 84 84 82 85 86] % temperature data temps = 82 81 80 82 84 79 63 61 65 81 84 84 82 85 86 矩阵temps表示整个海水的温度分布。temps的列与下标depth或y-维相联系，行与下标width或x-维相联系。为了估计在中间点的温