《科学计算引论》实验课：第三章 数值微积分.pptVIP

数值微分 数值微分 数值微分 例：设 用不同的方法求函数f(x)的数值导数，并在同一个坐标系中做出f (x)的图像。 解：用三种方法： 1、用一个5次多项式p(x)拟合函数f(x)，并对p(x)求一般意义下的导数dp(x)，求出dp(x)在假设点的值； 2、直接求f(x)在假设点的数值导数； 3、求出g(x)=f (x)的表达式 然后求f ‘ (x)在假设点的导数值。 数值微分 程序如下： f=inline(sqrt(x.^3+2*x.^2-x+12)+(x+5).^(1/6)+5*x+2); g=inline((3*x.^2+4*x-1)./sqrt(x.^3+2*x.^2-x+12)/2+1/6./(x+5).^(5/6)+5); x=-3:0.01:3; p=polyfit(x,f(x),5); %1用5次多项式p拟合f(x) dp=polyder(p); %对拟合多项式p求导数dp dpx=polyval(dp,x); %求dp在假设点的函数值 dx=diff(f([x,3.01]))/0.01; %2直接对f(x)求数值导数 gx=g(x); %3求函数f的导函数g在假设点的导数 plot(x,dpx,g.,x,dx,