資源描述:
《復(fù)化辛普森公式和高斯求積公式方法計(jì)算積分,matlab程序.pdf》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、一、實(shí)驗(yàn)?zāi)康募邦}目實(shí)驗(yàn)?zāi)康模赫莆绽脧?fù)化辛普森公式和高斯求積公式方法計(jì)算積分,熟悉matlab的操作。題目:1.利用復(fù)化辛普森公式計(jì)算積分:11、?xlnxdx02.利用高斯求積公式計(jì)算積分:11、?xlnxdx02、dx實(shí)驗(yàn)步驟:1.利用復(fù)化辛普森公式計(jì)算積分:1.1.建立M文件functiony=f(x)y=sqrt(x)*log(x);1.2.建立M文件functionT_n=F_H_T(a,b,n)h=(b-a)/n;fork=0:nx(k+1)=a+k*h;ifx(k+1)==0x(k+1)=10^(-10);endendT_1=h/2*(f(x(1))+
2、f(x(n+1)));fori=2:nF(i)=h*f(x(i));endT_2=sum(F);T_n=T_1+T_2;1.3.在命令窗口輸入T_n=F_H_T(0,1,20)輸出結(jié)果:2.利用高斯求積公式計(jì)算積分:2.1.建立M文件functions=guassl(a,b,n)h=(b-a)/n;s=0.0;form=0:(1*n/2-1)s=s+h*(guassf(a+h*((1-1/sqrt(3))+2*m))+guassf(a+h*((1+1/sqrt(3))+2*m)));ends;I=int('sin(x)',0,1);c=(I-s)/I;d=vpa(c
3、,10);2.2.1.建立M文件functiony=guassf(x)y=sqrt(x)*log(x);2.2.2.建立M文件functiony=guassf(x)y=sinx/(1+x*x);2.3.運(yùn)行結(jié)果2.3.1.在命令窗口輸入s=guassl(0,1,20)2.3.2.在命令窗口輸入s=guassl(0,1,20)