ex072800.m

matlab写的数字信号处理源代码,里面包括最基本的fft计算

% 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