1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) INRIA -
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
function [phi,db]=phasemag(z,mod)
//
mod_def='c' //continuous representation
//mod_def='m' //representation modulo 360 degrees
[lhs,rhs]=argn(0)
if lhs==2 then
db=20*log(abs(z))/log(10),
end
if rhs<>2 then
mod=mod_def
end
//compute first phase value in (-pi, pi]
phi1=atan(imag(z(:,1)),real(z(:,1)))
//compute phase increments in (-pi, pi]
z=z(:,2:$)./z(:,1:$-1)
dphi=atan(imag(z),real(z))
phi=cumsum([phi1 dphi],2)
if part(mod,1)<>'c' then // reset modulo 360
phi=modulo(phi,2*%pi)
end
phi=phi*180/%pi //transform in degree
endfunction
|
