线性微分方程组 矩阵怎么传参呀？？

jingjing122

请问一下，比如说：
dx/dt=x+2y
dy/dt=6y

我看各种教程都是把这两个式子分着像这样写的，但如果我的式子很多，有10个，这个时候用矩阵的形式就很方便，[x']=A[x] 这样， 可是不知道该怎么表示，是不是还用odeint函数呢？？

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

# function that returns dy/dt
def model(z,t):
    x=z[0]
    y=z[0]
    k = 0.3
    dxdt = x+2* y
    dydt = 6* y
    return [dxdt,dydt]

# initial condition
x0=1
y0=1

# time points
t = np.linspace(0,20)

# solve ODE
z = odeint(model,z0,t)


那如果不想这样一个一个式子写，想用矩阵写，该怎么写啊
def model(z,t):
        A=np.array([[1,2],[0,6])
        dxdt=np.dot(A,x)
        return dxdt
z0=[1,1]

不好意思哇。好像是可以的。。

jingjing122

import numpy as np
from scipy.integrate import odeint

def model(t,z):
    L = np.array([[1,1],[0,1]])
    dzdt = L*z
    return dzdt

z0=[1,1]

a=np.shape(z0)
print(a)
t= np.linspace(0,5,20)
z=odeint(z0,model,t)
#print (z)

会报错：error: The function and its Jacobian must be callable functions.
这里的Jacobian matrix是什么啊？如何求呢？

jingjing122

解决啦解决啦