clear all
% 初始化
x=-1:0.1:1;
n=length(x);
y=0.5+0.4*sin(2*pi*x);

for i=1:n
    noise=sqrt(0.5)*randn(1);
   yn(i)=y(i)+noise;
end

% M階多項試擬合
M=5;

F=zeros(n,M+1);
for k=1:M+1
    F(:,k)=x'.^(k-1);
end

B=F'*yn';


% 代入測試集進行測試
yc=0;

    for m=1:M+1
yc=yc+B(m).*(x'.^(m-1));
    end


%顯示結果
plot(x,y,'b-o',x,yn,'g-o',x,yc,'r--*');