看板 MATLAB 關於我們 聯絡資訊
小弟是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
iamapigtoo: 回一下老文章,應該是要用interp1去輔助 05/31 17:21