matlab Subfunctions-functions handles 题目。帮看下

写一个function deriv(f, a, h) 使用以下的有限差分公式:
Df(a) = (f(a+h)-f(a-h))/2*h
求导 f ' (x) f(x) 在 x=a

然后用function testDeriv 去完成求导
f(x) = cos(x/2)+sin(2*x) 在区间 [-pi, pi]
然后输出plots f(x)和Df(x)在同一张图上。

1
将下面的代码保存为M文件：
function testDeriv
syms x
f = cos(x/2)+sin(2*x);
h = pi/50;
a = -pi:h:pi;
Df = deriv(f, a, h); % 完成求导
plot(a,Df); % Df(x)图
hold on
plot(a,subs(f,'x',a),'.-'); % f(x)图
function Df = deriv(f, a, h) % Subfunction
fah1 = subs(f,'x',a+h);
fah2 = subs(f,'x',a-h);
Df = (fah1-fah2)/(2*h);
2 将下面的代码复制到命令窗口，按回车键，运行即可
testDeriv

可以私聊我~