資源描述：
《復化積分法(復化梯形求積-復化simpson公式-變步長求積法)matlab編程實驗報告》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在行業(yè)資料-天天文庫。
1、復化積分法（復化梯形求積，復化Simpson公式，變步長求積法）MATLAB編程實驗報告一、問題描述：編寫函數(shù)實現(xiàn)復化積分法。二、實驗步驟（過程）：（一）復化求積法（1）復化梯形求積：用復化梯形求積公式求解function[f]=Tn(a,b,n,y)symst;h=(b-a)/n;f=0;fork=1:n+1x(k)=a+(k-1)*hz(k)=subs(y,t,x(k));endfori=2:nf=f+z(i);endq=subs(y,t,a);ify=='sin(t)/t'&&a==0q=1;endp=subs(y,t,b);T=h/2*(q+p+2*f);T=v
2、pa(T,7)clc,clear;symst;a=0;b=1;y=sin(t)/t;n=8;Tn(a,b,n,y);（2）復化Simpson公式：用復化Simpson公式求解function[f]=simpson(a,b,n,y)symst;h=(b-a)/n;f=0;l=0;fork=1:n+1x(k)=a+(k-1)*hw(k)=0.5*h+x(k)z(k)=subs(y,t,x(k));endfori=2:nf=f+z(i);endfori=1:nl=l+w(i);endq=subs(y,t,a);ify=='sin(t)/t'&&a==0q=1;endp=sub
3、s(y,t,b);T=h/2*(q+p+2*f);T=vpa(T,7)clc,clear;symst;a=1;b=2;y=exp(1/t);n=5;simpson(a,b,n,y);（1）變步長求積法：以書本例4.5為例function[f]=TN(a,b,y,R0)symst;T=[];f=0;q=subs(y,t,a);ify=='sin(t)/t'&&a==0q=1;endp=subs(y,t,b);T(1)=(b-a)/2*(q+p);i=2;n=i-1;h=(b-a)/n;z1=a+h/2;z2=subs(y,t,z1);T(2)=T(1)/2+h/2*z2;
4、while((T(i)-T(i-1))/3)>R0i=i+1n=i-1;n=2^(n-1)h=(b-a)/n;f=0;fork=1:nx(k)=a+h*(k-1);w(k)=x(k)+h/2;z(k)=subs(y,t,w(k));f=f+z(k);endT(i)=T(i-1)/2+h/2*fif((T(i)-T(i-1))/3)<=R0break;endendtl=T(i)clc,clear;formatlong;symst;y=sin(t)/t;a=0;b=1;R0=0.5*10^(-3);TN(a,b,y,R0);結論如下：一、結論：復化梯形求積：復化Simpso
5、n公式變步長求積法
