移动最小二乘法MLS程序
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
x0=-1:0.1:5;
n=length(x0);
x=[0,1,2,3,4];
y=[0,1,1,2,2]';
w=zeros(1,5);
d=3;
f=zeros(1,n);
for j=1:n
for i=1:5
w(i)=exp(-(x0(j)-x(i))^2/d^2);
end
A=zeros(2,2);
for i=1:5
A=w(i)*[1;x(i)]*[1,x(i)]+A;
end
B=zeros(2,5);
for i=1:5
B(:,i)=w(i)*[1;x(i)];
end
f(j)=[1,x0(j)]*inv(A)*B*y;
end
plot(x0,f);
hold on
plot(x,y,'*');