[工学]计物05_matlab05_微积分与微分方程数值解.pptVIP

第五讲 数值微积分方法 ——以微分（差分）、积分的数值法以及微分方程数值解为例介绍数值方法 §5.1 差分 数值微分非常困难。积分描述了一个函数的整体或宏观性质，而微分则描述一个函数在一点处的斜率，这是函数的微观性质。因此积分对函数的形状在小范围内的改变不敏感。而微分却很敏感。一个函数小的变化，容易产生相邻点的斜率的大的改变。 由于微分这个固有的困难，所以尽可能避免数值微分，特别是对实验获得的数据进行微分。在这种情况下，最好用最小二乘曲线拟合这种数据，然后对所得到的多项式进行微分。 *diff*diff(x)——对向量,[x(2:n)-x(1:(n-1))];diff(X)——对矩阵列的差分[X(2:n,:)-X(1:(n-1),:)];diff(X,n)——n阶差分;diff(X,n,DIM)——沿DIM维n阶差分; 给定一些描述某函数的数据，MATLAB提供了一个计算其非常粗略的微分的函数。这个函数命名为diff，它计算数组中元素间的差分。 y=f(x)的微分可近似为： 例1 x=[0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1] y=[-.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; % data dy=diff(y) ./ di