数学实验习题概论.docx

题目1：导弹跟踪问题某军一导弹基地发现正北方向120km处海面上有敌艇一艘以90km/h的速度向正东方向行驶，该基地立即发射导弹追击敌艇，导弹速度为450km/h，自动导航系统使导弹在任一时刻都能对准敌艇。试问导弹在何时何处击中敌艇？如果当基地发射导弹的同时，敌艇立即由仪器发觉，假定敌艇为高速快艇，它即刻以135千米/小时的速度向与导弹方向垂直的方向逃逸，问导弹何时何地击中敌艇？敌艇与导弹方向成何夹角逃逸才好？从结论中你能得到些什么启示？思考：敌艇的逃跑策略是什么？在什么样的情形下可能逃脱？1问题分析：此问题涉及到得因素主要有：时间、敌艇逃逸速度、导弹追击速度、敌艇行驶方向、敌艇的位置等等因素。2假设及模型记导弹击中敌艇的时间为t,敌艇逃逸的速度为m,导弹追击速度为n,敌艇行驶方向与导弹的方向的夹角为θ，敌艇与基地的距离为H。先分析第一问，根据题意建立坐标系，设基地的为坐标原点（0，0），敌艇的坐标为（0，120），x轴指向正东方向，y轴指向正北方向。假设经过时间t后，导弹的位置为P（x(t),y(t)）,则由题意可得：+=; ① 其中 n=450km/h另外在t时刻，敌艇的位置为（mt,H）,其中m为90km/h，由于导弹轨迹的切线方向必须指向敌艇，即直线PM的方向就是导弹轨迹上的P点的切线方向，所以有=②公式①和公式②可得到一个微分方程组：又已知x(0)=0,y(0)=0。模型求解采用数值解，这里用差分方法：dx=x(k+1)-x(k),dy=y(k+1)-y(k),dt=t(k+1)-t(k)=h;其中h为时间的步长。Matlab程序如下：画出的函数图形如下：结论分析第一个问中，导弹会在经过0.277秒后，在（25.0058，120.0031）处击中敌艇。题目2：数据拟合Malthus人口指数增长模型中的参数1970年到1980年间美国每隔10年的人口记录如下表：年份1790180018101820183018401850人口（×）3.95.37.29.612.917.123.2年份1860187018801890190019101920人口31.438.650.2625年份193019401950196019701980人口123.2131.7150.7179.3204226.5问题：用以上数据检验马尔萨斯人口指数增长模型，根据检验结果进一步讨论马尔萨斯人口模型的改进。在matlab命令窗口中输入程序如下：画出的拟合曲线和差值曲线如下：其中，标有‘+’的是已知数据点，连接数据点的红色实线是线性插值函数曲线，光滑的蓝色函数曲线是最佳拟合曲线。题目3：现收集了华北地区不同行业的职工工资水平及年平均收入，如下表：地区工业建筑业交通运输、邮电通讯业教育、文化艺术和广播电视事业北京1327173715291244天津1161163115031162河北1113142812451026山西1210149312971108内蒙1246149913321199试问：1）不同行业的职工工资水平是否有显著差异？ 2）不同地区、不同行业的职工工资水平是否有显著差异？解：1）先提出零假设H0：各行业因素总体均值相等。在matlab下计算程序和结果如下： clear fish=[1327 1737 1529 1244; 1161 1631 1503 1162; 1113 1428 1245 1026; 1210 1493 1297 1108; 1246 1499 1332 1199]; save liu fish load liu fishfish = 1327 1737 1529 1244 1161 1631 1503 1162 1113 1428 1245 1026 1210 1493 1297 1108 1246 1499 1332 1199 anoval(fish)??? Undefined command/function anoval. anova1(fish)ans = 6.8562e-005单因素方差分析表：因为p= 6.8562e-0050.05=ɑ,故拒绝H0故可认为不同行业的职工工资水平有显著差异。解2）：首先提出零假设，H1：不同行业因素的效应值为0；H2：不同地区因素的效应值为0；H12：行业因素和地区因素交互效应值为0。双因素方差分析表：检验结