From 9fb41cd15692824bf47ff815908dc9ee33546e18 Mon Sep 17 00:00:00 2001
From: OverLordGoldDragon
 <16495490+OverLordGoldDragon@users.noreply.github.com>
Date: Wed, 2 Dec 2020 06:28:41 +0400
Subject: [PATCH] Include Nyquist bin, halve it for even lengths

---
 jWavelet/morsewave.m | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/jWavelet/morsewave.m b/jWavelet/morsewave.m
index 6ffa24f..fd14526 100644
--- a/jWavelet/morsewave.m
+++ b/jWavelet/morsewave.m
@@ -256,7 +256,10 @@
 X=vswap(X,inf,0);
 
 ommat=vrep(vrep(om,size(X,3),3),size(X,2),2);
-Xr=X.*rot(ommat.*(N+1)/2*fact); %ensures wavelets are centered 
+Xr=X.*(-1).^[0:N-1]'; %ensures wavelets are centered
+if ~mod(N, 2)
+  Xr(N/2+1, :)/=2; %ensures proper wavelet decay and analyticity
+endif
 %figure,plot(vrep(om,size(X,2),2).*(N+1)/2*fact)
 
 x=ifft(Xr);
@@ -267,7 +270,7 @@
 r=(2*be+1)./ga;
 c=r-1;
 L=0*om;
-index=(1:round(N/2));
+index=(1:floor(N/2)+1);
 psif=zeros(length(psizero),1,K);
 
 for k=0:K-1