homew2.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

N=10000;
n=1:N-1;
p=0.1;
u=rand(1,N);
u1=u*sqrt(p);
Mlag=50;
ru=xcorr(u1,Mlag,'biased');
A=10^0.5;
s=A*sin(pi/16*n);
x=s+u(1:N-1);
rx=xcorr(x,Mlag,'biased');
figure;
subplot(221);
plot(u1(1:N-1),'.');grid on;
title('white noise');xlabel('n');ylabel('u');
subplot(222);
plot((-Mlag:Mlag),ru,'.');grid on;
title('ru');xlabel('m');ylabel('ru');
subplot(223);
plot(n,x,'.');grid on;
title('signal&noise');xlabel('n');ylabel('x');
subplot(224);
plot((-Mlag:Mlag),rx,'.');grid on;
title('rx');xlabel('m');ylabel('ru');