afdscase.m

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

function h = afdscase(dstable,speca)

% afdscase.m  AFDS basic design cases
%
% 17:13  17/2/99
%
%   Authors: Dejan V. Tosic, Miroslav D. Lutovac, 1999.02.16
%                 [e-mail]
%                 [e-mail]
%
%   Copyright (c) 1999 by Tosic & Lutovac
%   $Revision: 1.0 $  $Date: 1999/02/16 21:27:42 $
%
%   References:
%        Miroslav D. Lutovac, Dejan V. Tosic, Brian L. Evans
%           Advanced Filter Design for Signal Processing
%                   Using MATLAB and Mathematica

%  calls:
%% afdswire.m

Fp = speca(1);
Fs = speca(2);
Kp = afda2k(speca(3));
Ks = afda2k(speca(4));

% n = afdnell(Fp,Fs,Kp,Ks);
% ksimin = afdxmin(n,Kp,Ks);
% ksimax = afdxmax(n,Fp,Fs,Kp,Ks);
% epsilonmin = Ks/afdl(n,ksimax);
% epsilonmax = Kp;
% fpmin = Fs/ksimax;
% fpmax = Fs/ksimin;

n = dstable(1,1);
Xmin = dstable(1,2);
Xmax = dstable(1,3);
Ymin = dstable(1,4);
Ymax = dstable(1,5);
Zmin = dstable(1,6);
Zmax = dstable(1,7);

ksimin = Xmin;
ksimax = Xmax;
epsilonmin = Ymin;
epsilonmax = Ymax;
fpmin = Zmin;
fpmax = Zmax;

afdswire(Xmin,Xmax,Ymin,Ymax,Zmin,Zmax,'b');
hold on;

xlabel('x', 'FontName','Symbol', 'FontSize',18);
ylabel('e', 'FontName','Symbol', 'FontSize',18);
zlabel('fp', 'FontName','Times', 'FontSize',18);
title('basic design cases');

plot3(...
      Fs/Fp,  epsilonmax,       Fp,    'go', ...
      Fs/Fp,  Ks/afdl(n,Fs/Fp), Fp,    'co', ...
      ksimin, epsilonmax,       Fp,    'rx', ...
      ksimin, epsilonmax,       fpmax, 'rx', ...
      ksimax, epsilonmax,       Fp,    'y*', ...
      ksimax, epsilonmin,       fpmin, 'y*'  ...
);

text( Fs/Fp,  epsilonmax,       Fp,    ' D1',  'VerticalAlignment', 'bottom' );
text( Fs/Fp,  Ks/afdl(n,Fs/Fp), Fp,    ' D2',  'VerticalAlignment', 'bottom' );
text( ksimin, epsilonmax,       Fp,    ' D3a', 'VerticalAlignment', 'bottom' );
text( ksimin, epsilonmax,       fpmax, ' D3b', 'VerticalAlignment', 'bottom' );
text( ksimax, epsilonmax,       Fp,    ' D4a', 'VerticalAlignment', 'bottom' );
text( ksimax, epsilonmin,       fpmin, ' D4b', 'VerticalAlignment', 'bottom' );

hold off;