资源描述：
《数学实验报告(同名9509).docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、数学实验报告实验序号：日期：2016年月日班级姓名学号实验名称问题背景描述：根据侦察,发现离我军大炮阵地水平距离10km的前方有一敌军的坦克群正以每小时50km向我军阵地驶来,现欲发射炮弹摧毁敌军坦克群.为在最短时间内有效摧毁敌军坦克,要求每门大炮都能进行精确射击,这样问题就可简化为单门大炮对移动坦克的精确射击问题.假设炮弹发射速度可控制在0.2km/s至0.6km/s之间,问应选择怎样的炮弹发射速度和怎样的发射角度可以最有效摧毁敌军坦克.实验目的：通过运用mathematica建立数学模型解决实
2、际问题来掌握数学建模的一些方法，加强对mathematica的灵活运用。实验原理与数学模型：该模型为斜抛模型。炮弹拥有初速度且与水平面成一定夹角，可将炮弹的运动分解成沿x轴方向与沿y轴方向解决。实验所用软件及版本：WolframMathematica9.0主要内容（要点）：1.选取一个初始速度和发射角，利用Mathematica画出炮弹运行的轨迹2.假定坦克在大炮前方10km处静止不动，炮弹发射的初速度为0.32km/s,应选择什么样的发射角才能击中坦克？画出炮弹运行的轨迹图，通过实验数据和图形来
3、说明你的结论的合理性.3.假定坦克在大炮前方10km处静止不动,探索降低或调高炮弹发射的初速度的情况下,应如何选择炮弹的发射角？从上述讨论中总结出最合理有效的发射速度和发射角.4.在上题结论的基础上,继续探索,假定坦克在大炮前方10km处以每小时50km向大炮方向前进,此时应如何制定迅速摧毁敌军坦克的方案？实验过程记录（含基本步骤、主要程序清单及异常情况记录等）：第一题选择初速度v=0.6km/s,发射角a=45°X轴方向运动为x=cosa*v*tY轴方向运动为y=sina*v*t-1/2*g*t
4、2统一单位将0.6km/s化为600m/s将数据代入利用函数ParametricPlot做出运动轨迹，函数式为ParametricPlot[{Cos[Pi4]*600*t,Sin[Pi4]*600*t-4.9*t^2},{t,0,90},AspectRatio→Automatic]第二题确定速度为320m/s，求最佳角度使得轨迹与X轴交点为(10000,0)先假定发射角为?/4作图ParametricPlot[{Cos[Pi/4]*320*t,Sin[Pi/4]*320*t-4.9*t^2},{t
5、,0,47},AspectRatio®Automatic]进行调整角度调整为?/3.5作图ParametricPlot[{Cos[Pi/3.5]*320*t,Sin[Pi/3.5]*320*t-4.9*t^2},{t,0,52},AspectRatio®Automatic]继续进行不断地调整，发现当发射角度为π/3.7时，落点十分接近(10000,0)点作图如下ParametricPlot[{Cos[Pi3.37]*320*t,Sin[Pi3.37]*320*t-4.9*t^2},{t,0,54}
6、,AspectRatio→Automatic]因此可以确定最合适的发射角就在π/3.7附近，此时可以利用FindRoot函数找出准确值首先需要对已知式做等量变换：∵X=cosa*v*t∴t=x/(cosa*v)将上式代入y=sina*v*t-1/2*g*t2中可得到Y=tana*x-1/2*g*(x/(cosa*v))2将y=0,x=10000,g=9.8,v=320代入利用FindRoot函数求解a的范围在π/3.7附近的a的值：FindRoot[Tan[a]*10000-4.9*(10000(
7、Cos[a]*320))^2==0,{a,Pi3.7}]得出{a→0.9325036430933912}将这个值由弧度制化为360度制a=53.4285°∴最佳发射角为53.4285°第三题由第二题的320m/s起步进行研究1.首先研究速度增大运用与第二题相似的研究方法，先大致计算符合要求的角度(1)V=350m/s时，最佳发射角为π/6.8:ParametricPlot[{Cos[Pi6.8]*350*t,Sin[Pi6.8]*350*t-4.9*t^2},{t,0,32},AspectRati
8、o→Automatic](2)V=400m/s时，最佳发射角为π/9.5:arametricPlot[{Cos[Pi9.5]*400*t,Sin[Pi9.5]*400*t-4.9*t^2},{t,0,26.5},AspectRatio→Automatic](3)V=450m/s时，最佳发射角为π/12.3:ParametricPlot[{Cos[Pi12.3]*450*t,Sin[Pi12.3]*450*t-4.9*t^2},{t,0,23.5},AspectRatio→Autom
