Diff of /wavelets/wfilt_lemarie.m [77a22a] .. [95fe2b]  Maximize  Restore

Switch to side-by-side view

--- a/wavelets/wfilt_lemarie.m
+++ b/wavelets/wfilt_lemarie.m
@@ -1,44 +1,44 @@
-function [h,g,a]=wfilt_lemarie(N)
-%WFILT_LEMARIE  Battle and Lemarie filters
-%   Usage: [h,g,a]=wfilt_lemarie(N)
-%
-%   Input parameters:
-%         N     : Filter length.
-%
-%   `[h,g,a]=wfilt_lemarie(N)` calculates coeficients of orthonormal
-%   Battle-Lemarie wavelets. Filter coefficients are obtainded by
-%   frequency domain sampling and trunctating the impulse response.
-%
-%   References: mallat89atheory
-%   
-%   %   Examples:
-%   ---------
-%
-%   Frequency responses of the analysis filters:::  
-%
-%      w = fwtinit({'lemarie',40});
-%      wtfftfreqz(w.h);
-%
-% Original copyright goes to:
-% Copyright (C) 1994, 1995, 1996, by Universidad de Vigo 
-% Author: Jose Martin Garcia
-% e-mail: Uvi_Wave@tsc.uvigo.es
-
-num_coefs = N;
-L = 1024;
-H = wfreq_lemarie(L);
-hh=real(ifft(H{1},L));
-hh=[ hh(L-floor(num_coefs/2)+1:L) hh(1:ceil(num_coefs/2))];
-hh=sqrt(2)/sum(hh)*hh;
-
-g{1} = fliplr(hh);
-g{2} = -(-1).^(1:length(hh)).*g{1}(end:-1:1);
- 
-h{1}=g{1}(length(g{1}):-1:1);
-h{2}=g{2}(length(g{2}):-1:1);
-
-a= [2;2];
-
-
-
-
+function [h,g,a]=wfilt_lemarie(N)
+%WFILT_LEMARIE  Battle and Lemarie filters
+%   Usage: [h,g,a]=wfilt_lemarie(N)
+%
+%   Input parameters:
+%         N     : Filter length.
+%
+%   `[h,g,a]=wfilt_lemarie(N)` calculates coeficients of orthonormal
+%   Battle-Lemarie wavelets. Filter coefficients are obtainded by
+%   frequency domain sampling and trunctating the impulse response.
+%
+%   References: mallat89atheory
+%   
+%   %   Examples:
+%   ---------
+%
+%   Frequency responses of the analysis filters:::  
+%
+%      w = fwtinit({'lemarie',40});
+%      wtfftfreqz(w.h);
+%
+% Original copyright goes to:
+% Copyright (C) 1994, 1995, 1996, by Universidad de Vigo 
+% Author: Jose Martin Garcia
+% e-mail: Uvi_Wave@tsc.uvigo.es
+
+num_coefs = N;
+L = 1024;
+H = wfreq_lemarie(L);
+hh=real(ifft(H{1},L));
+hh=[ hh(L-floor(num_coefs/2)+1:L) hh(1:ceil(num_coefs/2))];
+hh=sqrt(2)/sum(hh)*hh;
+
+g{1} = fliplr(hh);
+g{2} = -(-1).^(1:length(hh)).*g{1}(end:-1:1);
+ 
+h{1}=g{1}(length(g{1}):-1:1);
+h{2}=g{2}(length(g{2}):-1:1);
+
+a= [2;2];
+
+
+
+

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

Sign up for the SourceForge newsletter:

JavaScript is required for this form.





No, thanks