function [ys,W]=SNR_Max(x)
   for k=1:1000;
    s1=sin(0.011*k);
    s2=sin(0.011*k);
    zc=[s1,s2];
    A=[1.5 1;1 1.4];
    x=zc*A;
   [n,T]=size(x);
   x=x-mean(x')'*ones(1,T);
   [F,D]=eig((x*x')/T);
   v=F*(D^(-0.5))*F';
   x=v*x;
   p=80;
   a=ones(1,p)/p;
   x=x';
   S=filter(a,1,x);
   U=cov(S-x,1);
   V=cov(x,1);
   [W d]=eig(V,U);
   ys=(x*W)';
   mplot(zc);
   mplot(x);
   figure(2);
   mplot(ys);
function mplot(s)
 [dim,sample]=size(s);
for k=1:dim
    suplot(dim,1,k);
    plot(s(k,:));
    plot(s(:,k));
    axis('tight');
    set(gca,'XTick',[]);
end
运行时总是显示画图函数有问题  ：Illegal use of reserved keyword "function".请各位大侠指教
谢谢