Diff of /filterbank/filterbankresponse.m [34f3ea] .. [931198]  Maximize  Restore

Switch to unified view

a/filterbank/filterbankresponse.m b/filterbank/filterbankresponse.m
...
...
8
%   usefull tool to investigate the behaviour of the windows, as peaks
8
%   usefull tool to investigate the behaviour of the windows, as peaks
9
%   indicate that a frequency is overrepresented in the filterbank, while
9
%   indicate that a frequency is overrepresented in the filterbank, while
10
%   a dip indicates that it is not well represented.
10
%   a dip indicates that it is not well represented.
11
%
11
%
12
%   In mathematical terms, this function computes the diagonal of the
12
%   In mathematical terms, this function computes the diagonal of the
13
%   Fourier transform of the frame operator.
13
%   Fourier transform of the frame operator when the filterbank is painless.
14
%
14
%
15
%   `filterbankresponse(g,a,L,'real')` does the same for a filterbank
15
%   `filterbankresponse(g,a,L,'real')` does the same for a filterbank
16
%   intended for positive-only filterbank.
16
%   intended for positive-only filterbank.
17
%
17
%
18
%   See also: filterbank, filterbankbounds
18
%   See also: filterbank, filterbankbounds
...
...
22
[flags,kv]=ltfatarghelper({},definput,varargin);
22
[flags,kv]=ltfatarghelper({},definput,varargin);
23
23
24
24
25
% G1 is done this way just so that we can determine the data type.
25
% G1 is done this way just so that we can determine the data type.
26
G1=comp_transferfunction(g{1},L);
26
G1=comp_transferfunction(g{1},L);
27
gf=abs(G1).^2;
27
gf=abs(G1).^2*(L/a(1));
28
28
29
M=numel(g);
29
M=numel(g);
30
  
30
  
31
for m=2:M
31
for m=2:M
32
  gf=gf+abs(comp_transferfunction(g{m},L)).^2;
32
  gf=gf+abs(comp_transferfunction(g{m},L)).^2*(L/a(m));
33
end;
33
end;
34
  
34
  
35
if flags.do_real
35
if flags.do_real
36
  gf=gf+involute(gf);   
36
  gf=gf+involute(gf);   
37
end;
37
end;
38
38
39
gf=gf/L;
40
39
if flags.do_plot
41
if flags.do_plot
40
  plotfft(gf,'lin');
42
  plotfft(gf,'lin');
41
end;
43
end;

Get latest updates about Open Source Projects, Conferences and News.

Sign up for the SourceForge newsletter:

JavaScript is required for this form.





No, thanks