代碼：

M = 51; alpha = (M-1)/2; Dw = 2*pi/M; l = 0:M-1; wl = Dw*l; T1 = j*0.39; k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1;Hrs = [0, -j*ones(1, (M-3)/2), -T1, T1, j*ones(1, (M-3)/2)]; angH = [-alpha*Dw*k1, alpha*Dw*(M-k2)]; H = Hrs.*exp(j*angH); h = real(ifft(H, M));[db, mag, pha, grd, w] = freqz_m(h, 1); [Hr, ww, c, L] = Hr_Type3(h);%Plot figure('NumberTitle', 'off', 'Name', 'Exameple 7.20') set(gcf,'Color','white'); subplot(2,2,1); plot(wl(1:26)/pi, j*(Hrs(1:26)), 'o', ww/pi, -Hr); axis([0, 1, -1.0, 1.5]); grid on; xlabel('frequency in \pi nuits'); ylabel('Amplitude Response'); title('Hilbert Transformer, frequency sampleing design: M=51');subplot(2,2,2); stem(l, h); axis([-1, M, -1.1, 1.1]); grid on; xlabel('n'); ylabel('h(n)'); title('Impulse Response');subplot(2,2,3); plot(w/pi, pha/pi); axis([0, 1, -1.5, 2]); grid on; xlabel('frequency in \pi units'); ylabel('pha(w) in \pi units'); title('Phase Response');subplot(2,2,4); plot(w/pi, db); axis([0, 1, -60, 10]); grid on; xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response');

　　運行結果：

? ? ??過渡帶中有采樣點（左上圖）的效果。

?

轉載于:https://www.cnblogs.com/ky027wh-sx/p/6649462.html

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎勵來咯，堅持創(chuàng)作打卡瓜分現(xiàn)金大獎

總結

以上是生活随笔為你收集整理的《DSP using MATLAB》示例Example7.20的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。