Python 微分方程 #包括微分方程和参数 def型号（x，dydx，p）： s=p[0]#第一个参数 a=p[1]#第二个参数 dydx[0]=-2*（s+a）*y[0]+2*s*y[1]+s/2*y[2] dydx[1]=+2*（s+a）*y[1]-2*s*y[0]-s/2*y[2] dydx[2]=-（s+a）*y[2] 返回np.vstack（dydx[0]，dydx[1]，dydx[2]） #边界条件 def bc（ya、yb、yc、p）： s=p[0] a=p[1] I0=1 返回np.array（（[ya[0]，yb[0]，yc[0]，a，s]）） #x值 x=np.数组（[1,2,3,4,5,6,7,8,9]） y=np.零（（3，x.尺寸）） #p=np.零（（3，d.尺寸）） #y初始值 y[0]=3 y[1]=3 y[2]=2 I0=1 sol=solve_bvp（模型，bc，x，y，p=[1,1]）

我不知道如何写边界条件来解这三个微分方程

---

我想解方程，得到y值和参数值，另外一个问题，你有固定的参数和3个边界条件。这需要编码为

#边界条件
def bc（y0，yd，I0）：
返回np.数组（[y0[0]，y0[2]-I0，yd[1]]）

然后初始的

x

数组需要由边界点来限定

x=np.linspace（0，d，9）
y=np.零（（3，x.尺寸））
#y初始值
y[0]=3
y[1]=3
y[2]=2

需要在不使用可变参数的情况下调用解算器，通过包装器/部分求值将它们全部固定为常量值

sol=solve_bvp（λt，y:model（t，y[s，a]），λy0，yd:bc（y0，yd，I0），x，y）

---

什么是边界条件

solve_bvp

不处理两个区间端点以外的条件。边值问题是我的问题。我有这三个微分方程，我想一起求解，得到y0 y1 y2，参数s和pyu应该有5个条件或方程，在值

y_0（1），y_1（1），y_2（1）

和

y_0（9）、y_1（9）、y_2（9）

。在最简单的情况下，这会固定这6个数字中5个的值。然后，您可以应用

solve_bvp

（同时修复

模型

函数声明之后）。这里是一个matlab代码，我正试图将其转换为pythonThanks用于解决方案，但如果参数是可变的[s，a]，边界条件和解算器将如何？Thanks in advanceAs之前说过，您有#（边界条件）=#（状态变量）+#（未知参数），因此需要指定5个边界条件。ValueError:

bc

return应该有形状（3），但实际有（5，）.当我想将参数a和s估计为非常数的依赖参数时，我面临着这个错误。我在这里上传了代码。谢谢你的帮助