python scipy.odeint:如何放置任意边界条件

时间:2017-08-21 08:41:34

标签: python scipy odeint

我正在尝试使用scipy.odeint求解二阶常微分方程。

问题如下: f'(x)+ g(f(x))= 0.其中g是f(x)的已知函数。 我想在某个区间[-d,0]知道f(x),最重要的是我想要定义f(0)和f(-d)的边界条件!

然而,使用scipy.odeint我只找到了一种通过定义一组第一阶微分方程来解决这个方程的方法,基本上是: y1 = f(x) y2 = f'(x)

因此: Y1' = y2 Y2' = -g(y1)

据我所知,我形成边界条件的唯一方法是给出以下值: y1(0)= initial_condition_01 y2(0)= initial_condition_02

所以原则上(回到二阶微分方程)我必须为f(0)和f'(0)定义边界条件。

无论如何通过以下方式定义边界条件: f(0)= x0 f(-d)= xd? 那么不要在边界条件中包含一阶导数吗?

0 个答案:

没有答案