《数值逼近逼近》课程设计.docVIP

课程设计报告 课程名称 数 值 逼 近 专 业 信息与计算科学 班 级 计算091 姓 名 学 号 指导教师 秦新强、胡钢 日 期 2011-07-01 理学院应用数学系 一、目的意义 (1) 进一步掌握复化梯形积分公式，提高运用能力。 (2) 进一步掌握复化抛物线积分公式，提高运用能力。 (3) 进一步掌握c语言，提高编程能力。 二、内容要求 积分计算问题：分别用复化梯形和复化Simpson求积公式计算积分，并比较计算量（精度为10-8）。 三、问题解决的方法与算法 3.1解决方法： 利用复化梯形和复化Simpsom积分公式。 （1）根据两个方程式的误差估计式，求出各自满足精度的等分份数n。 （2）带入所编写的c语言程序，检验并计算出结果。 3.2算法设计： 输入所求的等分数n。 调用复化积分公式，算出n+1点和n点的积分估计数值。 用n+1点和n点的数值检验，比较数值是否达到精度要求。 若达不到精度要求，输出提示，结束程序。 若达到精度要求，输出积分近似结果，程序结束。 计算程序 4.1复化梯形积分程序源代码： #includestdio.h #includemath.h #define f(x)(13*((x)-(x)*(x))*exp(-1.5*(x))) #define epsilon 0double Tixing(double aa,double bb,int n){ int i; double fz; double h=(bb-aa)/n; fz=(f(aa)+f(bb))/2; for(i=1;in;i++) {fz+=f(aa+h*i);} fz*=h; return fz; } int main(){ double a=0.0; double b=4.0; int n; printf(请输入等分值：); scanf(%d,n); printf(\n); double T1; T1=Tixing(a,b,n); double T2; T2=Tixing(a,b,(n+1)); if(fabs(T2-T1)0{ printf(达不到精度要求，请检查等分值是否求解正确！\n); } else printf(用复化梯形积分方法所计算的积分的数值为：%.8lf\n,T1); return 0; } 4.2复化Simpson积分程序源代码: #includestdio.h #includemath.h #define f(x)(13*((x)-(x)*(x))*exp(-1.5*(x))) #define epsilon 0double Simp(double aa,double bb,int n){ int i; double S_n; double h=(bb-aa)/(2*n); double fz,fm,fk; fz=f(aa)+f(bb); fm=0; fk=0; for(i=1;in+1;i++) {fm+=f(aa+(2*i-1)*h);} for(i=1;in;i++) {fk+=f(aa+(2*i)*h);} S_n=h*(fz+4*fm+2*fk)/3; return S_n; } void main() { int n; double a=0.0; double b=4.0; printf(请输入等分值：); scanf(%d,n); double T1; T1=Simp(a,b,n); double T2; T2=Simp(a,b,(n+1)); if(fabs(T2-T1)0{ printf(达不到精度要求，请检查等分值是否求解正确！\n); } else printf(用复化Simpson积分方法所计算的积分的数值为：%.8lf\n,T1); } 五、计算结果与分析 5.1复化梯形积分计算结果与分析： 利用复化梯形积分的误差估计式，求出的值，带入公式计算，是满足精度要求的。 5.2复化Simpson积分计算结果与