MATLAB之怎样利用MATLAB中值差分法求一阶二阶导数

**

MATLAB初学之怎么利用中值差分法求一阶二阶导数

**

我们最近在学习MATLAB。在MATLAB中怎么求导数？
MATLAB中有专门求导的函数
针对f（x）类的函数：
diff(f,x) ：求f关于x的导数
diff（diff(f,x)，x）：求f二阶导数
针对f（x，y）类的函数求偏导：
diff(f,x) ：求f关于x的偏导数
diff（diff(f,x)，y）：求f关于x的偏导数再对y进行二阶偏导
但是以上方法求导不够精确，diff是基于向前欧拉差分或者向后欧拉差分实现的，在求导的方法中，中心差分法是精度比较高的。今天先大家分享如何利用中值差分法求一阶二阶导数。

##原理：

##程序源代码，以方程f=cos(x).*exp(-x.^2/2)为例：
clear all
clc
h=0.1;
x=-4:h:4;
f=cos(x).*exp(-x.^2/2);
for j=2:(length(x)-1)
%(从第二点开始，因为第一个点求不出来)
df(j)=(f(j+1)-f(j-1))/2/