新发现恒星S4716的Kepler轨道拟合：卡方值异常排查请求

问题概述 #

我正在撰写天体物理学论文，尝试为新发现的恒星S4716拟合Kepler轨道。拟合时严格采用了Peißker等人使用的约束条件和物理常数，但计算得到的约化卡方值（reduced chisq）高达206——尽管拟合结果的可视化图表看起来完全正常，从图的匹配程度判断，约化卡方值应该远小于2。

数据细节 #

• 总共有21个观测数据点，其中20个为有效数据点
• 这20个有效点里，有4个点的误差被设为0.5，这4个点主导了整个拟合过程和卡方值的计算

求助方向 #

• 这种异常是否是因为有效数据点数量不足导致的？
• 我的代码里是否存在未被发现的卡方值计算错误？

我已经编写了完整的拟合代码，其中包含拟合结果正常的S2恒星示例作为对照，希望有人能帮忙排查卡方值计算环节的问题。

内容的提问来源于stack exchange，提问作者Frederik Heyde-Petersen