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

img

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

你的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)

您好,我是有问必答小助手,您的问题已经有小伙伴帮您解答,感谢您对有问必答的支持与关注!
PS:问答VIP年卡 【限时加赠:IT技术图书免费领】,了解详情>>> https://vip.csdn.net/askvip?utm_source=1146287632