常微分方程78839new.docVIP

福建农林大学计算机与信息学院 数学系 课程实习报告 课程名称： 常微分方程 实习题目： 课程实习报告 姓 名： 张俊豪 系： 计算机与信息学院 专 业： 信息与计算科学 年 级： 2007级 学 号： 061152020 指导教师： 陈永雪 职 称： 助教 2009年 6 月 11日  福建农林大学计算机与信息学院数学系 课程实习报告结果评定 评语： 成绩： 指导教师签字： 评定日期： 目 录 1.实习目的和任务………………………………………………………………1 2.实习要求………………………………………………………………………1 3.实习地点………………………………………………………………………1 4.主要仪器设备…………………………………………………………………1 5.实习内容………………………………………………………………………1 5.1 欧拉格式方法……………………………………………………………1 5.2 改进欧拉格式方法………………………………………………………3 5.3四阶-龙格库塔方法………………………………………………………5 5.4 matlab软件画图…………………………………………………………7 6.问题讨论与分析………………………………………………………………9 7.结束语…………………………………………………………………………10 参考文献…………………………………………………………………………10 常微分方程 1.实习的目的和任务 目的：通过课程实习，了解，掌握并可以熟练的去运用matlab和c++数学软件，去解决在理论和实际上的问题。 任务：通过具体的案例实际的问题，利用matlab和c++数学软件来计算问题的结果，并通过图形图像的形式表现，并分析问题的结论。 2.实习要求 要求：能够熟练使用matlab和c++软件。 3.实习地点：数学实验室 4.主要仪器设备 计算机、 Microsoft Windows XP、 Matlab 6.5 Turbo CC++ 3 5.实习内容 运用欧拉格式方法、改进欧拉格式方法和龙格库塔格式方法计算y=y+exp(2*x) （0=x=1） 在x＝0时y＝1的初值解。步长设为h=0.05，区间进行20等分，进行求解，并比较各种方法计算的精确程度。 5.1 欧拉格式方法 #includeiostream.h #includestdio.h #includemath.h /*举例方程：y=y+exp(2*x) （0=x=1） y(0) = 1*/ double f(double x, double y, double h)//方程形式 { double re; if(x==0)re=1; else re=y+exp(2*x); return re; } int main() {double x0,y0,x1,y1,h;//h为步长 double y; int N; scanf(%.2f %.6f %.6f\n, x0, y0, h); while(cinx0y0hN) {int n=0; for(;nN;n++) { x1=x0+h; y=exp(2*x1); //精确值 y1=y0+h*f(x0,y0,h); printf(%.2f %.6f %.6f\n, x1, y1, y); x0=x1; y0=y1; } } return 0; } 在命令窗口中输入 0 1 0.05 20时（按回车键） 结果显示如下： 图5.1 5.2 改进欧拉格式方法 #include iostream.h #include math.h void main() { float x0,y0,h,**result; int N; float f(float x,float y); void euler(float (*f)(float,float),float,float,int,float,float **result); coutInputs the initial spot:; cinx0y0; coutInput solution integer and length of stride:; cinNh; result=new float *[N]; for(int i=0;iN;i++) result[i]=new float[2]; euler(f,x0,y0,N,h,result); coutThe operation result is:endl