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

作者: 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