唐汉 发表于 2023-7-12 04:30
我只是很好奇，如果按照你说的等式左边的五个量都是变量，其余都是常数。那么公式（1）的Qb不也是常数么？ ...

变量是Qb,Qc,Pb,Pd,wm.同样的在等号右边也是变量，除了这5个变量，其他的符号都是常数。

我只是很好奇，如果按照你说的等式左边的五个量都是变量，其余都是常数。那么公式（1）的Qb不也是常数么？如果Qb是常数，那公式（2）中Pb不就是个关于时间的线性函数么（假设微元是t）？而且这样一说公式（3），（4）也是个常数？
我觉得你可能要解决的问题：
（1）搞清楚这个方程的变量是什么；
（2）把这几个方程转化成一个线性独立的系统，按照你上面这几个方程，我觉得这个系统很有可能就是一个二元或者三元的常微分方程组；
（3）如果这个微分方程组的右边没有“奇奇怪怪”的非线性项，这个方程有可能是有解析解或者是有非常接近方程解的解析式的。你找找看，因为这对于你确定这个方程的初值/边界条件/是否收敛是有帮助的；
（3）最后，去找一个常微分方程的求解器（无论是不是Fortran的）求数值解就行了。你在第三步中找到的条件会帮你加快求解器迭代的次数。