ex072800.m

来自「matlab的一些函数,是初学习者的好资料,」· M 代码 · 共 15 行

% 第七章 例 7.28
% 用PM算法进行的希尔伯特变换器设计 
%
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('脉冲响应');
ylabel('h(n)'); axis([0,N,-0.8,0.8]);text(N+1,-0.8,'n')
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('幅度响应')
xlabel('频率（单位：pi）'); ylabel('|H|')
set(gca,'XTickMode','manual','XTick',[0,f,1])
set(gca,'YTickMode','manual','YTick',[0,1]);grid