日月山河 发表于 2023-7-11 16:10:33

求解方程组

大家好:      五个方程组中除了等号左侧的五个变量,其他符号都是常数。这个怎么在fortran中数值求解啊。

唐汉 发表于 2023-7-12 04:30:19

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

日月山河 发表于 2023-7-26 08:35:57

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

变量是Qb,Qc,Pb,Pd,wm.同样的在等号右边也是变量,除了这5个变量,其他的符号都是常数。
页: [1]
查看完整版本: 求解方程组