Table of Contents
建模和模拟工程系统是工程分析中必不可少的任务. 利用NumPy和SciPy等Python库通过提供强大的数字计算和系统分析工具来简化这些过程. 本条引入了与这些库一起建模和模拟工程系统的基本方法.
建立环境
要开始, 需要的库如果还没有的话, 则安装它们。 使用 pip 命令安装 NumPy 和 SciPy :
pip 安装 numby scipy
建模工程系统
工程系统往往可以使用微分方程或传输函数等数学模型来表示. NumPy提供了创建矩阵和向量的工具,而SciPy则提供了解方程和模拟系统响应的功能.
例如,一个简单的质量-泉水-雨量系统可以用微分方程来建模:
m * x' + c * x' + k * x = 0
模拟系统
SciPy 的集成模块允许在数字上解析微分方程。 通常用于此目的的“ solve ivp” 函数。 以下是模拟质量- 泉水- 遮挡系统的例子 :
将系统定义为函数 :
进口数字作为np]
从 scipy.intrate 导入解析 ivp ]
def assum spring damper(t,y,c,k):]
返回 [y[1], -(c/m) * y[1] - (k/m) * y[0]
设定初始条件和参数:
初始条件=[x0,v0]]
] t span = [0, 10]]
唤解者曰:
溶解=溶解 ivp(质量 泉水 damper,t span,初始 条件,args=(c,k))]
分析结果
模拟后, 分析结果, 绘制系统随时间推移的反应。 使用 Matplotlib 等库可视化数据 :
将 matplotlib.pyplot 导入为 plt
plt.plot(溶液.t,溶液.y[0])]
plt.xlabel('Time (s)') ] 互联网档案馆的存檔,存档日期2013-09-02.
plt.ylabel('Displacement')]] 互联网档案馆的存檔,存档日期2013-09-02.
plt.show()]