机器学习 最小二乘学习法 matlab

rand('state',0);randn('state',0);n=50;N=1000;x=linspace(-3,3,n)';X=linspace(-3,3,N)';pix=pi*x;y=sin(pix)./(pix)+0.1*x+0.05*randn(n,1);//要拟合的函数p(:,1)=ones(n,1);P(:,1)=ones(N,1);for j=1:15 P(:,2*j)=sin(j/2*X);P(:,2*j+1)=cos(j/2*X); p(:,2*j)=sin(j/2*x);p(:,2*j+1)=cos(j/2*x);endt=p\y;F=P*t;figure(1);clf;hold on;axis([-5 5 -2 2]);plot(X,F,'g+');plot(x,y,'bo');z