计算机仿真技术演示与验证性实验报告三

《计算机仿真技术演示与验证性实验报告三》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、重庆交通大学上机实验报告学院：信息科学与工程学院实验课程：计算机系统仿实验日期：45姓名陈彦陶学号631306020407专业电子信息工程实验名称连续系统线性多步法仿真实验成绩□实验目的明确□设计方案合理□实验分析全面□操作步骤正确□实验结果合理□实验报告规范评定日期一、实验目的(1)熟悉连续系统线性多步法的核心思想；(2)掌握连续系统线性多步法仿真方法与步骤。二、实验主要内容(1)对于线性微分方程3y+l/2y-2y=0,y(0)=1,y・(0)=0取步长h二0・1,0<丸<二1,采用线性多步法4阶显式公式对该系统进行仿真分

2、析，并将结果同RG法所得解相比较三、实验设计方案利用C语言编写#include#includevoidmain(){inti;floath=0.1;floatylL10J,y2[10J,kll,kl2,kl3,kl4,k21,k22,k23,k24;yl[0]=l；y2[0]=0；for(i=0;i<3;i++)kll=y2[i];k21=2*yl[i]/3-y2[i]/6;kl2=y2[i]+h*k21/2;k22=2*(yl[i]+h*kll/2)/3-l*(y2[i]+h*k21/2)

3、/6;kl3=y2[i]+h*k22/2;k23=2*(yl[i]+h*kl2/2)/3-(y2[i]+h*k22/2)/6;kl4=y2[i]+h*k23;k24=2*(yl[i]+h»kl3)/3-(y2[i]+h*k23)/6;yl[i+l]=yl[i]+h*(kll+2*kl2+2*kl3+kl4)/6;y2[i+l]=y2[i]+h*(k21+2*k22+2*k23+k24)/6;printf(n%fH,yl[i]);}for(i=3;i<10;i++){yl[i+l]=yl[i]+h*(55*y2[i]-59*

4、y2[i-l]+37*y2[i-2]-9*y2[i-3])/24;y2[i+l]=y2[i]+h*(55*(2*yl[i]/3-y2[i]/6)-59*(2*yl[i-l]/3-y2[i-l]/6)+37*(2*yl[i-2]/3-y2[i-2]/6)-9*(2*yl[i-3]/3-y2[i-3]/6))/24;printf(H%f,yl[i]);四、实验结果分析及说明实验结果1.0000001.0033171.0132161.0296541.0526311.0821921.1184231.1614531.2114561

5、.268646DressanykeytocontinueRG法结果如下.0.0000001.000000-0.0597501.641344-0.1461262.692802-0.2709944.4168614.2900964.3113634.2948174.3600204.4393854.5553170.2947640.5456891.0590141.4791241.9262752.372238Pressanykeytocontinue由线性多步法和RG法可知，线性多步法德精度更高单步法：在数值方程中，在求解n+1时刻数值y(

6、n+l)时，仅仅需要第n时刻yn和fn的值;多步法：在数值方程中，在求解(n+1)时刻数值y(n+1)时，不仅需要第n时刻yn和fn,还需要第n~l,n~2等过去时刻y和f信息。