5
square(-1.1,-1.1,1.1,1.1);
6
xtitle( ['loci with constant damping and constant frequencies';...
7
'in discrete plane'],' ',' ');
11
//roots of s + 2*xsi*w0*s +w0
12
//given by : w0*(-xsi+-%i*sqrt(1-sxi*xsi))
13
raci=((0:0.05:1)*%pi)'*(-xsi+%i*sqrt(ones(xsi)-xsi.*xsi))
14
// continuous --> discrete
15
raci=exp(raci);[mr,nr]=size(raci);
17
xstring(real(raci(mr-10,l)),-imag(raci(mr-10,l)),' '+string(xsi(l)),0,0);
19
plot2d(real(raci),imag(raci),1*ones(1,nr),"001");
20
plot2d(real(raci),-imag(raci),1*ones(1,nr),"001");
23
e_itheta=exp(%i*(%pi/2:0.05:%pi)')
24
zw=exp(e_itheta*w0);[mz,nz]=size(zw)
26
xstring(real(zw(1,l)),imag(zw(1,l)),' '+string(w0(l)/(2*%pi)),0,0);
29
plot2d(real(zw),imag(zw),1*ones(1,nr),"001");
30
plot2d(real(zw),-imag(zw),1*ones(1,nr),"001");