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
|
function test_failed=test_pbspline
Lr=[15,16,18,20];
ar=[ 3, 4, 6, 5];
or=[1, 1.5, 2,3];
%-*- texinfo -*-
%@deftypefn {Function} test_pbspline
%@verbatim
%btypes={'ed','xd','stard','ec','xc','starc'};
%@end verbatim
%@strong{Url}: @url{http://ltfat.github.io/doc/testing/test_pbspline.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/>.
btypes={'ed','xd','stard'};
centtypes={'wp','hp'};
test_failed=0;
disp(' =============== TEST_PBSPLINE ============');
for ii=1:length(Lr)
L=Lr(ii);
a=ar(ii);
N=L/a;
for jj=1:length(or)
order=or(jj);
for kk=1:numel(btypes)
btype=btypes{kk};
for ll=1:2
centstring=centtypes{ll};
[g,nlen]=pbspline(L,order,a,btype,centstring);
A=zeros(L,1);
for n=0:N-1
A=A+circshift(g,n*a);
end;
res=max(abs(A-1/sqrt(a)));
[test_failed,fail]=ltfatdiditfail(res,test_failed);
s=sprintf('PBSPLINE PU %2s %s L:%3i a:%3i o:%3.5g %0.5g %s', ...
btype,centstring,L,a,order,res,fail);
disp(s);
gcutextend=middlepad(middlepad(g,nlen,centstring),L,centstring);
res=norm(g-gcutextend);
[test_failed,fail]=ltfatdiditfail(res,test_failed);
s=sprintf('PBSPLINE NLEN %2s %s L:%3i a:%3i o:%3.5g %0.5g %s', ...
btype,centstring,L,a,order,res,fail);
disp(s);
end;
end;
end;
end;
|
