example1_13.m

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

%example1-13
clear
p=0.1;
N=50000;
u=randn(1,N);
a=sqrt(p);
u=u*a;
power_u=var(u);
subplot(2,1,1);
plot(u(1:100));grid on;
ylabel('u(n)');
subplot(2,1,2);
hist(u,50);grid on;
ylabel('histogram of u(n)')