c7_flicker111111111.m

通信系统仿真原理与无线应用书中matlab源程序

f = 0:1000;						% frequency points
fn = 1000;						% Nyquist rate
F = f/fn;						% frequency vector
M = (1-cos(0.001*pi*f)).*1000;  % normalized fequency response
M = [zeros(1,6),M(6:1000)];		% bound from zero frequency
[b1,a1] = yulewalk(3,F,M);		% generate order=3 filter
[b2,a2] = yulewalk(20,F,M);		% generate order=20 filter
[h1,w1] = freqz(b1,a1);			% generate 3-rd order H(f)
[h2,w2] = freqz(b2,a2);			% generate 20-th order H(f)
subplot(2,1,1)
plot(F,M,':',w1/pi,abs(h1))
xlabel('Normalized Frequency')
ylabel('Magnitude Response')
subplot(2,1,2)
plot(F,M,':',w2/pi,abs(h2))
xlabel('Normalized Frequency')
ylabel('Magnitude Response')