用ODE45求解微分方程组出现这样的报错应该怎么修改

Question

毕业论文卡在这一步了，有人指点的话，万分感谢。
大概就是一个耦合方程，好像是语句问题，我第一次学习ODE45，不太清楚应该怎么修改。

Answer 1

你的L是向量，这是不允许的，你可以把L定义成和时间相关的量，或者你的L就是一个数，这样的ODE才是可以求解的，比方说

L = @(t) t;
K0 = @(t) pi/2*L(t);
K1 = @(t)2*sqrt(2) *K0(t);
K2 = K1;
f = @(t,x)[i*K1(t)*x(2); i*K1(t)*x(1)+i*K2(t)*x(3); i*K2(t)*x(2)];
[t,y] = ode45(f, [0,1], [1;0;0]);
plot(t,y)

或者你就让L为常数1

L = 1;
K0 = pi/2*L;
K1 = 2*sqrt(2) *K0;
K2 = K1;
f = @(t,x)[i*K1*x(2); i*K1*x(1)+i*K2*x(3); i*K2*x(2)];
[t,y] = ode45(f, [0,1], [1;0;0]);
plot(t,y)