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
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) ????-2008 - INRIA
// Copyright (C) 2010 - DIGITEO - Allan CORNET
//
// This file is released under the 3-clause BSD license. See COPYING-BSD.
function demo_bode()
//////////////////
//BODE PLOTS//////
//////////////////
//First plot using state space representation
//
// xd=-2*%pi*x+u
// y=18*%pi*x+u
//
// h(s)=18*%pi/(s+2*%pi)
//
a = -2*%pi;
b = 1;
c = 18*%pi;
d = 1;
sl = syslin("c",a,b,c,d);
my_handle = scf(100001);
if is_handle_valid(my_handle) then
clf(my_handle,"reset");
bode(sl,.1,100)
halt(_("Press Return to continue ... \n"));
end
//Second plot using poly
s=poly(0,"s");
h1=1/real((s+2*%pi*(15+100*%i))*(s+2*%pi*(15-100*%i)));
h1=syslin("c",h1)
if is_handle_valid(my_handle) then
clf(my_handle,"reset");
bode(h1,10,1000,.01),
halt(_("Press Return to continue ... \n"));
end
//Third plot which combines the first two plots
if is_handle_valid(my_handle) then
clf(my_handle,"reset");
bode([h1;sl],.1,1000,.01);
halt(_("Press Return to continue ... \n"));
end
// series connection
if is_handle_valid(my_handle) then
clf(my_handle,"reset");
demo_viewCode("bode.dem.sce");
bode(h1*sl,.1,1000,.01);
end
endfunction
demo_bode();
clear demo_bode;
|
