-G2/2*cos(theta2)]?X = A\B; %用左除求解线性方程组Р程序exn712运行结果Р这样求解的方法不仅适用于全部静力学题目,而且可用于材料力学和结构力学中的超静定问题。因为那里只多了几个形变变量和变形协调方程,通常也是线性的,所以只不过是把矩阵方程扩大了几阶,解法没有什么差别。Р例7-1-3 质点运动学Р设导弹M速度为vm=800m/s,其速度向量始终对准速度为vt=300m/s的直线飞行目标T,发射点在目标运动方向的左(4000m)前(3000m)方,试求导弹轨迹及其加速度。?解:◆建模: 在与目标固连的等速直线运动坐标中(因而是惯性坐标系)列写动点M的方程。因动坐标与目标T固连,牵连速度。动点为M,它的绝对速度。由速度合成定理,相对速度,列出它在x,y两个方向的投影,得Р主程序exn713Рglobal vt vm?vt=input('vt=');vm=input('vm=');?%输入共用的参数?z0=input('[x0;y0]='); %输入数值积分需要的参数?tspan=input('tspan=[t0,tfinal]=');?%?[t,z] = ODE23('ex713f',tspan,z0);?%进行数值积分?plot(z(:,1),z(:,2)); %绘图?% M点位置的导数是相对速度,二次导数则为绝对加速度?dt=diff(t); Ldt=length(dt); %为了求导数先求t的增量?x=z(:,1);y=z(:,2); % 把z写成x,y两个分量形式?vx=diff(z(:,1))./dt;vy=diff(z(:,2))./dt; %?wx=diff(vx)./dt(1:Ldt-1);wy=diff(vy)./dt(1:Ldt-1); %二次导数?[t(2:Ldt),x(2:Ldt),y(2:Ldt),wx,wy] %显示数据
全文预览
