## 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;
```