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

Switch to side-by-side view

--- a/wavelets/wfilt_dtree.m
+++ b/wavelets/wfilt_dtree.m
@@ -1,85 +1,85 @@
-function [h,g,a] = wfilt_dtree(N)
-%WFILT_DTREE  Dual-TREE complex wavelet transform filters
-%   Usage: [h,g,a] = wfilt_dtree(N);
-%
-%   `[h,g,a]=wfilt_dtree(N)` computes filters used in the dual-tree complex wavelet transform. 
-%
-%   Examples:
-%   ---------
-%
-%   Frequency responses of the analysis filters:::  
-%
-%      w = fwtinit({'dtree',1});
-%      wtfftfreqz(w.h);
-%
-%   References: selesnick2005dtree
-%
-
-
-switch(N)
-case 1
-harr = [
-                  0                  0
-  -0.08838834764832  -0.01122679215254
-   0.08838834764832   0.01122679215254
-   0.69587998903400   0.08838834764832
-   0.69587998903400   0.08838834764832
-   0.08838834764832  -0.69587998903400
-  -0.08838834764832   0.69587998903400
-   0.01122679215254  -0.08838834764832
-   0.01122679215254  -0.08838834764832
-                  0                  0
-];
-case 2
-harr = [
-  0.01122679215254                   0
-   0.01122679215254                  0
-  -0.08838834764832  -0.08838834764832
-   0.08838834764832  -0.08838834764832
-   0.69587998903400   0.69587998903400
-   0.69587998903400  -0.69587998903400
-   0.08838834764832   0.08838834764832
-  -0.08838834764832   0.08838834764832
-                  0   0.01122679215254
-                  0  -0.01122679215254
-];
-
-case 3
-harr = [
-   0.03516384000000                  0
-                  0                  0
-  -0.08832942000000  -0.11430184000000
-   0.23389032000000                  0
-   0.76027237000000   0.58751830000000
-   0.58751830000000  -0.76027237000000
-                  0   0.23389032000000
-  -0.11430184000000   0.08832942000000
-                  0                  0
-                  0  -0.03516384000000
-];
-
-case 4
-harr = [
-                  0  -0.03516384000000
-                  0                  0
-  -0.11430184000000   0.08832942000000
-                  0   0.23389032000000
-   0.58751830000000  -0.76027237000000
-   0.76027237000000   0.58751830000000
-   0.23389032000000                  0
-  -0.08832942000000  -0.11430184000000
-                  0                  0
-   0.03516384000000                  0
-];
-
- otherwise
-        error('%s: No such Dual-Tree Complex Wavelet Transform Filters..',upper(mfilename));
-end
-a= [2;2];
-
-h=mat2cell(harr.',[1,1],length(harr));
-
-if(nargout>1)
-    garr = harr(end:-1:1, :);
-    g=mat2cell(garr.',[1,1],length(garr));
+function [h,g,a] = wfilt_dtree(N)
+%WFILT_DTREE  Dual-TREE complex wavelet transform filters
+%   Usage: [h,g,a] = wfilt_dtree(N);
+%
+%   `[h,g,a]=wfilt_dtree(N)` computes filters used in the dual-tree complex wavelet transform. 
+%
+%   Examples:
+%   ---------
+%
+%   Frequency responses of the analysis filters:::  
+%
+%      w = fwtinit({'dtree',1});
+%      wtfftfreqz(w.h);
+%
+%   References: selesnick2005dtree
+%
+
+
+switch(N)
+case 1
+harr = [
+                  0                  0
+  -0.08838834764832  -0.01122679215254
+   0.08838834764832   0.01122679215254
+   0.69587998903400   0.08838834764832
+   0.69587998903400   0.08838834764832
+   0.08838834764832  -0.69587998903400
+  -0.08838834764832   0.69587998903400
+   0.01122679215254  -0.08838834764832
+   0.01122679215254  -0.08838834764832
+                  0                  0
+];
+case 2
+harr = [
+  0.01122679215254                   0
+   0.01122679215254                  0
+  -0.08838834764832  -0.08838834764832
+   0.08838834764832  -0.08838834764832
+   0.69587998903400   0.69587998903400
+   0.69587998903400  -0.69587998903400
+   0.08838834764832   0.08838834764832
+  -0.08838834764832   0.08838834764832
+                  0   0.01122679215254
+                  0  -0.01122679215254
+];
+
+case 3
+harr = [
+   0.03516384000000                  0
+                  0                  0
+  -0.08832942000000  -0.11430184000000
+   0.23389032000000                  0
+   0.76027237000000   0.58751830000000
+   0.58751830000000  -0.76027237000000
+                  0   0.23389032000000
+  -0.11430184000000   0.08832942000000
+                  0                  0
+                  0  -0.03516384000000
+];
+
+case 4
+harr = [
+                  0  -0.03516384000000
+                  0                  0
+  -0.11430184000000   0.08832942000000
+                  0   0.23389032000000
+   0.58751830000000  -0.76027237000000
+   0.76027237000000   0.58751830000000
+   0.23389032000000                  0
+  -0.08832942000000  -0.11430184000000
+                  0                  0
+   0.03516384000000                  0
+];
+
+ otherwise
+        error('%s: No such Dual-Tree Complex Wavelet Transform Filters..',upper(mfilename));
+end
+a= [2;2];
+
+h=mat2cell(harr.',[1,1],length(harr));
+
+if(nargout>1)
+    garr = harr(end:-1:1, :);
+    g=mat2cell(garr.',[1,1],length(garr));
 end

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

Sign up for the SourceForge newsletter:





No, thanks