最近尝试把之前写过的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、或者是在写矩阵的时候出现问题,
但检查不出到底哪里有错误。
麻烦大家解惑了