[问题] 多项式回归曲线系数

楼主: pencilBen (铅笔)   2021-06-15 14:07:55
最近尝试把之前写过的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、或者是在写矩阵的时候出现问题,
但检查不出到底哪里有错误。
麻烦大家解惑了
作者: lycantrope (阿宽)   2021-06-15 20:08:00
测试很久,是input type问题 x改成Float就OK了在计算poly矩阵,用int64会overflow
楼主: pencilBen (铅笔)   2021-06-16 00:56:00
完全没注意到这个细节…感谢你的帮助!

Links booklink

Contact Us: admin [ a t ] ucptt.com