把我写的权函数给你参考一下吧,写的是多元的情况,用的是高斯核函数。附件不支持m文件的格式,我就直接复制粘贴了。
function W=Weight(N,x,xstar,h)
W=zeros(N,N);
for i=1:N
W(i,i)=exp(-(x(1,i)-xstar(1))^2/(2*h(1)^2)-(x(2,i)-xstar(2))^2/(2*h(2)^2))/(2*pi*h(1)*h(2));
end
N表示样本点的个数,x表示样本数据点,xstar表示在哪一个点展开,h表示你给出的带宽。
根据你的情况,我推荐你看看范建清写的局部多项式回归一书,当然其他非参数的方法也有,比较偏向样条方面的。Local Polynomial Modelling and Its Applications
Nonparametric and Semiparametric Models-2004
文件上传不了,说是论坛已存在相关附件,那就搜索之后下吧,没有论坛币我也可以传给你。