dfdspace.m

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

% dfdspace.m  Advanced Digital Filter Design Space (main script)
%
% 17:13  17/8/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:
%% dfdsbutt.m, dfdstabl.m

clear all; close all; clc; disp('Advanced Filter Design Space');

% demo lowpass attenuation limits specification
%        Fpass, Fstop, Apass, Astop
%                       (dB)   (dB)

speca = [0.2,  0.212,  0.2,   40];

% demo design space

dstable = ...
[8   1.04285  1.09245  0.0625785    0.217091  0.19406    0.203289
 9   1.022    1.11016  0.0143176    0.217091  0.190964   0.207436
10   1.01135  1.13668  0.002622     0.217091  0.186508   0.20962
11   1.00587  1.17518  0.000365604  0.217091  0.180398   0.210762
12   1.00304  1.23116  0.000036102  0.217091  0.172196   0.211357
13   1.00158  1.31502  2.23768e-6   0.217091  0.161215   0.211666
14   1.00082  1.44883  6.96606e-8   0.217091  0.146325   0.211827
15   1.00042  1.69028  6.68433e-10  0.217091  0.125423   0.21191
16   1.00022  2.269    4.3985e-13   0.217091  0.0934332  0.211953];

Nmin = min(dstable(:,1));
Nmax = max(dstable(:,1));
Nlength = length(dstable(:,1));
Nindex = 0;
Xmin = min(dstable(:,2));
Xmax = max(dstable(:,3));
Emin = min(dstable(:,4));
Emax = max(dstable(:,5));
FPmin = min(dstable(:,6));
FPmax = max(dstable(:,7));

Xmin1 = dstable(1,2);
Xmax1 = dstable(1,3);
Emin1 = dstable(1,4);
Emax1 = dstable(1,5);
FPmin1 = dstable(1,6);
FPmax1 = dstable(1,7);

cubeaxis = [1, Xmax*1.01, 0, Emax*1.01, FPmin*0.99, FPmax*1.01];
caseaxis = [Xmin1*0.99, Xmax1, 0, Emax1*1.01, FPmin1, FPmax1];

initview = [-40,30];
cubeview = initview;
caseview = [-20 40];

disp('Lowpass attenuation limits specification')
disp(['Fp = ', num2str(speca(1))])
disp(['Fs = ', num2str(speca(2))])
disp(['Ap = ', num2str(speca(3)), ' dB'])
disp(['As = ', num2str(speca(4)), ' dB'])

dfdstabl(dstable);

%-----------------------------------------------------------------%
fig1 = figure(1); 
axis off;
set(fig1, 'Name', 'Advanced Digital Filter Design Space', ...
          'NumberTitle', 'off')

initsize = get(fig1,'Position');
svgasize = [0 0 800 562];

dfdsbutt;

text(0,1.05,'Advanced Filter Design for Signal Processing', 'FontWeight', 'bold')
text(0,1.0,'Miroslav D. Lutovac, Dejan V. Tosic, Brian L. Evans', 'FontWeight', 'bold')
text(0.0,0.92,'Demo = load demo attenuation limits specification')
text(0.0,0.86,'User = load attenuation limits specification from userspec.dat')
text(0.0,0.80,'View = view current attenuation limits specification')
text(0.0,0.72,'Space = compute design space')
text(0.0,0.66,'Table = display design space in tabular form')
text(0.0,0.58,'plane1 = selectivity-ripple design space planes')
text(0.0,0.52,'plane2 = selectivity-frequency design space planes')
text(0.0,0.46,'plane3 = ripple-frequency design space planes')
text(0.0,0.38,'cases = basic design cases')
text(0.0,0.30,'cubes = selectivity-ripple-frequency design space cuboids')
text(0.0,0.24,'cube = selectivity-ripple-frequency design subspace cuboid')
text(0.0,0.18,'       cube+/cube- increment/decrement order by one')
text(0.0,0.10,'grid = toggle grid on/off')
text(0.0,0.04,'full/small = full/small screen display')
text(0.0,-0.02,'AZ/EL/2D/3D = set view angle azimuth/elevation/2D/3D')