应用MATLAB建模实例---吊门过程仿真.ppt

应用MATLAB建模实例——数学也是一门技术 数学是一门技术 一个例子——足球比赛中的吊门问题 谈谈数学建模竞赛及培训 数学是一门技术 技术的定义 《辞海》：泛指根据生产实践经验和自然科学原理而发展成的各种工艺操作方法与技能；除操作技能外， 广义的还包括相应的生产工具和其他物质设备，以及生产的工艺过程或作业程序、方法。 《科学学辞典》和《科技辞典》：是为社会生产和人类物质文化生活需要服务的，供人类利用和改造自然的物质手段、智能手段和信息手段的总和。 数学及其应用的特征 是一种智能形态的技术 在数学软件的平台上，又表现为一般的物化形态 数学建模技术 指数学及其应用于解决实际问题的整个过程 多样性、合理性、具体问题具体分析 艺术性 足球比赛中的吊门问题 考虑如下的因素：球与球门的距离为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-1 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) 程序1-1之抛射体轨迹函数 function [y,t]=paosheti1(x,alpha,v,g) % 函数返回抛射体轨迹及第一第一次落地飞行时间 y=x*tan(alpha)-x.^2*g/(2*v^2*(cos(alpha))^2); % 不考虑是直接进门还是落地后进门的轨迹 t=2*v*sin(alpha)/g; % 第一次落地前最大飞行时间 xmax=v*cos(alpha)*t; %第一次落地前最大飞行X向距离 n=length(x); for i=1:n if y(i)0 %判断是落地后进门 xx=x(i)-xmax; y(i)=xx*tan(alpha)-xx.^2*g/(2*v^2*(cos(alpha))^2); %计算落地后进门的后半段轨迹 end end 初步结果 对于第一组数据，吊门成功的最小角度1.53697（为弧度，下同），对应的时间大约在4.9281秒，最大角度1.53787，对应的时间是5.0627秒； 对于第二组数据，吊门成功的最小角度1.51437，对应的时间大约在4.1374秒，最大角度1.51587，对应的时间大约在4.2503秒； 对于第三组数据，吊门成功的最小角度1.4571