声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1525|回复: 3

[编程技巧] 我在解微分方程时遇到个郁闷的问题

[复制链接]
发表于 2007-5-27 13:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

方程如下: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 编辑 ]
回复
分享到:

使用道具 举报

发表于 2007-5-27 15:13 | 显示全部楼层
大小x X有区别吗?程序写得很乱
发表于 2007-5-27 15:38 | 显示全部楼层
请把你的调用代码也给出,以便他人调试.
发表于 2007-5-27 22:00 | 显示全部楼层
原帖由 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. 赋值时,保证左右两个变量是同样大小的。对此,请自己调试一下
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-12 14:32 , Processed in 0.073260 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表