计算方法结课大作业.docVIP

计算方法结课大作业 第一题：求数值积分，精确到。 解：Romberg方法基本思想： 根据 分别算出T0（0）， T0（1），T1（0）…….的值，若对角线上的值的差足够小则停止，否则继续计算，直到足够小为止。 T0（0）= = = 其中sin1= 0.841471 sin(-1)= -0.841471 ==0.438098 T0（1）= T0（0）+ = 0.469049 T1（0）== =0.479366 T0（2）=T0（1） =0.476729 T1（1） ===0.479289 T2（0）===0.479284 T0（3）=T0（2）=0.478645 T1（2）== =0.478645 T2（1）=T0（1）=0.438098 T3（0）===0.479284 化简成缩略图，根据公式可一目了然： T0（0）=, T0（1）=0.469049, T1（0）= 0.479366 T0（2）=0.476729, T1（1） =0.479289 T2（0）=0.479284 T0（3）= 0.478645, T1（2）=0.47, T2（）=,T3（0）= T3（0）=0.479284 第二题：设实验数据如下 1 2 3 4 5 6 7 8 9 10 5.8191 10.3626 18.4715 26.0003 34.9121 48.1684 60.2300 73.9750 92.2391 109.3331 建立形如最小二乘拟合。 解：已知曲线，将表格中所有点依次代入得到方程： 改写成矩阵形式： 根据公式进行转置： 上式整理得: = 整理为方程如下： 得解为： 于是所求拟合曲线为 第三题 ：给定初值问题 经典四阶R-K方法,分别取,计算并打印各点的值，并于准确值相比较。 解：四阶龙格-库塔算法公式如下： 当步长h=0.1时，运行结果如下： 当步长h=0.025时，运行结果如下： 当步长h=0.01时，运行结果如下： 与准确值相比较： 注：（因数据庞大，此处只取前十一个数作比较） 用编程实现如下： 注：此处只举例出当步长h=0.1时执行的程序，当h=0.025和h=0.01时只需稍加修改即可。 #includestdlib.h #includestdio.h /*n表示几等分，n+1表示输出的个数*/ int RungeKutta(double y0,double a,double b,int n,double *x,double *y,double (*function)(double,double)) { double h=(b-a)/n,k1,k2,k3,k4; int i; x[0]=a; y[0]=y0; for(i=0;in;i++) { x[i+1]=x[i]+h; k1=function(x[i],y[i]); k2=function(x[i]+h/2,y[i]+h*k1/2); k3=function(x[i]+h/2,y[i]+h*k2/2); k4=function(x[i]+h,y[i]+h*k3); y[i+1]=y[i]+h*(k1+2*k2+2*k3+k4)/6; } return 1; } double function(double x,double y) { return -30*y+30*x*x+2*x; } //例子求y=-30y+30*x*x+2*x(0x1);y0=1; int main() { double x[11],y[11]; printf(用四阶龙格-库塔方法\n); RungeKutta(1,0,1,10,x,y,function); for(int i=0;i11;i++) printf(x[%d]=%f,y[%d]=%f\n,i,x[i],i,y[i]); return 1; } 第1页 共8页