數(shù)值積分的matlab實(shí)現(xiàn)(編輯修改稿)
2024-09-19 01:55 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 %a,b表示被積區(qū)間[a,b]的端點(diǎn)%eps表示精度%quad是用Romberg加速算法求得的積分值%R為Romberg表%err表示誤差的估計(jì) 具體函數(shù)編程如下：（1）； function y=f1(x) 然后編寫(xiě)如下程序，%clear。a=0。b=1。T=chen_trap(@f1,a,b)。S=chen_simpson(@f1,a,b)。FT8=trap(@f1,a,b,8)。FS4=simpson(@f1,a,b,8)。quad=Romberg(@f1,a,b,10e6)。syms x f=x^2*sqrt(2*x^2+3)。S1=int(f, x, 0,1)。 format long, Sj=double(S1)。disp(39。積分準(zhǔn)確值為:39。),Sjdisp(39。積分近似值分別為:39。),[T,S,FT8,FS4,quad]disp(39。與精確值比較的誤差為:39。),wuchi=[double(abs(SjT)),double(abs(SjS)),double(abs(SjFT8)),double(abs(SjFS4)),double(abs(Sjquad))]。function y=f2(x)y=x*sin(x)。然后編寫(xiě)如下程序,%clear。a=0。b=2*pi。T=chen_trap(@f2,a,b)。S=chen_simpson(@f2,a,b)。FT8=trap(@f2,a,b,8)。FS4=simpson(@f2,a,b,8)。quad=Romberg(@f2,a,b,10e6)。syms x f=x*sin(x)。S1=int(f, x, 0, 2*pi)。 format long, Sj=double(S1)。disp(39。積分準(zhǔn)確值為:39。),Sjdisp(39。積分近似值分別為:39。),[T,S,FT8,FS4,quad]disp(39。與精確值比較的誤差為:39。),wuchi=[double(abs(SjT)),double(abs(SjS)),double(abs(SjFT8)),double(abs(SjFS4)),double(abs(Sjquad))] function y=f3(x)y=x*sqrt(1+x^2)。然后編寫(xiě)如下程序,%clear。a=0。b=3。T=c
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1