问题描述
我正在使用描述分子荧光发射的动力学模型。
我能够通过实验测量四个参数：两个寿命（τ1，τ2），荧光量子产率（ϕf）和辐射率（kr）。
我的模型包含三个未知比率，kMR，kRM和knr。我有一组包含所有这些值的三个方程，我想使用sympy求解三个未知数。
代码如下：
from sympy import *
kr,k1,k2,phi,kMR,kRM,knr = symbols('kr k1 k2 phi kMR kRM knr',real=True)
#kr = 0.00014
#k1 = 1/9
#k2 = 1/49
#phi= 0.005
Phi = kr/(kr+kMR-kMR*kRM/(kRM+knr))
X = kr + kMR
Y = kRM + knr
K1 = (X+Y+sqrt(X**2-2*X*Y+Y**2+4*kMR*kRM))/2
K2 = (X+Y-sqrt(X**2-2*X*Y+Y**2+4*kMR*kRM))/2
solutions = solve([K1-k1,K2-k2,Phi-phi],(kMR,knr))
print(solutions)
如果我不注释测量值，则将在几秒钟内找到一个数值解。但是，我的测量容易出错，因此我有兴趣探索解决方案的整个领域，以了解模型对每个参数的敏感程度。我也有来自不同实验的不同测量值。因此，我想根据测量值获得kMR，kRM和knr的符号表达式。不幸的是，如果运行此命令，它将无法收敛。
您能帮助我获得我正在寻找的符号解决方案吗？
谢谢！
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）