%%To creat the white noise E(k) and n(k)
subplot(2,1,1)
y=randn(1,3000);
y=y/std(y);
y=y-mean(y);
a=0;
b=sqrt(1);
y=a+b*y;
plot(1:1:3000,y);
title('White noise \xi(k)');

subplot(2,1,2)
z=randn(1,3000);
z=z/std(z);
z=z-mean(z);
c=0;
d=sqrt(1);
z=c+d*z;
plot(1:1:3000,z);
title('White noise \eta(k)');