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

格式: DiffSecondOrder(f,x)

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

说明:
1. 执行成功, 返回函数在x处对应的二阶偏导数值
2. 函数f里的变量个数必须为1
3. 本函数采用数值方法计算f在x处的导数值
4. 本函数主要用于无法符号求导的表达式,进行数值求导
5. 本函数采用6阶方法求解二阶偏导
$$
f(x)'' \approx \dfrac{2f(x+3h)-27f(x+2h)+270f(x+h)-490f(x)+270f(x-h)-27f(x-2h)+2f(x-3h)}{180h^2}
$$
6. 还有如下4阶方法求解二阶偏导
$$
f(x)'' \approx \dfrac{-f(x+2h)+16f(x+h)-30f(x)+16f(x-h)-f(x-2h)}{12h^2}
$$
例子:
f = @sin(x);
t = 0:0.2:1;
df = DiffSecondOrder(f,t)//回车得到如下结果,可以对于 -sin(t)
df =
[ 0.00000000000000   -0.19866933079461   -0.38941834230722   -0.56464247339416   -0.71735609090024   -0.84147098480636 ]