Matlab的ode45，In an assignment A(I) = B, the number of elements in B and

ukman

ode45 解一运动微分方程，
ic=[0,0]; tspan=[0,20]; [T,X]=ode45(@sigs,tspan,ic);
运行后出现In an assignment  A(I) = B, the number of elements in B and I must be the same.
Error in ==> sigs at 11
dq(2)=(k*y+c*dy-k*q(1)-c*q(2))/m;

估计问题出在y 和dy，因为把chirp换成sin和cos后就可以解，那用chirp的问题出在哪里，或是其他什么原因？
目标函数如下：
function dq=sigs(t,q)
%function to calculate time history of response
m=60;k=44130;c=1485;
% generate input
t = 0:0.001:2;
y = chirp(t,0,1,15); dy=[diff(y),1.001];
%q(1)=displacement;q(2)=velocity
dq=zeros(2,1); dq(1)=q(2); dq(2)=(k*y+c*dy-k*q(1)-c*q(2))/m;
end

[ 本帖最后由 ChaChing 于 2009-7-5 19:19 编辑 ]

无水1324

function dq=sigs(t,q)
%function to calculate time history of response
m=60;k=44130;c=1485;
% generate input
t = 0:0.001:2;           
y = chirp(t,0,1,15);
dy=[diff(y),1.001];
%q(1)=displacement;q(2)=velocity
dq=zeros(2,1);
dq(1)=q(2);
dq(2)=(k*y+c*dy-k*q(1)-c*q(2))/m;
end


这里面不能再设置t了，因为ode默认的就是时间t，而且你这两个时间t都是不同的区域内的，我不知道你里面的这个y是什么意图，这样子写是不是本身就有问题呢

ChaChing

Ref
13F of 常见的程序出错问题整理 (eight)
http://forum.vibunion.com/forum/thread-46001-1-1.html

还有勿重覆发帖!

[ 本帖最后由 ChaChing 于 2009-5-3 22:52 编辑 ]

ukman

这个能看懂，可是我这个问题就想不通，难道说右边是一个时间序列，左边只是一个值吗
恩，我把t = 0:0.001:2去掉就没有错误，当然结果对还是错再说。这里的y是位移激励，dy是对应的速度
为什么chirp产生的时间只有两秒，怎么设置更长的时间？

[ 本帖最后由 ChaChing 于 2009-7-5 19:21 编辑 ]

ChaChing

建议仔细看下chirp帮助文献!
t = 0:0.001:2; 当然只有两秒

ChaChing

原帖由 ukman 于 2009-5-3 23:08 发表

                               
登录/注册后可看大图

...难道说右边是一个时间序列，左边只是一个值吗

这个无需问别人, 自个儿设断点或display结果看看即可!
dq(1)=q(2);
zzz=(k*y+c*dy-k*q(1)-c*q(2))/m; size(zzz)
dq(2)=(k*y+c*dy-k*q(1)-c*q(2))/m;

无水1324

将t = 0:0.001:2去掉之后，y = chirp(t,0,1,15);这里面的t就是仿真的时间t了，所以我觉得你有必要把你的方程给出来，感觉你编程有点问题

ukman

方程很简单，如下。
mq''+k*(q-y)+c*(q'-y')=0

无水1324

那么y是一个怎么样的函数呢？

ukman

y是激励，自己定义的。