MCMA题让标枪飞

1、五一数学建模竞赛承诺书我们仔细阅读了五一数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网 上咨询等）与本队以外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其它公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参 考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规 则的行为，我们愿意承担由此引起的一切后果。我们授权五一数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包 括进行网上公示，在书籍、期刊和其

2、他媒体进行正式或非正式发表等）。参赛题号（从 A/B/C中选择一项填写）： A参赛组别（研究生、本科、专科、高中）： 本科所属学校（学校全称）参赛队员（打印并签名）：1. 2. 3. 日期：年月日获奖证书邮寄地址： 邮政编码： 收件人姓名： 联系电话： 五一数学建模竞赛题目 计标枪飞关键词：多元回归统计 K-means聚类 距离相关牛顿欧拉方程 刚体微分学 摘要:标枪的技术参数包括，标枪的长度、重量、几何形状、重心的位置、形心的位置等，标枪 的形心位置由枪身、枪头、把手三者综合判定而得，题目给出了三者长轴和直径的相关数据， 根据k-means聚类进行分析求解此问题。考虑到初始攻角受到了出手角的

3、影响，我们猜测初始 攻角和出手角之间存在联系，可以寻找初始攻角和出手角之间的关系，简化初始攻角、出手角、 投掷距离、出手速度之间的关系，利用多元回归统计进行分析。用前面已经拟合的公式，发现 初始攻角存在某一合法置信区间，出手速度固定，通过逐渐改变初始攻角的取值，可以得到多 组投掷距离，取最优解。阅读题意，顺风逆风改变了出手速度，前述问题，已经解决了出手角、 初始攻角对投掷距离的影响，而初始俯仰角速度是指标枪出手瞬间，标枪枪头半轴（枪尖至重 心部分）在单位时间内围绕重心旋转的角度，建立标枪飞行的数学模型，利用牛顿欧拉方程求 解。最后，根据上述问题得到的结论，对运动员出手速度、出手角、初始攻角、初

4、始俯仰角速 度、风向及风速等要素，对标枪投掷距离影响的相对重要性进行分析总结，用自然语言进行详 尽的描述。问题重述标枪的投掷是一项历史悠久的田径比赛项目。 标枪投掷距离的远近受到运动员水平 （出手速度、出手角、初始攻角、出手高度、出手时标枪的初始俯仰角速度等），标枪的技术参数（标枪的长度、重量、几何形状、重心的位置、形心的位置等）和比赛环境（空气的密度与粘度、风力、风向等）三方面因素的影响。1 .现有某型标枪的测量尺寸数据。请估算该型标枪沿标枪中轴线剖面面积、标枪表面 积和标枪形心的位置。2 .现有某标枪比赛中24名运动员使用同型号标枪投掷的实测数据。 请对这些数据进行 分析，建立合适的数学模

5、型，找出标枪飞行过程中的运动规律。3 .假设某型标枪的几何参数同问题1,重心位于形心前10cm （重心在枪尖与形心之间）, 出手时标枪的初始俯仰角速度为 00在无风的前提下，对标枪投掷出手瞬间、出手 后的受力及运动情况进行分析，建立标枪飞行的数学模型并求解如下问题：（1）假设某运动员投掷出手速度为29.70m/%'、出手角为36.6、初始攻角为-0书, 请估算出标枪的投掷距离。（2）假设某运动员投掷出手速度为30m/s,请给出最佳的出手角和初始攻角使得投 掷距离最大并估算出标枪的投掷距离。4 .假设标枪技术参数同问题3 ,风向分别为顺风和逆风，风速分别为 3m/s, 和运动员投掷出手速

6、度为31.70m/so请建立标枪飞行的数学模型，给出最佳出手角、最佳初始攻角、最佳初始俯仰角速度使得投掷距离最大。5 .标枪技术参数同问题3,请分析运动员出手速度、出手角、初始攻角、初始俯仰角 速度、风向及风速等要素对标枪投掷距离影响的相对重要性。题分析问题一的分析标枪的技术参数包括，标枪的长度、重量、几何形状、重心的位置、形心的位置等，标枪的形心位置由枪身、枪头、把手三者综合判定而得，题目给出了三者长轴和直径的相关数据，可以根据k-means 聚类进行分析，来求解此问题。问题二的分析题目给出了 24 位运动员的标枪比赛实测数据， 考虑到初始攻角受到了出手角的影响， 我们猜测初始攻角和出手角之

7、间存在联系，可以寻找初始攻角和出手角之间的关系，简化初始攻角、出手角、投掷距离、出手速度之间的关系，利用多元回归统计进行分析。问题三的分析第一小问，阅读题意，根据第二问得出的结论进行分析；第二小问，用前面已经拟合的公式，发现初始攻角存在某一合法置信区间，在区间中设置步长为 1 （题目给出的数据均为整数） ，猜想，初始攻角结果也为一个合法的整数，故出手速度固定，通过逐渐改变初始攻角的取值，可以得到多组投掷距离，取最优解问题四的分析阅读题意， 顺风逆风改变了出手速度， 前述问题， 已经解决了出手角、 初始攻角对投掷距离的影响， ，而初始俯仰角速度是指标枪出手瞬间，标枪枪头半轴（枪尖至重心部分）在单

8、位时间内围绕重心旋转的角度，建立标枪飞行的数学模型，给出最佳出手角、最佳初始攻角、最佳初始俯仰角速度使得投掷距离最大，利用牛顿欧拉方程求解。问题五的分析根据上述问题得到的结论，对运动员出手速度、出手角、初始攻角、初始俯仰角速度、风向及风速等要素，对标枪投掷距离影响的相对重要性进行分析总结，用自然语言进行详尽的描述。三、符号说明我在底下标注，自助总结下吧四、模型假设假设运动员：(1) 出手高度为2m(2) 标枪重量为800g(3) 不考虑标枪在飞行过程中的进动影响(4) 空气密度为L1H4X10(5) 空气粘度为L84 X 10 ”以?叱帕秒)(6) 如果考虑有风，空气阻力(推力)水平(7) 考

9、虑地心引力(8) 除了出手角度外，考虑初始俯仰角和初始攻角(初始攻角3=持枪角丫？出手角”)五、模型建立与求解5.1 问题一模型建立与求解5.1.1 模型建立标枪长轴与直径数据在不同的部位存在一定的关系逼近，利用题中已给的数据近似聚类，并进行统计描述，以完成对形心位置的准确预测。即利用k-means聚类研究它们之间的关系。聚类是一个将数据集中在某些方面相似的数据成员进行分类组织的过程，发现它们的内在结我们建立枪身、枪头、把手三者相关的聚类模型(1)选择3个点作为初始质心。(2)将每个点指派到最近的质心，形成 3个簇，重新计算每个簇的质心，直到质心不发生变化。(Result表示结果，v1v2v3

10、表示枪身、枪头、把手)resultkmeans(V1 ,V2 ,V3)D(x)表示沿枪轴方向的直径分布，我们将D(x)的度量为标枪长度分割为多个刻度，一共 70份，利用matlab进行拟合得到标枪沿枪轴方向的直径分布模型：D(x)表示沿枪轴方向的直径分布在matlab中，标枪长度刻度为 y轴，沿枪轴方向的直径分布为x轴，根据拟合优度选用以e为底的指数函数拟合，拟合得到如下图像拟合优度(x 2.506)2(x 0.9716)2(x 1.795)2D 1(x) 40.46e 230.36e 0.771618.04e反函数运算得到 D(x)表达式查阅相关资料，结合分析，发现，该型标枪沿标枪中轴线剖面

11、面积以及标枪表面积满足下列表达式(S剖表示标枪中轴线剖面面积，x表示从枪尖为起点的标枪长度，D(x)表示沿枪轴方向的直径分布，L表示标枪长度，800g标枪取)L° D(x)dxLL面心2(D'(x)/2)2dx° xD(x)dx5.1.2 模型求解根据上述模型，可以得到，三者样本量集群分布如下:集科大小由结合、分割轮廓可知三者品质如下，形心位于枪身（品质：劣）处:集群品第1 .砧育及分隔的轴廊渊触经过k-means聚类得到下述结果:口避123VARQ0001VARDDD92345.6517.871 408,8931.132348.7017.42最终经过计算：形心位于

12、距离枪尖1050mm处利用matlab软件对剖面面积以及标枪表面积进行辛普森梯形数值积分，得到表面积为145074（mm2）,剖面积 67336（mm2）5.2 问题二模型建立与求解5.2.1 模型的建立利用SPSS软件对变量进行距离相关性分析，发现运动员 出手速度v （米/秒），出手角（度），初始攻角（度），投掷距离dist （米）4者之间存在相关性，结果如下：但向盘之加的和偶MVAR0D0D2VAR30003VAR00004VAR000D5VAR000021 900137-.564979VAR000D31371.000-457251VAR00004-.564J571 oac-.610VAR

13、000D597 g251-6101 000我们以v,为自变量，dist为因变量进行多元回归分析103.00-9D.OO-03.00-7D.OO-6100-SD.QQ-50 DOEDOD7D 0D93 3D1 JU JijVARWVO道是相似性矩阵拟合准确度符合预期0%25%50%75%100%碎;蘸度查阅文献我们知道，出手角度，初始俯仰角，则，初始攻角，我们经过分析及查阅相关文档发现，考虑初始攻角由出手角以及初始俯仰角联合影响，故，只考虑初始攻角即可，拟合发现，初始攻角、出手速度V二者之间对投掷距离dist的影响满足三元三次关系：dist 2115 94.97V 221.10.5152v26.

14、849v7.788 2 0.003476v30.01733V20.1229V 2 0.08907 3函数图像如下：傅里叶函数6次拟合出手角和初始攻角之间的关系：33.73 0.7881cos0.72671.525sin 0.72671.04cos1.45340.1092sin1.45341.196cos2.18012.07sin 2.18010.51cos2.90680.07145sin 2.90680.6295cos3.63350.07075sin 3.63351.13cos4.36021.743sin 4.36025.2.2模型的求解标枪飞行过程中的运动规律如下：(1)在只考虑出手速度V,

15、出手角，高度h ,情况下则满足斜抛运动方程，距离s为SSiS2v2sin2q h g S2iaT(2)考虑初始攻角，出手角 ，出手速度v时，忽略高度的影响因素，距离 distdist 2115 94.97v 221.10.5152v26.849v7.788 2 0.003476v30.01733V20.1229v 2 0.08907 35.3 问题三模型建立与求解5.3.1 模型的建立假设某型标枪的几彳s参数同问题1,重心位于形心前10cm (重心在枪尖与形心之间)，出手时标枪的初始俯仰角速度为 0。在无风的前提下，对标枪投掷出手瞬间、出手后的受力及运 动情况进行分析：简化模型，忽略初始俯仰角

16、速度和重心对标枪投掷的影响；标枪出手瞬间，受到手的推力和重力；标枪出手后，只受重力。由第二问的分析以及第三问题意知，标枪投掷满足模型25.3.2 模型的求解(1)假设某运动员投掷出手速度为加.7。市/5、出手角为36.“、初始攻角为一。4 ,由第二问的结论dist 2115 94.97v 221.10.5152v26.849v7.788 20.003476v32230.01733v20.1229v 2 0.08907 3可以估算出标枪的投掷距离为2） 将初始攻角的区间设为 （-20,10 ） ， 步长设为 1， 建立程序计算得到， （程序见附录），初始攻角为 -7 度时，投掷距离较大，投掷距离

17、为，而用拟合的另一个公式，33.73 0.7881cos0.72671.525sin 0.72671.04cos1.45340.1092sin1.45341.196cos2.18012.07sin 2.18010.51cos2.90680.07145sin 2.90680.6295cos3.63350.07075sin 3.63351.13cos4.36021.743sin 4.3602由初始攻角得到，出手角为度5.4 问题四模型建立与求解5.4.1 模型的建立查阅相关文献，我们知道出手角度为 32 °，（ 1 ） 初始攻角为 0° 6°时,如有逆风 ,标枪飞行远

18、度均呈现出下降的趋势;（ 2 ） 当 初始攻角 >0 °时 ,逆风增大只会引起飞行远度大幅度地下降。（3）在顺风风速为012 m/s,初始攻角为-6° 6°时,标枪的飞行远度随风速增加有所增加 ,尤其是在初始攻角为6°时增加幅度更大。（ 4 ） 但 当风速超过12 m/s 时，飞行远度会急剧下降。出手角度为40 °时，（ 1 ） 初始攻角为-6° 6 时 ,如有逆风,标枪飞行远度均会减小,（ 2 ） 初 始攻角 > 0 °时下降的幅度更大。（ 3） 顺风时,在上述条件下,标枪的飞行远度会有较大幅度增加 ,但风速一

19、旦超过9 m/s12m/s,飞行远度会大大减小通过上面的论述， 我们可以知道， 题目中给出的顺风逆风速度没有让实验结果出现指数级别的误差，并且为我们后面的研究提供了理论和数值基础。在这里，尝试建立人体-标枪运动动力学方程：将标枪视为纯净刚体，利用牛顿欧拉方程建立刚体的动力学方程，作用在标枪上的力有重力以及空气阻力（反作用力）x x(t)x' x'(t)x'' x''(t)x'(t)v出v顺 cos 出x'(t)v出v逆 cos 出dxx x'(t) x''(t)tdtMx''(t) XJ M

20、 Xl sin 攻J Mr2x(t) h l cos 出符号说明：J标枪转动惯量r质点距离转轴，默认 15mm x(t)标枪位移x'(t)标枪速度x''(t)标枪加速度v出出手速度v顺顺风速度v逆逆风速度 出出手角 攻初始攻角M标枪质量l标枪长度(力矩)初始俯仰角速度h运动员高度X风向阻力5.4.2 模型的求解将初速代入上述*II型求解可得 ：顺风逆风3m/s6m/s9m/s|3m/sbm/s最佳出手角(度)最佳初始攻角(度)最佳初始俯仰角速度(度/秒)-8-8-9-9-9-8最大投掷距离(米)5.5问题五模型建立与求解通过对掷标枪运动的研究，主要将得出以下两个方面的结

21、论 一、人体投掷动作对标枪初始条件的影响本文先对飞行系统进行了统计预测，在对人体-标枪系统建立刚体力学模型基础上，通过动力学普遍定理推导了该系统运动的动力学方程，最终得出了各个刚体之间的相互作用力矩 与刚体之间相对角度之间的关系，从而间接建立了各刚体之间相互作用力及力矩与标枪飞行 的初始条件之间的关系。二、初始条件和风对标枪飞行远度的影响在对标枪飞行的计算机模拟以及模拟得到的数据进行分析后，将得出以下几个结论：(1)在其它初始条件不变的情况下，随着初速度的不断增加，标枪的飞行远度也不断增加，而月增加的幅度有加大的趋势。在本文研究的范围内，它们基本上呈线性关系。(2)在其它初始条件不变的情况下，

22、标枪在出手角度为0° 42。的范围内，随着出手角度的增加,标枪的飞行远度基本上是不断增加的。当出手角度为45。时，在空气阻力的影响下 ，并不能使标枪的飞行远度达到最大。当出手角度超过42。以后，标枪的飞行远度有所下降。本文认为该标枪的最佳出手角度大约在 42度左右。(3)在其它初始条件-定的情况下，初始俯仰角速度小于0时标枪要飞的更远一些。认为标枪的最佳初始俯仰角速度大约在 s s范围内。(4)在本文的研究中，初始攻角的变化和出手角存在一定的关系。六、模型评价问题一采用的是 K-means聚类，虽然直观，但是不便于数据的具体统计。问题二利用多元回归统计进行分析，采取的是忽略次要因素，

23、突出主要因素，可能会在结论 上出现些许偏差。问题三在问题二的基础上修缮，在程序上采用步长推断法，导出合适的结论。问题四更加复杂化，在数学模型的基础上，结合物理模型，建立牛顿欧拉方程，求解问题， 结果较为准确。问题五对整个标枪问题进行系统地总结概括。七、参考文献标枪的气动力特性和几何 物理参数测量 蔡国华标枪几何物理参数与空气动力学性能王倩标枪的仰角与阻力王吉寿世界优秀男子标枪运动员最后用力阶段运动学参数分析郑红岩风对标枪飞行初始条件及远度的影响王倩(一般力学与力学基础专业论文)人体一标枪系统对标枪飞行初始条件的影响及标枪飞行的计算机模拟八、附录附录：#include<bits/stdc+.h>using namespace std;int main()double x100;double y100;double z100;double p00 =2115;double p10=;double p01 =;double p20 =;double p11=;do