《2022数值计算第三次大作业参考答案宋宝瑞》由会员分享，可在线阅读，更多相关《2022数值计算第三次大作业参考答案宋宝瑞（5页珍藏版）》请在装配图网上搜索。
1、数值计算第三次大作业1. Euler法functiony=euler(x)%实现euler法%输入值x代表x=x0时，y0旳值x1=x(1);y1=x(2);h=1E-2;n=1;u(1)=y1;forx=x1+h:h:2*pin=n+1;y2=y1+h*fx(x);u(n)=y2;y1=y2;endy=y2;%绘图x=x1:h:2*pi;plot(x,u);gridonxlabel(角度/);ylabel(函数值);title(Euler法);成果为：y=5.00002. 梯形法functionu=tixingfa(x)%实现euler法%输入值x代表x=x0时，y0旳值x1=x(1);y1
2、=x(2);h=1E-2;n=1;y(1)=y1;forx=x1+h:h:2*piy(n+1)=y(n)+h/2*(fx(x)+fx(x+h);n=n+1;endu=y(n);x=x1:h:2*pi;plot(x,y);gridonxlabel(角度/);ylabel(函数值);title(梯形法);成果为：y=4.99993. 龙格库塔法functionu=runger_kuta(x)%实现龙格库塔法%采用四阶龙格库塔x1=x(1);y1=x(2);h=1E-2;n=1;y(1)=y1;%forx=x1+h:h:2*pik1=fx(x,y(n);k2=fx(x+h/2,y(n)+h*k1/2);k3=fx(x+h/2,y(n)+h*k1/2);k4=fx(x+h/2,y(n)+h*k3);y(n+1)=y(n)+h/6*(k1+2*k2+2*k3+k4);n=n+1;endu=y(n);x=x1:h:2*pi;plot(x,y);gridonxlabel(角度/);ylabel(函数值);title(龙格库塔法);成果为：y=5.00004、fx函数%函数fxfunctiony=fx(x,y)x=mod(x,pi);ifx=pi/2y=sin(2*x);elsey=-sin(2*x);end
2022数值计算第三次大作业参考答案宋宝瑞
