welcome 发表于 2008-5-21 19:34

请教怎样用MATLAB解一个方程

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

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

ch_j1985 发表于 2008-5-21 20:01

原帖由 welcome 于 2008-5-21 19:34 发表 http://www.chinavib.com/forum/images/common/back.gif
我们做毕业设计，需要解一个微分方程得出运动曲线，哪位帮忙解一下啊      

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

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

大混子 发表于 2008-5-21 20:09

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

yufeng 发表于 2008-5-22 08:31

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

sogooda 发表于 2008-5-22 09:08

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

welcome 发表于 2008-5-22 09:50

回复 2楼 的帖子

这个就是一个汽车二自由度振动系统的运动方程
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 发表于 2008-5-22 10:37

回复 6楼 的帖子

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

另： 强烈建议你搜索论坛

welcome 发表于 2008-5-22 18:57

求助：MATLAB求解二阶方程组

这个是一个汽车二自由度振动系统的运动方程
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=1500kgc=980N.s/m   k1=160000N/mk2=16000N/m q=1

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

ch_j1985 发表于 2008-5-22 20:15

原帖由 welcome 于 2008-5-22 18:57 发表 http://www.chinavib.com/forum/images/common/back.gif
这个是一个汽车二自由度振动系统的运动方程
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

%主程序
clc
clear
close all
global m1 m2 c k1 k2 q
m1=30;
m2=1500;
c=980;
k1=160000;
k2=16000;
q=1;
=ode45(@PendulumAbsorber,,);
plot(t,z(:,1),'b');
hold on
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=;
end

[ 本帖最后由 ch_j1985 于 2008-5-22 20:18 编辑 ]
页: [1]
查看完整版本: 请教怎样用MATLAB解一个方程