有六百组数据,需要进行最小二乘法拟合,1-7、2-8、3-9这样的七组进行一次拟合,请教!!
clear all;
close all;
clc;
N=10; %设置拟合阶数
x=1:0.5:10;
y=cos(x); %生成待拟合点
p=polyfit(x,y,N); %使用matlab函数拟合数据
xx=min(x):0.01:max(x);
yy=polyval(p,xx);
plot(xx,yy); %画出拟合结果
hold on;
plot(x,y,'r.')
%下面是使用公式来做最小二乘多项式拟合
F=zeros(N+1,length(x));
F(1,:)=1;
for i=2:N+1
for j=1:length(x)
F(i,j) = x(j)^(i-1);
end
end
F=F*F';
[m ~]=size(F);
Y=zeros(m,1);
Y(1) = sum(y);
for i=2:m
for j=1:length(y)
Y(i) = Y(i)+y(j)*x(j)^(i-1);
end
end
Re = F\Y;
Re=Re(end:-1:1)'; %数组反序
figure;
plot(x,y,'r.')
hold on;
yyy=polyval(Re,xx);
plot(xx,yyy,'g')
p
Re