中点公式法与五点公式法求数值微分.docVIP

《MATLAB程序设计实践》  1、编程实现以下科学计算算法，并举一例应用之。 “中点公式法和五点公式法求数值微分” 解：中点公式法和五点公式法求数值微分如下： ：中点公式法求导数应用实例。采用中点公式法求函数f=在x=4处的导数。 解：在MATLAB命令窗口中输入： df=MidPoint(sqrt(x),4) 输出结果为： df= 0.2500 采用中点公式法求函数f=在x=4处的导数为0.25，而导数的精确值也是0.25 . 详见以下： 中点公式法流程图： 开始 开始 判断输入参数个数 判断输入参数个数 If nargin=2,h=0.1 if nargin=3,h=0 输出：h不能为0结束用中点公式求数值微分 输出：h不能为0 结束 用中点公式求数值微分 源代码： function df=MidPoint(func,x0,h) if nargin == 2 h = 0.1; else if (nargin == 3 h == 0.0) disp(h2??ü?a0￡?); return; end end y1 = subs(sym(func), findsym(sym(func)),x0+h); y2 = subs(sym(func), findsym(sym(func)),x0-h); df = (y1-y2)/(2*h); 运行结果如下： ：五点公式法求导数应用实例。采用五点公式法求函数f=sin(x)在x=2处的导数。 解：在MATLAB命令窗口中输入： df1=FivePoint(sin(x),2,1); df2=FivePoint(sin(x),2,2); df3=FivePoint(sin(x),2,3); df4=FivePoint(sin(x),2,4); df5=FivePoint(sin(x),2,5); 用五种方法得到的结果为： df1= -0.4161 df2= -0.4161 df3= -0.4161 df4= -0.4161 df5= -0.4161 而函数在f=sin(x)在x=2的导数为cos(2)=-0.4161,从上面的结果来看，五点公式的精度是很高的。 详见以下： 五点公式法流程图： 开始 开始 判断输入参数个数 判断输入参数个数 If nargin=3,h=0.1 if nargin =4,h=0 输出: h不能为0 输出: h不能为0 用公式三求导用公式五求导用公式四求导用公式二求导用公式一求导 用公式三求导 用公式五求导 用公式四求导 用公式二求导 用公式一求导 结束 结束 源代码： function df=FivePoint(func,x0,type,h) %??μ?1?ê?·¨,?óè?oˉêyfunc?úx0′|μ?êy ×￡íòê?è?òa %oˉêy??￡ofunc %?óμ?μ?￡ox0 %1?ê?μ?D?ê?￡otype￡¨è?1,2,3,4,5,￡? %à?é￠2?3¤￡oh %μ?êy?μ￡odf if nargin ==3 h=0.1; else if (nargin ==4h==0.0) disp(h2??ü?a0); return; end end y0 = subs(sym(func),findsym(sym(func)),x0); y1 = subs(sym(func),findsym(sym(func)),x0+h); y2 = subs(sym(func),findsym(sym(func)),x0+2*h); y3 = subs(sym(func),findsym(sym(func)),x0+3*h); y4 = subs(sym(func),findsym(sym(func)),x0+4*h); y_1 = subs(sym(func),findsym(sym(func)),x0-h); y_2 = subs(sym(func),findsym(sym(func)),x0-2*h); y_3 = subs(sym(func),findsym(sym(func)),x0-3*h); y_4 = subs(sym(func),findsym(sym(func)),x0-4*h); switch type case 1, df=(-25*y0+48*y1-36*y2+16*y3-