1 clear all;
 2 close all;
 3 clc;
 4 
 5 x=(1:0.1:10)';
 6 y=x.^2+x+3 +rand(length(x),1)*6;
 7 plot(x,y,'.')
 8 
 9 sigma=0.1;              %设置局部窗口，越大越使用全局数据，越小越使用局部数据
10 W=zeros(length(x));
11 C=[];
12 for i=1:length(x)    
13     for j=1:length(x) 
14         W(j,j)=exp(-((x(i)-x(j))^2)/(2*sigma^2));   %权重矩阵
15     end
16 
17     XX=[x ones(length(x),1)];
18     YY=y;
19     C=[C inv(XX'*W*XX)*XX'*W*YY];           %加权最小二乘，计算求得局部线性函数的系数
20     
21 end
22 
23 re=diag(XX*C);
24 hold on;
25 plot(x,re);