数值分析实习报告.doc

数值分析实习报告 姓名:*** 学号：*** 班级：*** x =0 x =0 x =0.1000 x =0.2000 x =0.3000 x =0.4000 x = 0.5000 x = 0.6000 x = 0.7000 x = 0.8000 x =0.9000x =1.0000 前欧拉法 题目取h二0. 1,利用Euler公式求解 dy/dx=y-2*x/y (0=x= 1) y(o)=i 思路 利用左矩形公式得到的公式yn+l=yn+hf(xn, yn)进行迭代，在进行迭代的过程使用for循 环，让n从1収到10,每迭代一次输出一个x与y,同时x加0.1再进行下一次迭代，一直 到循环结束。 程序 clear; y=l, x=0, %初始化 for n=l: 10 y=l.l *y-0.2*x/y, x=x+0.1, end 运行结果 y=i y=1.1000 y=1.1918 y =1.2774 y =1.3582 y=1.4351 y =1.5090 y= 1.5803 y= 1.6498 y=1.7!78 y =1.7848 后欧拉法 题目取h=0.1,利用Euler公式求解 dy/dx=y-2*x/y (0=x= 1) y(o)=i 2思路 利用右矩形公式得到的公式yn+1二yn+hf(xn+l,yn+l)进行迭代，在进行迭代的过程使用for 循环，让n从1収到10,每迭代一次输出一个x与y,同时x力U 0. 1再进行下一次迭代，一 直到循环结束。 3.程序 clear; yO=l, x0=0, %初始化 for n=l:10 yp=y0+0.1 *(y0-2*x0/y0); x=x0+0.1; xO=x yp=y0+0.1 *(yp-2*x/yp) yo=yp end 4.运行结果 y 二 y 二1 y=1.0918 y=1.1763 y = 1.2546 y = 1.3278 y = 1.3964 y =1.4609 y=1.5216 y =1.5786 y = 1.6321 y = 1.6819 x=0 x =0.1000 x =0.2000 x =0.3000 x =0.4000 x =0.5000 x =0.6000 x =0.7000 x =0.8000 x =0.9000x =1.0000 梯形法 题目 取h=0丄利用梯形法求解 dy/dx=y-2*x/y (0=x= I) y(o)=i 思路 利用右矩形公式得到的公式yn+l=yn+hf(xn+l,yn+l)进行迭代，在进行迭代的过程使用for 循环，让n从1収到10,每迭代一次输出一个x与y,同时x加0. 1再进行下一次迭代，一 直到循环结束。 程序 clear; yO=l, x0=0, for n=l: 10 yp=y0+0.1 *(yO-2*xO/yO); x=x0+0.1; yc=y0+0.1 *(yp-2*x/yp); x0=x y0=( yp+yc)/2 end 4 ?运行结果 y=i x =0 y =1.0959 x=0」000 y=1.1841 x =0.2000 y =1.2662 x =0.3000 y =1.3434 x =0.4000 y =1.4164 x = 0.5000 y =1.4860 x = 0.6000 y =1.5525 x = 0.7000 y=1.6165 x = 0.8000 y =1.6782 x =0.9000 y=1.7379 x =1.0000 改进欧拉法 题目 取h=0丄利用改进欧拉法求解 dy/dx=y-2*x/y (0=x=l) y(o)=i 思路 首先利用欧拉公式公式yn+l=yn+hf (xn, yn)算出预测值yn+1 ,再用公式 yn+l=yn+h[f (xn, yn)+f (xn+1, yn+1)来进行校正,然后进行迭代。迭代的过程使用for循环, 让n从1取到10,每迭代一次输出一个x与y,同时x加0.1再进行下一次迭代，一直到 循环结束。 程序 clear; y0= 1, x0=0, for n=l: 10 yp=y0+0.1 *(y0?2*x0/y0); x=x0+0.1; yc=y0+0.1*(yp-2*x/yp); x0=x y0=(yp+yc)/2 end 运行结果 y=iy =1.0959 y=1.1841 y=i y =1.0959 y=1.1841 y =1.2662 y =1.3434 y=1.4164 y =1.4860 y =1.5525 y=1.6165 y =1.6782y =1.7379 x=0 x =01000 x =0.2000 x =0.3000 x =0.4000 x = 0.5000 x = 0.6000 x = 0.7000 x = 0.