變步長復(fù)化辛普森公式計算積分.doc

《變步長復(fù)化辛普森公式計算積分.doc》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、2.編寫用變步長復(fù)化辛普森公式計算積分的程序。用上面編寫的程序計算下列積分并分析計算結(jié)果（1）（2）（3）程序：functionS=bianfuhuasimpson(fx,a,b,eps,M)%變步長復(fù)合simpson求積公式%調(diào)用方式：S=fuhuasimpson(@fx,a,b,epsilon)%fx--求積函數(shù)（函數(shù)文件）%a,b--求積區(qū)間%eps--計算精度%M--最大允許輸出劃分?jǐn)?shù)n=1;h=(b-a)/n;T1=h*(feval(fx,a)-feval(fx,b))/2;Hn=h*feval(fx,(a+b)/

2、2);S1=(T1+2*Hn)/3;n=2*n;%最好與倒數(shù)第三行保持一致（變步長）whilen<=MT2=(T1+Hn)/2;Hn=0;h=(b-a)/n;forj=1:nx(j)=a+(j-1/2)*h;y(j)=feval(fx,x(j));Hn=Hn+y(j);endHn=h*Hn;S2=(T2+2*Hn)/3;fprintf('n=%2dS2=%-12.9fS2-S1=%-12.9f',n,S2,abs(S2-S1));ifabs(S2-S1)

3、dendS=S2;%達(dá)到下列條件之一，則運(yùn)算終止:%(1).abs(S2-S1)M%輸入1：S=bianfuhuasimpson(inline('sqrt(x)*cos(x)'),0,pi,10e-6,2000)%輸入2：S=bianfuhuasimpson(inline('2*x^2*cos(x^2)'),0,sqrt(pi),10e-6,2000)%輸入3：S=bianfuhuasimpson(inline('sqrt(x)'),0,1,10e-6,2000)輸出結(jié)果：（1）S=bianfu

4、huasimpson(inline('sqrt(x)*cos(x)'),0,pi,10e-6,2000)n=2S2=-0.016369112S2-S1=0.944423778n=4S2=-0.450266122S2-S1=0.433897010n=8S2=-0.669839370S2-S1=0.219573248n=16S2=-0.781318443S2-S1=0.111479074n=32S2=-0.837710689S2-S1=0.056392245n=64S2=-0.866141900S2-S1=0.028431211

5、n=128S2=-0.880440980S2-S1=0.014299080n=256S2=-0.887620063S2-S1=0.007179083n=512S2=-0.891220052S2-S1=0.003599989n=1024S2=-0.893023740S2-S1=0.001803689S=-0.8930（2）S=bianfuhuasimpson(inline('2*x^2*cos(x^2)'),0,sqrt(pi),10e-6,2000)n=2S2=1.076354541S2-S1=2.092222287n=4S

6、2=0.039359358S2-S1=1.036995183n=8S2=-0.430456535S2-S1=0.469815894n=16S2=-0.662796649S2-S1=0.232340113n=32S2=-0.778823323S2-S1=0.116026674n=64S2=-0.836827971S2-S1=0.058004648n=128S2=-0.865829756S2-S1=0.029001785n=256S2=-0.880330615S2-S1=0.014500859n=512S2=-0.8875810

7、42S2-S1=0.007250427n=1024S2=-0.891206256S2-S1=0.003625214S=-0.8912（3）S=bianfuhuasimpson(inline('sqrt(x)'),0,1,10e-6,2000)n=2S2=0.489859598S2-S1=0.185121744n=4S2=0.579745947S2-S1=0.089886349n=8S2=0.623731522S2-S1=0.043985575n=16S2=0.645384849S2-S1=0.021653327n=32S2=

8、0.656091436S2-S1=0.010706587n=64S2=0.661402273S2-S1=0.005310836n=128S2=0.664042680S2-S1=0.002640407n=256S2=0.665357576S2-S1=0.001314896n=512S2=0.