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

  Switch to side-by-side view

--- a/wavelets/wfilt_hden.m
+++ b/wavelets/wfilt_hden.m
@@ -1,107 +1,107 @@
-function [h,g,a] = wfilt_hden(N)
-%WFILT_HDEN  Higher DENsity dwt filters (tight frame, frame)
-%   Usage: [h,g,a] = wfilt_hden(N);
-%
-%   `[h,g,a]=wfilt_hden(N)` computes Higher DENsity dwt filters (tight frame, frame).
-%
-%   References: selesnick2006higher
-%
-%   Examples:
-%   ---------
-%
-%   Frequency responses of the analysis filters:::  
-%
-%      w = fwtinit({'hden',3});
-%      wtfftfreqz(w.h);
-%
-
-a= [2;2;1];
-
-switch(N)
-case 1
-% from the paper Example 1.
-harr = [
-    0                   0                    0 
-    0.353553390593274   0.353553390593274    0.5
-    0.707106781186548   0                    -0.5
-    0.353553390593274   -0.353553390593274   0
-];
-
-case 2
-% from the paper Example 2.
-harr = [
-    0                 0                0 
-    0.189604909379    0.025752563665   0.010167956157
-    0.631450512121    0.075463998066   0.046750380120
-    0.655505518357   -0.064333341412  -0.009172584871
-    0.099615139800   -0.327704691428  -0.354664087684
-    -0.163756210215   0.228185687127   0.499004628714
-    -0.023958870736   0.252240693362  -0.192086292435
-    0.025752563665   -0.189604909379   0
-];
-
-case 3
-% from the paper Example 3.
-harr = [
-   0                 0                0
-   0.022033327573    0.048477254777   0.031294135831
-   0.015381522616    0.019991451948   0.013248398005
-  -0.088169084245   -0.304530024033  -0.311552292833
-   0.051120949834    0.165478923930   0.497594326648
-   0.574161374258    0.308884916012  -0.235117092484
-   0.717567366340   -0.214155508410  -0.020594576659
-   0.247558418377   -0.074865474330   0.015375249485
-  -0.076963057605    0.028685132531   0.009751852004
-  -0.048477254777    0.022033327573   0
-];
-
-case 4
-    % from the paper Example 5. Is not a tight frame!
-    harr = [
-       0          0           0 
-       0          0           0
-       0.027222   0.044889    0
-       0.011217   0.005671    0.027671
-      -0.112709  -0.286349    0.007159 
-       0.096078   0.235789   -0.277671 
-       0.685299   0.235789    0.485682 
-       0.685299   -0.286349  -0.277671 
-       0.096078   0.005671    0.007159 
-      -0.112709   0.044889    0.027671 
-       0.011217   0           0 
-       0.027222   0           0  
-    ];
-    h=mat2cell(harr.',[1,1,1],length(harr));
-    if(nargout>1)
-
-        garr = [
-            0          0          0
-           -0.039237   0.023794   0.011029
-           -0.073518   0.037784   0.019204
-            0.181733  -0.070538  -0.020024
-            0.638129  -0.253037  -0.269204
-            0.638129   0.261996   0.517991
-            0.181733   0.261996  -0.269204
-           -0.073518  -0.253037  -0.020024
-           -0.039237  -0.070538   0.019204
-            0          0.037784   0.011029
-            0          0.023794   0
-            0          0          0
-        ];   
-
-        g=mat2cell(garr.',[1,1,1],length(harr));
-    end;
-
-    return
-
-
-otherwise
-        error('%s: No such Higher Density Wavelet Transform Filters..',upper(mfilename));
-end
-
-h=mat2cell(harr.',[1,1,1],length(harr));
-
-if(nargout>1)
-    garr = harr(end:-1:1, :);
-    g=mat2cell(garr.',[1,1,1],length(harr));
-end
+function [h,g,a] = wfilt_hden(N)
+%WFILT_HDEN  Higher DENsity dwt filters (tight frame, frame)
+%   Usage: [h,g,a] = wfilt_hden(N);
+%
+%   `[h,g,a]=wfilt_hden(N)` computes Higher DENsity dwt filters (tight frame, frame).
+%
+%   References: selesnick2006higher
+%
+%   Examples:
+%   ---------
+%
+%   Frequency responses of the analysis filters:::  
+%
+%      w = fwtinit({'hden',3});
+%      wtfftfreqz(w.h);
+%
+
+a= [2;2;1];
+
+switch(N)
+case 1
+% from the paper Example 1.
+harr = [
+    0                   0                    0 
+    0.353553390593274   0.353553390593274    0.5
+    0.707106781186548   0                    -0.5
+    0.353553390593274   -0.353553390593274   0
+];
+
+case 2
+% from the paper Example 2.
+harr = [
+    0                 0                0 
+    0.189604909379    0.025752563665   0.010167956157
+    0.631450512121    0.075463998066   0.046750380120
+    0.655505518357   -0.064333341412  -0.009172584871
+    0.099615139800   -0.327704691428  -0.354664087684
+    -0.163756210215   0.228185687127   0.499004628714
+    -0.023958870736   0.252240693362  -0.192086292435
+    0.025752563665   -0.189604909379   0
+];
+
+case 3
+% from the paper Example 3.
+harr = [
+   0                 0                0
+   0.022033327573    0.048477254777   0.031294135831
+   0.015381522616    0.019991451948   0.013248398005
+  -0.088169084245   -0.304530024033  -0.311552292833
+   0.051120949834    0.165478923930   0.497594326648
+   0.574161374258    0.308884916012  -0.235117092484
+   0.717567366340   -0.214155508410  -0.020594576659
+   0.247558418377   -0.074865474330   0.015375249485
+  -0.076963057605    0.028685132531   0.009751852004
+  -0.048477254777    0.022033327573   0
+];
+
+case 4
+    % from the paper Example 5. Is not a tight frame!
+    harr = [
+       0          0           0 
+       0          0           0
+       0.027222   0.044889    0
+       0.011217   0.005671    0.027671
+      -0.112709  -0.286349    0.007159 
+       0.096078   0.235789   -0.277671 
+       0.685299   0.235789    0.485682 
+       0.685299   -0.286349  -0.277671 
+       0.096078   0.005671    0.007159 
+      -0.112709   0.044889    0.027671 
+       0.011217   0           0 
+       0.027222   0           0  
+    ];
+    h=mat2cell(harr.',[1,1,1],length(harr));
+    if(nargout>1)
+
+        garr = [
+            0          0          0
+           -0.039237   0.023794   0.011029
+           -0.073518   0.037784   0.019204
+            0.181733  -0.070538  -0.020024
+            0.638129  -0.253037  -0.269204
+            0.638129   0.261996   0.517991
+            0.181733   0.261996  -0.269204
+           -0.073518  -0.253037  -0.020024
+           -0.039237  -0.070538   0.019204
+            0          0.037784   0.011029
+            0          0.023794   0
+            0          0          0
+        ];   
+
+        g=mat2cell(garr.',[1,1,1],length(harr));
+    end;
+
+    return
+
+
+otherwise
+        error('%s: No such Higher Density Wavelet Transform Filters..',upper(mfilename));
+end
+
+h=mat2cell(harr.',[1,1,1],length(harr));
+
+if(nargout>1)
+    garr = harr(end:-1:1, :);
+    g=mat2cell(garr.',[1,1,1],length(harr));
+end