数 值 分 析 课程设计摸版

1、数值分析课程设计报 告 专业：x 学号： x 学生姓名： xxx 指导教师： xx 一、题目用阿当姆斯方法解常微分方程 (用龙格库塔公式给出出发值)1、2、第一题的解第二题的解 求微分方程的精确解 公式两边分别对y和t求积分，又因为，得到，所以 ，t的取值范围为0,1。得到精确解：当t=0.1时，y=1.1111当t=0.2时，y=1.2500当t=0.3时，y=1.4286当t=0.4时，y=1.6667当t=0.5时，y=2当t=0.6时，y=2.5当t=0.7时，y=3.3334当t=0.8时，y=5当t=0.9时，y=10当t=1时，y取不到值二、理论设微分方程的初值问题为：，将微分方

2、程化为 ，然后再小区间上对上式的两边进行积分，即 则微分方程就化为 。其中积分 可以用数值积分的方法求得。通常采用不同的数值积分公式可以得到微分方程数值解的不同公式。如果对于被积函数，取，作为插值结点，则得到插值多项式用插值多项式的积分来近似代替在区间上的积分，即最后得到如下公式，这是一个显示阿当姆斯公式。如果取，作为插值结点，则得到插值多项式代替在区间上的积分得到，这是隐式阿当姆斯公式。四阶阿当姆斯预报校正系统。预报显示公式：；校正隐式公式：；其中。阿当姆斯方法是线性多步法，当计算时，需要用到前面四个点上的解函数值，所以用经典四阶龙格-库塔方法，给出初始值，公式如下：根据题目给出的第一个值，

3、步长h,算出， ，代入阿当姆斯公式得到解。三、算法与程序设计#includestdio.h#includemath.h#includestdlib.hdouble tf(double x,double y) double dy=y*y; return dy;void rkt4(double t,double h,double y) int i; for(i=0;i3;i+) double k1,k2,k3,k4; k1=tf(t+i*h,yi); k2=tf(t+i*h+1.0/2*h,yi+1.0/2*h*k1); k3=tf(t+i*h+1.0/2*h,yi+1.0/2*h*k2); k4

4、=tf(t+i*h+1.0/2*h,yi+h*k3); yi+1=yi+1.0/6*h*(k1+2*k2+2*k3+k4); void admas(double t,double h,double y) int m; double x11; for(m=0;m11;m+) xm=t+m*h; rkt4(t,h,y); for(m=0;m7;m+)ym+4=ym+3+1.0/24*h*(55*tf(xm+3,ym+3)-59*tf(xm+2,ym+2)+37*tf(xm+1,ym+1)-9*tf(xm,ym);ym+4=ym+3+1.0/24*h*(9*tf(xm+4,ym+4)+19*tf(x

5、m+3,ym+3)-5*tf(xm+2,ym+2)+tf(xm+1,ym+1); void main() int j; double t,h,x11,y11; t=0.0;h=0.1;y0=1.0; admas(t,h,y); for(j=0;j11;j+) xj=t+j*h; printf(x=%e,xj); printf(tt); printf(y=%e,yj); printf(n); 其中函数tf的功能是得到f(x,y)的值，故在上述程序中只需修改函数tf中的内容，便可以得到不同的f(x,y)的值。题目第二题只需将函数tf改为： double tf(double x,double y)

6、double dy=0.1*(x*x*x+y*y) return dy;即可得到第二题的解四、第一题误差分析精确解与近似解的比较=0 =0.1=0.2=0.3 =0.4 =0.5 =0.6=0.7 =0.8 =0.9 可以看到精确值与我们算得的值非常接近，四阶龙格库塔公式的截断误差为O()，即具有五阶精度。阿当姆斯预测校正系统也具有五阶精度。五、实际案例生物种群数量问题，设某生物种群在其适应的环境下生存，试讨论该生物种群的数量变化。问题假设：1、假设该生物种群的自然生长率为常数。2、设在其适应的环境下只有该生物种群生存或其他的生物种群的生存不影响该生物种群的生存。3、假设时刻t生物种群数量为N(t) ，由于