[問題] 多項式回歸曲線係數

[pencilBen (鉛筆)]

最近嘗試把之前寫過的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://www.ptt.cc/bbs/Python/M.1623737277.A.841.html
推 lycantrope: 測試很久，是input type問題 x改成Float就OK了
推 lycantrope: 在計算poly矩陣，用int64會overflow
→ pencilBen: 完全沒注意到這個細節…感謝你的幫助！