數(shù)值分析上機(jī)報(bào)告c

《數(shù)值分析上機(jī)報(bào)告c》由會(huì)員上傳分享，免費(fèi)在線(xiàn)閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

1、數(shù)值分析上機(jī)報(bào)告C數(shù)值分析實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)3數(shù)值積分實(shí)驗(yàn)?zāi)康耐ㄟ^(guò)本實(shí)驗(yàn)理解數(shù)值積分與微分的基本原理。掌握數(shù)值積分中常見(jiàn)的復(fù)合求積公式的編程實(shí)現(xiàn)。掌握龍貝格算法的基本思路和迭代步驟；培養(yǎng)編程與上機(jī)調(diào)試能力。算法描述龍貝格算法基本思路先算出他TO,從而計(jì)算TO以此類(lèi)推，在判斷

2、Tn-Tn-l

3、龍貝格算法計(jì)算步驟步驟1：輸入?yún)^(qū)間端點(diǎn)a,b,精度控制值e,循環(huán)次數(shù)M,定義函數(shù)f(x),取n?l,h?b?a;Rl,1?h?f(a)?f(b)?/2步驟2：fork?2toM9■{Rk,l??Rk?9for2k?2?h?kl,1toli?ll?f?a?2?i?k9?k?h/2??j2if

4、}?Rk,j?RR?klk,?jl??Re?k,?Rjl??kl,jl?/?4j?l??19■R?k,k,?kl?退岀循環(huán)步驟3:數(shù)據(jù)積分近似值Rk,ko利用Romberg方法計(jì)算函數(shù)I?910sinxxdx實(shí)驗(yàn)內(nèi)容用龍貝格算法計(jì)算:I?109sinxx實(shí)驗(yàn)步驟代碼#includezz//#include〃〃#defineedoublef(doublex){doubley;if(x==0){returny=;}elsey=sin(x)/x;returny;voidromberg(doublea,doubleb){}intn=l,k=0;doubleh,T2,S2=0,C

5、2=0,R2=0,Tl,Cl,SI,R1;h=(b-a)/2;T2二h*(f(a)+f(b));while(fabs((R2-R1))>e){R1=R2;T1=T2;S1=S2;C1=C2;doublesum=0;inti;for(i=l;isum=sum+f(a+(2*iT)*h);}T2二Tl/2+sum*h;S2=(4*T2-T1)/3;C2二(16*S2-SI)/15;R2=(64*C2-Cl)/63;n=n*2;k++;h=h/2;}coutvoidmain(){doublea,b;cout>a>>b;coutcoutcout}實(shí)驗(yàn)結(jié)果實(shí)驗(yàn)體會(huì)XT^XT^XT

6、^XT^XT%?TXXlXXlXXlXxjxXTXxrsXT^XT^XT^XT^zTsxTx?Tx?Tx?TxxTxxTs?Tx^TX^TX?TxXTXz7x?TxXjXxTxxTxXTX?Txxjxxjx?TxxTxxj%>1^>1^>1^>1^>1^>1^>1^>1^>1^>1^%fz%fz>1^%fzxjvZjxZjxzjxZjxzjx?Jx?Jx?Jx?Jx?Jx數(shù)值分析上機(jī)報(bào)告姓名：學(xué)號(hào)：專(zhuān)業(yè)：2013年10月27日弟一早舍入誤差與有效數(shù)設(shè)SN??j?2N311?o,其精確值為l?????j?12?2NN?l?2111,計(jì)算SN的通用程序。？???2222?1

7、3?1N?111,計(jì)算s的通用程序。編制按從小到大的順序SN?1????NN2?1(N?1)2?122?1編制按從大到小的順序SN?按兩種順序分別計(jì)算S102,S104,S106,并指出有效位數(shù)。通過(guò)本上機(jī)題，你明白了什么？解：、題程序見(jiàn)電子版按從大到小順序：按從小到大順序：S102二有效位數(shù)6位S104=有效位數(shù)3位S106=有效位數(shù)3位S102二有效位數(shù)5位S104二有效位數(shù)6位S106二有效位數(shù)6位(4)通過(guò)上述實(shí)驗(yàn)數(shù)據(jù)可以看出此次算法使用從小到大的順序進(jìn)行得到的數(shù)據(jù)相對(duì)而言更精確，可以得到這樣的啟示：在計(jì)算數(shù)值時(shí),要先分析不同算法對(duì)結(jié)果的影響,避免大數(shù)吃小數(shù)的現(xiàn)

8、象，找出能得到更精確的結(jié)果的算法。第一早Newton迭代法給定初值xO及容許誤差？，編制Newton法解方程f(x)?0根的通用程序。3給定方程f(x)?x/3?x?0,易知其有三個(gè)根x1?,x2?0,x3?9■9■9■1.由Newton方法的局部收斂性可知存在？?0,當(dāng)xO?(??,?)時(shí)，Newton迭代序列收斂于根x2。試確定盡可能大的?。9■2?試取若干初始值，觀察當(dāng)x0?(??,?l),(?1,??),(??,?),(?,1),(1,?)時(shí)Newton序列是否收斂以及收斂于哪一個(gè)根。通過(guò)本上機(jī)題，你明白了什么？解：OxO?(-?,?)時(shí)，迭代序列收斂于根x2*

9、9■由上表可以看出，xO在內(nèi)收斂于xl,在(??,?)內(nèi)收斂于x2,在內(nèi)收斂于x3,但在內(nèi)(？1,??)和(?，1)均可能收斂于xl和x3。在用迭代法進(jìn)行計(jì)算的時(shí)候，初值的選擇很關(guān)鍵，所選擇的初值必須在一定的范圍內(nèi)，如果超出這個(gè)范圍，將得不出正確的結(jié)果。9*?第二早39?列主元三角分解法對(duì)于某電路的分析，歸結(jié)為求解線(xiàn)性方程組RI=V,編制解n階線(xiàn)性方程組Ax=b的列主元三角分解法的通用程序；用所編制的程序解線(xiàn)性方程組RI=V,并打印出解向量,保留五位有效數(shù)；本編程之中，你提高了哪些編程能力？解：見(jiàn)電子版方程的解為：xl=-,x2=,x3=