bug 6744 fix
[scilab.git] / scilab / modules / cacsd / tests / nonreg_tests / bug_6744.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2010 - INRIA - Serge Steer
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7
8 // <-- Non-regression test for bug 6744 -->
9 // <-- JVM NOT MANDATORY -->
10 // <-- Bugzilla URL -->
11 // http://bugzilla.scilab.org/show_bug.cgi?id=6744
12 //
13 // <-- Short Description -->
14 //  p_margin() returns an erroneous result
15 z=poly(0,'z');
16 num=(0.0373327-0.0405775*z-0.1413155*z^2+0.1783403*z^3);
17 den=(-0.0373327+0.0405775*z+0.1413155*z^2-0.1783403*z^3+0.04867*z^4-0.04007*z^5-0.17482*z^6+0.2*z^7);
18 H=syslin(1.0,num/den);
19 [m,fr]=p_margin(H);
20
21
22 r=repfreq(H,fr);
23 //check if fr give a response with modulus=1
24 if abs(abs(r)-1)>1d-10 then pause,end
25 if abs(m-58.03)>0.01 then pause,end
26 if abs(m-pmodulo(atand(imag(r),real(r))-180,360))>1d-10 then pause,end