資源描述：
《數(shù)值分析報(bào)告某實(shí)驗(yàn)——數(shù)值積分》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。
1、實(shí)用標(biāo)準(zhǔn)文案桂林電子科技大學(xué)數(shù)學(xué)與計(jì)算科學(xué)學(xué)院實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)室：06406實(shí)驗(yàn)日期：2014年11月21日院（系）數(shù)學(xué)與應(yīng)用數(shù)學(xué)年級(jí)、專業(yè)、班姓名成績(jī)課程名稱數(shù)值分析實(shí)驗(yàn)實(shí)驗(yàn)項(xiàng)目名稱實(shí)驗(yàn)積分指導(dǎo)教師李光云一、實(shí)驗(yàn)?zāi)康耐ㄟ^(guò)實(shí)驗(yàn)掌握利用Matlab進(jìn)行數(shù)值積分的操作，掌握Matlab中的幾種內(nèi)置求積分函數(shù)，進(jìn)一步理解復(fù)化梯形，復(fù)化辛普生公式，并編程實(shí)現(xiàn)求數(shù)值積分二、實(shí)驗(yàn)原理Matlab中，有內(nèi)置函數(shù)計(jì)算積分：>>z=trapz(x,y)其中，輸入x，y分別為已知數(shù)據(jù)的自變量和因變量構(gòu)成的向量，輸出為積分值。>>z=quad(
2、fun,a,b)這個(gè)命令是使用自適應(yīng)求積的方法計(jì)算積分的命令。其中，fun為被積函數(shù)，a，b為積分區(qū)間。我們還可以利用復(fù)化梯形公式三、使用儀器，材料電腦MATLAB四、實(shí)驗(yàn)內(nèi)容與步驟1.編寫(xiě)復(fù)化辛普生公式的Matlab的程序。2.利用復(fù)化梯形法程序計(jì)算，記錄下計(jì)算結(jié)果隨著n增加的變化情況，畫(huà)圖與復(fù)化梯形公式的情況比較收斂速度。3.積分的原函數(shù)無(wú)法用初等函數(shù)表達(dá)，結(jié)合Matlab復(fù)化梯形程序，用描點(diǎn)法繪制其原函數(shù)在區(qū)間的圖形。五、實(shí)驗(yàn)過(guò)程原始記錄（數(shù)據(jù)，圖表，計(jì)算等）精彩文檔實(shí)用標(biāo)準(zhǔn)文案一、復(fù)化Simpson公式程序：f
3、unctions=Simpson(a,b,n)%輸出s為積分的數(shù)值解，輸入(a,b)為積分區(qū)間，n為等分區(qū)間的個(gè)數(shù).h=(b-a)/(n*2);s1=0;s2=0;s=h*(f(a)+f(b))/3;%先計(jì)算特殊兩點(diǎn)相加.fork=1:nx1=a+h*(2*k-1);%利用循環(huán)計(jì)算其他點(diǎn)的相加.s1=s1+f(x1);endfork=1:(n-1)x2=a+h*2*k;%利用循環(huán)計(jì)算其他點(diǎn)的相加.s2=s2+f(x2);ends=s+h*(4*s1+2*s2)/3;畫(huà)圖程序formatlong;%k為等分區(qū)間個(gè)數(shù)，t存
4、儲(chǔ)積分值.k=2:1:40;fori=1:length(k)t(i)=Simpson(0,1,k(i));disp([k(i),t(i)]);endplot(k,t,'.','MarkerSize',20)二、Simpson(0,1,26)ans=3.14159265358779>>Untitled62.000000000000003.141568627450983.000000000000003.141591780936044.000000000000003.141592502458715.0000000000000
5、03.141592613939226.000000000000003.141592640305387.000000000000003.141592648320658.000000000000003.14159265122482精彩文檔實(shí)用標(biāo)準(zhǔn)文案9.000000000000003.1415926524231710.000000000000003.1415926529697911.000000000000003.1415926532398112.000000000000003.1415926533821513.00000
6、0000000003.1415926534613414.000000000000003.1415926535074515.000000000000003.1415926535353616.000000000000003.1415926535528417.000000000000003.1415926535641118.000000000000003.1415926535715619.000000000000003.1415926535766120.000000000000003.1415926535801121.000
7、000000000003.1415926535825622.000000000000003.1415926535843223.000000000000003.1415926535856024.000000000000003.1415926535865525.000000000000003.1415926535872526.000000000000003.1415926535877927.000000000000003.1415926535881928.000000000000003.1415926535885129.0
8、00000000000003.1415926535887530.000000000000003.14159265358894精彩文檔實(shí)用標(biāo)準(zhǔn)文案31.000000000000003.1415926535890932.000000000000003.1415926535892233.000000000000003.141592653
