最近嘗試把之前寫過的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、或者是在寫矩陣的時候出現問題,
但檢查不出到底哪裡有錯誤。
麻煩大家解惑了