首页 > 代码库 > 《DSP using MATLAB》示例Example 8.8
《DSP using MATLAB》示例Example 8.8
%% ------------------------------------------------------------------------
%% Output Info about this m-file
fprintf(‘\n***********************************************************\n‘);
fprintf(‘ <DSP using MATLAB> Exameple 8.8 \n\n‘);
time_stamp = datestr(now, 31);
[wkd1, wkd2] = weekday(today, ‘long‘);
fprintf(‘ Now is %20s, and it is %9s \n‘, time_stamp, wkd2);
%% ------------------------------------------------------------------------
Wp = 0.2*pi; Ws = 0.3*pi; Rp = 1; As =16;
Ripple = 10 ^ (-Rp/20); Attn = 10 ^ (-As/20);
% Analog filter design:
[b, a] = afd_elip(Wp, Ws, Rp, As);
% Calculation of second-order sections:
[C, B, A] = sdir2cas(b, a)
% Calculation of Frequency Response:
[db, mag, pha, ww] = freqs_m(b, a, 0.5*pi);
% Calculation of Impulse Response:
[ha, x, t] = impulse(b, a);
%% -----------------------------------------------------------------
%% Plot
%% -----------------------------------------------------------------
figure(‘NumberTitle‘, ‘off‘, ‘Name‘, ‘Exameple 8.8‘)
set(gcf,‘Color‘,‘white‘);
M = 0.5; % Omega max
subplot(2,2,1); plot(ww/pi, mag); axis([0, M, 0, 1.2]); grid on;
xlabel(‘ Analog frequency in \pi units‘); ylabel(‘|H|‘); title(‘Magnitude Response‘);
set(gca, ‘XTickMode‘, ‘manual‘, ‘XTick‘, [0, 0.2, 0.3, M]);
set(gca, ‘YTickMode‘, ‘manual‘, ‘YTick‘, [0, 0.1585, 0.8913, 1]);
subplot(2,2,2); plot(ww/pi, db); axis([0, M, -30, 10]); grid on;
xlabel(‘Analog frequency in \pi units‘); ylabel(‘Decibels‘); title(‘Magnitude in dB ‘);
set(gca, ‘XTickMode‘, ‘manual‘, ‘XTick‘, [0, 0.2, 0.3, M]);
set(gca, ‘YTickMode‘, ‘manual‘, ‘YTick‘, [-30, -16, -1, 0]);
subplot(2,2,3); plot(ww/pi, pha/pi); axis([0, M, -1.1, 1.1]); grid on;
xlabel(‘Analog frequency in \pi nuits‘); ylabel(‘radians in \pi units‘); title(‘Phase Response‘);
set(gca, ‘XTickMode‘, ‘manual‘, ‘XTick‘, [0, 0.2, 0.3, M]);
set(gca, ‘YTickMode‘, ‘manual‘, ‘YTick‘, [-1:0.5:1]);
subplot(2,2,4); plot(t, ha); axis([0, 50, -0.1, 0.3]); grid on;
xlabel(‘time in seconds‘); ylabel(‘ha(t)‘); title(‘Impulse Response‘);
运行结果:
《DSP using MATLAB》示例Example 8.8
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。