数学实验“插值型数值微分,Taylor展开式数值微分”实验报告(内含matlab程序)..docVIP

西京学院数学软件实验任务书 课程名称 数学软件实验 班级 数0901 学号 0912020107 姓名 李亚强 实验课题 插值型数值微分，Taylor展开式数值微分 实验目的 熟悉插值型数值微分，Taylor展开式数值微分 实验要求 运用Matlab/C/C++/Java/Maple/Mathematica等其中一种语言完成 实验内容 插值型数值微分，Taylor展开式数值微分 成绩 教师  实验二十三实验报告 实验名称：插值型数值微分，Taylor展开式数值微分。 实验目的：进一步熟悉插值型数值微分，Taylor展开式数值微分。 实验要求：运用Matlab/C/C++/Java/Maple/Mathematica等其中一种语言完成程序设计。 实验原理： 插值型数值微分： 由得： 当为插值节点时，上式简化为： 故一般限于对节点上的导数值采用插值多项式的相应导数值进行近似计算，以便估计误差。 一般的 Taylor展开式数值微分： 理论基础：Taylor展开式 我们借助Taylor展开式，可以构造函数在点的一阶导数和二阶导数的数值微分公式。取步长则： （1） 所以 （2） 同理 （3） （4） 式(2)和式(4)是计算的数值微分公式，其截断误差为，为提高精度，将Taylor展开式多写几项 两式相减得 （5） 上式为计算的微分公式，其截断误差为O(h2),比式（2）和(4)精度高。 两式相加，如果，则有 （6） 其中, 式(6)是计算的数值微分公式,其截断误差为。 实验内容： %牛顿插值型微分 function [df,A,P]=diffnew(X,Y) n=length(X); A=Y; for j=2:n for i=n:-1:j A(i)=(A(i)- A(i-1))/(X(i)-X(i-j+1)); end end x0=X(1);df=A(2);chsh=1;m=length(A)-1; for k=2:m chsh=chsh*(x0-X(k)); df=df+chsh*(A(k+1)); end P=poly2sym(A); %拉格朗日插值型微分 function [C,L,dyk,k]=ndaolag(X,Y,n) m=length(X); n=m; L=ones(m,m); for k=1:m V=1; for i=1:m if k~=i V=conv(V,poly(X(i)))/(X(k)-X(i)); end end L1(k,:)=V; l(k,:)=poly2sym (V); end C=Y*L1;L=Y* l; syms x dyk for k=1:n k; dyk=diff(L,x,k) end - 0 -