数值计算实习.doc

数值计算 课程实习报告 学 校 西南林业大学 实习时间 2012-12-22 专业 信息与计算科学 班级 信息与计算科学2010级 分组成员 学 号 实习题目名称 足球比赛中的吊门初射角问题 实习条件 硬件环境 Windows XP 软件环境 MATLAB 7.1 参考文献 【1】赵静，但琦，数学建模与数学实验（第2版）。高等教育出版社，2008.1 【2】冉启康 张振宇 张立柱，常用数学软件教程。人民邮电出版社，2008.10 【3】张德丰，数值分析与应用。国防工业出版社，2007.1 【4】郑汉鼎，刁在筠，数学规划[M]。山东，山东教育出版社，1997.12 【5】马正飞，数学计算方法与软件的工程应用。化学工业出版社，2002.12 【6】黄有谦，计算方法。北京，高等教育出版社，1981 实习内容 前言、问题描述、具体理论知识点、具体实例、程序清单、程序实现的核心代码、小组成员分工合作清单。 实习过程 数值结果 出现的问题，解决方法及体会 评价 写实验报告，占 20 ％ 按照教学计划的实验：现场编程序，演示计算结果，占 60 ％ 将计算结果用数学软件作图和分析，占 10 ％ 与教学内容相关、自由选题或参与教师科研的选题，占 10 ％ 2012年 12 月 日 足球比赛中的吊门初射角问题 ——实习内容以及过程 前言：吊射是球员在距离球门很远的一种射门方式，一般是看到对方门将站位靠前或者注意力有所分散时做吊射选择，由于球在空中运行时间长，可以给守门员反应的时间较长，需要精确的脚法才能成功，因此在职业比赛中出现吊的情况很罕见。 h，守门员最大接球高度H，球在空中飞行时间t，球出脚的初速度v，与水平方向的夹角alpha（初射角）。 问题的提出与分析 考虑如下的因素： 球与球门的距离为a，守门员与球门的距离为b，球门高h，守门员最大摸高H，球出脚的初速度为v，与水平方向的夹角为alpha（称为初射角）．给定，h=2.44m，H=3.20m，v=30m/s，重力加速度g=10m/s2 提出问题： 针对下列几组数据分别给出能吊门成功的相应初射角范围，画出最小及最大角度的运动轨迹。要求精度在小数点后至少第3位。(所给数据见计算结果)。 a=6m，b=1m； a=10m，b=3m； a=20m，b=5m。 问题分析: 1.先考虑最简单情形，即不考虑空气阻力等…，此时，球的运动轨迹是抛物线，如果守门员不动，总有合适的角度使吊门成功。 2.这不是求一个角度值，而是求一个范围！通常的思路是把问题整理成两个方程求根问题：一个方程是求吊门成功的最小角度，一个方程是求吊门成功的最大角度。 3.有可能落地弹入球门，要考虑反弹入门的情况。 直观分析: 1.最简单情形，抛射体的运动轨迹为抛物线方程如下 2.借助于使用方便的数学软件，可直观地看到各种初射角对应的抛射体运动的轨迹图形。 最简情形 程序： v=30;g=10; h=2.44;H=3.2;a=6;b=1; l=a-b;L=a*1.01; %1.01表示进入门里 x=0:0.01:L; %对x采样 i1=floor(a/0.01)+1; i2=floor((a-b)/0.01)+1; %取整数 alpha=1.5368:0.00001:1.538; n=length(alpha);%弧度 for i=1:n; [y,tfinal]=paosheti1(x,alpha(i),v,g); tH=l/(v*cos(alpha(i))); %从射门到球到守门员位置的时间 if y(i1)=h | y(i2)=H; success(i)=0; else success(i)=1; end plot(l,H,r+,a,h,r+), hold on, plot(x,y), grid, hold off title([足球比赛中的吊门 ,初射角=,num2str(alpha(i),6) ,... 守门员的移动时间=,num2str(tH)]),pause end figure(2) plot(alpha,success) 程序抛射体轨迹函数（算法）： function [y,t]=paosheti1(x,alpha,v,g) % 函数返回抛射体轨迹及第一第一次落地飞行时间 y=x*tan(alpha)-x.^2*g/(2*v^2*(co