數(shù)值分析復(fù)化梯形公式,復(fù)化Simpson公式MATLAB程序.doc

《數(shù)值分析復(fù)化梯形公式,復(fù)化Simpson公式MATLAB程序.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫(kù)。

1、分別用復(fù)化梯形公式、復(fù)化Simpson公式計(jì)算定積分，取n=2,4,8,16分別驗(yàn)證結(jié)果（精確值I=4.）。復(fù)化梯形公式求定積分:functionI=tquad(x,y)%復(fù)化梯形求積公式，其中，%x為向量，被積函數(shù)自變量的等距結(jié)點(diǎn)；%y為向量，被積函數(shù)在結(jié)點(diǎn)處的函數(shù)值；n=length(x);m=length(y);%積分自變量的結(jié)點(diǎn)數(shù)應(yīng)與它的函數(shù)值的個(gè)數(shù)相同h=(x(n)-x(1))/(n-1);a=[12*ones(1,n-2)1];I=h/2*sum(a.*y);復(fù)化Simpson公式求定積分：functionI=squad(x,y)%復(fù)化Simpson求積公式，其中

2、，%x為向量，被積函數(shù)自變量的等距結(jié)點(diǎn)；%y為向量，被積函數(shù)在結(jié)點(diǎn)處的函數(shù)值；n=length(x);m=length(y);%積分自變量的結(jié)點(diǎn)數(shù)應(yīng)與它的函數(shù)值的個(gè)數(shù)相同ifrem(n-1,2)~=0I=tquad(x,y);return;endN=(n-1)/2;h=(x(n)-x(1))/N;a=zeros(1,n);fork=1:Na(2*k-1)=a(2*k-1)+1;a(2*k)=a(2*k)+4;a(2*k+1)=a(2*k+1)+1;endI=h/6*sum(a.*y);