exn534.m

MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5－3－4程序
% 多项式函数的泰勒级数展开
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月

clear,clf
a=input('输入多项式系数向量a=[an,...,a0]=');		% 输入参数
x0=input('输入展开点的坐标值x0= ');
[xm]=input('输入展开的坐标区间[xmin,xman]= ');
x=linspace(xm(1),xm(2));					% 设定自变量数组
y=polyval(a,x);y0=polyval(a,x0);		% 求y的准确值和y在x0点的值y(x0)
Da=polyder(a),Dy0=polyval(Da,x0);		% 求x0点的一阶导数的值
D2a=polyder(Da),D2y0=polyval(D2a,x0); % 求x0点的二阶导数的值
D3a=polyder(D2a),D3y0=polyval(D3a,x0); % 求x0点的三阶导数的值
yt(1,:)=y0+Dy0*(x-x0);						% 求y的一阶泰勒展开
yt(2,:)=yt(1,:)+D2y0*(x-x0).^2/prod(1:2); 	% 求y的二阶泰勒展开
yt(3,:)=yt(2,:)+D3y0*(x-x0).^3/prod(1:3); 	% 求y的三阶泰勒展开
plot(x,y,'.',x,yt(1:3,:)),grid		% 绘图, 准确值用点线表示