请教怎样用MATLAB解一个方程

welcome

我们做毕业设计，需要解一个微分方程得出运动曲线，哪位帮忙解一下啊       

m2z2＂+c(z2'-z1')+k(z2-z1)=0
m1z1＂+c(z1'-z2')+k(z1-z2)+k(z1-q)=0

ch_j1985

原帖由 welcome 于 2008-5-21 19:34 发表

                               
登录/注册后可看大图

我们做毕业设计，需要解一个微分方程得出运动曲线，哪位帮忙解一下啊        

m2z2＂+c(z2'-z1')+k(z2-z1)=0
m1z1＂+c(z1'-z2')+k(z1-z2)+k(z1-q)=0

这个有点像车身与车轮两个自由度振动系统的运动方程
还有没有其它的条件？

大混子

不赞同这种求解方式，对自己也没什么好处的！
多学习点论坛版规吧！:@L

yufeng

你应该是只要他们的曲线吧，初始条件是什么？我的qq604443022

sogooda

版面搜索“微分方程”或者到置顶帖里找找微分方程相关的帖子吧

welcome

这个就是一个汽车二自由度振动系统的运动方程
m2z2＂+c(z2'-z1')+k1(z2-z1)=0
m1z1＂+c(z1'-z2')+k1(z1-z2)+k2(z1-q)=0
m1为车轮质量 m2为车身质量 c为减振器阻尼系数 k1为弹簧刚度 k2为轮胎刚度

现在希望能通过MATLAB编程把这个解出来以得到运动曲线，以便进行系统仿真
哪位高手能帮忙啊？

[ 本帖最后由 welcome 于 2008-5-22 09:52 编辑 ]

sigma665

啥都没有
初始条件也没有
各个参数值也没有
让别人怎么帮

另： 强烈建议你搜索论坛

welcome

这个是一个汽车二自由度振动系统的运动方程
m2z2＂+c(z2'-z1')+k1(z2-z1)=0
m1z1＂+c(z1'-z2')+k1(z1-z2)+k2(z1-q)=0
m1为车轮质量 m2为车身质量 c为减振器阻尼系数 k1为弹簧刚度 k2为轮胎刚度
m1=30kg   m2=1500kg  c=980N.s/m   k1=160000N/m  k2=16000N/m q=1

哪位高手能帮忙解出来啊？

ch_j1985

原帖由 welcome 于 2008-5-22 18:57 发表

                               
登录/注册后可看大图

这个是一个汽车二自由度振动系统的运动方程
m2z2＂+c(z2'-z1')+k1(z2-z1)=0
m1z1＂+c(z1'-z2')+k1(z1-z2)+k2(z1-q)=0
m1为车轮质量 m2为车身质量 c为减振器阻尼系数 k1为弹簧刚度 k2为轮胎刚度
m1=30kg   m2=150 ...

我试着编了一段，不知道正确不正确？可能还要一些初始条件吧？
另外，可以参考此帖：
http://www.chinavib.com/forum/viewthread.php?tid=42367&extra=page%3D2%26amp%3Bfilter%3Ddigest

1. 
   
2. %主程序
   
3. clc
   
4. clear
   
5. close all
   
6. global m1 m2 c k1 k2 q
   
7. m1=30;
   
8. m2=1500;
   
9. c=980;
   
10. k1=160000;
    
11. k2=16000;
    
12. q=1;
    
13. [t z]=ode45(@PendulumAbsorber,[0 50],[0 0 0 0]);
    
14. plot(t,z(:,1),'b');
    
15. hold on
    
16. plot(t,z(:,3),'r');

复制代码

%函数
function Q=PendulumAbsorber(t,w)
global m1 m2 c k1 k2 q
A=-c*(w(2)-w(4))-k1*(w(1)-w(3))-k2*(w(1)-q);
B=-c*(w(4)-w(2))-k1*(w(3)-w(1));
x2dot=A/m1;
x4dot=B/m2;
Q=[w(2);x2dot;w(4);x4dot;];
end

[ 本帖最后由 ch_j1985 于 2008-5-22 20:18 编辑 ]