。
“散”
意味着在临界点位置计算出的结果越来越偏离真实值或理论值,甚至无法收敛(得出稳定答案)。
通俗来说,就好像你用地图软件规划路线,结果gps算法崩了,路线越来越偏,导航错得离谱——虽然你车没真的开偏,但系统给出的路线已经完全不可信了。
如果这样的模拟结果被用作控制依据,那火箭后续真的就有可能会飞偏。
顾远舟看着屏幕,语气凝重:“你的意思是,飞行轨迹本身没问题,真正导致偏差的是数值求解在关键时间点附近生了不稳定——结果不是物理误差,而是计算方式的问题?”
姜蕴宁点点头,眼神专注地盯着屏幕上的误差曲线:
“对。
你看这里,37到45秒是推进剂二次燃烧段,系统状态剧烈变化,刚好落在模型的临界过渡区。
我们现在用的是固定时间步长的r-k四阶法,遇到刚性区间,就很容易跑偏。”
她快在键盘上敲入几行代码,将原始解法切换为具备误差控制的自适应步长法。
“我在想,是否需要引入隐式法或者多步法,像g或b这类,处理刚性更稳定。”
顾远舟眼神一亮:“你考虑得很对。
这段区间变化太快,用普通方法根本不够稳定。
不如这样,我们加个局部刚性判断,让程序自动切换计算方法,怎么样?”
姜蕴宁点头,“可以的,顾老师。
而且我建议,把那段误差控制模块做成并行处理。
不然每次局部重算,整个解算器都得卡顿好几秒。”
她顿了顿,补充道:“就连我这种能心算多项式逼近的人,也觉得太慢了,说明它真的拖性能。”
在科研领域,能够快用心算多项式逼近,意味着对函数行为、近似公式非常熟悉,计算度快且准确,有利于快做出判断和调整。
顾远舟闻言挑了挑眉,嘴角浮起一抹若有若无的笑:“能把你逼得喊慢,那确实得改。”
