sxgyjq
2011-05-03, 17:03
在编写带有事件的微分方程组的积分问题时,
构造一个事件函数如下
function [value,isterminal,direction] = event1(t,qdq,ddq)
if qdq(1)==1
same_state=0;
else
same_state=1;
end
value=same_state;
isterminal=1;
direction=0
end
调试过程中发现same_state=0,为什么程序没有终止数值积分?
请前辈们不吝赐教,感激不尽!
构造一个事件函数如下
function [value,isterminal,direction] = event1(t,qdq,ddq)
if qdq(1)==1
same_state=0;
else
same_state=1;
end
value=same_state;
isterminal=1;
direction=0
end
调试过程中发现same_state=0,为什么程序没有终止数值积分?
请前辈们不吝赐教,感激不尽!