d85610b7bd26ba1d77824fdfa6c65308a936e431
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / remez.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2009 - DIGITEO
4 // Copyright (C) 2010 - DIGITEO
5 //
6 //  This file is distributed under the same license as the Scilab package.
7 // =============================================================================
8
9 iext = linspace(1,5,5);
10 ds = [ones(1:5),zeros(1:5)];
11 fg = linspace(0,0.5,10);
12 wt = ones(fg);
13 y = remez(iext,ds,fg,wt);
14
15 if max(abs(y - [0.500 0.852 0.000 0.000])) > 0.001 then pause, end
16
17
18 nc=21;
19 ngrid=nc*250;
20 fgrid=.5*(0:(ngrid-1))/(ngrid-1);
21 mag(1:ngrid/2)=ones(1:ngrid/2);
22 mag(ngrid/2+1:ngrid)=0*ones(1:ngrid/2);
23 weight=ones(fgrid);
24 guess=round(1:ngrid/nc:ngrid);
25 guess(nc+1)=ngrid;
26 guess(nc+2)=ngrid;
27 an=remez(guess,mag,fgrid,weight);