《DSP using MATLAB》Problem 7.16

最新推荐文章于 2024-11-06 10:36:48 发布
最新推荐文章于 2024-11-06 10:36:48 发布
阅读量181 收藏
点赞数
文章标签： matlab
原文链接：http://www.cnblogs.com/ky027wh-sx/p/10555840.html
版权

        使用一种固定窗函数法设计带通滤波器。

代码：

%% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%%            Output Info about this m-file
fprintf('\n***********************************************************\n');
fprintf('        <DSP using MATLAB> Problem 7.16 \n\n');

banner();
%% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

% bandpass
ws1 = 0.3*pi; wp1 = 0.4*pi; wp2 = 0.5*pi; ws2 = 0.6*pi; As = 40; Rp = 0.5;
tr_width = min((wp1-ws1), (ws2-wp2));
[delta1_ori, delta2_ori] = db2delta(Rp, As)

%% ---------------------------------------------------
%%         1  Rectangular Window
%% ---------------------------------------------------
M = ceil(1.8*pi/tr_width) + 1;                 % Rectangular Window
fprintf('\n\n#1.Rectangular Window, Filter Length M=%d.\n', M);

n = [0:1:M-1]; wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2;

hd = ideal_lp(wc2, M) - ideal_lp(wc1, M); 
w_rect = (boxcar(M))';  h = hd .* w_rect;
[db, mag, pha, grd, w] = freqz_m(h, [1]);  delta_w = 2*pi/1000;
[Hr,ww,P,L] = ampl_res(h);

Rp = -(min(db(wp1/delta_w+1 :1: floor(wp2/delta_w)+1)));     % Actual Passband Ripple
fprintf('\nActual Passband Ripple is %.4f dB.\n', Rp);

As = -round(max(db(ws2/delta_w+1 : 1 : 501)));        % Min Stopband attenuation
fprintf('\nMin Stopband attenuation is %.4f dB.\n', As);

[delta1_rect, delta2_rect] = db2delta(Rp, As)


%%  ----------------------------
%%             Plot
%% -----------------------------

figure('NumberTitle', 'off', 'Name', 'Problem 7.16.1 ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); stem(n, hd); axis([0 M-1 -0.3 0.3]); grid on;
xlabel('n'); ylabel('hd(n)'); title('Ideal Impulse Response');
subplot(2,2,2); stem(n, w_rect); axis([0 M-1 0 1.1]); grid on;
xlabel('n'); ylabel('w(n)'); title('Rectangular Window, M=19');
subplot(2,2,3); stem(n, h); axis([0 M-1 -0.3 0.3]); grid on;
xlabel('n'); ylabel('h(n)'); title('Actual Impulse Response');

subplot(2,2,4); plot(w/pi, db); axis([0 1 -100 10]); grid on;
set(gca,'YTickMode','manual','YTick',[-90,-26,0]);
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'26';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1]);
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.1 h(n) ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); plot(w/pi, db); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');
set(gca,'YTickMode','manual','YTick',[-90,-26,0])
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'26';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);

subplot(2,2,3); plot(w/pi, mag); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Absolute'); title('Magnitude Response in absolute');
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);
set(gca,'YTickMode','manual','YTick',[0.0,0.5,1.0])

subplot(2,2,2); plot(w/pi, pha); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Phase Response in Radians');
subplot(2,2,4); plot(w/pi, grd*pi/180);  grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Group Delay');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.1 h(n)')
set(gcf,'Color','white'); 

plot(ww/pi, Hr); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Hr'); title('Amplitude Response');
set(gca,'YTickMode','manual','YTick',[-delta2_rect,0,delta2_rect,1 - delta1_rect,1, 1 + delta1_rect])


%% ---------------------------------------------------
%%             2  Bartlett Window
%% ---------------------------------------------------
M = ceil(6.1*pi/tr_width) + 1;                 % Bartlett Window
fprintf('\n\n#2.Bartlett Window, Filter Length M=%d.\n', M);

n = [0:1:M-1]; wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2;

%wc = (ws + wp)/2,                    % ideal LPF cutoff frequency

hd = ideal_lp(wc2, M) - ideal_lp(wc1, M); 
w_bart = (bartlett(M))';  h = hd .* w_bart;
[db, mag, pha, grd, w] = freqz_m(h, [1]);  delta_w = 2*pi/1000;
[Hr,ww,P,L] = ampl_res(h);

Rp = -(min(db(wp1/delta_w+1 :1: floor(wp2/delta_w)+1)));     % Actual Passband Ripple
fprintf('\nActual Passband Ripple is %.4f dB.\n', Rp);

As = -round(max(db(ws2/delta_w+1 : 1 : 501)));               % Min Stopband attenuation
fprintf('\nMin Stopband attenuation is %.4f dB.\n', As);

[delta1_bart, delta2_bart] = db2delta(Rp, As)

%% --------------------------
%%            Plot
%% --------------------------

figure('NumberTitle', 'off', 'Name', 'Problem 7.16.2 ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); stem(n, hd); axis([0 M-1 -0.2 0.2]); grid on;
xlabel('n'); ylabel('hd(n)'); title('Ideal Impulse Response');

subplot(2,2,2); stem(n, w_bart); axis([0 M-1 0 1.1]); grid on;
xlabel('n'); ylabel('w(n)'); title('Bartlett Window, M=62');

subplot(2,2,3); stem(n, h); axis([0 M-1 -0.2 0.2]); grid on;
xlabel('n'); ylabel('h(n)'); title('Actual Impulse Response');

subplot(2,2,4); plot(w/pi, db); axis([0 1 -100 10]); grid on;
set(gca,'YTickMode','manual','YTick',[-90,-27,0]);
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'27';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1]);
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.2 h(n) ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); plot(w/pi, db); grid on; axis([0 2 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');
set(gca,'YTickMode','manual','YTick',[-90,-27,0])
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'27';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);

subplot(2,2,3); plot(w/pi, mag); grid on; %axis([0 2 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Absolute'); title('Magnitude Response in absolute');
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);
set(gca,'YTickMode','manual','YTick',[0.0,0.5,1.0])


subplot(2,2,2); plot(w/pi, pha); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Phase Response in Radians');
subplot(2,2,4); plot(w/pi, grd*pi/180);  grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Group Delay');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.2 h(n)')
set(gcf,'Color','white'); 

plot(ww/pi, Hr); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Hr'); title('Amplitude Response');
set(gca,'YTickMode','manual','YTick',[-delta2_bart,0,delta2_bart,1-delta1_bart,1, 1+delta1_bart])


%% ---------------------------------------------------
%%              3  Hann Window
%% ---------------------------------------------------
M = ceil(6.2*pi/tr_width) + 1;                 % Hann Window
fprintf('\n\n#3.Hann Window, Filter Length M=%4d.\n', M);

n = [0:1:M-1]; wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2;

%wc = (ws + wp)/2,                    % ideal LPF cutoff frequency

hd = ideal_lp(wc2, M) - ideal_lp(wc1, M); 
w_hann = (hann(M))';  h = hd .* w_hann;
[db, mag, pha, grd, w] = freqz_m(h, [1]);  delta_w = 2*pi/1000;
[Hr,ww,P,L] = ampl_res(h);

Rp = -(min(db(wp1/delta_w+1 :1: floor(wp2/delta_w)+1)));     % Actual Passband Ripple
fprintf('\nActual Passband Ripple is %.4f dB.\n', Rp);

As = -round(max(db(ws2/delta_w+1 : 1 : 501)));               % Min Stopband attenuation
fprintf('\nMin Stopband attenuation is %.4f dB.\n', As);

[delta1_hann, delta2_hann] = db2delta(Rp, As)

%% --------------------------
%%            Plot
%% --------------------------

figure('NumberTitle', 'off', 'Name', 'Problem 7.16.3 ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); stem(n, hd); axis([0 M-1 -0.2 0.3]); grid on;
xlabel('n'); ylabel('hd(n)'); title('Ideal Impulse Response');

subplot(2,2,2); stem(n, w_hann); axis([0 M-1 0 1.1]); grid on;
xlabel('n'); ylabel('w(n)'); title('Hann Window, M=63');

subplot(2,2,3); stem(n, h); axis([0 M-1 -0.2 0.3]); grid on;
xlabel('n'); ylabel('h(n)'); title('Actual Impulse Response');

subplot(2,2,4); plot(w/pi, db); axis([0 1 -100 10]); grid on;
set(gca,'YTickMode','manual','YTick',[-90,-43,0]);
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'43';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1]);
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.3 h(n) ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); plot(w/pi, db); grid on; axis([0 2 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');
set(gca,'YTickMode','manual','YTick',[-90,-43,0])
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'43';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);

subplot(2,2,3); plot(w/pi, mag); grid on; %axis([0 2 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Absolute'); title('Magnitude Response in absolute');
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);
set(gca,'YTickMode','manual','YTick',[0.0,0.5,1.0])

subplot(2,2,2); plot(w/pi, pha); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Phase Response in Radians');
subplot(2,2,4); plot(w/pi, grd*pi/180);  grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Group Delay');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.3 h(n)')
set(gcf,'Color','white'); 

plot(ww/pi, Hr); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Hr'); title('Amplitude Response');
set(gca,'YTickMode','manual','YTick',[-delta2_hann,0,delta2_hann,1 - delta1_hann,1, 1 + delta1_hann])



%% ---------------------------------------------------
%%             4  Hamming Window
%% ---------------------------------------------------
M = ceil(6.6*pi/tr_width) + 1;                 % Hamming Window
fprintf('\n\n#4.Hamming Window, Filter Length M=%4d.\n', M);

n = [0:1:M-1]; wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2;

%wc = (ws + wp)/2,                    % ideal LPF cutoff frequency

hd = ideal_lp(wc2, M) - ideal_lp(wc1, M); 
w_hamm = (hamming(M))';  h = hd .* w_hamm;
[db, mag, pha, grd, w] = freqz_m(h, [1]);  delta_w = 2*pi/1000;
[Hr,ww,P,L] = ampl_res(h);

Rp = -(min(db(wp1/delta_w+1 :1: floor(wp2/delta_w)+1)));     % Actual Passband Ripple
fprintf('\nActual Passband Ripple is %.4f dB.\n', Rp);

As = -round(max(db(ws2/delta_w+1 : 1 : 501)));               % Min Stopband attenuation
fprintf('\nMin Stopband attenuation is %.4f dB.\n', As);

[delta1_hamm, delta2_hamm] = db2delta(Rp, As)

%% --------------------------
%%            Plot
%% --------------------------

figure('NumberTitle', 'off', 'Name', 'Problem 7.16.4 ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); stem(n, hd); axis([0 M-1 -0.2 0.3]); grid on;
xlabel('n'); ylabel('hd(n)'); title('Ideal Impulse Response');

subplot(2,2,2); stem(n, w_hamm); axis([0 M-1 0 1.1]); grid on;
xlabel('n'); ylabel('w(n)'); title('Hamming Window, M=67');

subplot(2,2,3); stem(n, h); axis([0 M-1 -0.2 0.3]); grid on;
xlabel('n'); ylabel('h(n)'); title('Actual Impulse Response');

subplot(2,2,4); plot(w/pi, db); axis([0 1 -100 10]); grid on;
set(gca,'YTickMode','manual','YTick',[-90,-51,0]);
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'51';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1]);
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.4 h(n) ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); plot(w/pi, db); grid on; axis([0 2 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');
set(gca,'YTickMode','manual','YTick',[-90,-51,0])
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'51';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);

subplot(2,2,3); plot(w/pi, mag); grid on; %axis([0 2 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Absolute'); title('Magnitude Response in absolute');
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);
set(gca,'YTickMode','manual','YTick',[0.0,0.5,1.0])


subplot(2,2,2); plot(w/pi, pha); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Phase Response in Radians');
subplot(2,2,4); plot(w/pi, grd*pi/180);  grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Group Delay');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.4 h(n)')
set(gcf,'Color','white'); 

plot(ww/pi, Hr); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Hr'); title('Amplitude Response');
set(gca,'YTickMode','manual','YTick',[-delta2_hamm,0,delta2_hamm,1 - delta1_hamm,1, 1 + delta1_hamm])



%% ---------------------------------------------------
%%             5  Blackman Window
%% ---------------------------------------------------
M = ceil(11*pi/tr_width) + 1;                 % Blackman Window
fprintf('\n\n#5.Blackman Window, Filter Length M=%d.\n', M);

n = [0:1:M-1]; wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2;

%wc = (ws + wp)/2,                    % ideal LPF cutoff frequency

hd = ideal_lp(wc2, M) - ideal_lp(wc1, M); 
w_bla = (blackman(M))';  h = hd .* w_bla;
[db, mag, pha, grd, w] = freqz_m(h, [1]);  delta_w = 2*pi/1000;
[Hr,ww,P,L] = ampl_res(h);

Rp = -(min(db(wp1/delta_w+1 :1: floor(wp2/delta_w)+1)));     % Actual Passband Ripple
fprintf('\nActual Passband Ripple is %.4f dB.\n', Rp);

As = -round(max(db(ws2/delta_w+1 : 1 : 501)));               % Min Stopband attenuation
fprintf('\nMin Stopband attenuation is %.4f dB.\n', As);

[delta1_bla, delta2_bla] = db2delta(Rp, As)

%% --------------------------
%%            Plot
%% --------------------------

figure('NumberTitle', 'off', 'Name', 'Problem 7.16.5 ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); stem(n, hd); axis([0 M-1 -0.2 0.3]); grid on;
xlabel('n'); ylabel('hd(n)'); title('Ideal Impulse Response');

subplot(2,2,2); stem(n, w_bla); axis([0 M-1 0 1.1]); grid on;
xlabel('n'); ylabel('w(n)'); title('Blackman Window, M=111');

subplot(2,2,3); stem(n, h); axis([0 M-1 -0.2 0.3]); grid on;
xlabel('n'); ylabel('h(n)'); title('Actual Impulse Response');

subplot(2,2,4); plot(w/pi, db); axis([0 1 -120 10]); grid on;
set(gca,'YTickMode','manual','YTick',[-90,-73,0]);
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'73';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1]);
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.5 h(n) ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); plot(w/pi, db); grid on; axis([0 2 -120 10]); 
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');
set(gca,'YTickMode','manual','YTick',[-90,-73,0])
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'73';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);

subplot(2,2,3); plot(w/pi, mag); grid on; %axis([0 2 -120 10]); 
xlabel('frequency in \pi units'); ylabel('Absolute'); title('Magnitude Response in absolute');
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);
set(gca,'YTickMode','manual','YTick',[0.0,0.5,1.0])


subplot(2,2,2); plot(w/pi, pha); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Phase Response in Radians');
subplot(2,2,4); plot(w/pi, grd*pi/180);  grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Group Delay');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.5 h(n)')
set(gcf,'Color','white'); 

plot(ww/pi, Hr); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Hr'); title('Amplitude Response');
set(gca,'YTickMode','manual','YTick',[-delta2_bla,0,delta2_bla,1-delta1_bla,1, 1+delta1_bla])


%% ---------------------------------------------------
%%             6  Kaiser Window
%% ---------------------------------------------------
M = ceil((As-7.95)/(2.285*tr_width)) + 1;                 % Kaiser Window
if As > 21 || As < 50 
	beta = 0.5842*(As-21)^0.4 + 0.07886*(As-21);
else
	beta = 0.1102*(As-8.7);
end
fprintf('\n\n#6.Kaiser Window, Filter Length M=%d, beta=%.4f\n', M,beta);

n = [0:1:M-1]; wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2;

%wc = (ws + wp)/2,                    % ideal LPF cutoff frequency

hd = ideal_lp(wc2, M) - ideal_lp(wc1, M); 
w_kai = (kaiser(M,beta))';  h = hd .* w_kai;
[db, mag, pha, grd, w] = freqz_m(h, [1]);  delta_w = 2*pi/1000;
[Hr,ww,P,L] = ampl_res(h);

Rp = -(min(db(wp1/delta_w+1 :1: floor(wp2/delta_w)+1)));     % Actual Passband Ripple
fprintf('\nActual Passband Ripple is %.4f dB.\n', Rp);

As = -round(max(db(ws2/delta_w+1 : 1 : 501)));               % Min Stopband attenuation
fprintf('\nMin Stopband attenuation is %.4f dB.\n', As);

[delta1_kai, delta2_kai] = db2delta(Rp, As)

%% --------------------------
%%            Plot
%% --------------------------

figure('NumberTitle', 'off', 'Name', 'Problem 7.16.6 ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); stem(n, hd); axis([0 M-1 -0.2 0.2]); grid on;
xlabel('n'); ylabel('hd(n)'); title('Ideal Impulse Response');

subplot(2,2,2); stem(n, w_kai); axis([0 M-1 0 1.1]); grid on;
xlabel('n'); ylabel('w(n)'); title('Kaiser Window, M=92');

subplot(2,2,3); stem(n, h); axis([0 M-1 -0.2 0.2]); grid on;
xlabel('n'); ylabel('h(n)'); title('Actual Impulse Response');

subplot(2,2,4); plot(w/pi, db); axis([0 1 -120 10]); grid on;
set(gca,'YTickMode','manual','YTick',[-90,-72,0]);
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'72';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1]);
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.6 h(n) ideal_lp Method')
set(gcf,'Color','white'); 

subplot(2,2,1); plot(w/pi, db); grid on; axis([0 2 -120 10]); 
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');
set(gca,'YTickMode','manual','YTick',[-90,-72,0])
set(gca,'YTickLabelMode','manual','YTickLabel',['90';'72';' 0']);
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);

subplot(2,2,3); plot(w/pi, mag); grid on; %axis([0 2 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Absolute'); title('Magnitude Response in absolute');
set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,0.5,0.6,1,1.4,1.5,1.6,1.7,2]);
set(gca,'YTickMode','manual','YTick',[0.0,0.5,1.0])


subplot(2,2,2); plot(w/pi, pha); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Phase Response in Radians');
subplot(2,2,4); plot(w/pi, grd*pi/180);  grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Rad'); title('Group Delay');


figure('NumberTitle', 'off', 'Name', 'Problem 7.16.6 h(n)')
set(gcf,'Color','white'); 

plot(ww/pi, Hr); grid on; %axis([0 1 -100 10]); 
xlabel('frequency in \pi units'); ylabel('Hr'); title('Amplitude Response');
set(gca,'YTickMode','manual','YTick',[-delta2_kai,0,delta2_kai,1-delta1_kai,1, 1+delta1_kai])

　　运行结果：

        设计指标，As=40dB，Rp=0.5dB；换算成绝对指标，为δ1=0.0288，δ2=0.0103。

        由书中表7.1可知，矩形窗（Rectangular）和三角窗（Bartlett）不满足设计要求，这里我们也进行贴图验证。

        上图可知，rectangualr窗和bartlett窗不符合设计要求。

        下面将以上几种窗函数设计结果，总结如下

序号名  称长度MAsRp
1Rectangular19261.918
2Bartlett　　62270.1059
3Hann　　63430.1242
4Hamming67510.0488
5Blackman　　111730.0027
6Kaiser92720.0049

 

 

 

 

 

        由上表得知，满足设计要求的是用长M=63的Hann窗截断得到的滤波器。

 

转载于:https://www.cnblogs.com/ky027wh-sx/p/10555840.html

吃素的小动物
关注
elasticsearch-analysis-ik-7.16.2.zip
10-02
《Elasticsearch分词器：elasticsearch-analysis-ik-7.16.2深度解析》 在信息爆炸的时代，搜索引擎的效能成为了数据检索的关键。Elasticsearch作为一款强大的开源搜索引擎，其灵活性和可扩展性备受青睐。而在中文...
最新版windows kibana-7.16.2-windows-x86_64.zip
12-20
最新版的Windows Kibana 7.16.2 是针对64位（x86_64）架构设计的，这意味着它可以在运行Windows操作系统的64位计算机上高效运行。7.16.2版本可能包含了一系列性能改进、新功能和错误修复，旨在提升用户体验和数据...
《DSP using MATLAB》Problem 7.6
weixin_30654583的博客
02-26 130
代码： 子函数ampl_res function [Hr,w,P,L] = ampl_res(h); % % function [Hr,w,P,L] = Ampl_res(h) % Computes Amplitude response Hr(w) and its polynomial P of order L, % given a linear-phase...
《DSP using MATLAB》Problem 7.13
weixin_30701521的博客
03-07 95
代码： %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output Info about this m-file fprintf('\n**************************************************...
《DSP using MATLAB》Problem 6.18
weixin_30892037的博客
11-18 78
代码： %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output Info about this m-file fprintf('\n******************************************...
《DSP using MATLAB》Problem 7.25
weixin_30566063的博客
04-10 78
代码： %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output Info about this m-file fprintf('\n************************************************...
《DSP using MATLAB》Problem 7.36
weixin_30568591的博客
05-21 254
代码： %% ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ %% Output Info about this m-file fprintf('\n************************************************...
《DSP using MATLAB》示例Example7.16
weixin_30687811的博客
03-30 75
代码： M = 60; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l; Hrs = [ones(1, 7), 0.5925, 0.1099, zeros(1, 43), 0.1099, 0.5925, ones(1, 6)]; % Ideal Amp Res sampled Hdr = [1, 1, 0, 0]; wdl = [0, ...
7.16
weixin_42645474的博客
07-16 202
学习日志                     姓名：  尤嘉慧           日期： 2017.07.16   今日学习任务 单片机概念原理介绍及仿真软件的使用   今日任务完成情况 （详细说明本日任务是否按计划完成，开发的代码量。）   本日任务按计划完成 1、仿真环境操作流程。 ...
JxBrowser 7.16 Crack by Xacker
控件与插件之大全
07-08 4250
我们很高兴向您介绍 JxBrowser 7.16！在此版本中，我们使用新功能扩展了 API，并提高了库的稳定性。 什么是新的 现在您可以像在 Google Chrome 中一样查看加载的网页或框架的来源： browser.mainFrame().ifPresent(Frame::viewSource);????578867473 上面的代码告诉 Chromium 使用主框架的 HTML 源创建并打开一个弹出窗口。它将是一个具有以下外观的常规弹出窗口： 默认背景颜色 当 Chromium 不知道
最新matlab版本号,如何查看matlab版本信息？matlab版本信息怎么看
weixin_42516600的博客
03-16 954
查看MATLAB的版本信息，帮助我们更好的完成工作。学习使用ver函数，获取matlab的版本信息以及图像处理工具箱的版本信息。方法/步骤双击matlab快捷方式，运行matlab，程序启动比较慢，请等待！之后将电脑的输入法设置成英文状态，matlab的语言基本要求包括大小写，英文输入等，中文输入法是运行不了程序的，函数字母的大小写也有影响。ver函数的基本格式：ver()，单击回车表示运行。输入...
MATLAB 图像分割
qq_69604995的博客
12-08 4015
7.1点检测代码： 实验代码：7.2检测指定方向的线代码：clcclearf=imread('Fig1004(a).tif')subplot(321);imshow(f)w=[2 -1 -1;-1 2 -1;-1 -1 2];g=imfilter(tofloat(f),w);subplot(322);imshow(g,[])gtop=g(1:120,1:120);%Top,left section.gtop=pixeldup(gtop,4);%Enlarge by pixel duplication.
最新版linux logstash-7.16.3-linux-x86_64.tar.gz
01-14
最新版的 Logstash 7.16.3 为用户提供了一系列改进和新特性，以提升性能和稳定性。 首先，Logstash 包含三个主要组件：输入（Inputs）、过滤（Filters）和输出（Outputs）。输入插件用于从不同源获取数据，如文件、...
最新版linux logstash-7.16.1-linux-x86_64.tar.gz
12-16
这个最新的Linux版本，logstash-7.16.1-linux-x86_64.tar.gz，是专为64位Linux系统设计的，旨在提供高效的数据处理能力。 一、Logstash简介 Logstash主要由三个组件构成：输入（Inputs）、过滤（Filters）和输出...
最新版linux kibana-7.16.3-linux-x86_64.tar.gz
01-14
《Linux Kibana 7.16.3：探索数据可视化的新境界》 Kibana是Elastic Stack（以前称为ELK Stack）的关键组件之一，它是一个强大的数据分析和可视化工具，专为与Elasticsearch数据库配合使用而设计。Elasticsearch是...
MATLAB 将fig格式另存为可编辑的eps格式，但乱码问题解决
最新发布
qq_47885795的博客
11-06 230
fig格式图像正常，但通过手动导出后的eps格式图像导入到AI中会乱码，如下图所示。
MATLAB-数学建模-无约束规划求解方法（非线性规划）
一个大一博客！
11-02 493
MATLAB-数学建模-无约束规划求解方法（非线性规划）
（57）MATLAB使用迫零均衡器和MMSE均衡器的BPSK调制系统仿真
weixin_xxxxx的博客
11-03 532
本文给出仿真模型与MATLAB代码，分别使用具有ISI的三个不同传输特性的信道，仿真测试了使用迫零均衡器和MMSE均衡器的基带BPSK传输系统的性能，绘制出Eb/N0-SER性能曲线。
数据转换 | Matlab基于SP符号递归图（Symbolic recurrence plots）一维数据转二维图像方法
关注并私信文章链接，获取对应文章源码和数据，机器学习之心的博客。
11-02 290
数据转换 | Matlab基于SP符号递归图（Symbolic recurrence plots）一维数据转二维图像方法
Elasticsearch 7.16.1 for Linux系统发布
资源摘要信息:"最新版linux elasticsearch-7.16.1-linux-x86_64.tar.gz" Linux Elasticsearch 7.16.1版本是基于Linux操作系统的 Elasticsearch搜索服务器的最新稳定版本。Elasticsearch 是一个基于 Lucene 构建的...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值