code from digital signal processing using matlab, 3rd edition

Question

0 个投票

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 = '(hann(M))';
h = hd .* w_han.';
[Hr,w,c,L] = Hr_Type3(h);
subplot(2,2,1); stem(n,hd); title('Ideal Impulse Response')
axis([-1 M -1.2 1.2]); xlabel('n'); ylabel('hd(n)')
subplot(2,2,2);
stem(n,w_han);                
title('Hann Window');
axis([-1 M 0 1.2]); xlabel ' n '); ylabel('w(n)')
subplot(2,2,3); stem(n,h);title('Actual Impulse Response')
axis([-1 M -1.2 1.2]); xlabel('n'); ylabel('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');grid;
xlabel('frequency in pi units'); ylabel('Hr'); axis([-1 1 -1.1 1.1]);

2 个评论
显示无隐藏无

Image Analyst 2020-12-27

编辑：Image Analyst 2020-12-27

What about it? Do you have a question about it? Doesn't run for me because you forgot to include Hr_Type3(). Please supply the function Hr_Type3() if you want us to try the code.

muhammad nur adri nawi 2020-12-27

function [Hr,w,c,L] = Hr_Type3(h)

M=length(h);

L=(M-1)/2;

c=[2*h(L+1:-1:1)];

n=[0:1:L];

w=[0:1:500]'*pi/500;

Hr=sin(w*n)*c';

end

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question

Answer 1

KALYAN ACHARJYA 2020-12-27

编辑：KALYAN ACHARJYA 2020-12-27

在 MATLAB Online 中打开

0 个投票

Note: Here This

w_han = '(hann(M))';

To this > May be Hanning Window here

w_han =hann(M)';

Complete Code:

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 =hann(M)';
h = hd .* w_han.';
[Hr,w,c,L] = Hr_Type3(h);
subplot(2,2,1); stem(n,hd); title('Ideal Impulse Response')
axis([-1 M -1.2 1.2]); xlabel('n'); ylabel('hd(n)')
subplot(2,2,2);
stem(n,w_han);                
title('Hann Window');
axis([-1 M 0 1.2]); xlabel ' n '); ylabel('w(n)')
subplot(2,2,3); stem(n,h);title('Actual Impulse Response')
axis([-1 M -1.2 1.2]); xlabel('n'); ylabel('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');grid;
xlabel('frequency in pi units'); ylabel('Hr'); axis([-1 1 -1.1 1.1]);