File Coverage

blib/lib/Math/DWT/Wavelet/Coiflet.pm
Criterion Covered Total %
statement 0 43 0.0
branch 0 12 0.0
condition n/a
subroutine 0 7 0.0
pod 7 7 100.0
total 7 69 10.1


line stmt bran cond sub pod time code
1             package Math::DWT::Wavelet::Coiflet;
2              
3              
4             =head1 NAME
5              
6             Math::DWT::Wavelet::Coiflet - FIR lo- & hi-pass filter coefficients for the Coiflet wavelet.
7              
8             =head1 VERSION
9              
10             Version 0.01
11              
12             =cut
13              
14             our $VERSION = '0.01';
15              
16             =head1 SYNOPSIS
17              
18             This module provides the lo- and hi-pass decomposition and reconstruction filter coefficients for the Coiflet wavelet. It is meant to be used with other Math::DWT modules:
19              
20             use Math::DWT;
21             use Math::DWT::UDWT;
22              
23             my $dwt = Math::DWT->new('Coiflet');
24             my $udwt = Math::DWT::UDWT->new('Coiflet');
25              
26              
27             =cut
28              
29             =head1 SUBROUTINES/METHODS
30              
31             =head2 new(VAR)
32              
33             Takes a single argument, var, which determines the particular filter set which the object will provide. You can get a list of available filter sets with the vars() method.
34              
35             my $coiflet4 = Math::DWT::Wavelet::Coiflet->new(4);
36              
37             This method returns a Math::DWT::Wavelet::Coiflet object;
38              
39             =head2 vars()
40              
41             This method returns a list of possible choices for VAR when creating a new object, e.g.:
42              
43             my @v = Math::DWT::Wavelet::Coiflet->vars();
44             print scalar(@v); # 5
45              
46              
47             =head2 filters()
48              
49             Depending on the context in which it is called, returns an array or an arrayref containing (lo_d, hi_d, lo_r, hi_r) - the set of filters which are defined with the instantiated object.
50              
51              
52             =head2 lo_d()
53              
54             =head2 hi_d()
55              
56             =head2 lo_r()
57              
58             =head2 hi_r()
59              
60             Returns the requested set of filter coefficients as either an array or arrayref, depending on calling context.
61              
62             =head1 SEE ALSO
63              
64             Math::DWT(3pm), Math::DWT::UDWT(3pm), Math::DWT::Wavelet::Haar(3pm), Math::DWT::Wavelet::Daubechies(3pm), Math::DWT::Wavelet::Symlet(3pm), Math::DWT::Wavelet::Biorthogonal(3pm), Math::DWT::Wavelet::ReverseBiorthogonal(3pm), Math::DWT::Wavelet::DiscreteMeyer(3pm), perl(1)
65              
66              
67             =head1 AUTHOR
68              
69              
70             Mike Kroh, C<< >>
71              
72             =head1 BUGS
73              
74              
75             Please report any bugs or feature requests to C, or through the web interface at L. I will be notified, and then you'll automatically be notified of progress on your bug as I make changes.
76              
77              
78              
79             =head1 ACKNOWLEDGEMENTS
80              
81             These wavelet filter coefficients were scraped from this site: L.
82              
83             =head1 LICENSE AND COPYRIGHT
84              
85              
86             Copyright 2016 Mike Kroh.
87              
88             This program is free software; you can redistribute it and/or modify it
89             under the terms of the the Artistic License (2.0). You may obtain a
90             copy of the full license at:
91              
92             L
93              
94             Any use, modification, and distribution of the Standard or Modified
95             Versions is governed by this Artistic License. By using, modifying or
96             distributing the Package, you accept this license. Do not use, modify,
97             or distribute the Package, if you do not accept this license.
98              
99             If your Modified Version has been derived from a Modified Version made
100             by someone other than you, you are nevertheless required to ensure that
101             your Modified Version complies with the requirements of this license.
102              
103             This license does not grant you the right to use any trademark, service
104             mark, tradename, or logo of the Copyright Holder.
105              
106             This license includes the non-exclusive, worldwide, free-of-charge
107             patent license to make, have made, use, offer to sell, sell, import and
108             otherwise transfer the Package with respect to any patent claims
109             licensable by the Copyright Holder that are necessarily infringed by the
110             Package. If you institute patent litigation (including a cross-claim or
111             counterclaim) against any party alleging that the Package constitutes
112             direct or contributory patent infringement, then this Artistic License
113             to you shall terminate on the date that such litigation is filed.
114              
115             Disclaimer of Warranty: THE PACKAGE IS PROVIDED BY THE COPYRIGHT HOLDER
116             AND CONTRIBUTORS "AS IS' AND WITHOUT ANY EXPRESS OR IMPLIED WARRANTIES.
117             THE IMPLIED WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR
118             PURPOSE, OR NON-INFRINGEMENT ARE DISCLAIMED TO THE EXTENT PERMITTED BY
119             YOUR LOCAL LAW. UNLESS REQUIRED BY LAW, NO COPYRIGHT HOLDER OR
120             CONTRIBUTOR WILL BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, OR
121             CONSEQUENTIAL DAMAGES ARISING IN ANY WAY OUT OF THE USE OF THE PACKAGE,
122             EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
123              
124              
125             =cut
126              
127              
128              
129             my @vars=qw/1 2 3 4 5/;
130              
131             my %lo_d; my %hi_d; my %lo_r; my %hi_r;
132              
133              
134             $lo_d{1}=[qw/-0.01565572813546454 -0.0727326195128539 0.38486484686420286 0.8525720202122554 0.3378976624578092 -0.0727326195128539/];
135             $hi_d{1}=[qw/0.0727326195128539 0.3378976624578092 -0.8525720202122554 0.38486484686420286 0.0727326195128539 -0.01565572813546454/];
136             $lo_r{1}=[qw/-0.0727326195128539 0.3378976624578092 0.8525720202122554 0.38486484686420286 -0.0727326195128539 -0.01565572813546454/];
137             $hi_r{1}=[qw/-0.01565572813546454 0.0727326195128539 0.38486484686420286 -0.8525720202122554 0.3378976624578092 0.0727326195128539/];
138              
139             $lo_d{2}=[qw/-0.0007205494453645122 -0.0018232088707029932 0.0056114348193944995 0.023680171946334084 -0.0594344186464569 -0.0764885990783064 0.41700518442169254 0.8127236354455423 0.3861100668211622 -0.06737255472196302 -0.04146493678175915 0.016387336463522112/];
140             $hi_d{2}=[qw/-0.016387336463522112 -0.04146493678175915 0.06737255472196302 0.3861100668211622 -0.8127236354455423 0.41700518442169254 0.0764885990783064 -0.0594344186464569 -0.023680171946334084 0.0056114348193944995 0.0018232088707029932 -0.0007205494453645122/];
141             $lo_r{2}=[qw/0.016387336463522112 -0.04146493678175915 -0.06737255472196302 0.3861100668211622 0.8127236354455423 0.41700518442169254 -0.0764885990783064 -0.0594344186464569 0.023680171946334084 0.0056114348193944995 -0.0018232088707029932 -0.0007205494453645122/];
142             $hi_r{2}=[qw/-0.0007205494453645122 0.0018232088707029932 0.0056114348193944995 -0.023680171946334084 -0.0594344186464569 0.0764885990783064 0.41700518442169254 -0.8127236354455423 0.3861100668211622 0.06737255472196302 -0.04146493678175915 -0.016387336463522112/];
143              
144             $lo_d{3}=[qw/-3.459977283621256e-05 -7.098330313814125e-05 0.0004662169601128863 0.0011175187708906016 -0.0025745176887502236 -0.00900797613666158 0.015880544863615904 0.03455502757306163 -0.08230192710688598 -0.07179982161931202 0.42848347637761874 0.7937772226256206 0.4051769024096169 -0.06112339000267287 -0.0657719112818555 0.023452696141836267 0.007782596427325418 -0.003793512864491014/];
145             $hi_d{3}=[qw/0.003793512864491014 0.007782596427325418 -0.023452696141836267 -0.0657719112818555 0.06112339000267287 0.4051769024096169 -0.7937772226256206 0.42848347637761874 0.07179982161931202 -0.08230192710688598 -0.03455502757306163 0.015880544863615904 0.00900797613666158 -0.0025745176887502236 -0.0011175187708906016 0.0004662169601128863 7.098330313814125e-05 -3.459977283621256e-05/];
146             $lo_r{3}=[qw/-0.003793512864491014 0.007782596427325418 0.023452696141836267 -0.0657719112818555 -0.06112339000267287 0.4051769024096169 0.7937772226256206 0.42848347637761874 -0.07179982161931202 -0.08230192710688598 0.03455502757306163 0.015880544863615904 -0.00900797613666158 -0.0025745176887502236 0.0011175187708906016 0.0004662169601128863 -7.098330313814125e-05 -3.459977283621256e-05/];
147             $hi_r{3}=[qw/-3.459977283621256e-05 7.098330313814125e-05 0.0004662169601128863 -0.0011175187708906016 -0.0025745176887502236 0.00900797613666158 0.015880544863615904 -0.03455502757306163 -0.08230192710688598 0.07179982161931202 0.42848347637761874 -0.7937772226256206 0.4051769024096169 0.06112339000267287 -0.0657719112818555 -0.023452696141836267 0.007782596427325418 0.003793512864491014/];
148              
149             $lo_d{4}=[qw/-1.7849850030882614e-06 -3.2596802368833675e-06 3.1229875865345646e-05 6.233903446100713e-05 -0.00025997455248771324 -0.0005890207562443383 0.0012665619292989445 0.003751436157278457 -0.00565828668661072 -0.015211731527946259 0.025082261844864097 0.03933442712333749 -0.09622044203398798 -0.06662747426342504 0.4343860564914685 0.782238930920499 0.41530840703043026 -0.05607731331675481 -0.08126669968087875 0.026682300156053072 0.016068943964776348 -0.0073461663276420935 -0.0016294920126017326 0.0008923136685823146/];
150             $hi_d{4}=[qw/-0.0008923136685823146 -0.0016294920126017326 0.0073461663276420935 0.016068943964776348 -0.026682300156053072 -0.08126669968087875 0.05607731331675481 0.41530840703043026 -0.782238930920499 0.4343860564914685 0.06662747426342504 -0.09622044203398798 -0.03933442712333749 0.025082261844864097 0.015211731527946259 -0.00565828668661072 -0.003751436157278457 0.0012665619292989445 0.0005890207562443383 -0.00025997455248771324 -6.233903446100713e-05 3.1229875865345646e-05 3.2596802368833675e-06 -1.7849850030882614e-06/];
151             $lo_r{4}=[qw/0.0008923136685823146 -0.0016294920126017326 -0.0073461663276420935 0.016068943964776348 0.026682300156053072 -0.08126669968087875 -0.05607731331675481 0.41530840703043026 0.782238930920499 0.4343860564914685 -0.06662747426342504 -0.09622044203398798 0.03933442712333749 0.025082261844864097 -0.015211731527946259 -0.00565828668661072 0.003751436157278457 0.0012665619292989445 -0.0005890207562443383 -0.00025997455248771324 6.233903446100713e-05 3.1229875865345646e-05 -3.2596802368833675e-06 -1.7849850030882614e-06/];
152             $hi_r{4}=[qw/-1.7849850030882614e-06 3.2596802368833675e-06 3.1229875865345646e-05 -6.233903446100713e-05 -0.00025997455248771324 0.0005890207562443383 0.0012665619292989445 -0.003751436157278457 -0.00565828668661072 0.015211731527946259 0.025082261844864097 -0.03933442712333749 -0.09622044203398798 0.06662747426342504 0.4343860564914685 -0.782238930920499 0.41530840703043026 0.05607731331675481 -0.08126669968087875 -0.026682300156053072 0.016068943964776348 0.0073461663276420935 -0.0016294920126017326 -0.0008923136685823146/];
153              
154             $lo_d{5}=[qw/-9.517657273819165e-08 -1.6744288576823017e-07 2.0637618513646814e-06 3.7346551751414047e-06 -2.1315026809955787e-05 -4.134043227251251e-05 0.00014054114970203437 0.00030225958181306315 -0.0006381313430451114 -0.0016628637020130838 0.0024333732126576722 0.006764185448053083 -0.009164231162481846 -0.01976177894257264 0.03268357426711183 0.0412892087501817 -0.10557420870333893 -0.06203596396290357 0.4379916261718371 0.7742896036529562 0.4215662066908515 -0.05204316317624377 -0.09192001055969624 0.02816802897093635 0.023408156785839195 -0.010131117519849788 -0.004159358781386048 0.0021782363581090178 0.00035858968789573785 -0.00021208083980379827/];
155             $hi_d{5}=[qw/0.00021208083980379827 0.00035858968789573785 -0.0021782363581090178 -0.004159358781386048 0.010131117519849788 0.023408156785839195 -0.02816802897093635 -0.09192001055969624 0.05204316317624377 0.4215662066908515 -0.7742896036529562 0.4379916261718371 0.06203596396290357 -0.10557420870333893 -0.0412892087501817 0.03268357426711183 0.01976177894257264 -0.009164231162481846 -0.006764185448053083 0.0024333732126576722 0.0016628637020130838 -0.0006381313430451114 -0.00030225958181306315 0.00014054114970203437 4.134043227251251e-05 -2.1315026809955787e-05 -3.7346551751414047e-06 2.0637618513646814e-06 1.6744288576823017e-07 -9.517657273819165e-08/];
156             $lo_r{5}=[qw/-0.00021208083980379827 0.00035858968789573785 0.0021782363581090178 -0.004159358781386048 -0.010131117519849788 0.023408156785839195 0.02816802897093635 -0.09192001055969624 -0.05204316317624377 0.4215662066908515 0.7742896036529562 0.4379916261718371 -0.06203596396290357 -0.10557420870333893 0.0412892087501817 0.03268357426711183 -0.01976177894257264 -0.009164231162481846 0.006764185448053083 0.0024333732126576722 -0.0016628637020130838 -0.0006381313430451114 0.00030225958181306315 0.00014054114970203437 -4.134043227251251e-05 -2.1315026809955787e-05 3.7346551751414047e-06 2.0637618513646814e-06 -1.6744288576823017e-07 -9.517657273819165e-08/];
157             $hi_r{5}=[qw/-9.517657273819165e-08 1.6744288576823017e-07 2.0637618513646814e-06 -3.7346551751414047e-06 -2.1315026809955787e-05 4.134043227251251e-05 0.00014054114970203437 -0.00030225958181306315 -0.0006381313430451114 0.0016628637020130838 0.0024333732126576722 -0.006764185448053083 -0.009164231162481846 0.01976177894257264 0.03268357426711183 -0.0412892087501817 -0.10557420870333893 0.06203596396290357 0.4379916261718371 -0.7742896036529562 0.4215662066908515 0.05204316317624377 -0.09192001055969624 -0.02816802897093635 0.023408156785839195 0.010131117519849788 -0.004159358781386048 -0.0021782363581090178 0.00035858968789573785 0.00021208083980379827/];
158             ;
159              
160             sub new {
161 0     0 1   my $class=shift;
162 0           my $self={};
163 0           my $var=shift;
164              
165             $self={lo_d=>$lo_d{$var},
166             hi_d=>$hi_d{$var},
167             lo_r=>$lo_r{$var},
168 0           hi_r=>$hi_r{$var}
169             };
170            
171 0           bless $self,$class;
172 0           return $self;
173             };
174              
175             sub vars {
176 0     0 1   my $self=shift;
177 0 0         if (wantarray()) {
178 0           return (@vars);
179             };
180 0           return \@vars;
181             };
182              
183             sub filters {
184 0     0 1   my $self=shift;
185 0           my $lo_d=$self->lo_d;
186 0           my $hi_d=$self->hi_d;
187 0           my $lo_r=$self->lo_r;
188 0           my $hi_r=$self->hi_r;
189 0           my @a=( $lo_d,$hi_d,$lo_r,$hi_r);
190 0 0         if (wantarray()) {
191 0           return (@a);
192             };
193 0           return \@a;
194             };
195              
196             sub lo_d {
197 0     0 1   my $self=shift;
198 0           my $a=$self->{lo_d};
199 0 0         if (wantarray()) {
200 0           return (@{$a});
  0            
201             };
202 0           return $a;
203             };
204             sub hi_d {
205 0     0 1   my $self=shift;
206 0           my $a=$self->{hi_d};
207 0 0         if (wantarray()) {
208 0           return (@{$a});
  0            
209             };
210 0           return $a;
211             };
212             sub lo_r {
213 0     0 1   my $self=shift;
214 0           my $a=$self->{lo_r};
215 0 0         if (wantarray()) {
216 0           return (@{$a});
  0            
217             };
218 0           return $a;
219             };
220             sub hi_r {
221 0     0 1   my $self=shift;
222 0           my $a=$self->{hi_r};
223 0 0         if (wantarray()) {
224 0           return (@{$a});
  0            
225             };
226 0           return $a;
227             };
228              
229            
230             1;