matlab基础数值运算基础参考.ppt

第二章 数组的运算基础 第四章 Matlab数值运算基础 多项式运算 数据插值 数据分析 第一节 多项式运算 一、多项式创建 二、多项式求值和求根 三、多项式乘法和除法 四、多项式的微分和积分 五、多项式曲线拟合 一、多项式创建 1. 系数矢量直接输入法 【例4-1】在Matlab中，将多项式x5-5x4-4x3+3x2-2x+1用 行向量表示如下： p=[1 -5 -4 3 -2 1] p= 1 -5 -4 3 -2 1 y= poly2sym(p) y= x^5-5*x^4-4*x^3+3*x^2-2*x+1 【例4-2】某多项式的行矢量 p=[ 1 0 3 2 -1]，创建与之 相对应的多项式。 p=[ 1 0 3 2 -1]; y=poly2sym( p) y= x^4+3*x^2+2*x-1 【例4-4】根矢量 r =[ 2 -1+i -1-i]，创建对应此根 的多项式。 r =[ 2 -1+i -1-i]; p=poly(r) p= 1 0 -2 -4 y =poly2sym (pr) y= x^3-2*x-4 二、多项式求值和求根 1. 多项式求值，调用格式如下： y=polyval( p ,x) y=polyvalm(p ,x) 通过示例观察polyval和polyvalm的用法。 【例4-5】求多项式 2x3+2x+1在0、1、2处的值。 p=[2 0 2 1]; pv =polyval( p, [0 1 2]) pv= 1 5 21 【例4-6】求多项式2x3+2x+1 对于矩阵［1 0 ; 0 1］及 标量5的值。 p =[ 2 0 2 1]; pv1= polyvalm( p , [ 1 0 ;0 1]) pv1 = 5 0 0 5 pv2=polyvalm(p,5) pv2= 261 pv3=polyval(p,5) pv3= 261 可见，多项式在一点的值，既可以用polyval，也可以用polyvalm。 三、多项式乘法和除法 多项式的乘、除法分别用函数conv和deconv来实现。 【例4-8】求多项式a(x) =x3+x+1和b(x)= x2-1的乘积。 a= [1 0 1 1]; b=[ 1 0 -1]; c =conv (a ,b) c = 1 0 0 1 -1 -1 cx=poly2sym(c) 四、多项式的微分和积分 多项式的微分和积分分别用polyder和polyint实现。 【例4-10】 求多项式4x3 +3x2-2x+1的微分和积分。 p=[ 4 3 -2 1]; polyder(p) ans = 12 6 -2 poly2sym(ans) 五、多项式曲线拟合 【例4-11】 对数据 x=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1]和y=[-0.232 0.647 1.877 3.565 5.134 7.443 9.221 10.011 11.678 12.566 13.788]做二次曲线拟合，并图示拟合曲线和原来的数据。 x=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1]; y=[-0.232 0.647 1.877 3.565 5.134 7.443 9.221 10.011 11.678 12.566 13.788]; p=polyfit(x ,y ,2);