作者ivs86 (泪雨)
看板MATLAB
标题[绘图]求助 matlab画时变系统
时间Thu Nov 28 23:30:13 2013
小弟是matlab的初学者，
希望功力深的大大可以教一下matlab该如何画时变系统
x1'=-10*x1+e^3t*x2
x2'=-2*x2
我的做法如下：
function dxdt=f(t,x)
dxdt(1)=-10*x(1)+exp(3*t)*x(2);
dxdt(2)=-2*x(2);
dxdt=dxdt';
执行的function
[t,x]=ode45('f',[0 500],[0 10]);
x1=x(:,1);
x2=x(:,2);
plot(t,x1,'b',t,x2,'g');
可是我跑出来的图是收敛，但答案因该是发散才对，
原因因该是我忘记对"t"做定义，但该如何定义t呢?
有听说可以用for回圈写，不需用function~~~不知高手大大可否指导一下
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.34.67.209
1F：推 iamapigtoo: 回一下老文章，应该是要用interp1去辅助 05/31 17:21