﻿功能: 求函数一阶偏导值

格式: DiffFirstOrder(f,x)

f  : 符号函数对象
x  : 数值或者矩阵变量

说明:
1. 执行成功, 返回函数在x处对应的一阶偏导数值
2. 函数f里的变量个数必须为1
3. 本函数采用数值方法计算f在x处的导数值
4. 本函数主要用于无法符号求导的表达式,进行数值求导
5. 本函数采用4阶方法求解一阶偏导
$$
f(x)'\approx \dfrac{-f(x+2h)+8f(x+h)-8f(x-h)+f(x-2h)}{12h}
$$

参考: CRC standard mathematical tables and formulate[M]. CRC press, 2002.

例子:
f = @sin(x);
t = 0:0.2:1;
df = DiffFirstOrder(f,t)//回车得到如下结果,可以对于 df与理论解cos(t)
df =
[ 1.00000000000000    0.98006657784121    0.92106099400269    0.82533561490936    0.69670670934704    0.54030230586797 ]