为什么我用MATLAB多项式拟合得到的函数,再用函数去作图时两者的曲线会不一样呢?

为什么我用MATLAB多项式拟合得到的函数,再用函数去作图时两者的曲线会不一样呢?
clear all
x0=[0 0.5 1 1.5 2 2.5 3 7 8 9 10 11.82 13 14 15 16 17.73 19 20 21 22 23.64 26 27 28 29 29.55];
y0=[1.892e9 1.89e9 1.87e9 1.86e9 1.85e9 1.84e9 1.83e9 3e5 3e5 3.1e5 3e5 2.9e5 2.4e5 2.3e5 2e5 2.2e5 2e5 1.85e5 1.8e5 1.78e5 1.75e5 1.7e5 1.62e5 1.65e5 1.6e5 1.55e5 1.5e5];
p3=polyfit(x0,y0,18);
vpa(poly2sym(p3),5)
x=0:0.1:29.6;
y1=polyval(p3,x);
>> plot(x,y1,x0,y0,'o')
grid on
得到如下函数：
y=1.604*10^(-9)*x^18 - 4.4275*10^(-7)*x^17 + 0.000056159*x^16 - 0.0043377*x^15 + 0.22791*x^14 - 8.6173*x^13 + 241.87*x^12 - 5123.9*x^11 + 82462.0*x^10 - 1.0064*10^6*x^9 + 9.2252*10^6*x^8 - 6.2331*10^7*x^7 + 3.0125*10^8*x^6 - 9.9516*10^8*x^5 + 2.1031*10^9*x^4 - 2.5932*10^9*x^3 + 1.6075*10^9*x^2 - 3.7944*10^8*x + 1.8927*10^9
我再输入
clear all
x=0:0.1:30;
y=1.604e-9*x.^18 - 4.4275e-7*x.^17 + 0.000056159*x.^16 - 0.0043377*x.^15 + 0.22791*x.^14 - 8.6173*x.^13 + 241.87*x.^12 - 5123.9*x.^11 + 82462.0*x.^10 - 1.0064e6*x.^9 + 9.2252e6*x.^8 - 6.2331e7*x.^7 + 3.0125e8*x.^6 - 9.9516e8*x.^5 + 2.1031e9*x.^4 - 2.5932e9*x.^3 + 1.6075e9*x.^2 - 3.7944e8*x+ 1.8927e9;
plot(x,y)
为什么曲线会不一样,我仔细看过只有曲线前一部分会是相同的大概在x小于5.5内是相同的,之后就完全不一样了.请教给位了.实在不明白为什么了.
主要是这个得到的曲线，一部分小于0，我是希望把常数项再变大一些，使得曲线完全在x轴上方。哪知道画出来的曲线完全不是那么回事。

clear all
x0=[0 0.5 1 1.5 2 2.5 3 7 8 9 10 11.82 13 14 15 16 17.73 19 20 21 22 23.64 26 27 28 29 29.55];
y0=[1.892e9 1.89e9 1.87e9 1.86e9 1.85e9 1.84e9 1.83e9 3e5 3e5 3.1e5 3e5 2.9e5 2.4e5 2.3e5 2e5 2.2e5 2e5 1.85e5 1.8e5 1.78e5 1.75e5 1.7e5 1.62e5 1.65e5 1.6e5 1.55e5 1.5e5];
p3=polyfit(x0,y0,18);
x1=0:0.1:29.6;
y1=polyval(p3,x1);
figure(1)
plot(x1,y1,x0,y0,'o')
grid on
figure(2)
x=0:0.1:30;
yy=poly2sym(p3);
y2=eval(yy);
plot(x,y2)
grid on
楼主你的问题可能出在精度问题上,因为你的18阶拟合本身拟合效果就很不好,而你用这句
vpa(poly2sym(p3),5)
是他的系数精度变低,也就导致后面的图形不一样,用我上面的代码,结果是一样的,
matlab的精度问题很重要,尤其在高阶拟合中