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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
function test_failed=test_fbwarped_framebounds
test_failed = 0;
Ls = 44100;
fs = 44100;
fmax = fs/2;
bins = 1;
fac = [1,7/8,3/4,5/8,1/2];
eigstol = 1e-4;
eigsmaxit = 100;
pcgmaxit = 150;
pcgtol = 1e-4;
warpfun = cell(4,1);
invfun = cell(4,1);
%-*- texinfo -*-
%@deftypefn {Function} test_fbwarped_framebounds
%@verbatim
% ERBlet warping
%@end verbatim
%@strong{Url}: @url{http://ltfat.github.io/doc/testing/test_fbwarped_framebounds.html}
%@end deftypefn
% Copyright (C) 2005-2016 Peter L. Soendergaard <peter@sonderport.dk>.
% This file is part of LTFAT version 2.2.0
%
% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program. If not, see <http://www.gnu.org/licenses/>.
warpfun{3} = @freqtoerb;
invfun{3} = @erbtofreq;
% constant-Q warping
warpfun{4} = @(x) 10*log(x);
invfun{4} = @(x) exp(x/10);
% sqrt-warping
warpfun{2} = @(x) sign(x).*((1+abs(x)).^(1/2)-1);
invfun{2} = @(x) sign(x).*((1+abs(x)).^2-1);
% Linear warping
warpfun{1} = @(x) x/100;
invfun{1} = @(x) 100*x;
fmin = [0,0,0,50];
A = zeros(4,length(fac));
B = ones(4,length(fac));
red = A;
for jj = 1:4
[g,a,fc,L]=warpedfilters(warpfun{jj},invfun{jj},fs,fmin(jj),fmax,bins,Ls,'bwmul',1.5,'fractional','complex');
gf=filterbankresponse(g,a,Ls); framebound_ratio = max(gf)/min(gf);
disp(['Painless system frame bound ratio: ', num2str(framebound_ratio)]);
for kk = 1:length(fac)
%[jj,kk]
atemp = a;
idx = [2:length(fc)/2,length(fc)/2+2:length(fc)];
atemp(idx,2) = ceil(atemp(idx,2).*fac(kk));
red(jj,kk) = sum(atemp(:,2)./atemp(:,1));
[g,asan]=filterbankwin(g,atemp,L,'normal');
gtemp=comp_filterbank_pre(g,asan,L,10);
gtemp{1}.H = gtemp{1}.H.*sqrt(fac(kk));
gtemp{length(fc)/2+1}.H = gtemp{length(fc)/2+1}.H.*sqrt(fac(kk));
F = frame('filterbank',gtemp,asan,numel(gtemp));
[A(jj,kk),B(jj,kk)] = framebounds(F,Ls,'tol',eigstol,'pcgtol',pcgtol,'maxit',eigsmaxit,'pcgmaxit',pcgmaxit);
end
end
disp('This is a ratio B/A. Rows - warping sunction, Cols - redundancy compared to te minimal painless case')
B./A
|
