【正文】 都 能 準(zhǔn) 確 成 立 但 對( ) 1 , , , mf x x x?1() mf x x ??不能準(zhǔn)確成立 . 顯然，一個求積公式的代數(shù)精度越高， 它就能對更多的被積函數(shù) f(x)準(zhǔn)確成立， 從而具有更好的實際計算意義。0( 1)1[ ] ( ) [ ]( ) ()1[ ] ( ) ( )( 1 ) !插 值 型 求 積 公 式bnann k kkbkkabnnaI f f x dx I R fI A f xA l x dxR f f x dxn??????? ? ???????? ???????????( 1 )0, [ ] 0( ) ( )nnnbkkakf n f x R ff x dx A f x????? ??若 為 次 數(shù) 的 多 項 式 則 ( )= 0, 從 而 此 時 也就是說，當(dāng)被積函數(shù) f為次數(shù)不超過 n 的多項式時，其相應(yīng)的插值型求積公式不是近似公式，而是準(zhǔn)確公式。 需要做的工作： 1. 確定節(jié)點和求積系數(shù)； 2. 估計余項； 3. 討論公式的算法設(shè)計及其數(shù)值穩(wěn)定性。 (3) 記數(shù)值積分公式為 0,ninnniiI A fI I R????? 即 特點： 把求積過程（極限過程）轉(zhuǎn)化為有限次的乘法與加法的代數(shù)運算。 0 1 1 2 11010()2 2 22bnnannni i iiif f f f f ff x d x h h hffh h f A f????? ? ?? ? ? ??? ? ????0 1 2 111,22 nnA h A A A h A h?? ? ? ? ? ?(2) 數(shù)值積分的一般形式 數(shù)值積分的一般形式是： 0()nbi i naif x d x A f R?????其中， fi 是函數(shù) f(x)在節(jié)點 xi 上的函數(shù)值，它可能以列表形式給出，也可以是由函數(shù)的解析式計算出的函數(shù)值； Ai 稱為節(jié)點 xi 上的權(quán)系數(shù)，也稱求積系數(shù)。 從 幾何上 看，就是計算 曲邊梯形面積 的近似值。?? ba dxxffI )()(對于積分 公式有則由的原函數(shù)如果知道 L e i b n i zN e w t o nxFxf ?),()(?ba dxxf )( )()()( aFbFxF ba ???但是在工程技術(shù)和科學(xué)研究中 ,常會見到以下現(xiàn)象 : 的一些數(shù)值只給出了的解析式根本不存在 )(,)()1( xfxf不是初等函數(shù)如求不出來的原函數(shù) )(,)()()2( xFxFxf求原函數(shù)較困難的表達(dá)式結(jié)構(gòu)復(fù)雜 ,)()3( xf傳統(tǒng)方法的困境 以上這些現(xiàn)象 ,NewtonLeibniz很難發(fā)揮作用 ! 只能建立積分的近似計算方法 數(shù)值積分 正是為解決這樣的困難而提出來的， 不僅如此，數(shù)值積分也是微分方程數(shù)值解法的工具之一。第七章 微積分的數(shù)值計算方法 ? 傳統(tǒng)方法的困境 ? 數(shù)值積分的基本思想 ? 數(shù)值積分的一般形式 ? 代數(shù)精度問題 求函數(shù) f(x) 在區(qū)間 [a,b] 上的定積分 ()baI f x d x? ?是微積分學(xué)中的基本問題。 返回章 基本概念 167。 數(shù)值積分的基本思想 數(shù)值積分 是計算定積分的具有一定精度的近似值的各種計算方法。 最簡單的辦法，是用許多 小矩形之和 近似曲邊梯形 的面積，如 圖 70所示，這就是 矩形公式 ： 圖 70 矩形規(guī)則 y x a=x0 x1 x2 xi xi+1 xn1 xn =b f0 f1 f2 fi fi+1 fn1 fn f(x) 1000 1 2 1( ) ( )(),0nnbi i iaiinnf x d x f x h A fbaA A A A h An???????? ? ? ? ? ? ????(1) 圖 71 梯形規(guī)則 x a=x0 x1 x2 xi xi+1 xn1 xn =b y f0 f1 f2 fi fi+1 fn1 fn f(x) 如果改用許多小梯形之和近似曲邊梯形的面積，如圖 71，就會更精確些，這就是 梯形公式 。 正是由于權(quán)系數(shù)的構(gòu)造方法不同，從而決定了數(shù)值積分的不同方法。 xi為節(jié)點 ， Ai 為求積系數(shù)。 最常用的一種方法是利用插值多項式來構(gòu)造數(shù)值求積公式 , 具體步驟如下 : 上取一組節(jié)點在積分區(qū)間 ],[ babxxxa n ????? ?10次插值多項式的作 nxf )(???nkkkn xlxfxL0)()()(為插值基函數(shù)其中： ),1,0)(( nkxl k ??不同的 插值方法 有不同的 基函數(shù) ,不同的表示形式 插值型求積公式 有的近似作為被積函數(shù)用 ,)()( xfxL n?ba dxxf )( ?? ba n dxxL )( ? ??? bankkk dxxlxf0)()(? ???nkba kk dxxlxf0 )()(則，若記 ?? ba kk dxxlA )(?? ba dxxffI )(][ n0( ) = I ( 1 )nkkkA f x?? ? (1)式為數(shù)值求積公式 . Ak為求積系數(shù) , 且僅與積分區(qū)間和求積節(jié)點 xk 有關(guān) . 0 [ ] [ ] ( ) = [ ] ( 2 )nk k nkR f I f A f x I f I?? ? ??稱 為 求 積 余 項 。 ——— 當(dāng)然期望公式能對越多的被積函數(shù)精確成立，并與此作為 判斷求積公式“好”與“差”的一個標(biāo)準(zhǔn)。 結(jié)論： 含有 n+1個節(jié)點的插值型求積公式 的代數(shù)精度至少為 n. 例 1. 試確定下面積分公式中的參數(shù)使其代數(shù)精確度盡量高 . )()]()0([)]()0([2)()( 120 fIhffahhffhdxxffI h ???????? ??? h dxxI 0 0解 : 221hI ?]20[2 231 hahhI ???0)( xxf ?對于 hI ?1h??? h dxxI 0 11)( xxf ?對于 22h??? h dxxI 0 22)( xxf ?對于33h? 3)221( ha??1II ?令 121?a]30[2 2241 hahhI ????? h dxxI033)( xxf ?對于44h? 44h?]40[2 3251 hahhI ????? h dxxI044)( xxf ?對于55h?65h?3,2,1,0)()( 1 ?? jxIxI jj)()( 414 xIxI ?因此 所以該積分公式具有 3次代數(shù)精確度 ? 1． NewtonCotes公式 ? 2． 常用的 NC公式 ? 3. NewtonCotes公式的穩(wěn)定性 NewtonCotes求積公式 NewtonCotes數(shù)值求積公式 NewtonCotes公式是指等距節(jié)點下使用 Lagrange插值 多項式建立的數(shù)值求積公式 ],[)( baCxf ?設(shè)函數(shù)為插值多項式及余項分別的 L a g r a n g exf )(等份分割為將積分區(qū)間 nba ],[nkkhax k ,1,0, ????為步長其中 n abh ??各節(jié)點為 ???nkkkn xlxfxL0)()()( )()!1( )()( 1)1(xnfxR nnn ???? ??10( ) ( ) , [ , ]nniix x x a b????? ? ??其中 139。 所以，人們一般常用前兩個方法。 響的舍入誤差對公式的影只需討論 )( kxf)()()(,)( 計算值的近似值作為而以為精確值假設(shè) kkk xfxfxf為誤差)()( kkk xfxf ???)( fIn ????nkknk xfCab0)( )()(記 )( 計算值的近似值為 nI而理論值為 )( fIn ????nkknk xfCab0)( )()(的誤差為與 nn II)()( fIfI nn ? ?????nkkknk xfxfCab0)( )]()([)(nn II ? ????nkknkCab0)()( ?????nkknkCab0)()( ?nn II ?????nknkCab0)()( ? |}m a x {| k?? ?有若 ,0, )( ??? nkCnknn II ? ????nknkCab0)()( ??)( ab ??()01nnkkC???性 質(zhì) ：即 nn II ? ?)( ab ??NewtonCotes公式的舍入誤差只是函數(shù)值誤差的 倍)( ab ?時，公式都是穩(wěn)定的當(dāng)事實上 8, ?n公式是穩(wěn)定的時即 C o t e sN e w t o nCnk nk ???? ,0, )(???nknkCab0)()( ?有有正有負(fù)若 ,)( nkC????nknkCab0)()( ? ?)( ab ??此時 ,公式的穩(wěn)定性將無法保證 因此 ,在實際應(yīng)用中一般不使用高階 NewtonCotes公式 而是采用低階復(fù)合求積法 (下節(jié) ) 復(fù)化求積法 167。 為了提高公式的精度 ,又使算法簡單易行 ,往往使用復(fù)合方法 [ , ]ab即 將 積 分 區(qū) 間 分 成 若 干 個 子 區(qū) 間 ,然后在每個小區(qū)間上使用低階 NewtonCotes公式， 最后將每個小區(qū)間上的積分的近似值相加。 計量經(jīng)濟(jì)學(xué)模型的應(yīng)用 一、 結(jié)構(gòu)分析 二、 經(jīng)濟(jì)預(yù)測 三、 政策評價 四、 理論檢驗與發(fā)展 60 一、結(jié)構(gòu)分析 ? 經(jīng)濟(jì)學(xué)中的結(jié)構(gòu)分析是對經(jīng)濟(jì)現(xiàn)象中變量之間相互關(guān)系的研究。 ? 計量經(jīng)濟(jì)學(xué)模型的功能是揭示經(jīng)濟(jì)現(xiàn)象中變量之間的相互關(guān)系，即通過模型得到彈性、乘數(shù)等。 ? 計量經(jīng)濟(jì)學(xué)模型是以模擬歷史、從已經(jīng)發(fā)生的經(jīng)濟(jì)活動中找出變化規(guī)律為主要技術(shù)手段。 ? 模型理論方法的發(fā)展以適應(yīng)預(yù)測的需要。 ? 經(jīng)濟(jì)政策的不可試驗性。 64 四、理論檢驗與發(fā)展 ? 實踐是檢驗真理的唯一標(biāo)準(zhǔn)。 ? 計量經(jīng)濟(jì)學(xué)模型提供了一種檢驗經(jīng)濟(jì)