请教如何向ode中传递变量值
一般来说,要求解微分方程,首先定义一个方程,然后调用ode(例如ode45)就可以了。例如,我要求解Duffing方程,首先定义微分方程:function du=duffin(t,x,p)
wd=1;
f=1;
dy(1)=y(2);
dy(2)=-k*y(2)+y(1)-y(1)^3+f*cos(wd*t)+p;
然后调用
=ode45('duffin',,)
就可以得到不同时刻t对应的y的解。
现在我的问题是:在每一次迭代计算中,需要p是一个变量,而不是常量,那么应该如何传递参数p呢?也就是说,在每一次迭代过程中p的取值都要更新一次,请问如何解决,谢谢解答啊~
[ 本帖最后由 wannete 于 2009-8-20 23:00 编辑 ] 就是duffing方程添加高斯白噪声的问题,怎么没人回答啊!:@( :'(
页:
[1]