求解车桥振动方程用matlab如何写程序？急求，大侠们帮帮忙，明天就要交出稿了

Question

大侠帮我一下好吗？做毕业设计卡住了，不会matlab，能帮我列出下面方程的具体程序吗？急需....
车辆运动方程：M1*z''+c(z'-y'*sin(pi*v*t/l)+k(z-ysin(pi*v*t/l)+u*k(z-y*sin(pi*v*t/l）=0
桥梁振动方程：y''+(EI*pi^4/m*l^4)*y=2/m*l[(M1+M2)g-M1*z''-M2*y''*sin(pi*v*t/l)]*sin(pi*v*t/l)
初始条件t=0;y'0=0;y''0=0;z'0=0;z''0=0;z0=0;y0=0;
参数：m=1.2*(10^4);EI=1.275*10^11;M1=4.69*10^4;M2=1.69*10^4;
l=50;g=9.8;c=3.14*10^5;k=4.87*10^6;u=1.
做出（y。t）和M1,M2与荷载位置vt/l的响应图。
我的邮箱：673860305@qq.com
谢谢！

Answer 1

你的 v是什么?
补充一下条件吧

还有你的式子括号根本就不匹配,怎么解?
同志,认真点吧!

这是大概解法 自己再修改修改吧

m=1.2*(10^4);
EI=1.275*10^11;
M1=4.69*10^4;
M2=1.69*10^4;
l=50;
g=9.8;
c=3.14*10^5;
k=4.87*10^6;
u=1;
v=1;
dsolve('M1*D2z+c*(Dz-Dy*sin(pi*v*t/l))=-k*(z-y*sin(pi*v*t/l))-u*k*(z-y*sin(pi*v*t/l))',...
'D2y+2/m*l*(M1*D2z+M2*D2y*sin(pi*v*t/l))*sin(pi*v*t/l)=2/m*l*(M1+M2)*g*sin(pi*v*t/l)-(EI*pi^4/m*l^4)*y',...
'Dy(0)=0','D2y(0)=0','Dz(0)=0','D2z(0)=0','z(0)=0','y(0)=0')

祝君好运