1 files changed, 44 insertions, 44 deletions
diff --git a/macros/rceps.sci b/macros/rceps.sci
index f0433a3..831611f 100644
--- a/macros/rceps.sci
+++ b/macros/rceps.sci
@@ -1,48 +1,48 @@
+ /*Description
+        Real cepstrum and minimum-phase reconstruction
+        If  called with two output arguments,
+        the minimum phase reconstruction of the signal x is returned in ym.Calling Sequence
+        [y,ym] = rceps(x)
+    Parameters
+        x : A vector or a Matirx
+        y : Real cepstrum
+        ym : Minimum-phase reconstruction
+    Dependencies fft1 , ifft1
+    Example:
+        // create a 45 Hz sine wave sampled at 100 Hz. 
+        t = 0:0.01:1.27;
+        s1 = sin(2*%pi*45*t);
+        s2 = s1 + 0.5*[zeros(1,20) s1(1:108)]; //Add an echo of the signal, with half the amplitude, 0.2 seconds after the beginning of the signal.
+        c = rceps(s2); // real cepstrum of signal 
+        plot(t,c)
+*/
 function [y, xm]= rceps(x)
-//Produce the cepstrum of the signal x, and if desired, the minimum phase reconstruction of the signal x.
-//Calling Sequence
-//[y, xm] = rceps(x)
-//Parameters 
-//x: real or complex vector input
-//Produce the cepstrum of the signal x, and if desired, the minimum phase reconstruction of the signal x. If x is a matrix, do so for each column of the matrix.
-//Examples
-// f0 = 70; Fs = 10000;                   # 100 Hz fundamental, 10kHz sampling rate
-// a = poly (0.985 * exp (1i*pi*[0.1, -0.1, 0.3, -0.3])); # two formants
-// s = 0.005 * randn (1024, 1);           # Noise excitation signal
-// s(1:Fs/f0:length(s)) = 1;              # Impulse glottal wave
-// x = filter (1, a, s);                  # Speech signal in x
-// [y, xm] = rceps (x .* hanning (1024)); # cepstrum and min phase reconstruction
-
-if(argn(2)~= 1 )
-error("Wrong number of Input Arguments");
-end
-
-if(argn(1)>2)
-error("Wrong number of Output Arguments")
-end
-  f = abs(fft1(x));
-  if (or(f == 0))
-    error ("The spectrum of x contains zeros, unable to compute real cepstrum");
-  end
-   
-  y = real(ifft1(log(f)));
-  
-  if argn(1) == 2 then
-    n=length(x);
-    if size(x,1)==1 then
-      if (n-fix(n./2).*2) ==1 then
-        xm = [y(1), 2*y(2:n/2+1), zeros(1,n/2)];
-      else
-        xm = [y(1), 2*y(2:n/2), y(n/2+1), zeros(1,n/2-1)];
+    if(argn(2)~= 1 )
+    error("Wrong number of Input Arguments");
+    end
+    if(argn(1)>2)
+    error("Wrong number of Output Arguments")
+    end
+      f = abs(fft1(x));
+      if (or(f == 0))
+        error ("The spectrum of x contains zeros, unable to compute real cepstrum");
       end
-    else
-      if (n-fix(n./2).*2)==1
-        xm = [y(1,:); 2*y(2:n/2+1,:); zeros(n/2,size(y,2))];
-      else
-        xm = [y(1,:); 2*y(2:n/2,:); y(n/2+1,:); zeros(n/2-1,size(y,2))];
+      y = real(ifft1(log(f)));
+      if argn(1) == 2 then
+        n=max(size(x));
+        if size(x,1)==1 then
+          if (n-fix(n./2).*2) ==1 then
+            xm = [y(1), 2*y(2:n/2+1), zeros(1,n/2)];
+          else
+            xm = [y(1), 2*y(2:n/2), y(n/2+1), zeros(1,n/2-1)];
+          end
+        else
+          if (n-fix(n./2).*2)==1
+            xm = [y(1,:); 2*y(2:n/2+1,:); zeros(n/2,size(y,2))];
+          else
+            xm = [y(1,:); 2*y(2:n/2,:); y(n/2+1,:); zeros(n/2-1,size(y,2))];
+          end
+        end
+        xm = real(ifft1(exp(fft1(xm))));
       end
-    end
-    xm = real(ifft1(exp(fft1(xm)')));
-  end
-  
 endfunction