+关注
已关注

分类  

暂无分类

标签  

暂无标签

日期归档  

暂无数据

使用complex_ode求解微分方程

发布于2021-03-05 09:37     阅读(93)     评论(0)     点赞(12)     收藏(2)


0

1

2

3

4

5

我有一个尝试使用complex_ode求解的微分方程,但是我的代码未产生与理论结果相符的结果谁能发现我在做什么错?

from scipy.integrate import complex_ode
import numpy as np
import matplotlib.pyplot as plt

def equation(x,t): 
    return 1j*x

sol = complex_ode(equation)
sol.set_initial_value(1)
t_max = 10
dt = 0.1

sol_array = np.array([])
t_array = np.array([])

while sol.successful() and sol.t < t_max:
    sol_array = np.append(sol_array, sol.y)
    t_array = np.append(t_array, sol.t)
    sol.integrate(sol.t+dt)

plt.plot(t_array, np.abs(sol_array))

解决方案


暂无回答

0

1

2

3

4

5

6



所属网站分类: 技术文章 > 问答

作者:黑洞官方问答小能手

链接: https://www.pythonheidong.com/blog/article/879502/c9c925636e5bb298d56d/

来源: python黑洞网

任何形式的转载都请注明出处,如有侵权 一经发现 必将追究其法律责任

12 0
收藏该文
已收藏

评论内容:(最多支持255个字符)