1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
function dgt_test_mat( filepath )
% Generate some signals and the associated dgt and save them into file
if nargin < 1
filepath = BaseFolder.get_BasePath();
end;
Lr = [24, 35, 35, 24,144,108,144,135,77,77];
ar = [ 6, 5, 5, 4, 9, 9, 12, 9, 7, 7];
Mr = [ 8, 7, 7, 6, 16, 12, 24, 9,11,11];
glr = [16, 14, 21, 12, 48, 12, 24, 18,22,11];
filename=[filepath 'dgt_signal_ex.mat'];
disp(filename)
ind = 1;
for phase=0:1
pt = 'timeinv';
if(phase == 0)
pt = 'freqinv';
end;
for ii=1:length(Lr)
L=Lr(ii);
M=Mr(ii);
a=ar(ii);
gl=glr(ii);
fprintf('L=%d, M=%d, a=%d, gl=%d\n',L,M,a,gl);
for rtype=1:2
for W=1:3
if rtype==1
rname='REAL';
f=tester_rand(L,W);
g=tester_rand(gl,1);
cc = dgtreal(f,g,a,M,L,pt);
gd = gabdual(g,a,M);
frec = idgtreal(cc,gd,a,M,L,pt);
else
rname='CMPLX';
f=tester_crand(L,W);
g=tester_crand(gl,1);
cc = dgt(f,g,a,M,L,pt);
gd = gabdual(g,a,M);
frec = idgt(cc,gd,a,L,pt);
end;
Data{ind}.phase = pt;
Data{ind}.L = L;
Data{ind}.M = M;
Data{ind}.a = a;
Data{ind}.gl = gl;
Data{ind}.rtype = rtype;
Data{ind}.W = W;
Data{ind}.rname = rname;
Data{ind}.f = f;
Data{ind}.g = g;
Data{ind}.cc = cc;
Data{ind}.gd = gd;
Data{ind}.freq = frec;
ind = ind + 1;
end;
end;
end;
end
save(filename,'Data');
end
|
