MATLAB实例分析.pptVIP

仿真模型设计（程序）:从数学模型中可见，小球在空中自由运动时刻与撞击时刻的动力方程不同。通过小球所处位置（位移）是否为零可判定小球处于何种状态。程序中采用if 语句来作出判断，以决定使用（1）还是使用（3）来计算。程序文件代码如下。 g=9.8; v0=0; y0=1; m=1; t0=0; N=5000; dt=0.001; v=v0; y=y0; K=1; for k=1:N if (y0)|(v0) v =v -g*dt; y =y +v*dt; else y =y-K.*v*dt; v =-K.*v-g*dt; end s(k)=y; end t=t0:dt:dt*(N-1); plot(t,s,--k); xlabel(时间t); ylabel(位移y(t)); hold on K=0.85; N=5000; dt=0.001; v=v0; y=y0; for k=1:N if (y0)|(v0) v =v -g*dt; y =y +v*dt; * 对乒乓球的弹跳过程进行仿真 0700103212 王剑 系统过程分析： 忽略空气对球的影响，乒乓球垂直下落，落点为光滑的水平面，乒乓球接触落点立即反弹。如果不考虑弹跳中的能量损耗，则反弹前后的瞬时速率不变，但方向相反。如果考虑撞击损耗，则反弹速率有所降低。 数学模型的建立 首先对乒乓球弹跳过程进行一些理想化假设。设球是刚性的，质量为m，垂直下落。碰击面为水平光滑平面。在理想情况下碰击无能量损耗。如果考虑碰击面损耗，则碰击前后速度方向相反，大小按比例系数K; 0 K ≤ 1 下降。在t 时刻的速度设为v = v(t)，位移设为y = y(t)，并以碰击点为坐标原点，水平方向为坐标横轴建立直角坐标系。球体的速度以竖直向上方向为正方向。重力加速度为 。 初始条件假设：设初始时刻t0 = 0 球体的初始速度为v0 = v(t0)，初始位移为y0 = y(t0)。 受力分析：在空中时小球受重力F = mg 作用,其中，g = -dv/dt。则在t + dt 时刻小球的速度为 在t + dt 时刻小球的位移为 在小球撞击水平面的瞬间，即y(t) = 0 的时刻，它的速度方向改变，大小按比例K 衰减。当K = 1 时，就是无损耗弹跳情况。因此，小球反弹瞬间（t + dt 时刻）的速度为 反弹瞬间的位移为 else y =y-K.*v*dt; v =-K.*v-g*dt; end s(k)=y; end t=t0:dt:dt*(N-1); plot(t,s,-g); xlabel(时间t); ylabel(位移y(t)); hold on K=0.5; N=5000; dt=0.001; v=v0; y=y0; for k=1:N if (y0)|(v0) v =v -g*dt; y =y +v*dt; else y =y-K.*v*dt; v =-K.*v-g*dt; end s(k)=y; end t=t0:dt:dt*(N-1); plot(t,s,:); legend(K=1,K=0.85,K=0.5); xlabel(时间t); ylabel(位移y(t)); axis([0 5 0 1.1]); 图中分别作出了是碰击衰减系数K = 1 ，K = 0.85和K=0.5 三种情况下的小球弹跳位移曲线。 *