計(jì)算方法公式總結(jié)_第1頁
計(jì)算方法公式總結(jié)_第2頁
計(jì)算方法公式總結(jié)_第3頁
計(jì)算方法公式總結(jié)_第4頁
計(jì)算方法公式總結(jié)_第5頁
已閱讀5頁,還剩34頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、x為準(zhǔn)確值,x為近似值。,e為正數(shù),稱為絕對(duì)誤差限F x x e通常用e 表示相對(duì)誤差x x計(jì)算方法公式總結(jié)緒論絕對(duì)誤差 ex x,絕對(duì)誤差限|e| |x x|r x x e相對(duì)誤差e - 一 x x相對(duì)誤差限|er |或|er|Word 資料有效數(shù)字一元函數(shù)y=f (x) '絕對(duì)誤差e(y) f (x)e(x)''e(y) f (x) e(x) xf (x)相對(duì)誤差er(y) y- 74(幻二元函數(shù) y=f (x1,x2)/ f (Xi,X2) f (Xi,X2) 絕對(duì)誤差 e(y)dX11dx2相對(duì)誤差er(y)f (Xi,X2)X1 / 、f (X1,X2)X2

2、/ 、- -er(Xi)- -er(X2)XiyX2ye(xi + x2) = (町)+ e(x2): 巴3 -=eJCi)-4, 巴(£"2)七 ! " 1 I + j;ie(x2)tF (巴)=- ?收工2), 5/ T2 喘S(E1 + J:2)七一H勺(12),1+£211 +2巴丁(血 -X2)右一外(11)%(12),X 一七2工1 一12立(工"2)七今(11)+外(6人機(jī)器數(shù)系1)數(shù)的浮點(diǎn)表小:/ = i(0*Qir>2JL 中三 Qi < /? (1 < i < n), L < p < U.

3、F(U.U) =。 U 12)機(jī)器中的數(shù)集是有限的:石=±(0.修如。計(jì))用,1 W因W 3 L 10 W 出 0 8 - L - 2. 3. ,% A S P 0 U J全人況差.3)實(shí)數(shù)£二 3)(機(jī)器數(shù)).注:1. B冷,且通常取2、4、6、82 . n為計(jì)算機(jī)字長(zhǎng)3 .指數(shù)p稱為階碼(指數(shù)),有固定上下限L、U4 .尾數(shù)部s0.a1a2L an,定位部5 .機(jī)器數(shù)個(gè)數(shù)1 2(1) n1(U L 1)機(jī)器數(shù)誤差限,1 n p截?cái)嘞鄬?duì)|x fl(x)l|x|舍入絕對(duì)|x fl(x)| - 截?cái)嘟^對(duì)| x fl (x) |x fl(x)|舍入相對(duì).一.|x|積分中值定理設(shè)

4、f屈iMhb上連紈g(幻在a可上保號(hào)(即II負(fù)或II I:). 在£ 舊陽,使得j f 8gs d舉=A0,(丁 jdw.JoJ a秦九韶算法方程求根f (x) (x x)mg(x), g(x) 0, x* 為 f (x) =0 的 m 重根。二分法一以| < «1) < 即7。-。).對(duì)于給定精度E,若取k使得則有迭代法f (x)0xk 1(xk)k=0、1、2*xk為迭代序列,(x)為迭代函數(shù),lim xk x (x )k定理1流尸.一在Lb內(nèi)存在 階連續(xù)導(dǎo)數(shù),且滿足事/) "J W 時(shí),E a. b ;2;存在止常數(shù)L < L當(dāng)工丘a, b

5、時(shí),爐(0| < L< L則1)工=奴工)在區(qū)b上有唯一實(shí)根,記為了*;2)時(shí)任意初值叼e 同迭代格式(5)收斂,且Um?。?)L«|工* 以I W r|i -網(wǎng)|: A 1,2.3 - - ;(7)1 _ L5)M -T)(8)Jr ;r* 八局部收斂定義1對(duì)方程工=/公在小的某個(gè)鄰域£ = 司工一"* < 匹內(nèi),對(duì)任意初值w S迭代格式以十1 =1(以)都收斂,則稱 迭代法在二的附近局部收斂.定理3法方程I =火)有根廠H一在工*的某個(gè)笫域S = xx - x*| < 6內(nèi)品工)一階連續(xù)可導(dǎo),則1)當(dāng)|")| < 1時(shí),

6、迭代格式局部收斂;2)當(dāng)|“(力| > 1時(shí),迭代格式發(fā)散.注:如果知道近似值,可以用近似值代替根應(yīng)用定理 3判斷是否局部收斂定義2設(shè)序列琛收斂于邕 并記年=xt - x如果存在常 數(shù)P之1及非零常數(shù)C使得limk boc”+1則稱序列磔是p階收斂的.定理4花奴在H附近的某個(gè)鄰域內(nèi)仃p(之1)階連續(xù)導(dǎo)數(shù),且114)(15),好(H J = 01卜=1,2p 1, d心)等o,則迭代格式住/附近是p階局部收斂的,11有X-A+1 - X* lim Atx(I* I*),0叫力p!(16)如果P= L要求牛頓迭代法f(x) f(Xk) f(Xk)(X Xk) 0xk 1xkf()(k 0,1

7、,2,L ) f (Xk)注:牛頓迭代對(duì)單根重根均局部收斂,只要初值足夠靠近真值。.當(dāng)二為方程單根時(shí),Newton迭代是二階局部收斂."舊為方程.皿m > 2;重根時(shí),Newlon迭代一階局部收斂.牛頓迭代法對(duì)初值要求很高,要保證初值在較大范圍內(nèi)也收斂,加如下四個(gè)條件定理5設(shè)函數(shù)/(©在區(qū)間&目?jī)?nèi)2階連續(xù)可導(dǎo) 且滿足:1)< 0;2)當(dāng)z E a山時(shí),開工)¥ 0;/而3)當(dāng)工丘45)時(shí)尸(切保號(hào);/仙)、>?一則對(duì)V$o G迭代格式/(必)i n 1隊(duì)+1 =4377= 0, 1, * * *收斂到方程穴r) = 0在6內(nèi)的唯一實(shí)根注:

8、證明牛頓迭代法大范圍收斂性,要構(gòu)造一個(gè)區(qū)間£,M( £),其中M()ff;,在這個(gè)區(qū)間內(nèi)驗(yàn)證這四個(gè)條件。如果知道根的位置,構(gòu)造e,M (e)時(shí)應(yīng)該包括根,即e+常數(shù)3.4割線法/(以) /£店十I =龍為-4不一石以-人一卜=1":,,jxk) - fxk-l)線性方程組求解有兩種方法:消去法 和迭代法高斯消去法利用線性代數(shù)中初等行變換將 增廣矩陣轉(zhuǎn)化為等價(jià)上三角矩陣。注意:第一行第一列為0,將第一列不為0的某一行與第一行交換位置,繼續(xù)初 等行變換LalnLa2nMLann對(duì)角占優(yōu)矩陣aii42a21a22AM Man1an2n|aij |(i1,2,L

9、 , n)則稱a為按行嚴(yán)格對(duì)角占優(yōu)矩陣j 1n|aij |(j1,2,L ,n)則稱a為按列嚴(yán)格對(duì)角占優(yōu)矩陣i 1au aji (i 1, j n) x Rn,x0,(x, Ax) 0則稱A是對(duì)稱正定的當(dāng)A是上面三種情況時(shí),用高斯消去法消元時(shí)akk0 ,不用換行。追趕法是高斯消元法的一種特例考慮二對(duì)角方桿組Ui其系數(shù)矩陣元素滿足1) |6i| > ci > 0;2) bi > |aj + |q|, g #0(,= 2,3,中1);3) bn > | > 0.易證(4)的系數(shù)矩陣指奇異.利用Gmi,消去法解方程組(4),每步 消元只要消一個(gè)元素.消元過程算法如下:B

10、i =瓦, 01 =九 對(duì)晨一 2.:1, * * ”1做li = 3-,氏=b$ - Lg-lA-i今回代算法:bi ci is b >w«381 5 0 02 0比 n = l/n / Br,對(duì),=71 17 71 2. * ,.,1 做麻=(柒-Q1+1)/片.列主元高斯消元法當(dāng)1a(k)sk即第k次消元把kn行第k列絕對(duì)值最設(shè)進(jìn)行1步消元c(1)J1)fl)r的1aV2口13"1限* "1"<f Ln+ 1U12)小冷a23電小-1門4k口上找1 L00'(3) “3上-1* *的上-(3)。而»“出翼+1«

11、;1 1U4<1)* «/】)/l)(島)1 iak-,k"-1最0000(A)的臉*)盯+1制IJ000*) S-+1, A-件)'" + 1,1 田* -4* *-*4 -»*«曾 鼻 *»9 -U000, a.1,71»的% "十10V004*9V0,(A- H心* 44渴* * , ann%+1i (k) imax©k |,大的行(s行)調(diào)到第k行,再進(jìn)行高斯消元2.1向量范數(shù)設(shè) 1 =(£1,石r)T G R ".定義1設(shè)/=|到是R"上的函數(shù),如果滿

12、足;非負(fù)性)Vh G R?,有|同| > 0,且|同| = 00宓=0;2"齊次性)Vi G R", AeR W | | = |N|磯3三角不等式)Mr, y G RV有|"十M| < h| + h|.則稱|為R"上的范數(shù).常用的三個(gè)向量范數(shù):Di范數(shù):|同h=£|闖;2 = 12)00范數(shù):I罔|g = max |/l<i<n3) 2范數(shù):|訓(xùn)2 =定理2設(shè)4 6 Rnxn,則1)Mlh=摩 |4班=lklli=1max1<.?<«nj=l2)Mb=躁|4唯=明£|%|.k II oo=

13、i,THII2 = mac 114M2 =I 112=1定義5設(shè)B Rf 入 小汴為B的辦個(gè)特征值.稱。但)=眄|%1l<i<n為衽陣B的滔半徑.定理3設(shè)| |是R"答上的任意一個(gè)矩陣范數(shù).A w IV.則有。 < MH定理4如果4 e R"印為對(duì)稱矩陣,則p(A) = |Al定義3設(shè)|是RR中的 個(gè)范數(shù),hS工.工,是R療中的 一個(gè)向量件列“E H"為常I可量,如果liin xk' c| 0, k ->0C則稱向量序列工0收斂于c記為lim a?"c.迭代序列構(gòu)造Ax b x Bx fx(k 1) Bx(k) f第三個(gè)等

14、式為迭代序列,B為迭代矩陣 迭代收斂判別1 .充分條件:迭代矩陣范數(shù)小于1, PBP 1結(jié)論:Ax=b有唯一解x對(duì)任意初始向量:E迭代格乃收斂4且有II”) c網(wǎng) |- 1 -1四I ".JA-1)JU12 3-IIII” " k2 .充要條件:迭代矩陣譜半徑小干 1, (B) 1Jacobi 迭代法A L D U其中L (low)為下三角,U為上三角,D為對(duì)角線元 素(k 1)1(k)1.迭代格式:x D (L U)x D b上,+ " =(_ a.*陽一。1313 一. 一 以舊七9)/。11K"" = (& -一 023% 一以加

15、或)/。22=('-。31式;?川一。32 叼一)/a33= (% 。北- 。汽萬-1 汽)3)/。加i迭代矩陣J D (L U)收斂性判據(jù):| I J | 0 |D 1 |?|L D U | 0 |L D U | 0求出 最大值小于1 (J的譜半徑小于1)即迭代格式收斂.Gauss-Seidel 迭代法迭代格式x(k1) D1( Lx(k1) Ux(k) b)x(k。(D l_) 1 Ux(k) (D l_) 1b1 :迭代矩陣:G (D L) U1.常數(shù)矩陣:g (D L) b工” =(61 -即濯-勾3追)如飛)/以1if + n = 口力工十口 儀2M:)法工華"&#

16、169;2二(加一 口同4&*"-"駝芯"1 -"時(shí)匚%)W,)*(ft+1)( I(卜+1)3+1)(*+】)、/琉 =("-口/1叼-通網(wǎng)2%用一 11計(jì)1 )/口.收斂性判據(jù):| I G | 0 |(D L) 1 |?| (D L) U | 0 | (D L) U | 0求出最大值小于1 (G的譜半徑小于1)即迭代格式收斂.結(jié)論:當(dāng)A是嚴(yán)格對(duì)角占優(yōu)的,則Jacobi和Gauss-Seidal迭代法均是文斂的宙值法用插值多項(xiàng)式p (x)代替被插函數(shù)f(x)插值多項(xiàng)式:P(X) 3o 3iX LanXn ,n+1 個(gè)點(diǎn) P(x。 y(

17、i 0: n)插值區(qū)間:a,b,插值點(diǎn)滿足a x0 x1 Lxn b求插值多項(xiàng)式P (x),即求多項(xiàng)式系數(shù)的過程為插值法帶入可知求系數(shù)的插值點(diǎn)行列式為范德蒙行列式,不為0,有唯一解。即n+1插值條件對(duì)應(yīng)的不超過 n次的插值函數(shù) P (x)只有一個(gè)。次線性插_ x x x xnP(x) y。y1 y°l°(x) y(x)x。 x1x x0lk(x)LagrangeLn(X)插值余項(xiàng)n(x xi)i 0i kn0(xk X)k插值多項(xiàng)式非插值節(jié)點(diǎn)上Q(x)(x Xi)k(xk xi)y*(x)k 0Lagrange 插值多項(xiàng)式為被插函數(shù)f(x) Ln(x)(a,b)帶導(dǎo)數(shù)插值條

18、件的余項(xiàng)估計(jì)(I)求一個(gè)二次多項(xiàng)式“(上)使褥它的余頊為- HU)=公廠 L求-個(gè)三次多項(xiàng)式H(*)使得(. i ikxkx4 )ykf(x)的近似值f (n 1)()(n 1)!n(xi 0xixi)H-,H"S(3.” 3 s它的余項(xiàng)為/(上)-H (工)(工-6), E W (口,/定理2沒*)在也可上連續(xù),尸計(jì)1)(工)在(&內(nèi)存在,如皿 *為互異節(jié)點(diǎn),££/)是滿足的插值多項(xiàng)式,則對(duì)W 凡6. 3 e (a, b) (e依賴于工),使得心=/(0 一 Ln(x) -3門+ !),(10)其中電計(jì)1(£)= J(£i=0注:推導(dǎo)

19、過程用羅爾中值定理構(gòu)造輔助函數(shù)(t) Rn(t) K(X)Wni(t)注2,乂依敕于二后即£ = i上仆,:1,一, , 1", IIKLX2 0. £代 L2)當(dāng)f(工)本身是一個(gè)次數(shù)不超過m的多項(xiàng)式吐f(T)4(工)= 0,因而乙Q) = /").特別當(dāng)/(£) = L則有£43 = 13)由于£一般不能牯確求出.因此只能估計(jì)誤差,沒 max |-+D | = A/n+1,則行IE(功二告嚴(yán)(功i -L J 第二條性質(zhì)用于可以證明階數(shù)不大于n的f(x)的插值余項(xiàng)為0.差商和Newton插值法定義2設(shè)知函數(shù)人£在

20、門+ 1個(gè)互異節(jié)點(diǎn)以也山的函數(shù)俏為力Mo),一,/(廝).稱/(勺)-/(石)/心叼=町雹為(力關(guān)于節(jié)點(diǎn)“:了的/階年兩f均差人稱/價(jià)差商/,勺和f 町叫的差商f4以=/勺以麻為1位)關(guān)于節(jié)點(diǎn)上的2階差商,一般電稱2個(gè)上一 1階的 差商為人階差商,即/10,陽,* ,,須-一1;1#_ / ,1,工2廠上*一 1 :一 /工5 4:1 : * * * : 42,1工k 工仆約定。階差商是函數(shù)值.計(jì)算函數(shù)的差商,可以通過列表法計(jì)算。記憶方法:先記分母,最后一個(gè)減去第一個(gè),對(duì)應(yīng)的分子第一項(xiàng)是最后一個(gè)臨近k元素的差商,第二項(xiàng)是第一個(gè)臨近 k個(gè)元素的差商。性質(zhì)3尼階差商和比階導(dǎo)數(shù)之間有如下關(guān)系:聲 1

21、 /其中7; E(1疝1H0,111一 4,max",勺,一以)牛頓插值多項(xiàng)式&=/30)+ /住5 詞(,-工0)+叼42(4-工1)+ * * . + /,0,1 /,%(遼一比0)(工一£)*(# n-1)*( I 5)(15)稱為幾次New ton插值多項(xiàng)式.通常記作Nn(X) 分段樣條插值給定在門+ 1個(gè)節(jié)點(diǎn)口 =如 < 也4=6上的函數(shù) 值:£于0_1 _ . . 于T>一于3 MH /山)小二)f心記兒=修+i H h = max hi,在每個(gè)小區(qū)間工七+i上作八力)的線性插位"-n 1£耳(1) = /(&#

22、169;) + f(/ - &), r E 出心+4其誤差為/-=""(S)(工一用(工一f 卜) 宣 £ (陽,41)從而有5f "(&)(£ - M(£ - g*i)max4 <T<i + i£"幻| <max應(yīng) Wt+i-hH ?屢,"31,分段二次樣條插值 討論n為奇偶情況時(shí)的三個(gè)點(diǎn)余項(xiàng)估計(jì)式三次樣條插值函數(shù)第一類邊界條件(端點(diǎn)一階導(dǎo)數(shù)已知)ZMq 十%=n.xj-琬,Mn_2Mn =-八工冕") 1TD0等于第一個(gè)式子,dn等于第二個(gè)式子自然邊界條件(

23、端點(diǎn)二階導(dǎo)數(shù)已知 二階導(dǎo)數(shù)和M0 , Mn=0 )/+%=16()"一力必一萬7、也+& k % )=6/%_力,毛川S(e) = !fj + /町,叼+1(必 + g%+1)%(1 巧)+$必(* 叼產(chǎn)+(必+1 必)("一叼)12ont e 叼,叼+ 1, j = 0,1,,泣一L (26)曲線擬合最小二乘原理函數(shù)關(guān)于n個(gè)點(diǎn)線性無關(guān)定義1設(shè)項(xiàng)物戶口是幾個(gè)互不相同幽點(diǎn),.壯0(工),。1(萬人,。是(m4 1)個(gè)已知的函數(shù).如果存在不全為 零的常勤加匕I:,為使力次(勺)+勺腐(叼)+ 卜編 =。,=12 ,小1)則稱強(qiáng)3M (頊岫凡(工*關(guān)廠點(diǎn)%孫嚴(yán)力是線性無關(guān)

24、的; 否則稱為線性相關(guān)的.機(jī)(1)次(11)血(打一 1)方程組(D可以寫為陶1-1!1110m (1)篇1(;誦)狐6) L:* * *腳北1 (中比一 1)4/l(刀一11* 0m l(£n)0/7?(77)方程組(2)只有零解的充分必要條件是方數(shù)矩陣列線性無關(guān).定義2給定函數(shù)y = /上的“個(gè)點(diǎn)(孫為)>,= m.|11 3叼工口ylyi英成K設(shè)仰M(叫,鼠(沖關(guān)廣點(diǎn)% *線性無關(guān).令rn?=')用僅0.心:心)二工(心卜)一機(jī)汽Jt=l求% %Gn,使得力(勾,門, .6)= inin 中(小 的? 一 . tzj. (3) 心如號(hào)eER稱 mP(工)=

25、3;埼&(工)f=0為數(shù)據(jù)的擬合函數(shù),如果如.=xkr則和為加次最小二乘卷頂式.23 n注:線性無關(guān)的函數(shù)為1,x, x ,x ,L ,x才是最小二乘多項(xiàng)式(00,00)他5血)(0(h 以“: : * B : 0 _ (。如 00)(。的曲), M'(y,。l ) 睢(y<M1r瓦2 21- 7« d> aattti(5)匕m _其中加 n,系數(shù)矩陣4的列向量線性無關(guān),方程(5)稱為超定方 程組.該方程組般沒仃精確解.記m / n4*(X1, X2f ,&)-E E*=1 J=1QjjXj bf求工工紇,:琮,使得噸攻")二min 4(

26、丁卜12,.工口).由駐點(diǎn)方程組的理論可知,若產(chǎn)力一,琮是卜.面方程組的解: AtAx = A1 b注:記住公式即可。數(shù)值積分和數(shù)值微分性質(zhì)7 (定積分中值定理)如果函數(shù)人工)在閉區(qū)間加上連續(xù), 則在積分區(qū)間14可上至少存在一個(gè)席, 使ff3 444方)蔗而值公式般的數(shù)值積分公式為:itf擊七£ 4跺), A-0xk為求積節(jié)點(diǎn), Ak 為求積系數(shù)。插值求積公式baf(3;)dx 七 / Lu(j:)dj:J Ct.k=o n=£>wm).A'=0a + th, t 0. nl,a + j兒石k = Q+kli定義2如果求積點(diǎn)“體=ojm)是等距的 即b - q

27、七卜=a+kh. h =A,= 0,1, * * * , zi,n則稱對(duì)應(yīng)的插值型求積公式為Nch付-C”處公式.n * fc!(ATI 也 TT(i j)dt,卜=0,1(3) 周! Jo ;=0?。〤%kf(£k, kO梯形公式1 ) 71 = 1, li=bay j;()工1 = bl由(3)可以求彳¥GJ得2個(gè)等距節(jié)點(diǎn)的插值型求枳公式: JT(f)=二 /(回 + 則.(明稱為梯形公式.Simpson 公式2) :工一 2,/丁 一 -;工口 一 口,1 I - h : 11 (3)目 木丁。2.0 191361 =不。22 =:得3個(gè)等距節(jié)點(diǎn)的插伯邛求枳公式636

28、I-jS=+) + /(“ .,5)U/B-S(5)稱為Simpson公式.Cotes 公式r、 一 b a3a + 6«+ fej) n = 4, h = *j;o = a. Xi =沙=4-_門h、由求儲(chǔ)329():32i爪得5個(gè)等距節(jié)點(diǎn)的抽值型求積公式一 b a r“:缶一方 ,口,十 b。 -7仙)+32/1>+ 12/( ?。?JIJ Lnt乙n + 3Z>+324彳當(dāng)+7”.(6)4(6)稱為Cbte醇公式.截?cái)嗾`差/?(/) = /(/) -稱為求積公式的截?cái)嗾`差.Zb" ,b以工法一£ gm一k=0 山I f I / L門(工)di;J

29、qJ afb-廣品(W&產(chǎn)"1伙)5+1)!口(工一 xdi, £二0£三(匹b).代數(shù)精度定義3紿定一個(gè)求枳分/(/)= r/(£)曲;的求積公式/七乙= £4/(以,A'=0當(dāng)f(x)為不超過m次多項(xiàng)式時(shí)上式成立,f (x)為m+1多項(xiàng)式時(shí)上式不成立。則稱為求積 公式有m次代數(shù)精度。由插值型求積公式的戕斷誤差(2)知.n + 1個(gè)節(jié)點(diǎn)的插值型求 枳公式的代數(shù)精度至少是m定理1求枳公式naj /) - £jt0至少具有八次代數(shù)粘度o該求積公式是插值型求積公工,即月卜=/上.,k = 0 J Ju定理2求枳公式n(/)

30、為/£/) £/(以)(8)Ar=O的代數(shù)粘度是它對(duì)g(工)=1必(1)=叫g(shù)a(£)= n,Sm(x)=工精確成立,而對(duì)外(工)=工職+1不精確成立.即%)-La), , fc - U J - -# A(.7m+1)-梯形公式代數(shù)精度為 1, Simpson公式代數(shù)精度為3, Cotes公式代數(shù)精度為 5一般,幾十1個(gè)節(jié)點(diǎn)的Newton-Cotes(等距節(jié)點(diǎn)插值型)公式的代 數(shù)精度_(畫 門是奇數(shù)二,一L/是偶數(shù)截?cái)嗾`差梯形公式Rr(f) = 1(f)- r(/)Simpson 公式用(/) =/(/) s(力r廣f w a=I ;X a)x - bjdx =儂

31、7().心=,9 f 、V G (a,b).rbrb=/(£)/£ / Hxdx=- /嚴(yán)(E)( #4)( /:4.J ab a (b _180 (Cotes 公式Rc(f) = A/) - c-Gauss求積公式求積公式代數(shù)精度為 2n+1-1,1上的兩點(diǎn)Gauss公式111f 3dxf(工卜1,1上的三/JGauss公式JaJ aH (.r)dj:,/ a + b2£ 一 )(£ ) x b)dxf. ( a + b f. Tx a i i 1 jj b)dx2 ) #書(力,£(a,b).216 a (b aG ,,945 4 ) f)

32、,刈(3次代數(shù)精度) 1)f(T3)(5次代數(shù)精度)f(x)dx 9f(4)8f lf(3)3.2 區(qū)間q.句上的Gauss公式 考慮區(qū)間同可上的積分,(/)= f也自可得作變換I =/=J 1b a/a + b b a 丁 f +1 dt.由-:1上的Gauss公式得®可上的Gauss公式“Jw)= Ea/k-O記住xktk, A A的關(guān)系,tk A查表即可29Amt4T k = 0,1,2,.匹則儲(chǔ)小 "上的Gjusb積分公式為n乙(/)= £/(八)A,-0定義4設(shè)有計(jì)算積分的復(fù)化求積公式八(力,如果存在正整 數(shù)p和非零常數(shù)G使Uni:- C,占hr則稱公式

33、4(,)是p階的.復(fù)化梯形公式 2 階,復(fù)化 Simpson 公式 4 階,復(fù)化 Cote 公式 6 階計(jì)算機(jī)通過不斷把區(qū)間二分,所得前后兩次積分差值滿足精度條件即可 給定精度21112nIn1時(shí),1,II(f) 12n(f)| 2II2n(f) In(f)|因而可以取I2n(f )為I (f)的近似值。is 乙冗;%(/) TM梯形3Simpson/-即定3匹-S“(川15(向前差商)(向后差商)數(shù)值微分f Q);>flMt 1 打必)一 I% 一”)J 3。) ,h(外)°數(shù)值微分截?cái)嗾`差rf( + ) f (珀h ff Mb= 一$/ (/o) + Oh ),£f1/(比。)一/(I。一 ) h . z-)/l2f M;= -/ (No)十 0(h)t f(j:o + h) - f(x() h) A2 ff/q尸(皿)- = -ffM +。(h3), 2fi0中點(diǎn)公式:D(h)f (Xo h) f (Xo h

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論