%%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)