ex072800.m

来自「主要是介绍DSP算法在matlab中的实现」· M 代码 · 共 29 行

% Chapter 7 Example 7.28

%           Hilbert transformer design using PM algorithm

%

f = [0.05,0.95]; m = [1 1]; M = 51; N = M-1;

h = remez(N,f,m,'hilbert');

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

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

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

xlabel('n'); ylabel('h(n)'); axis([0,N,-0.8,0.8])

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

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

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

xlabel('frequency in pi units'); ylabel('|H|')

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

set(gca,'YTickMode','manual','YTick',[0,1]);grid