基于MATLAB抛体运动的分析.docx

基于 MATLAB 抛体运动的分析作者：昊学长 yv y v0y O v0 x v0x 一、理论分析 我们研究的抛体运动为斜抛运动，设物体斜抛的初速度为v0 ，重力加速度为 g，在无阻尼抛体运动中，我们不计空气阻力， 在无阻尼抛体运动中，物体的斜抛运动是水平方向的匀速直线运动和竖直上抛运动的合运动，物体的运动方程为 0?x ? v cos?? t 0 ? ? y ? v sin?? t ? 1 gt2 消除 t 可得轨迹方程 ?? 0 2 0gx2 0 令 y ? 0 ，可得水平射程 y ? x tan?? 2 v2 cos2? v2 sin 2? X ? 0 g 在有阻尼抛体运动中，设阻尼系数为 b，物体质量为 m，我们近似认 为阻力 F 只与速度 v 的量值成正比，即 F ? ?bv 抛体在xy 平面内运动，其运动的微分方程为 ?m dvx ? ?bv ? dt ??? dvy m ? ? ? dt x ? ?mg  bvy 设当t ? 0 ， vx ? vx0 ， vy ? vy0 ，对上式进行一次积分后，得到 ??x? ? ? vx ? vx0 ?bt e m ? y? ? v ? ? mg ? v ?e?bt mg ? y ? ? ? b y0 ? m ? b 再进行一次积分，并设t ? 0， x ? 0 ， y ? 0 ，则得 ?m?x ? mvx0 ??1? e?bt ?? ? m ? b ? ? ? 2 ? mv ? bt ????? y ? ? m g ? y0 ???1? e? m ?? ? mg t ? ? ? ? ? ? b2 b ? ? b 消除 t，得到轨道方程为 ? mg vy0 ? m2 g ? mvx0 ? y ? ? ? bvx0 ? ?x ? vx0 ? b2 ln? ? mvx0 ? bx ? ? vy0 ? v0 sin?, vx0 ? v0 cos? ? vy0 vx0 ? v0 sin? ? tan? v0 cos? 即 ? mg  ? m2 g ?  mvx0 ? y ? ? ? bvx0 ? tan???x ? 2 b? b ln? ?mvx0 ? mv bx ? ?二、程序设计 ? （1）初速度v0 相同，θ不同，有阻尼与无阻尼的比较 clear v0=10; g=9.8; m=5; b=0.2; theta=15:15:75; th=theta*pi/180; x1=linspace(0,11,1000); x2=linspace(0,11,1000); [TH,X1]=meshgrid(th,x1); [TH2,X2]=meshgrid(th,x2); Y1=X1.*tan(TH)-g*X1.^2/v0^2/2./cos(TH).^2; Y2=(m.*g./b./v0./cos(TH2)+tan(TH2)).*X2-m.^2.*g./b.^2.*log(m.*v0.*cos(TH)./… (m.*v0.*cos(TH)-b.*X2)); plot(x1,Y1(:,1),x2,Y2(:,1)) hold on plot(x1,Y1(:,2),x2,Y2(:,2)) hold on plot(x1,Y1(:,3),x2,Y2(:,3)) hold on plot(x1,Y1(:,4),x2,Y2(:,4)) hold on plot(x1,Y1(:,5),x2,Y2(:,5)) hold on axis([0 11 0 5]) title(有阻力与无阻力的抛体运动，初速度相同θ不同) xlabel(x/m) ylabel(y/m) grid on 5 4.5 4 3.5 3 y/m2.5 y/m 2 1.5 1 0.5 0 有 阻 力 与 无 阻 力 的 抛 体 运 动 ， 初 速 度 相 同 θ不 同 0 1 2 3 4 5 6 7 8 9 10 11 x/m （2）初速度θ相同， v0 不同，有阻尼与无阻尼的比较 clear theta=45; g=9.8; m=5; b=0.2; v0=10:10:50; th=theta*pi/180; x1=linspace(0,260,1000); x2=linspace(0,220,1000); [V0,X1]=meshgrid(v0,x1); [V02,X2]=meshgrid(v0,x2); Y1=X1.*tan(th)-g*(X1.^2)./V0.^2/2/cos(th).^2; Y2=(m.*g./b./V02./cos(th)+tan(th)).*X2-m.^2.*g./b.^2.*log(m.*V02.*cos(th)./… (m.*V02.*co