数值微积分适合求解没有或很准求出微分或积分表达式的问题的计算。
1、数值微分
(1)数值差分与差商
任意函数f(x)在x0点的导数是通过极限定义的:
如果去掉极限定义中h趋向于0的极限过程,得到函数在x0点处以h为步长的向前差分。
当步长h充分小时得到函数在x0点处,以h为步长的向前差商
函数f(x)在x0点的微分接近于函数在该点的差分。而f(x)在x0点的导数接近于函数在该点的差商。
(2)数值微分的实现
MATLAB提供了求 向前差分的函数diff,其调用格式有三种:
(1)dx=diff(x): 计算向量x的一阶向前差分,dx(i)=x(i+1)-x(i),i=1,
2,…,n-1。n是向量x的元素的个数
(2)dx=diff(x,n): 计算向量x的n阶向前差分。例如,diff(x,2)=diff(diff(x)):x的二阶差分等于x的一阶差分再求一阶差分
(3)dx=diff(A,n,dim): 计算矩阵A的n阶差分,dim=1时(默认状态)按列计算差分;dim=2,按行计算差分。
注意:dff函数计算的是向量元素间的差分,故差分向量元素的个数比原向量少了一个;同样对于矩阵来说差分后的矩阵比原矩阵少了一行或一列。另外计算差分之后,可以用f(x)在某点处的差商作为其导数的近以值。
例1设f(x)=sin x,在[0,2π]范围内随机采样,计算f’(x)的近似值,并与
理论值f’(x)=cos x进行比较。
x = [0,sort(2*pi*rand(1,5000)),2*pi]; %确定x向量,其首尾元素分别是0和2π,中闻是0到2π开区间的随机数并按照从小到大排序
y = sin(x);
f1 = diff(y)./diff(x); %利用diff函数求一阶向前差分,y的差分除以x的差分,得到差商向量f1。差商作为导数的近似值
f2 = cos(x(1:end-