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
|
FREQ(G) Scilab Function FREQ(G)
NAME
freq - frequency response
CALLING SEQUENCE
[x]=freq(A,B,C [,D],f)
[x]=freq(NUM,DEN,f)
PARAMETERS
A, B, C, D : real matrices of respective dimensions nxn, nxp, mxn,
mxp.
NUM,DEN : polynomial matrices of dimension mxp
x : real or complex matrix
DESCRIPTION
x=freq(A,B,C [,D],f) returns a real or complex mxp*t matrix such that:
x(:,k*p:(k+1)*p)= C*inv(f(k)*eye-A)*B + D.
Thus, for f taking values along the imaginary axis or on the unit circle x
is the continuous or discrete time frequency response of (A,B,C,D).
x=freq(NUM,DEN,f) returns a real or complex matrix x such that columns
k*(p-1)+1 to k*p of x contain the matrix NUM(f(k))./DEN(f(k))
EXAMPLE
s=poly(0,'s');
sys=(s+1)/(s^3-5*s+4)
rep=freq(sys("num"),sys("den"),[0,0.9,1.1,2,3,10,20])
[horner(sys,0),horner(sys,20)]
//
Sys=tf2ss(sys);
[A,B,C,D]=abcd(Sys);
freq(A,B,C,[0,0.9,1.1,2,3,10,20])
SEE ALSO
repfreq, horner
|
