程序如下:附件中data.txt必须洞戚和程序放在同一文件夹中。
load data.txty=data(:,1)
x=1:1:length(y)
a1=polyfit(x',y,5)
syms X
f3=vpa(poly2sym(a1,X),4)%多项式5次旅岁拟合结果
x1=1:0.1:356
h1=polyval(a1,x1)
plot(x,y,'.',x1,h1,'r')
xlabel('x')
ylabel('y')
title('拟合纳镇陵拟合曲线')
legend('原始数据点','拟合')
结果:
f3 =
- 5.582e-10*X^5 + 6.199e-7*X^4 - 0.0002665*X^3 + 0.05546*X^2 - 5.605*X + 237.2
1.随便粗兆给一个b,c的值,拟合a和d
2,
根据拟合出的a和d和b拟合c
3.
根据a,c,d拟合b
4.
迭迟虚代循环,直到收敛
5.如果不收码凳燃敛,就试试把上边的过程反过来,先拟合b,c之类的
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)