MATLAB history on Proakis & Man Ch 14

%-- 30/10/09 11:35 --%
N=2^14;n=0:N-1;w=randn(size(n))*3+2;
mean(w)
var(w)
b=[1 -2 1];
x=filter(b,1,w;
x=filter(b,1,w);
mean(x)
var(x)
help xcorr
rxx=xcorr(w,4)
rxx=xcorr(w,4,'biased')
rxx=xcorr((w-mean(w)),4,'biased')
rxx=xcorr(x,4,'biased')
rxx=xcorr(x,8,'biased')
stem(rxx)
M=2^7;n=0:M-1;x=sin(2*pi*.1*n);
clear fax
X=fft(x,M*8);
N=length(X)
stem(fax,abs(X))
stem(fax,20*log10(abs(X)))
plot(fax,20*log10(abs(X)))
grid
ginput(2)
mea=ginput(2)
mea(2)-meaa(1)
mea(2)-mea(1)
M=100;n=0:M-1;x=sin(2*pi*.1*n);
X=fft(x,M*8);
N=length(X)
plot(fax,20*log10(abs(X)))
grid
mea=ginput(2)
mea(2,1)-mea(1,1)
M=90;n=0:M-1;x=sin(2*pi*.1*n);
X=fft(x,M*8);
X=fft(x,M*16);
plot(fax,20*log10(abs(X)))
N=length(X)
plot(fax,20*log10(abs(X)))
grid
mea=ginput(2)
mea(2,1)-mea(1,1)
[H,w]=freqz(b,a,512,Fs);
semilogx(w,20*log10(abs(H)))
Fs=1000;
10*log10(3/Pn)
Q=2/2^6;Pn=Q^2/12
10*log10(3/Pn)