Matlab：使用匿名函數求解二階ODE

Question

用matlab解決下面的ODE

t = 0:.01:20;
x0 = [0 0]';
xdot = @(t,x) [0 1; 0 0] * x;
ode45(@(t,x) xdot(x(1), x(2)),t,x0)

我收到此錯誤：

使用odearguments時出錯（第91行）@（T，X）XDOT（X（1），X（2））必須返回列向量。

Answer 1

您的x0是2 by 1矩陣，請在函數句柄中對其進行精確調整

xdot = @(t,x) [0 1; 0 0] * x; ---> xdot = @(t,x) [0 1; 0 0] * [x(1);x(2)];

因為xdot已經被定義為一個函數，所以您不需要@在ode45傳遞它

代碼如下

t = 0:.01:20;
x0 = [0 0]';
xdot = @(t,x) [0 1; 0 0] * [x(1);x(2)];
ode45(xdot,t,x0)