Question

我正在研究很长的微分方程列表，我最终需要将其与实验数据相匹配，但我在让NDSolve接受我的下标变量时遇到了问题。由于方程式非常长，我只是简单地给出了一个样本，但是有人可以指导我进一步帮助吗？我已经发现我需要使用Symbolize和Patten - 还有更多步骤吗？此外，我为糟糕的形式道歉，复制和粘贴采取输入形式而不是漂亮的形式。如果屏幕截图更有帮助，请告诉我。谢谢！

<< Notation`

Symbolize[ParsedBoxWrapper[SubscriptBox["h", "_"]]]
Symbolize[ParsedBoxWrapper[SubscriptBox["R", "_"]]]
Symbolize[ParsedBoxWrapper[SubscriptBox["T", "_"]]]
Symbolize[ParsedBoxWrapper[SubscriptBox["\[Lambda]", "_"]]] 

Pattern[Subscript[h, \[Placeholder]], _]
Pattern[Subscript[R, \[Placeholder]], _]
Pattern[Subscript[T, \[Placeholder]], _]
Pattern[Subscript[\[Lambda], \[Placeholder]], _]

model[Subscript[h, 1] _?NumericQ, ...] :=

NDSolve[{Derivative[1][g1][t] == 
Subscript[R, 1]/(
 1 + E^(Subscript[h, 1] + g13[t] Subscript[T, 1, 13] + 
   g15[t] Subscript[T, 1, 15] + g16[t] Subscript[T, 1, 16] + 
   g22[t] Subscript[T, 1, 22] + g24[t] Subscript[T, 1, 24] + 
   g25[t] Subscript[T, 1, 25])) - g1[t] Subscript[\[Lambda], 1], ...,

g1[0] == 532.031,...}, {g1, g2, g3, g4, g5, g6, g7,
g8, g9, g10, g11, g12, g13, g14, g15, g16, g17, g18, g19, g20, 
g21, g22, g23, g24, g25, g26, g27, g28, g29}, {t, 0, 4}]

Mathematica NDSolve中的订阅参数

0 个答案: