给出一组横纵坐标,绘制曲线逼近这些点,获得曲线方程的过程就是曲线拟合
代码实现:
clear all
clc
close all
x=[1.8 3.2 5 6 6.8 7.5 9 12 13.3 15];
[~,k]=size(x);
y=[-10 -6.9 -4.2 -2 0 2.1 3 5.2 6.4 4.5];
for n=1:6 %从1次方到6次方拟合
ANSS=polyfit(x,y,n); %用ployfit拟合曲线,ANSS是一个系数行向量
for i=1:n+1 %每次求得的方程系数按列存储到answer中
answer(i,n)=ANSS(i);
end
x0=0:0.01:20;
y0=ANSS(1) * x0.^n; %根据系数初始化并构造多项式方程
for num=2:n+1
y0=y0+ANSS(num)*x0.^(n+1-num);
end
subplot(2,3,n)
plot(x,y,'*')
hold on
plot(x0,y0)
end
suptitle('不同次数方程拟合结果,从1到6')
结果如下: