polyfit:给定n+1个点将可以唯一确定一个n阶多项式。利 用命令 polyfit可容易确定多项式的系数 例 >>X=[1.1,2.3,3.9,5.1 y=[3.8874.276,4.651,2.117]; >a=polyfit(x, y, length(x)-1) a -0.201514385-2.74775.4370 poly2sym(a) ans -403/2000*x^3+2877/2000*x^2-27477/10000*x+5437/1000 多项式为y=-02015x3+1.4385x2-2.7477x+54370 Polyfit的第三个参数是多项式的阶数
• polyfit:给定n+1个点将可以唯一确定一个n阶多项式。利 用命令polyfit可容易确定多项式的系数。 例： >> 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 >> poly2sym(a) ans = -403/2000*x^3+2877/2000*x^2-27477/10000*x+5437/1000 多项式为 Polyfit的第三个参数是多项式的阶数。 3 2 y x x x = − + − + 0.2015 1.4385 2.7477 5.4370
多项式积分 y=c1x+C2x+…+cnx+ n+1 Y Jdr Cx++=2x"+.+x Tom+l-c+Cnt n+1 功能:求多项式积分 调用格式:py= poly itg(p) p:被积多项式的系数 py:求积后多项式的系数 poly itg. m function py=poly itg(p) Length(p) py=[p.*[ni:-1:1]^(-1),0 不包括最后一项积分常数
多项式积分： 功能：求多项式积分 调用格式：py=poly_itg(p) p:被积多项式的系数 py：求积后多项式的系数 poly_itg.m function py=poly_itg(p) n=length(p); py=[p.*[n:-1:1].^(-1),0] 不包括最后一项积分常数1 2 1 2 1 2 1 2 n n n n n c c c Y ydx x x x c x c n n + = = + + + + + + + + 1 1 2 1 n n n n y c x c x c x c − = + + + + +
多项式微分: y=Cr +c2x+.''tCnxtcn+ y=nC1x+(n-1)c2x2+…+C Polder:求多项式一阶导数的系数 调用格式为:b= polymer(c) c为多项式y的系数,b是微分后的系数, 其值为: nC1,(n-1)c2,…,Cn
多项式微分： • Polyder: 求多项式一阶导数的系数。 调用格式为: b=polyder(c ) c为多项式y的系数，b是微分后的系数， 其值为： 1 1 2 1 n n n n y c x c x c x c − = + + + + + ' 1 2 1 2 ( 1) n n n y nc x n c x c − − = + − + + 1 2 [ ,( 1) , , ] nc n c c − n
两个多项式的和与差: 71 Va =ax tax ∴+a.x+a m+1 yb=bx+b2x+…+bnx+bn 命令 poly add:求两个多项式的和,其调用格式为 poly add(a, b) 多项式a减去b,可表示为 c- poly add(a-b)
两个多项式的和与差： 命令poly_add:求两个多项式的和,其调用格式为: c= poly_add(a,b) 多项式a减去b，可表示为: c= poly_add(a,-b) 1 1 2 1 n n b n n y b x b x b x b − = + + + + + 1 1 2 1 m m a m m y a x a x a x a − = + + + + +
功能:两个多项式相加 调用格式:b= poly add(plp2) b:求和后的系数数组 poly add. m function p3=poly add(pl, p2) nI=length(pl) n2=length(p2) ifnI==n2 p3-p1+p2; end if nI>n2 p3=p1+zeros(l, n1-n2),p2]: end if nI<n2 p3=zeros(l, n2-n1)p1+p2; end
功能：两个多项式相加 调用格式：b=poly_add(p1,p2) b:求和后的系数数组 • poly_add.m function p3=poly_add(p1,p2) n1=length(p1); n2=length(p2); if n1==n2 p3=p1+p2;end if n1>n2 p3=p1+[zeros(1,n1-n2),p2];end if n1<n2 p3=[zeros(1,n2-n1),p1]+p2;end