版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、數(shù)學(xué)軟件實驗任務(wù)書課程名稱數(shù)學(xué)軟件實驗班級數(shù) 0901Romberg積分法, Gauss 型積分法,高斯 - 勒讓德積分實驗課題法,高斯 - 切比雪夫積分法,高斯 - 拉蓋爾積分法,高斯- 埃爾米特積分法熟悉 Romberg 積分法,Gauss型積分法,高斯 - 勒讓德實驗?zāi)康姆e分法,高斯 - 切比雪夫積分法, 高斯 - 拉蓋爾積分法,高斯 - 埃爾米特積分法實驗要求運用 Matlab/C/C+/Java/Maple/Mathematica 等其中一種語言完成Romberg積分法, Gauss 型積分法,高斯 - 勒讓德積分實驗容法,高斯 - 切比雪夫積分法,高斯- 拉蓋爾積分法,高斯- 埃爾
2、米特積分法成績教師實驗 1Romberg積分法1 實驗原理Romberg 方法是實用性很強(qiáng)的一種數(shù)值積分方法,其收斂速度是很快的,這里給出 Romberg積分的計算方法。( 1)計算 R(0,0)1(ba) f ( a)f (b)2( 2)計算1hi 12 i 212(i1,0)2f (a(k2)hi 1)R(i,0)k1( 3)計算 R( m, j )4j 1R(m, j1)R( m1, j1)4 j112 實驗數(shù)據(jù)用 Romberg 積分方法計算:15x2 dxS4x03 實驗程序程序 1function s=rombg(a,b,TOL)n=1;h=b-a;delt=1;x=a;k=0;R
3、=zeros(4,4);R(1,1)=h*(rombg_f(a)+rombg_f(b)/2;while deltTOLk=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)=(4i)*R(k+1,i)-R(k,i)/(4i-1);enddelt=abs(R(k+1,k)-R(k+1,k+1);ends=R(k+1,k+1);程序 2function f=rombg_f(x)f=x/(4+x2);程序 3s=rombg(0,1.5,1
4、.e-6)%作出圖形x=0:0.02:1.5;y=x./(4+x.2);area(x,y)grid4 實驗結(jié)果s =0.2231實驗 2 高斯 - 勒讓德積分法1 實驗原理Gauss-Legendre 求積公式為1nf ( x)dxAk f ( xk )11k其中 xk 為 Legendre 多項式在1,1 區(qū)間上的零點。n 階Legendre 多項式定義為:Pn (t )1d n2nnn! dtn ( t1) 2Ak 為權(quán)系數(shù),Ak22(1 xk )22 Pn ( xn )2 n2 Pn ( xk )2(1 xk )對于一般的積分區(qū)間為a, b問題,可以做變換xabb a2t2bba na
5、bb aaf ( x)dxAk f (xk )2k1222 實驗數(shù)據(jù)用 Gauss-Legendre 積分方法計算定積分S2 x2 cosxdx03 實驗程序function s=gau_leg(a,b)%5階 Legendre 多項式結(jié)點node=-0.9061798459,-0.5384693101,0,0.5384693101,0.9061798459;%結(jié)點對應(yīng)的權(quán)quan=0.2369268851,0.4786286705,0.5688888889,0.478628670 5,0.2369268851;%t 為( 1, 5)的行向量,整個區(qū)間上的結(jié)點 t=(b+a)/2+(b-a)*
6、node/2; s=(b-a)/2)*sum(quan.*gau_leg_f(t);function f=gau_leg_f(x)f=(x.2).*cos(x);disp(計算結(jié)果為 :)s=gau_leg(0,pi/2)%畫出圖形x=0:0.01:pi/2;y=(x.2).*cos(x);bar(x,y)grid4 實驗結(jié)果計算結(jié)果為 : s =0.4674實驗 3 高斯 - 拉蓋爾積分法1 實驗原理n 個結(jié)點 Gauss-Laguerre 求積公為:nSAk f ( xk )k 1其中 xk 為零點, Ak 為權(quán)系數(shù)Akxk 2 Ln 1 ( xk ) 2(n1)Laguerre 多項式為
7、Ln (x) ex dnn ( xn e x ),0 xdx2 實驗數(shù)據(jù)計算反常積分Sxe xdx03 實驗程序function s=gau_lag()%多項式結(jié)點node=0.26355990,1.41340290,3.59624600,7.08580990,12.640800;%權(quán)重向量quan=0.6790941054,1.638487956,2.769426772,4.31594400,7.10489623;%求和s=sum(quan.*gau_lag_f(node)%以下為畫出積分示意圖clearx=0:0.1:20;y=x.*exp(-x);area(x,y)gridfunctio
8、n f=gau_lag_f(x)f=x.*exp(-x);4 實驗結(jié)果s =1.0000實驗 4 高斯 - 埃爾米特積分法1 實驗原理n 個結(jié)點點 Guass-Hermite 求積公式為nkkk 1其中 xk , Ak 分別為結(jié)點以及相應(yīng)的權(quán)系數(shù)。2 實驗數(shù)據(jù)采用 Gauss-Hermite 方法計算反常積分Sxe xdx3 實驗程序function s=gau_lag()%多項式結(jié)點node=-2.02018200-0.958571900.000000000.958571902.02018200;%權(quán)重向量quan=1.1814695990.98657914170.94530892370.98657914171.181469599;%求和s=su
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 順豐同城新人考試
- 業(yè)財一體化智慧樹知到答案2024年青島理工大學(xué)
- 30高中數(shù)學(xué)新教材課堂導(dǎo)學(xué)案(復(fù)合函數(shù)的導(dǎo)數(shù)及函數(shù)圖像的切線)及答案
- 第一章 有理數(shù) 復(fù)習(xí)課 教學(xué)設(shè)計-2024-2025學(xué)年青島版數(shù)學(xué)七年級上冊
- 第十一章三角形專題練習(xí)2024-2025學(xué)年人教版數(shù)學(xué)八年級上冊
- 分散型控制系統(tǒng)(DCS)相關(guān)行業(yè)投資方案范本
- 專題01 全等圖形與全等三角形壓軸題五種模型全攻略(原卷版)
- 金屬礦區(qū)生態(tài)修復(fù)與污染控制一體化
- 納米材料增強(qiáng)電熱封裝
- 云遷移策略與最佳實踐
- 第三單元名著導(dǎo)讀《紅星照耀中國》(公開課一等獎創(chuàng)新教學(xué)設(shè)計+說課稿)
- 數(shù)據(jù)中心儲能應(yīng)用需求技術(shù)報告2024
- 2024年上半年教師資格證《高中數(shù)學(xué)》真題及答案
- 2024年統(tǒng)編版新教材語文小學(xué)一年級上冊第二單元測試題(有答案)
- 風(fēng)光制氫一體化項目可行性研究報告
- 第一章 有理數(shù)單元測試卷2024-2025學(xué)年人教版(2024)七年級數(shù)學(xué)上冊
- 輸變電工程施工質(zhì)量驗收統(tǒng)一表式附件1:線路工程填寫示例
- 2024年云南楚雄州楚雄市事業(yè)單位選調(diào)30人歷年(高頻重點提升專題訓(xùn)練)共500題附帶答案詳解
- 2024年紅十字應(yīng)急救護(hù)知識競賽考試題庫500題(含答案)
- 剪刀式升降車的安全管理
- 簡易用工協(xié)議書范本
評論
0/150
提交評論