我在解微分方程时遇到个郁闷的问题

lanling01

我在解微分方程时遇到个郁闷的问题，怎么改也不行．

方程如下：m2*x''+c*x+k*x-c*x=k*x
               c*X'-c*X1'-N*K2*X1=-N*K2*X2
               m1*X2''+C*X1'-C*X'+(K1+K2)*X2-K2*X=K1*U

代码如下：function dydt=Ruzickafangcheng(t2,y2)
  global  dis col dcol v c k1 k2 m1 m2   N u2  
slong=t2*v;
[u2]=chazhi(slong,dis,col,dcol);
dy(1)=y2(2);
dy(2)=k2*((1+N)*y2(4)-N*y2(3)-y2(1))/m2;
dy(3)=N*k2(y2(3)-y2(4))-c*y2(2)/(-c);
dy(4)=y2(5);
dy(5)=(k1*u2-N*k2*y2(2)+k2*y2(1)+N*k2*y2(3)-(k1+k2)*y2(4))/m1;
dydt=[dy(1);dy(2);dy(3);dy(4);dy(5)];

错误提示如下：
???  In an assignment  A(I) = B, the number of elements in B and
I must be the same.
Error in ==> D:\毕业设计\haiming22\haiming\Ruzickafangcheng.m
On line 7  ==> dy(2)=k2*((1+N)*y2(4)-N*y2(3)-y2(1))/m2;

朋友们帮帮我啊，小弟急用！！

[ 本帖最后由 lanling01 于 2007-5-27 13:37 编辑 ]

无水1324

大小x X有区别吗?程序写得很乱

xjzuo

请把你的调用代码也给出,以便他人调试.

eight

原帖由 lanling01 于 2007-5-27 13:24 发表

                               
登录/注册后可看大图

我在解微分方程时遇到个郁闷的问题，怎么改也不行．

方程如下：m2*x''+c*x+k*x-c*x=k*x
               c*X'-c*X1'-N*K2*X1=-N*K2*X2
               m1*X2''+C*X1'-C*X'+(K1+K2)*X2-K2*X=K1*U

代码如下 ...

1. 建议预分配空间，对此，可看看 置顶贴：聚宝盆 里面提到的精华贴
2. 赋值时，保证左右两个变量是同样大小的。对此，请自己调试一下