d50ee1ab0cecdf23849a30c8324c24bf51abaf0a
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / lattp.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2013 - Scilab Enterprises - Charlotte Hecquet
4 // Copyright (C) 2013 - Scilab Enterprises - Sylvestre Ledru
5 //
6 //  This file is distributed under the same license as the Scilab package.
7 // =============================================================================
8
9 //Generate the process
10 t1=0:0.1:100;
11 y1=sin(2*%pi*t1)+sin(2*%pi*2*t1);
12
13 //Covariance of y1
14 nlag=128;
15 c1=corr(y1,nlag);
16 c1=c1';
17
18 //Compute the filter with maximum order=15 and p=5
19 n=5; p=2;
20 [la1,sig1]=lattp(n,p,c1);
21 poly(0,"z");
22 assert_checkequal(la1(1),  1-1.83287468880270521*%z+2.08814963103010243*%z^2-1.02130636786998474*%z^3+0.09534323012550860*%z^4+0.39938577578836609*%z^5);
23
24 assert_checkequal(la1(2), 1+5934927594313.27051*%z-10877978567494.7891*%z^2+12393016866258.4355*%z^3-6061379344921.20703*%z^4+565855167403.415283*%z^5+2370325661503.31836*%z^6);
25
26 assert_checkequal(sig1, list());
27