版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
1、數(shù)學(xué)軟件實(shí)驗(yàn)任務(wù)書課程名稱數(shù)學(xué)軟件實(shí)驗(yàn)班級(jí)數(shù) 0901實(shí)驗(yàn)課題Romberg積分法,Gauss型積分法,咼斯-勒讓德積分 法,咼斯-切比雪夫積分法,咼斯-拉蓋爾積分法,咼 斯-埃爾米特積分法實(shí)驗(yàn)?zāi)康氖煜omberg積分法,Gauss型積分法,咼斯-勒讓德 積分法,咼斯-切比雪夫積分法,咼斯-拉蓋爾積分法, 咼斯-埃爾米特積分法實(shí)驗(yàn)要求運(yùn)用 Matlab/C/C+/Java/Maple/Mathematica 等其中一種語言完成實(shí)驗(yàn)容Romberg積分法,Gauss型積分法,咼斯-勒讓德積分 法,咼斯-切比雪夫積分法,咼斯-拉蓋爾積分法,咼 斯-埃爾米特積分法成績教師實(shí)驗(yàn)IRomberg積分法
2、1實(shí)驗(yàn)原理Romberg方法是實(shí)用性很強(qiáng)的一種數(shù)值積分方法,其收斂速度是很快的,這里給出Rombergs分的計(jì)算方法。1(1) 計(jì)算 R(0,0)-(b a)f (a) f (b)2i 2(2) 計(jì)算 R(i,0)-(i 1,0) 山 f(a (k22 k 1知)(3)計(jì)算 R(m, j)4 (m, j 1) R(m 1, j 14j 112實(shí)驗(yàn)數(shù)據(jù)用Romberg積分方法計(jì)算:S3實(shí)驗(yàn)程序 程序1fun cti on s=rombg(a,b,TOL) n=1;h=b-a;delt=1;15 x0 4x2dxx=a;k=0;R=zeros(4,4);R(1,1)=h*(rombg_f(a)+r
3、ombg_f(b)/2; while delt>TOLk=k+1; h=h/2; s=0;for j=1: nx=a+h*(2*j-1); s=s+rombg_f(x);endR(k+1,1)= R(k, 1)/2+h*s; n=2* n;for i=1:kR(k+1,i+1)=(4Ai)*R(k+1,i)-R(k,i)/(4Ai-1); enddelt=abs(R(k+1,k)-R(k+1,k+1);ends=R(k+1,k+1);程序2fun cti on f=rombg_f(x) f=x/(4+xA2);程序3s=rombg(0,1.5,1.e-6) %乍出圖形x=0:0.02:1
4、.5;y=x./(4+x.A2); area(x,y) grid4實(shí)驗(yàn)結(jié)果s =0.22310.511.5實(shí)驗(yàn)2高斯-勒讓德積分法1實(shí)驗(yàn)原理Gauss-Legendre求積公式為(x)dxAJ(Xk)k 1其中xk為Legendre多項(xiàng)式在Legendre多項(xiàng)式定義為:1,1區(qū)間上的零點(diǎn)。n階i)npn(t)A為權(quán)系數(shù),A22(1 Xk)2Ak2222(1 Xk) Pn (Xn) n Pn(Xk)對于一般的積分區(qū)間為a, b問題,可以做變換abbaxt2 2f(x)dxa bAkf(-k 1Xk)2實(shí)驗(yàn)數(shù)據(jù)用Gauss-Legendre積分方法計(jì)算定積分S02 x2 cosxdx3實(shí)驗(yàn)程序fu
5、n cti on s=gau_leg(a,b)%5階Legendre多項(xiàng)式結(jié)點(diǎn)node=-0.9061798459,-0.5384693101,0,0.5384693101,0.90617 98459;% 吉點(diǎn)對應(yīng)的權(quán)quan=0.2369268851,0.4786286705,0.5688888889,0.4786286705,0.2369268851;%t為(1, 5)的行向量,整個(gè)區(qū)間上的結(jié)點(diǎn)t=(b+a)/2+(b-a)* no de/2;s=(b-a)/2)*sum(qua n.*gaueg_f(t);fun cti on f=gau_leg_f(x) f=(x.A2).*cos(x
6、);disp('計(jì)算結(jié)果為:') s=gau_leg(0,pi/2) %畫出圖形 x=0:0.01:pi/2; y=(x42).*cos(x); bar(x,y)grid4實(shí)驗(yàn)結(jié)果 計(jì)算結(jié)果為:s =0.4674實(shí)驗(yàn)3咼斯-拉蓋爾積分法 1實(shí)驗(yàn)原理n個(gè)結(jié)點(diǎn)Gauss-Laguerre求積公為:nSAk f (xQk 1其中Xk為零點(diǎn),Ak為權(quán)系數(shù)AkLaguerre多項(xiàng)式為Ln(X)x dn x、e 護(hù)xe ),02實(shí)驗(yàn)數(shù)據(jù)計(jì)算反常積分S xe Xdx03實(shí)驗(yàn)程序fun cti on s=gau_lag()%多項(xiàng)式結(jié)點(diǎn)node=0.26355990,1.41340290,3.5
7、9624600,7.08580990,12.640800;%權(quán)重向量quan=0.6790941054,1.638487956,2.769426772,4.31594400,7.10489623;%求和s=sum(qua n. *gau_lag_f( no de)%以下為畫出積分示意圖clearx=0:0.1:20;y=x.*exp(-x);area(x,y)grid fun cti on f=gau_lag_f(x) f=x.*exp(-x);4實(shí)驗(yàn)結(jié)果1.00000.402468101211 1E 1620實(shí)驗(yàn)4高斯-埃爾米特積分法1實(shí)驗(yàn)原理n個(gè)結(jié)點(diǎn)點(diǎn)Guass-Hermite求積公式為nSAk f ( xk )k 1其中Xk,Ak分別為結(jié)點(diǎn)以及相應(yīng)的權(quán)系數(shù)。2實(shí)驗(yàn)數(shù)據(jù)采用Gauss-Hermite方法計(jì)算反常積分S xe xdx3實(shí)驗(yàn)程序fun cti on s=gau_lag()%多項(xiàng)式結(jié)點(diǎn)node=-2.02018200-0.958571900.000000000.958571902.02018200;%權(quán)重向量quan=1.1814695990.98657914170.94530892370.98657914171.181469599;%求和s=sum(quan
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(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ǔ)空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 發(fā)貨倉庫合同模板
- 高校教師編制合同模板
- 原油進(jìn)口中文合同模板
- 裝修協(xié)會(huì)質(zhì)保合同模板
- 飯店簽合同模板
- 銅管買賣合同模板
- 賣方房子定金合同模板
- 瓷磚購銷付款合同模板
- 建設(shè)合伙合同模板
- 父母和孩子上學(xué)合同模板
- 2024年4月自考00322中國行政史試題及答案含評分標(biāo)準(zhǔn)
- 2022乳腺癌中西醫(yī)結(jié)合診療共識(shí)
- (正式版)SHT 3115-2024 石油化工管式爐輕質(zhì)澆注料襯里工程技術(shù)規(guī)范
- 武裝押運(yùn)人員崗前培訓(xùn)課件
- 幼兒英語演講與口才課程
- 腫瘤與冠心病 - 副本
- 高三攻略:家長與學(xué)生
- 2024年考研英語真題及答案(完整版)
- 食品質(zhì)量與安全專業(yè)職業(yè)生涯規(guī)劃
- 世界衛(wèi)生統(tǒng)計(jì)報(bào)告
- 值長口試試題
評論
0/150
提交評論