signal_processing plugged.
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / remez.dia.ref
1
2 // =============================================================================
3
4 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
5
6 // Copyright (C) 2009 - DIGITEO
7
8 // Copyright (C) 2010 - DIGITEO
9
10 //
11
12 //  This file is distributed under the same license as the Scilab package.
13
14 // =============================================================================
15
16 iext = linspace(1,5,5);
17
18 ds = [ones(1:5),zeros(1:5)];
19
20 fg = linspace(0,0.5,10);
21
22 wt = ones(fg);
23
24 y = remez(iext,ds,fg,wt);
25
26 if max(abs(y - [0.500 0.852 0.000 0.000])) > 0.001 then bugmes();quit;end
27
28 nc=21;
29
30 ngrid=nc*250;
31
32 fgrid=.5*(0:(ngrid-1))/(ngrid-1);
33
34 mag=ones(1:ngrid/2);
35
36 mag(ngrid/2+1:ngrid)=0*ones(1:ngrid/2);
37
38 weight=ones(fgrid);
39
40 guess=round(1:ngrid/nc:ngrid);
41
42 guess(nc+1)=ngrid;
43
44 guess(nc+2)=ngrid;
45
46 an=remez(guess,mag,fgrid,weight);
47 ************ failure to converge **********
48 probable cause is machine rounding error
49 - number of iterations =   5
50 if the number of iterations exceeds 3 the design
51 may be correct, but should be verified with an fft
52