|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我想画一个关于微分方程的图,但是感觉画出来的一直有问题,本来应该是周期函数的,但出来的完全不是,郁闷了好几天
代码在如下,请看看有没有什么问题,
function lsf
[t,y]=ode45(@f,[0,5],[0.3;0.7 ]);
plot(t,y(:,2));
% ------------------------------------------------------------
function dydt = f(t,y)
dydt = [ y(1)*(1+0.25*sin(t)-y(1))+0.5*(y(2)-y(1))
y(2)*(0.25*sin(t)-0.25-y(2))+0.5*(y(2)-y(1)) ];
[ 本帖最后由 xinyuxf 于 2007-7-2 15:37 编辑 ] |
|