使用 SymPy 的指数增长模型

Question

我正在尝试使用 SymPy 实现指数增长模型。

n = Function('N')
t = symbols('t')
r = 1.7

eq = Eq(n(t).diff(t), r*n(t))
sol = dsolve(eq, ics={n(0): 10})
plot(sol.rhs, t, title='Exponential growth', xlabel='Time', ylabel='Population')
t = np.linspace(0,50,5)

这是我的输出：

为什么人口似乎从 0 开始？ 我该如何解决这个问题，例如从 10 开始？

Answer 1

它确实从 10 开始，但指数增长很快变得很大，所以你的地块有一个1e9的垂直比例，这使得 10 看起来非常小。 这就是为什么呈指数增长的事物通常以对数标度绘制的原因。 尝试减少 xlimits 或使用yscale="log"进行绘图以获得更好的视图。