您所在位置：主题：求助 ：关于解一个微分方程组的问题
sccbjj
[专家分：0] 发布于 2007-09-18 10:33:00
dx1(t)/dt=x1+sin(t)*x2;
dx2(t)/dt=x2+cos(t)*x1;
我自己用的方法是这样的：
function F=myfunction(x,t)
F=[x(1)+sin(t)*x(2);x(2)+cos(t)*x(1)];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[t,x]=ode45('myfunction',[0,100],[0 0]);
plot(t,y);
不知道这种解法对吗？书上的例子sin（t）和cos（t）那个位置都是常数，可是我想解得方程组就是上面那种带时间常数t的形式，还能用这种解法吗？
多者大虾们有什么别的解法给小弟讲讲，多谢啦！
dx2(t)/dt=x2+cos(t)*x1;
我自己用的方法是这样的：
function F=myfunction(x,t)
F=[x(1)+sin(t)*x(2);x(2)+cos(t)*x(1)];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[t,x]=ode45('myfunction',[0,100],[0 0]);
plot(t,y);
不知道这种解法对吗？书上的例子sin（t）和cos（t）那个位置都是常数，可是我想解得方程组就是上面那种带时间常数t的形式，还能用这种解法吗？
多者大虾们有什么别的解法给小弟讲讲，多谢啦！
