zhengshuhe 发表于 2011-3-17 21:04

为什么位移全算出来都为0？

MATLAB里运行一下程序：
tspan=;
y0=zeros(6,1);
=ode45(@zcx3fun,tspan,y0);
plot(t,y);grid;
hold on;
%%
函数.m文件：
function dy=zcx3fun(t,y)
%质量矩阵
M=eye(3,3);
%刚度矩阵
K=;
%阻尼矩阵
C=;
%%
%激振力矩阵
F=
%函数定义
x=y(1:3);
dx=y(4:6);
dy=;
这样的代码，为什么运行结果出来后工作表中变量数值如下：
y(:,1)   y(:,2)       y(:,3)             y(:,4)                     y(:,5)                            y(:,6)
0        0        0        0.00027583        5.044e-006        -4.5144e-008
0        0        0        0.00060967        1.6558e-005        -2.5499e-007
0        0        0        0.0010649        3.8314e-005        -7.9677e-007
0        0        0        0.0016351        7.3123e-005        -1.8962e-006
0        0        0        0.0023141        0.00012347        -3.8438e-006
为什么前面应该是位移，怎么全为0？
后面速度为正余弦，跟解析解接近呢？
？？、、？？？感觉我自己是不是有点理解错误了呢，请高手帮我解答下，谢谢！

appleseed05 发表于 2011-3-18 03:26

你自己想想这句对吗？
dy=;

zhengshuhe 发表于 2011-3-18 08:55

回复 2 # appleseed05 的帖子

dy=;?

zhengshuhe 发表于 2011-3-18 09:08

非常感谢appleseed05在凌晨三点多给我回复的点拨，令人感动！我已经把上面一个例子结合5自由度一个振动微分方程组给算出来了！

meiyongyuandeze 发表于 2011-3-18 09:49

回复 2 # appleseed05 的帖子

哦，太敬业了{:{39}:}

zhouyang664 发表于 2011-3-18 12:16

能给出具体的解决方案就更好了！

appleseed05 发表于 2011-3-18 16:26

回复 3 # zhengshuhe 的帖子

这你是误会我了，我有7个小时的时差{:4_85:}

appleseed05 发表于 2011-3-18 16:29

zhouyang664 发表于 2011-3-18 12:16 static/image/common/back.gif
能给出具体的解决方案就更好了！

咱们不能光给个答案，要启发人家自己解决不是{:4_74:}

zhouyang664 发表于 2011-5-9 22:30

回复 8 # appleseed05 的帖子

我是说楼主将这个帖子总结一下，给出个解决方案，方便后头人！
页: [1]
查看完整版本: 为什么位移全算出来都为0？