作者pencilBen (铅笔)
看板Python
标题[问题] 多项式回归曲线系数
时间Tue Jun 15 14:07:55 2021
最近尝试把之前写过的matlab code转换成python,其中一个尝试如下所示:
https://github.com/Frankkao813/Numerical-Analysis/blob/main/PolynominalFit
里面的大致步骤为:
1)汇入资料
2)写出一个任意order的polynominalFit矩阵,矩阵的公式如下所示。
https://i.imgur.com/gyLpxvc.jpg
3)解完矩阵之後再将原始资料跟回归曲线叠在一起。
目前遇到的问题是看起来matlab解出来的回归曲线系数跟python似乎不太一样,且用n>5
时会出现奇怪的结果。怀疑是在用np.linalg.solve、或者是在写矩阵的时候出现问题,
但检查不出到底哪里有错误。
麻烦大家解惑了
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 223.139.55.148 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Python/M.1623737277.A.841.html
1F:推 lycantrope: 测试很久,是input type问题 x改成Float就OK了 06/15 20:08
2F:推 lycantrope: 在计算poly矩阵,用int64会overflow 06/15 20:17
3F:→ pencilBen: 完全没注意到这个细节…感谢你的帮助! 06/16 00:56