signal_processing plugged.
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / remez.dia.ref
index 603662c..7f98a90 100644 (file)
@@ -1,12 +1,52 @@
+
 // =============================================================================
+
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+
 // Copyright (C) 2009 - DIGITEO
+
+// Copyright (C) 2010 - DIGITEO
+
 //
+
 //  This file is distributed under the same license as the Scilab package.
+
 // =============================================================================
+
 iext = linspace(1,5,5);
+
 ds = [ones(1:5),zeros(1:5)];
+
 fg = linspace(0,0.5,10);
+
 wt = ones(fg);
+
 y = remez(iext,ds,fg,wt);
+
 if max(abs(y - [0.500 0.852 0.000 0.000])) > 0.001 then bugmes();quit;end
+
+nc=21;
+
+ngrid=nc*250;
+
+fgrid=.5*(0:(ngrid-1))/(ngrid-1);
+
+mag=ones(1:ngrid/2);
+
+mag(ngrid/2+1:ngrid)=0*ones(1:ngrid/2);
+
+weight=ones(fgrid);
+
+guess=round(1:ngrid/nc:ngrid);
+
+guess(nc+1)=ngrid;
+
+guess(nc+2)=ngrid;
+
+an=remez(guess,mag,fgrid,weight);
+************ failure to converge **********
+probable cause is machine rounding error
+- number of iterations =   5
+if the number of iterations exceeds 3 the design
+may be correct, but should be verified with an fft
+