ex3卫星轨道的周长.doc

数值分析实验报告专业 数学学院 班级 信计081班 姓名 陈姝 学号 2008062008 协作队员 实验日期 2010 年 11 月 18 日 星期 四 成绩评定 教师签名 批改日期 题目一、 问题提出教材136页12题用龙贝格求积算法计算其近似值（加速3次）地球卫星轨道是一个椭圆，椭圆周长的计算公式是这里a是椭圆的半长轴，c是地球中心与轨道中心（椭圆中心）的距离，记h为近地点距离，H为远地点距离，R=6371（km）为地球半径，则，我国第一颗人造地球卫星近地点距离h=439（km），远地点距离H=2384（km），试求卫星轨道的周长。二、 模型建立当加速3次时，龙贝格求积算法如下：三、 求解方法#ifndef A#define A#include math.hclass Tpublic:double operator()(const double& x)constreturn (double)sin(x)/(1+x);template class Rombegpublic :void Solve(double a,double b,int MAX_N,int DengFenShu,double e);private :T R;#endif/*方法实现*/#include iostream#include LongBeiGe.h#include math.h#define MAX_SIZE 255using namespace std;template void Rombeg:Solve(double a,double b,int MAX_XunHuan,int DengFenShu,double e)double tMAX_SIZEMAX_SIZE;int n=DengFenShu;double h=(b-a)/n;double s=0;int flag=0;int i;for(i=1;in;i+)s+=R(a+i*h);t01=h*(s+(R(a)+R(b)/2);n*=2;int m;for(m=1;mMAX_XunHuan;m+)for(i=0;im;i+)ti0=ti1;h=(b-a)/n;s=0;int j;for(j=1;jn;j+)s=s+R(a+j*h);t01=h*(s+(R(a)+R(b)/2.0);n*=2;for(i=1;i=m;i+)ti1=ti-11+(ti-11-ti-10)/(pow(double)2,(double)2*m)-1);if(fabs(tm1-tm-11)e)flag=1;couttm1endl;break;if(flag=0)cout无解;/无解/*主函数*#include iostream#include LongBeiGe.h#include LongBeiGe.cppint main()Rombeg r;double a,b,e;int DengFenShu;int Max_XunHuan;coutb;couta;coutDengFenShu;coutMax_XunHuan;coute;r.Solve(a,b,Max_XunHuan,DengFenShu,e);return 0;四、 输出结果由以