数学模型上机程序.doc

数学模型作业 姓名： 学号： 班级： 实验报告193.食饵捕食者模型程序：function xdot=shier(t,x)r=1;d=0.5;a=0.1;b=0.02;xdot=(r-a*x(2).*x(1);(-d+b*x(1).*x(2);ts=0:0.1:15;x0=25,2;t,x=ode45(shier,ts,x0);t,x,plot(t,x),grid,gtext(x(t),gtext(y(t),pause,plot(x(:,1),x(:,2),gridans = 0 25.0000 2.0000 0.1000 27.0818 2.0041 0.2000 29.3344 2.0170 0.3000 31.7689 2.0394 0.4000 34.3961 2.0726 0.5000 37.2258 2.1178 0.6000 40.2673 2.1767 0.7000 43.5012 2.2534 0.8000 46.9360 2.3503 0.9000 50.6072 2.4683 1.0000 54.5301 2.6106 1.1000 58.6999 2.7819 1.2000 63.0917 2.9891 1.3000 67.6604 3.2411 1.4000 72.3409 3.5484 1.5000 77.0479 3.9238 1.6000 81.6759 4.3819 1.7000 86.0996 4.9391 1.8000 90.1732 5.6140 1.9000 93.7311 6.4268 2.0000 96.5873 7.4000 2.1000 98.5360 8.5577 2.2000 99.3055 9.9234 2.3000 98.6143 11.5085 2.4000 96.2851 13.3067 2.5000 92.2472 15.2882 2.6000 86.5853 17.3947 2.7000 79.5349 19.5427 2.8000 71.5364 21.6225 2.9000 63.0848 23.5300 3.0000 54.6236 25.1819 3.1000 46.5441 26.5163 3.2000 39.1860 27.4921 3.3000 32.7932 28.0978 3.4000 27.3368 28.3766 3.5000 22.7375 28.3764 3.6000 18.9134 28.1426 3.7000 15.7771 27.7178 3.8000 13.2354 27.1426 3.9000 11.1873 26.4556 4.0000 9.5278 25.6911 4.1000 8.1758 24.8740 4.2000 7.0684 24.0245 4.3000 6.1592 23.1580 4.4000 5.4185 22.2853 4.5000 4.8129 21.4162 4.6000 4.3124 20.5582 4.7000 3.8967 19.7164 4.8000 3.5500 18.8948 4.9000 3.2601 18.0961 5.0000 3.0181 17.3219 5.1000 2.8154 16.5735 5.2000 2.6460 15.8515 5.3000 2.5047 15.1562 5.4000 2.3872 14.4875 5.5000 2.2898 13.8454 5.6000 2.2099 13.2295 5.7000 2.1455 12.6393 5.8000 2.0953 12.0740 5.9000 2.0578 11.5330 6.0000 2.0317 11.0154 6.1000 2.0162 10.5207 6.2000 2.0106 10.0479 6.3000 2.0143 9.5964 6.4000 2.0270 9.1653 6.5000 2.0483 8.7539 6.6000 2.0781 8.3614 6.7000 2.1164 7.9870 6.8000 2.1632 7.6300 6.9000 2.2188 7.2897 7.0000 2.2833 6.9654 7.1000 2.3573 6.6565 7.2000 2.4410 6.3622 7.3000 2.5351 6.0821 7.4000 2.6401 5.8155 7.5000 2.7566 5.5618 7.6000 2.8855 5.3204 7.7000 3.0276 5.0910 7.8000 3.1837 4.8729 7.9000 3.3549 4.6656 8.0000 3.5424 4.4688 8.1000 3.7478 4.2819 8.2000 3.9723 4.1046 8.3000 4.2174 3.9365 8.4000 4.4847 3.7773 8.5000 4.7761 3.6265 8.6000 5.0937 3.4838 8.7000 5.4398 3.3490 8.8000 5.8171 3.2217 8.9000 6.2283 3.1017 9.0000 6.6766 2.9888 9.1000 7.1653 2.8826 9.2000 7.6978 2.7831 9.3000 8.2781 2.6900 9.4000 8.9100 2.6032 9.5000 9.5980 2.5225 9.6000 10.3468 2.4479 9.7000 11.1620 2.3794 9.8000 12.0494 2.3166 9.9000 13.0153 2.2596 10.0000 14.0665 2.2084 10.1000 15.2102 2.1629 10.2000 16.4541 2.1233 10.3000 17.8063 2.0898 10.4000 19.2755 2.0627 10.5000 20.8708 2.0422 10.6000 22.6016 2.0287 10.7000 24.4779 2.0228 10.8000 26.5102 2.0249 10.9000 28.7093 2.0355 11.0000 31.0844 2.0558 11.1000 33.6416 2.0875 11.2000 36.3982 2.1308 11.3000 39.3680 2.1869 11.4000 42.5598 2.2573 11.5000 45.9776 2.3447 11.6000 49.6205 2.4525 11.7000 53.4829 2.5849 11.8000 57.5539 2.7468 11.9000 61.8180 2.9441 12.0000 66.2549 3.1834 12.1000 70.8428 3.4692 12.2000 75.5222 3.8129 12.3000 80.1900 4.2348 12.4000 84.7213 4.7545 12.5000 88.9686 5.3918 12.6000 92.7622 6.1658 12.7000 95.9102 7.0955 12.8000 98.1986 8.1995 12.9000 99.3910 9.4962 13.0000 99.2287 11.0036 13.1000 97.4360 12.7384 13.2000 93.8972 14.6791 13.3000 88.6885 16.7638 13.4000 82.0129 18.9151 13.5000 74.2204 21.0377 13.6000 65.7971 23.0201 13.7000 57.2540 24.7617 13.8000 49.0238 26.1948 13.9000 41.4342 27.2791 14.0000 34.7062 28.0016 14.1000 28.9351 28.3802 14.2000 24.0660 28.4610 14.3000 20.0105 28.2919 14.4000 16.6759 27.9181 14.5000 13.9650 27.3821 14.6000 11.7742 26.7239 14.7000 9.9982 25.9788 14.8000 8.5540 25.1736 14.9000 7.3762 24.3300 15.0000 6.4158 23.4645实验结果：139传染病模型模型3程序：function y=ill(t,x)a=1;b=0.3;y=a*x(1)*x(2)-b*x(1),-a*x(1)*x(2);ts=0:50; x0=0.02,0.98; t,x=ode45(ill,ts,x0);t,xans = 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