matlab数学建模有组实验数据为0 mm 8.7mv5 mm 8.9mv 10 mm 9.1mv 15 mm 9.3m

这个很难办哦,但是一般来说大部分都是先拟合多项式的,那我们试试多项式的结果
clc
clear
data=[0 8.7
5 8.9
10 9.1
15 9.3
20 9.5
25 9.7
30 9.9
35 10.1
40 10.3
45 10.6
50 10.9
55 11.1
60 11.4
65 11.7
70 12.0
75 12.4
80 12.8
85 13.2
90 13.7
95 14.2
100 15.1
105 15.8
110 16.8
115 18.3
120 20.2 ];
x=data(:,2);
y=data(:,1);
xi=linspace(min(x),max(x),100);
for n=1:5;
p{n}=polyfit(x,y,n);
yi(n,:)=polyval(p,xi);
end
plot(x,y,'o',xi',yi')
legend('n=1','n=2','n=3','n=4','n=5')
从图形中你可以看出几次多项式拟合比较好,我看了下,应该是n=4的效果做好
对应参数为
p{4}
>> vpa(poly2sym(ans),3)
ans =
.195e-3*x^4+.702e-1*x^3-4.41*x^2+88.9*x-488.