论述matlab数值微积分计算,第七章 基于MATLAB的科学计算—数值微积分1

科学计算—理论、方法

及其基于MATLAB 的程序实现与分析

数 值 微 积 分

§1 数值微分

对于给定的函数()x f y =，如果 １、()x f y =的函数关系式比较复杂时；

２、()x f y =未知，而仅仅知道()x f y =在1+n 个相异点k x ，n k ,,1,0 =处 的函数值k y ；

则希望能用相对简单的计算方法，求得()x f y =导数的(近似)值。

基于上述考虑，选择的方法之一是利用函数()x f y =的插值多项式的导数作为函数()x f y =导数的近似值，例如Lagrange 插值多项式，由于

()()()()∑==≈n

k k k n x f x l x L x f 0 (１)

因而有

()()x L x f n

'≈' (２) 这里需要说明一点的是，尽管()x f 和()x L n 的函数值可能相差不

多，但是它们的导数有可能相差很大，如下面的例子

例1： 考虑函数()2

2511

x x f +=

在区间[-1，1]的插值问题，取区间

[-1，1]的11个点k x k ?+-=2.01，10,,1,0 =k ，作函数()2

2511

x x f +=的

10次插值多项式：

()18552.163597.1234338.3819095.4949417.220246810+-+-+-=x x x x x x L n

函数()x f 和插值多项式()x L n 的导数分别为

()()

2225150x

x

dx x df +-=

()x x x x x dx

x dLn 7.334.4936.22883.39594.22093579-+-+-= 对函数()x f 和插值多项式()x L n 及其导数分别比较，结果如图所示：

Derivative_Runge