关于`DifferentialEquation.jl`计算变量是矩阵的常微分方程

我是一个Julia新手,最近我准备使用DifferentialEquation.jl来求解我遇到的Linbland Master Equation(林布拉德主方程),方程的形式如下:
d607f48cc3e1b2643b4651a5a1896f7
其中ρ是密度矩阵算符,H是哈密顿量,L(ρ)是刘维尔超算符。
正如你所看到的,主方程的变量都是矩阵或者高阶张量。
下面的代码是我网上看到的一段使用DifferentialEquation.jl求解方程的事例:

using DifferentialEquations 
using Plots; gr()
f(u,p,t) = 0.98u            #新建一个待求方程
u0 = 1.0                    #初值
tspan = (0.0,1.0)           #时间范围
#新建一个常微分方程问题
prob = ODEProblem(f,u0,tspan)
sol = solve(prob)
#画图
plot(sol)
savefig("ode_1.png")

它计算的是非常简单的一个常微分线性方程,变量都是一个数字。
我的问题是如果我建立一个如下的方程:

LME(rho,t)= -i*(H*rho-rho*H) + L_rho

DifferentialEquation.jl可以正确计算出结果吗?
如果不可以请问Julia中是否有其他库可以解决我的问题,请各位指教。

我最近在学这个 DifferentialEquations.jl 我们可以一起探讨,不过有个问题
你这个函数要怎么表示出来,

LME(rho,t)= -i*(Hrho-rhoH) + L_rho

这个不是个我能看懂的函数,请你重新写一遍

你好,这个函数就是按照林布拉德主方程的形式写出来,其中的变量都会对应定义一个张量,带入方程即可

https://docs.sciml.ai/DiffEqDocs/stable/getting_started/#Example-2:-Solving-Systems-of-Equations

image

是的,这个文档非常详细的介绍了DifferentialEquation.jl

是的,这个回答是一针见血的,我也意识到了计算速度的问题。我会先读一下这篇文章。 1312.0111.pdf (arxiv.org)