한 창함수를 사용하여 길이가 25인 디지털 힐버트 변환기를 설계하라.

ex7_13.cem

/*   Chapter 7: Example 7.13   Hilbert Transformer design - Hanning window
*/
 
 M = 25; alpha = (M-1)/2;
 n = 0:M-1;
 hd = (2/pi)*((sin((pi/2)*(n-alpha)).^2)./(n-alpha)); hd(alpha+1)=0;
 w_han = (hanning(M))';
 h = hd .* w_han;
 [Hr,w,P,L] = Hres_t3(h);
 /* plots */
 subplot(1,1,1);
 subplot(2,2,1); stem(n,hd); title("Ideal Impulse Response") xtitle("n"); ytitle("hd(n)") subplot(2,2,2); stem(n,w_han);title("Hanning Window") xtitle("n"); ytitle("w(n)") subplot(2,2,3); stem(n,h);title("Actual Impulse Response") xtitle("n"); ytitle("h(n)") w = w'; Hr = Hr';
 w = [-fliplr(w), w(2:501)]; Hr = [-fliplr(Hr), Hr(2:501)];
 subplot(2,2,4);plot(w/pi,Hr); title("Amplitude Response") xtitle("frequency in pi units"); ytitle("Hr")