ode数值求解中的,'NonNegative',使用问题。我仔细阅读完了帮助文件,并多次尝试,仍然搞不太明白。求教大家。
就拿matlab自带的例子来说。其中有个脚本m文件 rigidode。
如果在调用ode45前面加一行程序,“options=odeset('NonNegative',[1 2 3]);”强制y的每1,2,3个结果都非负。
但是运行结果仍然和设置前一样。y(2),y(3)的结果都有正有负。请问我错在哪儿里。
“,”后面的矩阵中的数字,代表分别将第几个分量 强制为非负吧。想让第2,3个分量为非负值。则那个矩阵应该是[2 3]吧。
同时odeset('event',@eventfunction)中,
function [value,istermina,direction]=eventfunction(t,y)
这个函数的三个参量都什么意思啊?怎么使用event函数啊?
有没有 比较好的例证?帮助中的那个弹球的例子,看不懂。 |