demods.m

来自「AFD - Advanced Filter Design using MATLA」· M 代码 · 共 87 行

clear all; close all
fig1 = figure;

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];

set(fig1, 'Name', 'Advanced Filter Design Digital Space', 'NumberTitle', 'off')

%echo on
Nindex=dfdscubp(dstable,Nindex,Nlength);cubeview=initview;view(cubeview);axis(cubeaxis);pause(3)
Nindex=dfdscubp(dstable,Nindex,Nlength);cubeview=initview;view(cubeview);axis(cubeaxis);pause(1)
Nindex=dfdscubp(dstable,Nindex,Nlength);cubeview=initview;view(cubeview);axis(cubeaxis);pause(1)
Nindex=dfdscubp(dstable,Nindex,Nlength);cubeview=initview;view(cubeview);axis(cubeaxis);pause(1)
grid;pause(3)
grid;pause(3)
Nindex=dfdscubm(dstable,Nindex,Nlength);cubeview=initview;view(cubeview);axis(cubeaxis);pause(1),clf
Nindex=dfdscubm(dstable,Nindex,Nlength);cubeview=initview;view(cubeview);axis(cubeaxis);pause(1),clf
Nindex=dfdscubm(dstable,Nindex,Nlength);cubeview=initview;view(cubeview);axis(cubeaxis);pause(1),clf
dfdscube(dstable); cubeview=initview; view(cubeview); axis(cubeaxis);pause(3), clf
Nindex=dfdscubm(dstable,Nindex,Nlength);cubeview=initview;view(cubeview);axis(cubeaxis);pause(3)
cubeview=cubeview+[10 0];  view(cubeview);pause(3)
cubeview=cubeview+[10 0];  view(cubeview);pause(1)
cubeview=cubeview+[10 0];  view(cubeview);pause(1)
cubeview=cubeview+[10 0];  view(cubeview);pause(3)
cubeview=cubeview+[-10 0]; view(cubeview);pause(1)
cubeview=cubeview+[-10 0]; view(cubeview);pause(1)
cubeview=cubeview+[-10 0]; view(cubeview);pause(1)
cubeview=cubeview+[-10 0]; view(cubeview);pause(3)
cubeview=cubeview+[0 -10]; view(cubeview);pause(1)
cubeview=cubeview+[0 -10]; view(cubeview);pause(1)
cubeview=cubeview+[0 -10]; view(cubeview);pause(3)
cubeview=cubeview+[0 10];  view(cubeview);pause(1)
cubeview=cubeview+[0 10];  view(cubeview);pause(1)
cubeview=cubeview+[0 10];  view(cubeview);pause(3)
dfdspla1;pause(3)
dfdspla2;pause(3)
dfdspla3;pause(3)
view(2);pause(3)
cubeview=initview; view(cubeview);pause(3)
dfdscase(dstable,speca);cubeview=caseview;view(cubeview);axis(caseaxis);pause(5), clf
echo off
disp(' ')
disp(' ')
disp(' Invoke this toolbox by executing')
disp('                         dfdspace')
disp('---------------------------------')
disp(' ')
disp(' ')
dfdspace