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
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 1984-2011 - INRIA - Serge STEER
//
// 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.1-en.txt
function [frq,rep,splitf]=repfreq(sys,fmin,fmax,pas)
pas_def="auto";
l10=log(10);
[lhs,rhs]=argn(0)
//discretization
if and(typeof(sys)<>[ "rational" "state-space" ]) then
error(msprintf(gettext("%s: Wrong type for input argument #%d: Linear dynamical system expected.\n"),"repfreq",1))
end
dom=sys.dt
if dom==[]|dom==0 then error(96,1),end
if dom=="d" then dom=1;end
select rhs
case 1 then
pas=pas_def
if dom=="c" then fmax=1.d3; else fmax=1/(2*dom),end
fmin=0
case 2 then
if type(fmin)==10 then
rhs=1
pas=pas_def
if dom=="c" then fmax=1.d3; else fmax=1/(2*dom),end
if fmin=="sym" then
fmin=-fmax
else
fmin=0
end
else
frq=fmin
end
case 3 then
pas=pas_def
case 4 then ,
else
error(msprintf(gettext("%s: Wrong number of input arguments: %d to %d expected.\n"), "repfreq",1,4))
end;
splitf=1
if rhs<>2 then
if fmin==[] then
fmin=0,
end
if fmax==[]|fmax==%inf then
if dom=="c" then
fmax=1.d3;
else
fmax=1/(2*dom);
end
end
if type(pas)==1 then
splitf=1
eps=1.e-14
if fmin<0&fmax>=0 then
frq=- [exp(l10*((log(eps)/l10):pas:(log(-fmin)/l10))) -fmin];
if fmax>eps then
frq1=[exp(l10*((log(eps)/l10):pas:(log(fmax)/l10))) fmax];
frq=[frq($:-1:1) frq1]
else
frq=frq($:-1:1);
end
elseif fmin<0&fmax<0 then
frq= [exp(l10*((log(-fmax)/l10):pas:(log(-fmin)/l10))) -fmin];
frq=-frq($:-1:1);
elseif fmin >= fmax then
error(msprintf(gettext("%s: Wrong value for input arguments #%d and #%d: %s < %s expected.\n"),..
"repfreq",2,3,"fmin","fmax"));
else
fmin=max(eps,fmin);
frq=[exp(l10*((log(fmin)/l10):pas:(log(fmax)/l10))) fmax];
end
else
[frq,bnds,splitf]=calfrq(sys,fmin,fmax)
end;
end
//
typ=sys(1)
select typ(1)
case "r" then
[n,d]=sys(["num","den"]),
[mn,nn]=size(n)
if nn<>1 then error(95,1),end
if dom=="c" then
rep=freq(n,d,2*%pi*%i*frq),
else
rep=freq(n,d,exp(2*%pi*%i*dom*frq)),
end;
case "lss" then
[a,b,c,d,x0]=sys(2:6),
[mn,nn]=size(b)
if nn<>1 then error(95,1),end
if dom=="c" then
rep=freq(a,b,c,d,2*%pi*%i*frq)
else
rep=freq(a,b,c,d,exp(2*%pi*%i*dom*frq))
end;
else error(97,1),
end;
//representation
if lhs==1 then frq=rep,end
endfunction
|
