50万买房半年难拿购房合 同6.docVIP

实验6 微分方程模型上机实验 （1） 一、实验目的和要求 掌握用Matlab软件求解微分方程模型的编程方法。 二、实验环境 Windows系列操作系统，Matlab软件。 三、实验内容 使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解)； 利用图形对解的特征作定性分析。 四、实验步骤 开启软件平台——Matlab，开启Matlab编辑窗口； 根据微分方程求解步骤编写M文件 保存文件并运行； 观察运行结果(数值或图形)； 根据观察到的结果和体会写出实验报告。 五、实验任务： 根据实验内容和步骤，完成以下实验，要求写出实验报告 1．编程练习课本P140页SIR模型（14），并画出图形7，8。 ts=0:50; x0=[0.02,0.98]; [t,x]=ode45(ill,ts,x0);[t,x] ans = 0 0.0200 0.9800 1.0000 0.0390 0.9525 2.0000 0.0732 0.9019 3.0000 0.1285 0.8169 4.0000 0.2033 0.6927 5.0000 0.2795 0.5438 6.0000 0.3312 0.3995 7.0000 0.3444 0.2839 8.0000 0.3247 0.2027 9.0000 0.2863 0.1493 10.0000 0.2418 0.1145 11.0000 0.1986 0.0917 12.0000 0.1599 0.0767 13.0000 0.1272 0.0665 14.0000 0.1004 0.0593 15.0000 0.0787 0.0543 16.0000 0.0614 0.0507 17.0000 0.0478 0.0480 18.0000 0.0371 0.0460 19.0000 0.0287 0.0445 20.0000 0.0223 0.0434 21.0000 0.0172 0.0426 22.0000 0.0133 0.0419 23.0000 0.0103 0.0415 24.0000 0.0079 0.0411 25.0000 0.0061 0.0408 26.0000 0.0047 0.0406 27.0000 0.0036 0.0404 28.0000 0.0028 0.0403 29.0000 0.0022 0.0402 30.0000 0.0017 0.0401 31.0000 0.0013 0.0400 32.0000 0.0010 0.0400 33.0000 0.0008 0.0400 34.0000 0.0006 0.0399 35.0000 0.0005 0.0399 36.0000 0.0004 0.0399 37.0000 0.0003 0.0399 38.0000 0.0002 0.0399 39.0000 0.0002 0.0399 40.0000 0.0001 0.0399 41.0000 0.0001 0.0399 42.0000 0.0001 0.0399 43.0000 0.0001 0.0399 44.0000 0.0000 0.0398 45.0000 0.0000 0.0398 46.0000 0.0000 0.0398 47.0000 0.0000 0.0398 48.0000 0.0000 0.0398 49.0000 0.0000 0.0398 50.0000 0.0000 0.0398 plot(t,x(:,1),t,x(:,2)),grid,pause pl