《数学模型实验四new.doc》由会员分享，可在线阅读，更多相关《数学模型实验四new.doc（5页珍藏版）》请在金锄头文库上搜索。
1、数学系数学模型课程实验报告姓 名学号班级10级指导教师实验地点逸夫楼212实验日期2012-05-14专业数学与应用数学实验组员一、 实验项目名称Matlab软件程序设计二、实验目的 用Matlab程序设计求解统计分析模型三、实验内容及要求 （1）Matlab程序设计； （2）用Matlab软件统计分析研究按年龄分组的种群增长； （3）编写Matlab画图程序。四、实验材料、工具Matlab软件 计算机 五、实验过程1.建立M文件，在文件中输入繁殖率，存活率，种群数量等数据，保存文件，2.输入数据：程序：b=0,0.2,1.8,0.8,0.2;s=diag(0.5,0.8,0.0.8,0.1)
2、;l=b;s,zeros(4,1);x(:,1)=100*ones(5,1);n=30;for k=1:n x(:,k+1)=l*x(:,k)endround(x)y=diag(1./sum(x);z=x*y,k=0:30;subplot(1,2,1);plot(k,x),gridsubplot(1,2,2);plot(k,x),grid运行结果如下:x = Columns 1 through 8 100.0000 300.0000 220.0000 154.8000 264.8800 251.3200 196.1840 256.8824 100.0000 50.0000 150.0000 1
3、10.0000 77.4000 132.4400 125.6600 98.0920 100.0000 80.0000 40.0000 120.0000 88.0000 61.9200 105.9520 100.5280 100.0000 80.0000 64.0000 32.0000 96.0000 70.4000 49.5360 84.7616 100.0000 10.0000 8.0000 6.4000 3.2000 9.6000 7.0400 4.9536 Columns 9 through 16 269.3688 232.9739 263.7238 284.2604 264.7160
4、279.6724 300.1434 293.0213 128.4412 134.6844 116.4869 131.8619 142.1302 132.3580 139.8362 150.0717 78.4736 102.7530 107.7475 93.1895 105.4895 113.7042 105.8864 111.8690 80.4224 62.7789 82.2024 86.1980 74.5516 84.3916 90.9633 84.7091 8.4762 8.0422 6.2779 8.2202 8.6198 7.4552 8.4392 9.0963 Columns 17
5、through 24 300.9651 318.6957 319.6984 325.4988 340.3531 346.2450 352.2769 365.0525 146.5106 150.4825 159.3478 159.8492 162.7494 170.1766 173.1225 176.1384 120.0574 117.2085 120.3860 127.4783 127.8794 130.1995 136.1413 138.4980 89.4952 96.0459 93.7668 96.3088 101.9826 102.3035 104.1596 108.9130 8.470
6、9 8.9495 9.6046 9.3767 9.6309 10.1983 10.2303 10.4160 Columns 25 through 31 373.7377 380.9616 392.6104 402.8953 411.5666 422.8871 434.1877 182.5263 186.8689 190.4808 196.3052 201.4477 205.7833 211.4436 140.9108 146.0210 149.4951 152.3846 157.0442 161.1581 164.6266 110.7984 112.7286 116.8168 119.5961
7、 121.9077 125.6353 128.9265 10.8913 11.0798 11.2729 11.6817 11.9596 12.1908 12.5635 x = Columns 1 through 7 100.0000 300.0000 220.0000 154.8000 264.8800 251.3200 196.1840 100.0000 50.0000 150.0000 110.0000 77.4000 132.4400 125.6600 100.0000 80.0000 40.0000 120.0000 88.0000 61.9200 105.9520 100.0000
8、80.0000 64.0000 32.0000 96.0000 70.4000 49.5360 100.0000 10.0000 8.0000 6.4000 3.2000 9.6000 7.0400 Columns 8 through 14 256.8824 269.3688 232.9739 263.7238 284.2604 264.7160 279.6724 98.0920 128.4412 134.6844 116.4869 131.8619 142.1302 132.3580 100.5280 78.4736 102.7530 107.7475 93.1895 105.4895 11
9、3.7042 84.7616 80.4224 62.7789 82.2024 86.1980 74.5516 84.3916 4.9536 8.4762 8.0422 6.2779 8.2202 8.6198 7.4552 Columns 15 through 21 300.1434 293.0213 300.9651 318.6957 319.6984 325.4988 340.3531 139.8362 150.0717 146.5106 150.4825 159.3478 159.8492 162.7494 105.8864 111.8690 120.0574 117.2085 120.
10、3860 127.4783 127.8794 90.9633 84.7091 89.4952 96.0459 93.7668 96.3088 101.9826 8.4392 9.0963 8.4709 8.9495 9.6046 9.3767 9.6309 Columns 22 through 28 346.2450 352.2769 365.0525 373.7377 380.9616 392.6104 402.8953 170.1766 173.1225 176.1384 182.5263 186.8689 190.4808 196.3052 130.1995 136.1413 138.4
11、980 140.9108 146.0210 149.4951 152.3846 102.3035 104.1596 108.9130 110.7984 112.7286 116.8168 119.5961 10.1983 10.2303 10.4160 10.8913 11.0798 11.2729 11.6817 Columns 29 through 31 411.5666 422.8871 434.1877 201.4477 205.7833 211.4436 157.0442 161.1581 164.6266 121.9077 125.6353 128.9265 11.9596 12.
12、1908 12.5635ans = Columns 1 through 13 100 300 220 155 265 251 196 257 269 233 264 284 265 100 50 150 110 77 132 126 98 128 135 116 132 142 100 80 40 120 88 62 106 101 78 103 108 93 105 100 80 64 32 96 70 50 85 80 63 82 86 75 100 10 8 6 3 10 7 5 8 8 6 8 9 Columns 14 through 26 280 300 293 301 319 320 325 340 346 352 365 374 381 132 140 150 147 150 159 160 163 170 173 176 183 187 114 106 112 120 117 120 127 128 130 136 138 141 146 84 91 85 89 96 94 96 102 102 104 109 111 113 7 8 9 8 9 10 9 10 10 10 10 11 11 Columns 27 through 31 393 403 412 423 434 190 196 20
