Question

matlab曲线拟合的问题，做不来，但是急着用！！求助！！！

clear
x=[1986 1996 2009 2011 ];
y=[0.0758 0.0754 0.0758 0.0762];
p=polyfit(x,y,2)
xi=1985:5:2015;
yi=polyval(p,xi);
plot(x,y,'*',xi,yi,'r')
用这个拟合出来有警告：
Warning: Polynomial is badly conditioned. Remove repeated data points
or try centering and scaling as described in HELP POLYFIT.
应该怎样改呢 而且那个y矩阵还是每个元素都扩大了14倍的，，因为如果用实际的话就更有问题了 对于这样的应该怎样做啊 ，可不可以把程序弄出来啊 ，今晚就要。。谢谢！！！

Answer 1

将年份改成自然数，即x=[1986 1996 2009 2011 ];改写成x=[1 10 23 25];

xi=1985:5:2015;改写成xi=1:31;

图形输出可以这样写，

xi=1985:2015;

x=[1986 1996 2009 2011];

plot(x,y,'*',xi,yi,'r'),xlabel('x'),ylabel('y(x)');legend('原始数据曲线','二阶拟合曲线')

追问

请问为什么要这样改啊 改了后我想拟合的曲线的自变量还是原来那个吗

？如果不是不就不对了吗？而且我还需要拟合出来函数的系数哦

追答

因为用年份作自变量，实际上是不合理的，如果9999的话，那算出的值你敢相信吗？

所以一般都以数字序列来作自变量。
在具体出图的时候，可以用年份号来表示。

追问

哦哦 有一点明白了 可是还是做不出来啊 怎样友谊时间序列为自变量 然后出图的时候还用年份号？ 哎 时间不够来了 麻烦你了谢谢了 我还是扩大y吧

追答

你可以把数据给我吗，可以帮你一下。

追问

恩，，，可以啊 太感谢你了！！
x=[1986 1996 2009 2011 ];%就是年份
y1=[0.0758 0.0754 0.0758 0.0762];%这里每一数据都扩大了十四倍，因为如果不扩大就会忽律小数编变成一样的了
y=y1/14;
我想拟合的是 时间关于y的

追答

y=4.15588381534824*10^(-6)*x² -9.54639595021616*10^(-5)*x+ 0.0759041848555914复制去Google翻译翻译结果

追问

恩恩 需要拟合函数和图像 图像叫做“古塔弯曲变形” 拜托了 时间有点紧了

追答

Answer 2

警告可以不用管的，主要是你数据太少，而提示说可以降维（因为高维变量多），但你的数据的趋势决定了没法降维的。
拟合出来基本反映趋势即可！
PS：Y数据扩大是不影响拟合效果的（你可以去尝试，缩小和扩大都是不影响拟合效果的）