数值分析计算机实验3和4.docVIP

实验三 复化辛卜生法，龙贝格法 计算机121班 杨朝晖 122479 一、 目的与要求： ? 通过实际计算体会各种方法的精确度； ? 会编写用复化辛卜生、龙贝格算法求定积分的程序。 二、 实验内容： ? 通过实际计算体会各种方法的精确度并且会编写用复化辛卜生、龙贝格算法求定积分的程序 三、 程序与实例 复化辛卜生公式 算法：复化辛卜生公式为Sn=h/6,计算过程为： 1．令 2．对 计算 3． 。 程序与实例 例 用复化辛卜生法计算积分 源程序： #include iostream.h int main() { double a=0.0,b=1.0; double x=a; double Sn=0; int n; coutinput a number:; cinn; int m=2*n+1; double h=0.5/n; double *f; f=new double[m]; for(int i=0;im;i++) { f[i]=1.0/(1+x*x); x=x+h; } for(int k=0;k=n-1;k++) { Sn=Sn+f[k*2]+4*f[2*k+1]+f[2*(k+1)]; } coutSn=((h/3)*Sn)endl; return 0; } 运行结果为 s(2)=0.785392 s(4)=0.785398 s(8)=0.785398 说明：本例运行了三次，当时，就与时有6位数字相同，若用复化梯形法计算，当n=512时有此结果。 ? 龙贝格算法计算 算法 用事后估计法控制精度 。 源程序：//龙贝格法计算积分，f(x)=1/(1+x^2) #include iostream.h #include math.h #define f(x) (1.0/(1.0+(x)*(x))) void main() { double e=2.71828; double a=0.0,b=1.0; double wu=5*e-6; cout误差精度为：wuendl; double T1=0.5*(f(0)+f(1)); double T2=0.5*T1+0.5*f(0.5); coutT1=T1endl; coutT2=T2endl; if(fabs(T2-T1)wu) {cout|T2-T1|满足精度要求！endl; } double S1=4.0*T2/3.0-1.0*T1/3.0; coutS1=S1endl; double T4=0.5*T2+0.25*(f(0.25)+f(0.75)); double S2=(1.0/3.0)*(4*T4-T2); if(fabs(S2-S1)wu) {cout|S2-S1|满足精度要求！endl; } double C1=(1.0/15.0)*(16*S2-S1); coutC1=C1endl; double T8=0.5*T4+(1.0/8.0)*(f(1.0/8.0)+f(3.0/8.0)+f(5.0/8.0)+f(7.0/8.0)); double S4=(1.0/3.0)*(4*T8-T4); double C2=(16.0/15.0)*S4-(1.0/15.0)*S2; if(fabs(C2-C1)wu) {cout|C2-C1|满足精度要求！endl; } double R1=(64.0/63.0)*C2-(1.0/63.0)*C1; coutR1=R1endl; double T16=0.5*T8+(1.0/16.0)*(f(1.0/16.0)+f(3.0/16.0)+f(5.0/16.0)+f(7.0/16.0)+f(9.0/16.0)+f(11.0/16.0)+f(13.0/16.0)+f(15.0/16.0)); double S8=(1.0/3.0)*(4*T16-T8); double C4=(16.0/15.0)*S8-(1.0/15.0)*S4; double R2=(64.0/63.0)*C4-(1.0/63.0)*C2; coutR2=R2endl; double rr=fabs((R2-R1)); if(rrwu) { cout|R2-R1|5e-6 满足精度要求！endl; } else {cout不满足精度要求！endl; } } 运行结果： 误差精度为：7.5914 T1=0.75 T2=0.775 |T2-T1|满足精度要求！ S1=0.783333 |S2-S1|满足精度要求！ C1=