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

《MATLAB程序设计实践》

1、编程实现以下科学计算算法，并举一例应用之。

“中点公式法和五点公式法求数值微分”

解：中点公式法和五点公式法求数值微分如下：

：中点公式法求导数应用实例。采用中点公式法求函数f=在x=4处的导数。

解：在MATLAB命令窗口中输入：

df=MidPoint(sqrt(x),4)

输出结果为：

df=

0.2500

采用中点公式法求函数f=在x=4处的导数为0.25，而导数的精确值也是0.25

.

详见以下：

中点公式法流程图：

开始

开始

判断输入参数个数

判断输入参数个数

Ifnargin=2,h=0.1ifnargin=3,h=0

输出：h不能为0结束用中点公式求数值微分

输出：h不能为0

结束

用中点公式求数值微分

源代码：

functiondf=MidPoint(func,x0,h)

ifnargin==2

h=0.1;

elseif(nargin==3h==0.0)

disp(h2??ü?a0￡?);

return;

end

end

y1=subs(sym(func),findsym(sym(func)),x0+h);

y2=subs(sym