syzeng2002的博客

原创 用scipy的solve_ivp 求解复数常微分方程组

注意代码中只对y1进行了求解，初始条件一定要写成a+bj（或bj）的形式，告诉它你的初值是复数，即使虚部可能为零。如果只填个实数，后面求解会把dydt中取实部。然后算出来就是一条直线。因科研需要，我要求复数odes，首先想python scipy 能不能求，例如dx/dt=ix ,应该能求出振荡解。