Update of the license file + update of headers
Sylvestre Ledru [Tue, 30 Sep 2008 11:27:37 +0000 (11:27 +0000)]
79 files changed:
scilab/modules/signal_processing/macros/%asn.sci
scilab/modules/signal_processing/macros/%k.sci
scilab/modules/signal_processing/macros/%sn.sci
scilab/modules/signal_processing/macros/analpf.sci
scilab/modules/signal_processing/macros/bilt.sci
scilab/modules/signal_processing/macros/buttmag.sci
scilab/modules/signal_processing/macros/casc.sci
scilab/modules/signal_processing/macros/cepstrum.sci
scilab/modules/signal_processing/macros/cheb1mag.sci
scilab/modules/signal_processing/macros/cheb2mag.sci
scilab/modules/signal_processing/macros/chepol.sci
scilab/modules/signal_processing/macros/convol.sci
scilab/modules/signal_processing/macros/cspect.sci
scilab/modules/signal_processing/macros/czt.sci
scilab/modules/signal_processing/macros/detrend.sci
scilab/modules/signal_processing/macros/dft.sci
scilab/modules/signal_processing/macros/ell1mag.sci
scilab/modules/signal_processing/macros/eqfir.sci
scilab/modules/signal_processing/macros/eqiir.sci
scilab/modules/signal_processing/macros/faurre.sci
scilab/modules/signal_processing/macros/fft2.sci
scilab/modules/signal_processing/macros/fftshift.sci
scilab/modules/signal_processing/macros/find_freq.sci
scilab/modules/signal_processing/macros/findm.sci
scilab/modules/signal_processing/macros/frfit.sci
scilab/modules/signal_processing/macros/frmag.sci
scilab/modules/signal_processing/macros/fsfirlin.sci
scilab/modules/signal_processing/macros/group.sci
scilab/modules/signal_processing/macros/hank.sci
scilab/modules/signal_processing/macros/hilb.sci
scilab/modules/signal_processing/macros/hilbert.sci
scilab/modules/signal_processing/macros/iir.sci
scilab/modules/signal_processing/macros/iirgroup.sci
scilab/modules/signal_processing/macros/iirlp.sci
scilab/modules/signal_processing/macros/iirmod.sci
scilab/modules/signal_processing/macros/intdec.sci
scilab/modules/signal_processing/macros/jmat.sci
scilab/modules/signal_processing/macros/kalm.sci
scilab/modules/signal_processing/macros/lattn.sci
scilab/modules/signal_processing/macros/lev.sci
scilab/modules/signal_processing/macros/levin.sci
scilab/modules/signal_processing/macros/lindquist.sci
scilab/modules/signal_processing/macros/mese.sci
scilab/modules/signal_processing/macros/mfft.sci
scilab/modules/signal_processing/macros/mrfit.sci
scilab/modules/signal_processing/macros/phc.sci
scilab/modules/signal_processing/macros/pspect.sci
scilab/modules/signal_processing/macros/remezb.sci
scilab/modules/signal_processing/macros/sincd.sci
scilab/modules/signal_processing/macros/srfaur.sci
scilab/modules/signal_processing/macros/srkf.sci
scilab/modules/signal_processing/macros/sskf.sci
scilab/modules/signal_processing/macros/system.sci
scilab/modules/signal_processing/macros/trans.sci
scilab/modules/signal_processing/macros/wfir.sci
scilab/modules/signal_processing/macros/wiener.sci
scilab/modules/signal_processing/macros/wigner.sci
scilab/modules/signal_processing/macros/window.sci
scilab/modules/signal_processing/macros/yulewalk.sci
scilab/modules/signal_processing/macros/zpbutt.sci
scilab/modules/signal_processing/macros/zpch1.sci
scilab/modules/signal_processing/macros/zpch2.sci
scilab/modules/signal_processing/macros/zpell.sci
scilab/modules/signal_processing/sci_gateway/fortran/Ex-corr.f
scilab/modules/signal_processing/sci_gateway/fortran/bgetx.f
scilab/modules/signal_processing/sci_gateway/fortran/bgety.f
scilab/modules/signal_processing/sci_gateway/fortran/sci_amell.f
scilab/modules/signal_processing/sci_gateway/fortran/sci_corr.f
scilab/modules/signal_processing/sci_gateway/fortran/sci_fft.f
scilab/modules/signal_processing/sci_gateway/fortran/sci_fiir.f
scilab/modules/signal_processing/sci_gateway/fortran/sci_rpem.f
scilab/modules/signal_processing/sci_gateway/fortran/sci_syredi.f
scilab/modules/signal_processing/src/fortran/amell.f
scilab/modules/signal_processing/src/fortran/belan.f
scilab/modules/signal_processing/src/fortran/desi24.f
scilab/modules/signal_processing/src/fortran/dspln.f
scilab/modules/signal_processing/src/fortran/nstabl.f
scilab/modules/signal_processing/src/fortran/sci_ffir.f
scilab/modules/signal_processing/src/fortran/trbipo.f

index e37ca96..f396749 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [y]=%asn(x,m)
 //Calculates the elliptic integral:
 //  y = integral from 0 to x of
@@ -8,14 +17,6 @@ function [y]=%asn(x,m)
 //  y :Value of the integral
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
   m=real(m);
 
   if m<0|m>1 then 
index f111204..b8c0b4a 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [K]=%k(m)
 //K=%k(m)
 //Calculates Jacobi's complete elliptic integral
@@ -12,14 +21,6 @@ function [K]=%k(m)
 //    :on the real axis.
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
    [n1,n2]=size(m);
    un=ones(n1,n2);
index 719fb66..4356769 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [y]=%sn(x,m)
 //Jacobi 's elliptic function with parameter m
 //which computes the inverse of the elliptic
@@ -11,14 +20,6 @@ function [y]=%sn(x,m)
 //  y :Result
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 [n1,n2]=size(x);
 n=n1*n2;
index e0f574c..4b5141a 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [hs,pols,zers,gain]=analpf(n,fdesign,rp,omega)
 //[hs,pols,zers,gain]=analpf(n,fdesign,rp,omega)
 //Creates analog low-pass filter with cut-off frequency at omega
@@ -22,14 +31,7 @@ function [hs,pols,zers,gain]=analpf(n,fdesign,rp,omega)
 //     hs=gain*poly(zers,'s')/poly(pols,'s')
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 
    select fdesign
    case 'butt' then
index 4da1fd6..a38caf2 100644 (file)
@@ -1,3 +1,13 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1989 - C. Bunks
+// Copyright (C) INRIA - 1997 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [npl,nzr,ngn]=bilt(pl,zr,gn,num,den)
 //[npl,nzr,ngn]=bilt(pl,zr,gn,num,den)
 //macro for calculating the gain poles and zeros
@@ -11,15 +21,6 @@ function [npl,nzr,ngn]=bilt(pl,zr,gn,num,den)
 //  num :numerator of transform
 //  den :denominator of transform
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1989 - C. Bunks
-// Copyright (C) INRIA - 1997 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
   if type(pl)<>1 then
     error(msprintf(gettext("%s: Wrong type for input argument #%d: Array of floating point numbers expected.\n") ,"bilt",1))
   end
index d528416..9a0a8b3 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [h]=buttmag(order,omegac,sample)
 //<h>=buttmag(order,omegac,sample)
 //Squared Magnitude response of a Butterworth filter
@@ -8,14 +17,6 @@ function [h]=buttmag(order,omegac,sample)
 //  h      :Butterworth filter values at sample points
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
    [n1,n2]=size(sample);
    un=ones(n1,n2);
index 461512d..15e97e6 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - F.D
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function cels=casc(x,z)
 //cels=casc(x,z)
 //Creates cascade realization of filter
@@ -24,14 +33,6 @@ function cels=casc(x,z)
 //  !              2               2               2 !
 //  !   7 + 10z + z     8 + 11z + z     9 + 12z + z  !
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - F.D
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 cels=[];
 for col=x,
index f4a861e..10625c2 100644 (file)
@@ -1,12 +1,5 @@
-function fresp = cepstrum(w,mag)
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 // Copyright (C) INRIA
-// 
-// Uses the complex-cepstrum (Oppenheim & Schafer, Digital
-//Signal Processing, p. 501) to generate, at the frequencies 
-//w, a complex frequency response fresp whose magnitude is 
-//equal to magnitude data mag and whose phase corresponds
-//to a stable, minimum phase transfer function. 
 //
 // This file must be used under the terms of the CeCILL.
 // This source file is licensed as described in the file COPYING, which
@@ -14,6 +7,13 @@ function fresp = cepstrum(w,mag)
 // are also available at    
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
+function fresp = cepstrum(w,mag)
+// Uses the complex-cepstrum (Oppenheim & Schafer, Digital
+//Signal Processing, p. 501) to generate, at the frequencies 
+//w, a complex frequency response fresp whose magnitude is 
+//equal to magnitude data mag and whose phase corresponds
+//to a stable, minimum phase transfer function. 
+
   if ~isreal(w) then
      error(msprintf(gettext("%s: Input argument #%d must be real.\n"),"cepstrum",1)); 
   end
index ad83b1f..93cf2b8 100644 (file)
@@ -1,3 +1,11 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 function [h2]=cheb1mag(n,omegac,epsilon,sample)
 //<h2>=cheb1mag(n,omegac,epsilon,sample)
 //Square magnitude response of a type 1 Chebyshev filter
@@ -12,14 +20,7 @@ function [h2]=cheb1mag(n,omegac,epsilon,sample)
 //  h2      :Chebyshev I filter values at sample points
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 
    [n1,n2]=size(sample);
    un=ones(n1,n2);
index d7eb782..bcfc2cf 100644 (file)
@@ -1,3 +1,11 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 function [h2]=cheb2mag(n,omegar,A,sample)
 //<h2>=cheb2mag(n,omegar,A,sample)
 //Square magnitude response of a type 1 Chebyshev filter
@@ -11,14 +19,7 @@ function [h2]=cheb2mag(n,omegar,A,sample)
 //  h2      :Chebyshev II filter values at sample points
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 
    [n1,n2]=size(sample);
    un=ones(n1,n2);
index fb86d83..be92be3 100644 (file)
@@ -1,10 +1,3 @@
-function Tn=chepol(n,var)
-//Chebychev polynomial
-//  n   :Polynomial order
-//  var :Polynomial variable (character string)
-//  Tn  :Polynomial in var
-//
-//!
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 // Copyright (C) INRIA - F.D
 // 
@@ -13,6 +6,14 @@ function Tn=chepol(n,var)
 // you should have received as part of this distribution.  The terms
 // are also available at    
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+function Tn=chepol(n,var)
+//Chebychev polynomial
+//  n   :Polynomial order
+//  var :Polynomial variable (character string)
+//  Tn  :Polynomial in var
+//
+//!
 if n==0 then
   Tn=poly(1,var,'coeff'),
 elseif n==1 then
index 42c1a8d..5dbf4c4 100644 (file)
@@ -1,3 +1,11 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 function [y,e1]=convol(h,x,e0)
 //  convol - convolution 
 //%CALLING SEQUENCE
@@ -18,14 +26,7 @@ function [y,e1]=convol(h,x,e0)
 //  Final call : [yN]=convol(h,xN,eNm1)
 //  Finally y=[y1,y2,...,yNm1,yN]
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 [lhs,rhs]=argn(0)
 n=prod(size(x))
 m=prod(size(h))
index c4efcbc..8b88fae 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [sm,cwp]=cspect(nlags,ntp,wtype,x,y,wpar)
 //<sm,cwp>=cspect(nlags,ntp,wtype,x,y,wpar)
 //Spectral estimation using the modified periodogram method.
@@ -17,14 +26,6 @@ function [sm,cwp]=cspect(nlags,ntp,wtype,x,y,wpar)
 //           :window parameter
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
   [lhs,rhs]=argn(0);
   cross=0;
index cd73685..5a86900 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [czx]=czt(x,m,w,phi,a,theta)
 //czx=czt(x,m,w,phi,a,theta)
 //chirp z-transform algorithm which calcultes the
@@ -11,14 +20,6 @@ function [czx]=czt(x,m,w,phi,a,theta)
 //  theta :Initial phase
 //  czx   :Chirp z-transform output
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 //get the size of x and find the maximum of (n,m)
  
index 196188b..4b3a889 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) Bruno Pincon
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+//
 function [y] = detrend(x, flag, bp)
    //
    // this function removes the constant or linear or
@@ -17,9 +26,7 @@ function [y] = detrend(x, flag, bp)
    // (m = length(x) if x is a vector and m = size(x,1) in case
    // x is a matrix). So bp must be reals in [0 m-1].
    //
-   // AUTHOR
-   //   Bruno Pincon
-   //
+
    rhs = argn(2)
    if rhs < 1 | rhs > 3 then
       error(msprintf(gettext("%s: Wrong number of input arguments: %d to %d expected.\n"),'detrend',1,3));
index b2ec264..efb2b1b 100644 (file)
@@ -1,10 +1,3 @@
-function xf=dft(x,flag);
-//xf=dft(x,flag);
-//macro which computes dft of vector x
-// x    :input vector
-// flag :indicates dft or idft
-// xf   :output vector
-//!
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 // Copyright (C) INRIA - 1988 - C. Bunks
 // 
@@ -14,6 +7,14 @@ function xf=dft(x,flag);
 // are also available at    
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
+function xf=dft(x,flag);
+//xf=dft(x,flag);
+//macro which computes dft of vector x
+// x    :input vector
+// flag :indicates dft or idft
+// xf   :output vector
+//!
+
 n=maxi(size(x));
 arg=(0:n-1);
 am=-2*%pi*%i*arg'*arg/n;
index e76dd79..7fc0177 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [v]=ell1mag(eps,m1,z)
 //Function used for squared magnitude of an elliptic filter
 //Usually m1=eps*eps/(a*a-1);
@@ -7,14 +16,6 @@ function [v]=ell1mag(eps,m1,z)
 //  v       :Elliptic filter values at sample points
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 s=%sn(z,m1);un=ones(z);
 v=real(un./(un+eps*eps*s.*s))
 
index 18a33bd..f585fd7 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [hn]=eqfir(nf,bedge,des,wate)
 //<hn>=eqfir(nf,bedge,des,wate)
 //Minimax approximation of multi-band, linear phase, FIR filter
@@ -7,14 +16,6 @@ function [hn]=eqfir(nf,bedge,des,wate)
 //  wate  :M-vector giving relative weight of error in each band
 //  hn    :Output of linear-phase FIR filter coefficients
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 //get number of cosines
  
index f477b66..f1eb72d 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [cells,fact,zzeros,zpoles]=eqiir(ftype,approx,om,deltap,deltas)
 //[cells,fact,zzeros,zpoles]=eqiir(ftype,approx,om,deltap,deltas)
 //Design of iir filter :interface with eqiir (syredi)
@@ -21,14 +30,6 @@ function [cells,fact,zzeros,zpoles]=eqiir(ftype,approx,om,deltap,deltas)
 //     hz=fact*prod(cells(2))./prod(cells(3))
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
   select part(approx,1);
   case 'b'
     iapro=1
index 489c529..26eb0c3 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1989 - G. Le Vey
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [P,R,T]=faurre(n,H,F,G,R0)
 //[P,R,T]=faurre(n,H,F,G,R0)
 //macro which computes iteratively the minimal solution of the algebraic
@@ -8,14 +17,6 @@ function [P,R,T]=faurre(n,H,F,G,R0)
 //   P    : solution of the Riccati equation after n iterations.
 //   R,T  : gain matrix of the filter.
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1989 - G. Le Vey
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 //initialization
 Pn=G*inv(R0)*G'
index da6712c..9b0409d 100644 (file)
@@ -1,4 +1,4 @@
-function x=fft2(varargin)
+
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 // Copyright (C) INRIA - F.B
 // 
@@ -8,6 +8,7 @@ function x=fft2(varargin)
 // are also available at    
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 //
+function x=fft2(varargin)
 // Two-dimensional fast Fourier transform
 // Syntax : y = fft2(x) or y = fft2(x,m,n)
 // Inputs :
index 2235af4..793985c 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function x = fftshift(x,job)
   if argn(2)<2 then job='all',end
   deff('sel=fun(sk)','c=ceil(sk/2);sel=[c+1:sk,1:c]')
index c750870..3a4f0cf 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [m]=find_freq(epsilon,A,n)
 //Search for m such that n=K(1-m1)K(m)/(K(m1)K(1-m))
 //with m1=(epsilon*epsilon)/(A*A-1);
@@ -11,14 +20,6 @@ function [m]=find_freq(epsilon,A,n)
 //          :elliptic filter
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
    m1=(epsilon*epsilon)/(A*A-1);
    chi1=%k(1-m1)/%k(m1);
index 19a9670..feb4aa6 100644 (file)
@@ -1,6 +1,3 @@
-function m=findm(chi)
-//Search for m such that chi = %k(1-m)/%k(m)
-//!
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 // Copyright (C) INRIA
 // 
@@ -10,6 +7,10 @@ function m=findm(chi)
 // are also available at    
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
+function m=findm(chi)
+//Search for m such that chi = %k(1-m)/%k(m)
+//!
+
 if chi < 1 then
   t=1;
   tn=2;
index 78203c7..0464f58 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [num,den]=frfit(w,fresp,order,weight)
 //Calling sequence:
 //[num,den]=frfit(w,fresp,r,weight)
@@ -14,14 +23,6 @@ function [num,den]=frfit(w,fresp,order,weight)
 //  freq(num,den,%i*w) should be close to fresp
 //
 // changing frequencies to rad/s 
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 w=2*%pi*w;
 [LHS,RHS]=argn(0);
 if RHS==3 
index bfe60a2..4f06eec 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [xm,fr]=frmag(num,den,npts)
 //[xm,fr]=frmag(num[,den],npts)
 //Calculates the magnitude of the frequency respones of
@@ -16,14 +25,7 @@ function [xm,fr]=frmag(num,den,npts)
 //  fr   :Points in the frequency domain where
 //       :magnitude is evaluated
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
   select argn(2)
   case 2 then //frmag(sys,npts)
     npts=den;
index c72844e..3c6c44b 100644 (file)
@@ -1,3 +1,11 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1989 - G. Le Vey
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 function [hst]=fsfirlin(hd,flag)
 //<hst>=fsfirlin(hd,flag)
 //macro for the design of FIR, linear phase filters
@@ -8,14 +16,7 @@ function [hst]=fsfirlin(hd,flag)
 //  hst  : vector giving the approximated continuous response
 //         on a dense grid of frequencies
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1989 - G. Le Vey
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 
    n1=prod(size(hd));//filter length
    if int(n1/2)==n1/2,n=2*n1;else,n=2*n1+1;end;//even or odd length
index 0ea7250..290de8e 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [tg,fr]=group(npts,a1i,a2i,b1i,b2i)
 //Calculate the group delay of a digital filter
 //with transfer function h(z).
@@ -21,14 +30,7 @@ function [tg,fr]=group(npts,a1i,a2i,b1i,b2i)
 //       h(z)=prod(a1i+a2i*z+z**2)/prod(b1i+b2i*z+z^2)
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 
 //get frequency grid values in [0,.5)
  
index 354a754..3b1f72c 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1989 - G. Le Vey 
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [hk]=hank(m,n,cov)
 //<hk>=hank(m,n,cov)
 //this macro builds the hankel matrix of size (m*d,n*d)
@@ -7,14 +16,6 @@ function [hk]=hank(m,n,cov)
 //   cov: sequence of covariances; it must be given as :[R0 R1 R2...Rk]
 //   hk : computed hankel matrix
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1989 - G. Le Vey 
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
    hk=[];
    d=mini(size(cov));
index 1ad2b7c..8ec1efb 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function xh=hilb(n,wtype,par)
 //xh=hilb(n [,wtype [,par]])
 //returns the first n points of the
@@ -13,11 +22,6 @@ function xh=hilb(n,wtype,par)
 // xh    :Hilbert transform
 //
 //!
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 // References:
 // http://ieeexplore.ieee.org/iel4/78/7823/00330385.pdf?tp=&arnumber=330385&isnumber=7823
 //  A. Reilly, G. Frazer, and B. Boashash, "Analytic signal generation
index 8ad36c0..4149cc3 100644 (file)
@@ -1,3 +1,14 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2008 - INRIA
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+//
+//
+
 function x = hilbert(xr)
 // Marple, S.L., "Computing the discrete-time analytic signal via FFT,"
 // IEEE Transactions on Signal Processing, Vol. 47, No.9 (September
index d490f4c..e17dcc0 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function hz=iir(n,ftype,fdesign,frq,delta)
 // hz=iir(n,ftype,fdesign,frq,delta)
 //macro which designs an iir digital filter
@@ -25,14 +34,6 @@ function hz=iir(n,ftype,fdesign,frq,delta)
 //           :                    0<ripple<delta(2)   in stopband
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 //select analog filter design for low-pass filter with fc=.25
   if type(n)<>1 then
index af055d2..c47e0a7 100644 (file)
@@ -1,10 +1,3 @@
-function [lt,grad]=iirgroup(p,r,theta,omega,wt,td)
-//    p===> Lp
-//    r=magnitude of poles and zeros of filters
-//theta=phase   "    "   "    "    "
-//omega=frequencies corresponding to filters specs.
-//wt=weighting function for group delay
-//!
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 // Copyright (C) INRIA
 // 
@@ -13,6 +6,14 @@ function [lt,grad]=iirgroup(p,r,theta,omega,wt,td)
 // you should have received as part of this distribution.  The terms
 // are also available at    
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+function [lt,grad]=iirgroup(p,r,theta,omega,wt,td)
+//    p===> Lp
+//    r=magnitude of poles and zeros of filters
+//theta=phase   "    "   "    "    "
+//omega=frequencies corresponding to filters specs.
+//wt=weighting function for group delay
+//!
 [m,n]=size(td);if m>n,td=td';end;
 [m,n]=size(omega);if m>n,omega=omega';end;
 [m,n]=size(r);if n>m,r=r';m=n;end;
index d112b03..0baff0b 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [cout,grad,ind]=iirlp(x,ind,p,flag,lambda,omega,ad,wa,td,wt)
 //
 //optimization of IIR filters IIR with  Lp criterium for magnitude
@@ -14,14 +23,6 @@ function [cout,grad,ind]=iirlp(x,ind,p,flag,lambda,omega,ad,wa,td,wt)
 //wa,wt=fonctions de ponderation pour l'amplitude et le
 //retard de groupe ad,td=amplitudes et retard de groupe desires
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 r=x(:,1);theta=x(:,2);
 [m,n]=size(ad);if m>n,ad=ad';end
 [m,n]=size(td);if m>n,td=td';end
index da0bed0..a01c220 100644 (file)
@@ -1,10 +1,3 @@
-function [la,grad,ind]=iirmod(x,ind)
-//    p===> critere Lp
-//    r=module des poles et zeros des filtres
-//theta=argument des  "    "   "    "    "
-//omega=frequences ou sont donnees les specifications des filtres
-//wa=fonction de ponderation pour l'amplitude
-//!
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 // Copyright (C) INRIA
 // 
@@ -13,6 +6,14 @@ function [la,grad,ind]=iirmod(x,ind)
 // you should have received as part of this distribution.  The terms
 // are also available at    
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+function [la,grad,ind]=iirmod(x,ind)
+//    p===> critere Lp
+//    r=module des poles et zeros des filtres
+//theta=argument des  "    "   "    "    "
+//omega=frequences ou sont donnees les specifications des filtres
+//wa=fonction de ponderation pour l'amplitude
+//!
 r=x(:,1);theta=x(:,2);
 [m,n]=size(ad);if m>n,ad=ad';end;
 [m,n]=size(omega);if m>n,omega=omega';end;
index ceb162c..0966d27 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1990 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [y]=intdec(x,lom)
 //y=intdec(x,lom)
 //Changes the sampling rate of a 1D or 2D signal by the rates in lom
@@ -7,14 +16,6 @@ function [y]=intdec(x,lom)
 //         :changes lom=(col rate change,row rate change)
 //  y      :Output sampled signal
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1990 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 //Get dimensions of vectors
  
index 26c4766..9a442be 100644 (file)
@@ -1,10 +1,3 @@
-function [j]=jmat(n,m)
-//j=jmat(n,m)
-//This macro permutes block rows or block columns of a matrix
-//
-//   n   : number of block rows or block columns of the matrix
-//   m   : size of the (square) blocks
-//!
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 // Copyright (C) INRIA
 // 
@@ -13,6 +6,14 @@ function [j]=jmat(n,m)
 // you should have received as part of this distribution.  The terms
 // are also available at    
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+function [j]=jmat(n,m)
+//j=jmat(n,m)
+//This macro permutes block rows or block columns of a matrix
+//
+//   n   : number of block rows or block columns of the matrix
+//   m   : size of the (square) blocks
+//!
 j=[];
 for k=1:n,j(k,n-k+1)=1;end;
 j=j.*.eye(m,m);
index 80dc128..d15c954 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [x1,p1,x,p]=kalm(y,x0,p0,f,g,h,q,r)
 //[x1,p1,x,p]=kalm(y,x0,p0,f,g,h,q,r)
 //macro which gives the Kalman update and error variance
@@ -14,14 +23,6 @@ function [x1,p1,x,p]=kalm(y,x0,p0,f,g,h,q,r)
 // x,p    :updated estimate and error covariance at t=0
 //        :based on data up to t=0
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
    k=p0*h'*(h*p0*h'+r)**(-1);
    p=(eye(p0)-k*h)*p0;
index 1ac1bdb..b4f191e 100644 (file)
@@ -1,3 +1,11 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1989 - G. Le Vey
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 function [la,lb]=lattn(n,p,cov)
 //[la,lb]=lattn(n,p,cov)
 //macro which solves recursively on n (p being fixed)
@@ -33,14 +41,7 @@ function [la,lb]=lattn(n,p,cov)
 //   la  : list-type variable, giving the successively calculated
 //       : polynomials (degree 1 to degree n),with coefficients Ak
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1989 - G. Le Vey
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 
    if argn(2)<>3 then 
      error(msprintf(gettext("%s: Wrong number of input argument(s): %d expected.\n"),'lattn',3)) 
index c7393e3..f8eb1aa 100644 (file)
@@ -1,3 +1,13 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// Copyright (C) INRIA - 1991 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [ar,sigma2,rc]=lev(r)
 //[ar,sigma2,rc]=lev(r)
 //Resolve the Yule-Walker equations:
@@ -14,15 +24,6 @@ function [ar,sigma2,rc]=lev(r)
 //  sigma2 :Scale constant
 //  rc     :Reflection coefficients
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// Copyright (C) INRIA - 1991 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 //get the size of the correlation vector
  
index e5db688..4b16535 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - G. Le Vey
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [la,sig,lb]=levin(n,cov)
 //[la,sig,lb]=levin(n,cov)
 //macro which solves recursively on n
@@ -33,14 +42,6 @@ function [la,sig,lb]=levin(n,cov)
 //  sig : list-type variable, giving the successive
 //      : mean-square errors.
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - G. Le Vey
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
    [lhs,rhs]=argn(0);
    if rhs<>2 then
index ba00a67..31edbbf 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1989 - G. Le Vey
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [P,R,T]=lindquist(n,H,F,G,R0)
 //[Pn,Rn,Tn]=lindquist(n,H,F,G,R0)
 //macro which computes iteratively the minimal solution of the algebraic
@@ -9,14 +18,6 @@ function [P,R,T]=lindquist(n,H,F,G,R0)
 //   Pn    : solution of the Riccati equation after n iterations.
 //   Rn,Tn : gain matrices of the filter.
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1989 - G. Le Vey
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 [d,m]=size(H);
 //initialization
 Gn=G;
index 76cc534..62f73f7 100644 (file)
@@ -1,12 +1,3 @@
-function [sm,fr]=mese(x,npts);
-//<sm,fr]=mese(x [,npts]);
-//Calculate the maximum entropy spectral estimate of x
-//  x    :Input sampled data sequence
-//  npts :Optional parameter giving number of points of fr and sm
-//        (default is 256)
-//  sm   :Samples of spectral estimate on the frequency grid fr
-//  fr   :npts equally spaced frequency samples in [0,.5)
-//!
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 // Copyright (C) INRIA - 1988 - C. Bunks
 // Copyright (C) INRIA - 1991 - C. Bunks
@@ -17,6 +8,16 @@ function [sm,fr]=mese(x,npts);
 // are also available at    
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
+function [sm,fr]=mese(x,npts);
+//<sm,fr]=mese(x [,npts]);
+//Calculate the maximum entropy spectral estimate of x
+//  x    :Input sampled data sequence
+//  npts :Optional parameter giving number of points of fr and sm
+//        (default is 256)
+//  sm   :Samples of spectral estimate on the frequency grid fr
+//  fr   :npts equally spaced frequency samples in [0,.5)
+//!
+
 //default evaluation
  
    [rhs,lhs]=argn(0);
index 92361f2..3daa1de 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [xk]=mfft(x,flag,dim)
 //<xk>=mfft(x,flag,dim)
 //macro which calculates the fft for a multi-dimensional signal
@@ -26,14 +35,6 @@ function [xk]=mfft(x,flag,dim)
 //      dim=<3,2,3>
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
    xk=x;
    dims=[1 matrix(dim,1,maxi(size(dim)))];
index 79c0b37..6433f02 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [num,den]=mrfit(w,mod,r)
 //Calling sequence:
 //[num,den]=mrfit(w,mod,r)
@@ -13,14 +22,6 @@ function [num,den]=mrfit(w,mod,r)
 //
 //  abs(freq(num,den,%i*w)) should be close to mod
 //
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 w=w(:);mod=mod(:);
 [LHS,RHS]=argn(0);
 if w(1)==0 then w(1)=%eps;end
index b03f0f2..e95074e 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1989 - G. Le Vey
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [h,f,g]=phc(hk,d,r)
 //[h,f,g]=phc(hk,d,r)
 //macro which computes the matrices h, f, g by the principal hankel
@@ -9,14 +18,7 @@ function [h,f,g]=phc(hk,d,r)
 //       : for the approximated model
 //   h,f,g  : relevant matrices of the state-space model
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1989 - G. Le Vey
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 
    [p,q]=size(hk);
    [u,s,v]=svd(hk);
index 49f620f..5358e4c 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [sm,cwp]=pspect(sec_step,sec_leng,wtype,x,y,wpar)
 //<sm,cwp>=pspect(sec_step,sec_leng,wtype,x,y,wpar)
 //Cross-spectral estimate between x and y if both are given
@@ -14,14 +23,6 @@ function [sm,cwp]=pspect(sec_step,sec_leng,wtype,x,y,wpar)
 // sm       :power spectral estimate in the interval [0,1]
 // cwp      :unspecified Chebyshev window parameter
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
   [lhs,rhs]=argn(0);
   cross=0;
index ff93fc1..1863127 100644 (file)
@@ -1,3 +1,11 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 function an=remezb(nc,fg,ds,wt)
 //an=remezb(nc,fg,ds,wt)
 //Minimax approximation of a frequency domain
@@ -20,14 +28,7 @@ function an=remezb(nc,fg,ds,wt)
 //  wt :Weighting function on error on grid fg
 //  an :Cosine filter coefficients
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 
 //get frequency grid size
   
index 4b88ec4..7211f85 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1989 - G. Le Vey
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [s]=sincd(n,flag)
 //<s>=sincd(n,flag)
 //macro which calculates the function Sin(N*x)/Sin(x)
@@ -7,14 +16,6 @@ function [s]=sincd(n,flag)
 //  s    : vector of values of the function on a dense
 //       : grid of frequencies
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1989 - G. Le Vey
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
    npt=4*n;
    pas=%pi/npt;
index 0e07e7d..14cb056 100644 (file)
@@ -1,3 +1,11 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 function [p,s,t,l,rt,tt]=srfaur(h,f,g,r0,n,p,s,t,l)
 //square-root algorithm for the algebraic Riccati equation.
 //
@@ -13,14 +21,7 @@ function [p,s,t,l,rt,tt]=srfaur(h,f,g,r0,n,p,s,t,l)
 //          : p,s,t,l may be given as input if more than one recursion
 //          : is desired (evaluation of intermediate values of p, e.g.).
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 [lhs,rhs]=argn(0);
 [d,m]=size(h);
 if rhs==5,
index d4f8665..297d2bd 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [x1,p1]=srkf(y,x0,p0,f,h,q,r)
 //square root kalman filter algorithm
 //Input to the macro is:
@@ -12,14 +21,6 @@ function [x1,p1]=srkf(y,x0,p0,f,h,q,r)
 // x1,p1  :updated estimate and error covariance
 //        :at t=1 based on data up to t=0
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
    n=maxi(size(x0));
    p=maxi(size(y));
index 26b8dcf..75a371a 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [xe]=sskf(y,f,h,q,r,x0)
 //<xe>=sskf(y,f,h,q,r,x0)
 //steady-state kalman filter
@@ -9,14 +18,6 @@ function [xe]=sskf(y,f,h,q,r,x0)
 //
 // xe  :estimated state
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 //get steady-state Kalman gain
  
index 929f4e3..33bf249 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 1988 - INRIA - Bunks
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+//
 function [x1,y]=system(x0,f,g,h,q,r)
 //<x1,y>=system(x0,f,g,h,q,r)
 //define system macro which generates the next
index 8189b43..a20eb2b 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [hzt]=trans(pd,zd,gd,tr_type,frq)
 //hzt=trans(pd,zd,gd,tr_type,frq)
 //macro for transforming standardized low-pass filter into
@@ -8,14 +17,6 @@ function [hzt]=trans(pd,zd,gd,tr_type,frq)
 // frq     :frequency values
 // hzt     :output polynomial
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
   select argn(2)
   case 3 then //trans(hz,tr_type,frq): filter given by a siso tranfer function
     hz=pd
index a0dea04..d1fef74 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [wft,wfm,fr]=wfir(ftype,forder,cfreq,wtype,fpar)
 //[wft,wfm,fr]=wfir(ftype,forder,cfreq,wtype,fpar)
 //Macro which makes linear-phase, FIR low-pass, band-pass,
@@ -17,14 +26,6 @@ function [wft,wfm,fr]=wfir(ftype,forder,cfreq,wtype,fpar)
 //  wfm    :Frequency domain filter response on the grid fr
 //  fr     :Frequency grid
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 wft=[];wfm=[];fr=[]
 //check arguments of macro call
  
index 04ba7da..4ca2660 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [xs,ps,xf,pf]=wiener(y,x0,p0,f,g,h,q,r)
 //<xs,ps,xf,pf>=wiener(y,x0,p0,f,g,h,q,r)
 //macro which gives the Wiener estimate using
@@ -22,14 +31,6 @@ function [xs,ps,xf,pf]=wiener(y,x0,p0,f,g,h,q,r)
 //form of xs and xf:  [x0,x1,...,xf], and xk is a column n-vector
 //form of ps and pf:  [p0,p1,...,pf], and pk is a nxn matrix
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 //obtain the dimensions of xk and yk.
 //Get the time interval [t0,tf].
index f364cc2..699924e 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function tab=wigner(x,h,deltat,zp)
 //macro which computes the 'time-frequency' wigner
 //spectrum of a signal.
@@ -11,14 +20,6 @@ function tab=wigner(x,h,deltat,zp)
 //
 //   Initializations
 //
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 l=prod(size(x));
 n=prod(size(h));
 npr=2*n;
index 2ce5bc4..21e4d71 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1988 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [win_l,cwp]=window(wtype,n,par)
 //[win_l,cwp]=window(wtype,n,par)
 //macro which calculates symmetric window
@@ -10,14 +19,6 @@ function [win_l,cwp]=window(wtype,n,par)
 // win   :window
 // cwp   :unspecified Chebyshev window parameter
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1988 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
   [lhs,rhs]=argn(0);
   cwp=-1;
index 9bbf2bf..20b916f 100644 (file)
@@ -1,3 +1,12 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [Nz,Dz]=yulewalk(Norder, frq, mag)
 //YULEWALK  filter design using a least-squares method.
 //     Hz = YULEWALK(N,frq,mag) finds the N-th order iir filter
@@ -22,14 +31,6 @@ function [Nz,Dz]=yulewalk(Norder, frq, mag)
 //xbasc(1);xset('window',1);plot2d(fs*frq',abs(repf'));
 //xtitle('Obtained Frequency Response')
 //
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
   [LHS,RHS]=argn(0);
   if RHS <>3
     error(msprintf(gettext("%s: Wrong number of input argument(s): %d expected.\n"),'yulewalk',3));
index 7b9a931..8a9fe83 100644 (file)
@@ -1,3 +1,13 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// Copyright (C) INRIA - 1996 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [pols,gain]=zpbutt(n,omegac)
 //<pols,gain>=zpbutt(n,omegac)
 //Computes the poles of a Butterworth analog
@@ -10,15 +20,6 @@ function [pols,gain]=zpbutt(n,omegac)
 //  gain   :Resulting gain of filter
 //
 //! 
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// Copyright (C) INRIA - 1996 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 angles=ones(1,n)*(%pi/2+%pi/(2*n))+(0:n-1)*%pi/n;
 pols=omegac*exp(%i*angles);
 gain=abs((-omegac)^n);
index 9e832ea..95227c0 100644 (file)
@@ -1,3 +1,13 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// Copyright (C) INRIA - 1996 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [pols,gain]=zpch1(n,epsilon,omegac)
 //Poles of a Type 1 Chebyshev analog filter
 //The transfer function is given by :
@@ -9,15 +19,6 @@ function [pols,gain]=zpch1(n,epsilon,omegac)
 //  gain    :Resulting filter gain
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// Copyright (C) INRIA - 1996 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 Gamma=((1+sqrt(1+epsilon**2))/epsilon)^(1/n);
 a=omegac*(Gamma-1/Gamma)/2;
 b=omegac*(Gamma+1/Gamma)/2;
index 7e22540..d43f10c 100644 (file)
@@ -1,3 +1,13 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - F.D
+// Copyright (C) INRIA - 1996 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [zers,pols,gain]=zpch2(n,A,omegar)
 //[zers,pols,gain]=zpch2(n,A,omegar)
 //Poles and zeros of a type 2 Chebyshev analog filter
@@ -11,15 +21,6 @@ function [zers,pols,gain]=zpch2(n,A,omegar)
 //  gain   :Resulting filter gain
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - F.D
-// Copyright (C) INRIA - 1996 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 un=ones(1,n);
 v=%pi/(2*n)*(1:2:2*n-1);
 w=exp(%i*v);
index 567d444..418c849 100644 (file)
@@ -1,3 +1,13 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA - 1989 - F.Delebecque
+// Copyright (C) INRIA - 1996 - C. Bunks
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function [ze,po,gain]=zpell(epsilon,A,omegac,omegar)
 //[ze,po,gain]=zpell(epsilon,A,omegac,omegar)
 //Poles and zeros of prototype lowpass elliptic filter
@@ -11,15 +21,7 @@ function [ze,po,gain]=zpell(epsilon,A,omegac,omegar)
 //  gain    :Resulting gain of filter
 //
 //!
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) INRIA - 1989 - F.Delebecque
-// Copyright (C) INRIA - 1996 - C. Bunks
-// 
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at    
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 m1=(epsilon*epsilon)/(A*A-1);
    K1=%asn(1,m1);
    K1t=imag(%asn(1/sqrt(m1),m1));
index c1ea02d..406b222 100644 (file)
@@ -1,5 +1,3 @@
-c     Exemple for corr
-c     ----------------
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
index fa3970b..86b345a 100644 (file)
@@ -1,8 +1,3 @@
-      subroutine bgetx(x, incr, istart)
-c ======================================================================
-c     macros or list externals for corr 
-c ======================================================================
-
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
@@ -12,6 +7,12 @@ c you should have received as part of this distribution.  The terms
 c are also available at    
 c http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
+      subroutine bgetx(x, incr, istart)
+c ======================================================================
+c     macros or list externals for corr 
+c ======================================================================
+
+
 
       INCLUDE 'stack.h'
       integer iadr,sadr
index c6e6407..9d58eb8 100644 (file)
@@ -1,8 +1,3 @@
-      subroutine bgety(y, incr, istart)
-c ======================================================================
-c     macros or list externals for corr 
-c ======================================================================
-
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
@@ -12,6 +7,12 @@ c you should have received as part of this distribution.  The terms
 c are also available at    
 c http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
+      subroutine bgety(y, incr, istart)
+c ======================================================================
+c     macros or list externals for corr 
+c ======================================================================
+
+
       INCLUDE 'stack.h'
       integer iadr,sadr
       double precision y(*)
index 380a76f..2af0e7a 100644 (file)
@@ -1,4 +1,3 @@
-c     =======================================
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
@@ -60,4 +59,4 @@ c
       return
       end
 c     =======================================
-      
\ No newline at end of file
+      
index 392b8ad..2ba8791 100644 (file)
@@ -1,4 +1,3 @@
-c     =======================================
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
index 4c26fe5..fe3f241 100644 (file)
@@ -1,4 +1,3 @@
-c     =======================================
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
index 70067dd..44144a9 100644 (file)
@@ -1,4 +1,3 @@
-c     =======================================
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
index 6dae034..69f5370 100644 (file)
@@ -1,4 +1,3 @@
-c     =======================================
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
@@ -146,4 +145,4 @@ c
       return
       end
 c     =======================================
-      
\ No newline at end of file
+      
index f953dd5..ff4906e 100644 (file)
@@ -1,4 +1,3 @@
-c     =======================================
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
@@ -267,4 +266,4 @@ c      endif
       return
       end
 c     =======================================
-      
\ No newline at end of file
+      
index 9d85419..5a7da6a 100644 (file)
@@ -1,3 +1,12 @@
+c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+c Copyright (C) INRIA
+c 
+c This file must be used under the terms of the CeCILL.
+c This source file is licensed as described in the file COPYING, which
+c you should have received as part of this distribution.  The terms
+c are also available at    
+c http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
       subroutine amell(du, dk, dsn2 ,n)
 c!purpose
 c calculation of the jacobi's elliptic function am(u,k)
index 85daa4e..bef3a36 100644 (file)
@@ -1,3 +1,12 @@
+c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+c Copyright (C) INRIA
+c 
+c This file must be used under the terms of the CeCILL.
+c This source file is licensed as described in the file COPYING, which
+c you should have received as part of this distribution.  The terms
+c are also available at    
+c http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 C/MEMBR ADD NAME=BELAN,SSI=0
       subroutine belan(ak1,ak2,ak3,ak4,eps,a,wc,wr,
      1                 zerr,zeri,polr,poli)
index bd38647..0980a34 100644 (file)
@@ -1,3 +1,12 @@
+c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+c Copyright (C) INRIA
+c 
+c This file must be used under the terms of the CeCILL.
+c This source file is licensed as described in the file COPYING, which
+c you should have received as part of this distribution.  The terms
+c are also available at    
+c http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
       subroutine desi24(ndeg,adelp,adels,adelta,dsk,ugc,ogc,ack,
      *nj,nh,dk,dks,dcap02,dcap04,acx,ac,rdelp,rdels,sfa,spr,spi)
 c!purpose
index 4e2ff25..5f37588 100644 (file)
@@ -1,3 +1,12 @@
+c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+c Copyright (C) INRIA
+c 
+c This file must be used under the terms of the CeCILL.
+c This source file is licensed as described in the file COPYING, which
+c you should have received as part of this distribution.  The terms
+c are also available at    
+c http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 C/MEMBR ADD NAME=DSPLN,SSI=0
       subroutine dspln(n,x,f,d,a,fail)
 c!but
index b8cb7f4..50e617d 100644 (file)
@@ -1,3 +1,12 @@
+c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+c Copyright (C) INRIA
+c 
+c This file must be used under the terms of the CeCILL.
+c This source file is licensed as described in the file COPYING, which
+c you should have received as part of this distribution.  The terms
+c are also available at    
+c http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 C/MEMBR ADD NAME=NSTABL,SSI=0
       subroutine nstabl(a,n,w,ist)
 c     test de stabilite
index 859e6ac..78e18b2 100644 (file)
@@ -1,4 +1,3 @@
-c     =======================================
 c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 c Copyright (C) INRIA
 c 
index de4d202..f1bdc4e 100644 (file)
@@ -1,3 +1,12 @@
+c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+c Copyright (C) INRIA
+c 
+c This file must be used under the terms of the CeCILL.
+c This source file is licensed as described in the file COPYING, which
+c you should have received as part of this distribution.  The terms
+c are also available at    
+c http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 C/MEMBR ADD NAME=TRBIPO,SSI=0
       subroutine trbipo(nmaxi,maxdeg,nj,sfa,sm,spr,spi,zfa,zpr,zpi)
 c!purpose