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
|
function []=m_circle(gain)
// Copyright INRIA
[lhr,rhs]=argn(0)
//
d36=11;
if rhs<1 then
gain=[-12 -8 -6 -5 -4 -3 -2 -1.4 -1 -.5 0.25 0.5 0.7 1 1.4 ..
2 2.3 3 4 5 6 8 12]
else
gain=matrix(gain,1,prod(size(gain)))
end
//
titre='isogain contours for y/(1+y)'
l10=log(10);
//
lambda=exp(l10*gain/20)
rayon=lambda./(lambda.*lambda-ones(lambda))
centre=-lambda.*rayon
rayon=abs(rayon)
rect=[mini(centre-rayon),mini(-rayon),maxi(centre+rayon),maxi(rayon)];
//
plot2d(0,0,1,"000"," ");
//xtitle(titre,'R','I');
llrect=xstringl(0,0,'1')
hx=llrect(3);
//
for i=1:prod(size(gain))
sgn=1;if 2*int(i/2)==i then sgn=-1,end
g=string(gain(i)),ng=length(g)
if gain(i)<0 then
w=0:0.03:%pi;
xx=centre(i)+rayon(i)*sin(%pi/3)+hx
xy=sgn*cos(%pi/3)*rayon(i)
else
xx=centre(i)-sin(%pi/3)*rayon(i)-hx*ng
xy=sgn*cos(-%pi/3)*rayon(i)
w=-%pi:0.03:0;
end;
n=prod(size(w))
rf=centre(i)*ones(w)+rayon(i)*exp(%i*w);
rf=rf';
plot2d(real(rf),imag(rf),[3],"000")
plot2d(real(rf),-imag(rf),[3],"000")
// xstring(xx,xy,g)
end;
|
