【多项式拟合】用polyfit+poly1d拟合水在高温下的物性参数

前几天导师给我派了个活：
对水在高温（600K-1000K）下的物性曲线（密度、导热系数、粘性、定压比热等）拟合出多项式曲线，并按照规定格式输出多项式系数，方便后续查找和调用。

问题出在水在650K左右的物性会有剧烈变化，导师又要求拟合残差与真实值的比值不高于0.001，拟合阶数不超过五阶，所以我们需要按照温度区间分小段再拟合。

本来用origin也能完成这项工作，但是origin里拟合曲线的参数复制到excel里居然会自动降低精度，是可忍孰不可忍！我怎么可能用手输入！于是只能自己写代码处理。。。

以下就是调用polyfit+poly1d进行拟合+残差分析的代码：

import matplotlib.pyplot as plt
import numpy as np
import pandas as pd

filename = ["fluid635-640", "fluid640-645", "fluid643-649", "fluid649-650",
            "fluid650-650p1", "fluid650p1-650p3", "fluid650p3-650p8", "fluid650p8-651p2",
            "fluid651p1-651p3", "fluid651p3-653",  "fluid653-660",
            "fluid660-680", "fluid680-750", "fluid750-900"]
p_arr = np.array([[0