实验一关于多项式的MATLAB命令.docVIP

第2章 实验一关于多项式的MATLAB命令 实验目的：熟悉并能灵活使用多项式的matlab命令。 2.1 多项式的幂系数形式 一个多项式的幂系数其形式可以表示为： 其中是多项式的阶数，是系数且。多项式也可表示为嵌套形式： 或因子形式： 其中是多项式根。例如，多项式 可以等价写为： 或 2.2 关于多项式的MATLAB命令： 常用的多项式的MATLAB命令有 roots; poly; polyval; polyfit 例2.1 已知 ， 在matlab中输入表示为： y=x^4-6*x^3+3*x^2-26*x-24; 用矩阵表示为 p=[1 -6 3 26 -24] roots：多项式零点可用命令roots求得。 roots(p) % 求多项式的零点的命令 ans = -2.0000 4.0000 3.0000 1.0000 输入的多项式、求多项式y=x^4-6*x^3+3*x^2-26*x-24的零点的命令，以及得出的结果如图2.1 图 2.1 多项式求根命令及结果 poly： 当所有零点已知时，是否可以恢复原多项式？答案在一定程度上是肯定的，命令是poly(r)。 poly给出的是经过标准化的，即最高项次数的系数为1。 如果重根存在，这种转化可能会降低精度，例如 y=(x-1)^6 polyval：polyval 可以计算多项式的值 p=[1 -6 3 26 -24] xi=2.5; yi=polyval(p,xi) % yi = 5.0625 polyfit：polyfit 给定n+1个点将可以唯一确定一个n阶多项式 x=[1.1 2.3 3.9 5.1]; y=[3.887 4.276 4.651 2.117]; a=polyfit(x,y,length(x)-1); % a = -0.2015 1.4385 -2.7477 5.4370 % 这是多项式的系数向量。参数length(x)-1是多项式的阶数。 x1=1:0.1:5.3; % 自变量取值范围 C=[-0.2015 1.4385 -2.7477 5.4370]; % 多项式的系数 y1=polyval(C,x1); % 计算多项式的值 plot(x1,y1,x,y,o) % 作图 图 2.2 .m文件 图 2.3 （通过4个数据点的三次多项式曲线图） 2.3 关于多项式的微分与积分 对多项式 p=[1 -6 3 26 -24]积分，可以用 poly_itg 计算，调用格式为： c = poly_itg(d) 程序2.1： function py=poly_itg(p) n=length(p) py=[p.*[n:-1:1].^(-1),0] 例如 p=[1 -6 3 26 -24]； c = poly_itg(p) 结果： 0.2000 -1.5000 1.0000 13.0000 -24.0000 0 一阶导数的系数可用polyder计算，调用格式为 b=polyder(p) b=polyder(p) 结果： b = 4 -18 6 26 2.4 关于多项式的加减法 程序：ploy_add.m function p3=poly_(p1,p2) n1=length(p1); n2=length(p2); if n1==n2 p3=p1+p2;end if n1n2 p3=p1+[zeros(1,n1-n2),p2];end if n1n2 p3=[zeros(1,n2-n1),p1]+p2;end 例如： p1=[1 2 3]; p2=[3 2 1]; p3=poly_(p1,p2) 结果： p3 = 4 4 4