[an error occurred while processing this directive]
|
%This function calculates
%Least-mean square n+1 degree polinomial approximation
%for data yi, so p(xi)~yi.
function a=lms(xi,yi,n);
xn=zeros(1,2*n+1);
x=zeros(n+1);
y=zeros(1,n+1);
for i=1:2*n, xn(i+1)=sum(xi.^i); end
xn(1)=length(xi);
for i=1:n+1, x(i,:)=xn(i:i+n); end;
for i=0:n, y(i+1)=sum(yi.*xi.^i); end;
a=x\y.';
E-mail: info@telesys.ru