这些日子鼻子一直不爽,早上起来后总感觉鼻子堵得慌,老多鼻屎了。。。额。。。
昨天老板指出我的数值计算结果的可疑性之后,今天努力了一天。
早上开始,先是调整昨天使用的积分程序的一些参数,比如区间长度以及区间等分的数目,但是结果基本没什么变化,曲线看着也挺奇怪的。
然后将数值积分的程序换了,用新的算法来计算,并且这次的算法可以控制精度(当精度选取不适当的时候,积分得到的结果就错的一塌糊涂)。换了新算法,发现对于那不通的4项积分,选择误差控制参数的要求也不一样。当我找到合适的误差控制参数后,发现,计算的结果稍有改观,但是和mathematica的结果还存在一些差异。不过有一点好的就是,今天计算得到的曲线,重子振荡的信号比较清楚了。
第一个图是mathematica给出的结果,,第二个图是今天重新计算得到的,与昨天的结果还是差不多,改进就是曲线更光滑了。
这最后两张图表示的是对红移空间功率的修正,分别是昨天和今天的结果
现在还有4个程序在跑,明天才能得到最后结果(4组参数曲面)。
评论