这是一个GUI,它在一些等间距的点上执行函数的数值微分。它还附带了一个代码,授予用于数值微分的系数。
图片和示例应该足够理解如何使用该文件。
例子:
npoints = 3;
= 1;
d = datnum (npoints顺序)
d =
- 1.5% 2 - 0.5%前进
-0.5 -0 0.5%中央
0.5% -2 1.5%向后
结果是一个由系数组成的矩阵
%用于数值上的微分,像这样:
x = 1;
内联(f = cos (x))
h = 1;
s = (d (1,1) * f (x) + d(1、2)* f (x + h) + d(1、3)* f (x + 2 * h)) / h ^秩序
s =
-0.8444
s = (d (2, 1) * f (x h) + d (2, 2) * f (x) + d(2、3)* f (x + h)) / h ^秩序
s =
-0.84007
s = (d (3,1) * f (x - 2 * h) + d (3 2) * f (x h) + d (3,3) * f (x)) / h ^秩序
s =
-0.84413
正确的答案是s = -0.84147
该代码使用符号数学工具箱获得真实值(以便计算误差)。如果您没有符号数学工具箱,那么您将无法享受这个好处(程序仍然有效)。
引用作为
Husam Aldahiyat(2022年)。数值微分(//www.tianjin-qmedu.com/matlabcentral/fileexchange/22807-numerical-differentiation), MATLAB中央文件交换。检索.