diff scripts/signal/spectral_xdf.m @ 3426:f8dde1807dee

[project @ 2000-01-13 08:40:00 by jwe]
author jwe
date Thu, 13 Jan 2000 08:40:53 +0000
parents 041ea33fbbf4
children 858695b3ed62
line wrap: on
line diff
--- a/scripts/signal/spectral_xdf.m
+++ b/scripts/signal/spectral_xdf.m
@@ -1,66 +1,66 @@
 ## Copyright (C) 1995, 1996, 1997  Friedrich Leisch
-## 
+##
 ## This program is free software; you can redistribute it and/or modify
 ## it under the terms of the GNU General Public License as published by
 ## the Free Software Foundation; either version 2, or (at your option)
 ## any later version.
-## 
+##
 ## This program is distributed in the hope that it will be useful, but
 ## WITHOUT ANY WARRANTY; without even the implied warranty of
 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
-## General Public License for more details. 
-## 
+## General Public License for more details.
+##
 ## You should have received a copy of the GNU General Public License
 ## along with this file.  If not, write to the Free Software Foundation,
 ## 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
 
 ## usage:  retval = spectral_xdf (X, [win, [b]])
-## 
+##
 ## Returns the spectral density estimator.
 ## X ....... data vector
 ## win ..... window name, eg. "triangle" or "rectangle"
 ##           spectral_adf searches for a function called win_sw ()
 ## b ....... bandwidth
-## 
+##
 ## If win is omitted, the triangle window is used as default.
 ## If b is omitted, 1 / sqrt (length (X)) is used as default.
-  
+
 ## Author:  FL <Friedrich.Leisch@ci.tuwien.ac.at>
 ## Description:  Spectral density estimation
-  
+
 function retval = spectral_xdf (X, win, b)
-  
+
   xr = length (X);
-  
+
   if (columns (X) > 1)
     X = X';
   endif
-  
+
   if (nargin < 3)
     b = 1 / ceil (sqrt (xr));
   endif
-  
+
   if (nargin == 1)
     w = triangle_sw (xr, b);
   else
     win = [win, "_sw"];
     w = feval (win, xr, b);
   endif
-  
+
   X = X - sum (X) / xr;
 
   retval = (abs (fft (X)) / xr).^2;
   retval = real (ifft (fft(retval) .* fft(w)));
-  
+
   retval = [(zeros (xr, 1)), retval];
   retval(:, 1) = (0 : xr-1)' / xr;
-  
+
 endfunction
-    
-
-  
 
 
 
 
 
+
+
+