ex072700.m

Digital Signal Processing Using MATLAB V.4随书代码

% Chapter 7 Example 7.27

%           Differentiator design using PM algorithm

%

f = [0 0.2 0.4 0.6 0.8 1];         % in w/pi unis

%m = [1 1 2 2 3 3];                 % in sam/cycle (old version)

m = [0,0.1,0.4,0.6,1.2,1.5];       % new Student Edition

h = remez(25,f,m,'differentiator');

[db,mag,pha,grd,w] = freqz_m(h,[1]);

figure(1); subplot(1,1,1)

subplot(2,1,1); stem([0:25],h); title('Impulse Response');

xlabel('n'); ylabel('h(n)'); axis([0,25,-0.6,0.6])

set(gca,'XTickMode','manual','XTick',[0,25])

set(gca,'YTickMode','manual','YTick',[-0.6:0.2:0.6]);

subplot(2,1,2); plot(w/(2*pi),mag); title('Magnitude Response')

xlabel('Normalized frequency f'); ylabel('|H|')

set(gca,'XTickMode','manual','XTick',f/2)

set(gca,'YTickMode','manual','YTick',m)

grid