afdtran.m

AFD - Advanced Filter Design using MATLABMiroslav D. Lutovac, Dejan V. Tosicversion 1.00 released 15

function twoaxis=afdtran(num,den,speca,filnumb,moreaxis)

% afdtran.m  AFD Plot attenuation in transition
% 1:12  11/2/99
%
%   Authors: Dejan V. Tosic, Miroslav D. Lutovac, 1999.02.08
%                 [e-mail]
%                 [e-mail]
%
%   Copyright (c) 1999 by Tosic & Lutovac
%   $Revision: 1.0 $  $Date: 1999/02/08 03:07:42 $
%
%   References:
%        Miroslav D. Lutovac, Dejan V. Tosic, Brian L. Evans
%           Advanced Filter Design for Signal Processing
%                   Using MATLAB and Mathematica

initaxis = [0.13 0.11 0.775 0.815];

twoaxis = 0;
if moreaxis ~= 0
 delete(gca);
 set(gca, 'Position', initaxis);
end

if filnumb == 1

twoaxis = 0;
if moreaxis ~= 0
 delete(gca);
 set(gca, 'Position', initaxis);
end
f1 = speca(1); f2 = speca(2);
f = f1:(f2-f1)/100:f2;
s = j*2*pi*f;
H = polyval(num,s)./polyval(den,s);
A = -20*log10(abs(H));
plot(f, A, [speca(1) speca(2)], [speca(3) speca(4)], 'm:');
xlabel('frequency (Hz)');
ylabel('attenuation (dB)');
title('Transition');
a = axis;
a(1) = f1;
a(2) = f2;
axis(a);
%set(gca,'Color',[0 0 0]);

elseif filnumb == 2

twoaxis = 0;
if moreaxis ~= 0
 delete(gca);
 set(gca, 'Position', initaxis);
end
f1 = speca(1); f2 = speca(2);
f = f1:(f2-f1)/100:f2;
s = j*2*pi*f;
H = polyval(num,s)./polyval(den,s);
A = -20*log10(abs(H));
plot(f, A, [speca(1) speca(2)], [speca(4) speca(3)], 'm:');
xlabel('frequency (Hz)');
ylabel('attenuation (dB)');
title('Transition');
a = axis;
a(1) = f1;
a(2) = f2;
axis(a);
%set(gca,'Color',[0 0 0]);

elseif filnumb == 3

twoaxis=1;
subplot(211);
f1 = speca(1); f2 = speca(2);
f = f1:(f2-f1)/100:f2;
s = j*2*pi*f;
H = polyval(num,s)./polyval(den,s);
A = -20*log10(abs(H));
plot(f,A, ...
 [speca(1) speca(2)],[speca(5) speca(6)],'m:');
xlabel('frequency (Hz)');
ylabel('attenuation (dB)');
title('Transition 1');
a = axis;
a(1) = f1;
a(2) = f2;
axis(a);
%set(gca,'Color',[0 0 0]);

subplot(212);
f1 = speca(3); f2 = speca(4);
f = f1:(f2-f1)/100:f2;
s = j*2*pi*f;
H = polyval(num,s)./polyval(den,s);
A = -20*log10(abs(H));
plot(f,A, ...
 [speca(3) speca(4)],[speca(6) speca(7)],'m:');
xlabel('frequency (Hz)');
ylabel('attenuation (dB)');
title('Transition 2');
a = axis;
a(1) = f1;
a(2) = f2;
axis(a);
%set(gca,'Color',[0 0 0]);

elseif filnumb == 4

twoaxis=1;
subplot(211);
f1 = speca(1); f2 = speca(2);
f = f1:(f2-f1)/100:f2;
s = j*2*pi*f;
H = polyval(num,s)./polyval(den,s);
A = -20*log10(abs(H));
plot(f,A, ...
 [speca(1) speca(2)],[speca(5) speca(6)],'m:');
xlabel('frequency (Hz)');
ylabel('attenuation (dB)');
title('Transition 1');
a = axis;
a(1) = f1;
a(2) = f2;
axis(a);
%set(gca,'Color',[0 0 0]);

subplot(212);
f1 = speca(3); f2 = speca(4);
f = f1:(f2-f1)/100:f2;
s = j*2*pi*f;
H = polyval(num,s)./polyval(den,s);
A = -20*log10(abs(H));
plot(f,A, ...
 [speca(3) speca(4)],[speca(6) speca(7)],'m:');
xlabel('frequency (Hz)');
ylabel('attenuation (dB)');
title('Transition 2');
a = axis;
a(1) = f1;
a(2) = f2;
axis(a);
%set(gca,'Color',[0 0 0]);

end