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
|
// 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 [] = sgrid(zeta,wn,col)
// sgrid()
// sgrid(Z,Wn)
// sgrid('new')
[lhs,rhs]=argn(0)
axes=gca();
select rhs
case 0 then
wn= 0:1:10;
zeta = linspace(0,1,10);
col=3;
case 1 then
if type(zeta)<>10 then
error(msprintf(_("%s: Wrong type for input argument #%d: String array expected.\n"),"sgrid",1))
end
if size(zeta,'*')<>1 then
error(msprintf(_("%s: Wrong size for input argument #%d: A string expected.\n"),"sgrid",1))
end
wn = 0:1:10;
zeta = [ 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1 ];
col=3;
clf()
axes.data_bounds=[-20,-20;20,20];axes.axes_visible='on';
case 2 then
col=3;
end
wmax = 10 .^(floor(log10(max(abs(axes.data_bounds)))));
// building a grid
zx = 0:.01:1;
[rx,cx]=size(wn);[ry,cy]=size(zx);
w=wn.*.ones(cy,1);z=zx'.*.ones(1,cx);
// plot : part I
re = - w .* z;
[zr,zc] = size(z);
im = w .* sqrt( 1 - z .* z );
plot2d(re,im,style=col*ones(1,zc));
plot2d(re,-im,style=col*ones(1,zc));
// info on curves
[rer,rec] = size(re)
axes.clip_state = "clipgrf";
xnumb(re(1,:),im(1,:),wn);
axes.clip_state = "off";
// building an other grid
wn = [0,wn,2*wmax];
[rx,cx]=size(wn);[ry,cy]=size(zeta);
w=wn.*.ones(cy,1);z=zeta'.*.ones(1,cx);
// plot part II
[zr,zc] = size(z);
re = -w .* z;
im = w .* sqrt( 1 - z .* z );
plot2d(re',im',style=col*ones(1,zr));
plot2d(re',-im',style=col*ones(1,zr));
// info on each curve ( straight lines )
[rer,rec] = size(re)
axes.clip_state = "clipgrf";
xnumb(re(:,$)',im(:,$)',zeta);
axes.clip_state = "off";
endfunction
|
