home *** CD-ROM | disk | FTP | other *** search
- ! Nice example to send 16-Bit gle's to nirwana
- ! by Axel Rohde
- ! size 25 25 ! too large for DIN A4 -> shrink to fit
- ! set all colours to black for b/w printers
- size 19 19
- set font pldr fontlwidth 0.03 cap round
- beginn = -0.15
- endd = 0.6
- stepp = 0.001
- pi = 3.14156294
-
- ! shrink to fit
- begin scale 0.75 0.75
-
- amove 0 14
- begin graph
- nobox
- size 25 10
- xaxis min 0 max 10
- yaxis min -0.3 max 2*pi+0.3 nticks 6
- yplaces 0 pi/2 pi 3/2*pi 2*pi
- ynames "0" "\pi/2" "\pi" "3\pi/2" "2\pi"
- xtitle "time/a.u." hei .60
- ytitle "phase/radian" hei .60
- xlabels hei .4
- ylabels hei .4
-
- let d1 = pi+2*atn(0.3-sqrt(1-0.3**2)*tan(25*x*sqrt(1-0.3**2)/2)) from beginn to endd step stepp
- d1 line color black lstyle 0 xmin beginn xmax endd
- let d2 = pi+2*atn(0.5-sqrt(1-0.5**2)*tan(25*x*sqrt(1-0.5**2)/2)) from beginn to endd step stepp
- d2 line color red lstyle 2 xmin beginn xmax endd
- let d3 = pi+2*atn(0.8-sqrt(1-0.8**2)*tan(25*x*sqrt(1-0.8**2)/2)) from beginn to endd step stepp
- d3 line color blue lstyle 4 xmin beginn xmax endd
- let d4 = pi+2*atn(0.99-sqrt(1-0.99**2)*tan(25*x*sqrt(1-0.99**2)/2)) from beginn to endd step stepp
- d4 line color green lstyle 9 xmin beginn xmax endd
- let d5 = pi+2*atn(0-sqrt(1-0**2)*tan(25*x*sqrt(1-0**2)/2)) from beginn to endd step stepp
- d5 line color cyan lstyle 3 xmin beginn xmax endd
- title "driven Van der Pol oscillator - Adler's equation" hei .7
- end graph
-
- begin key
- hei 0.25
- position tr
- line color green lstyle 9 text "\alpha=0.99"
- line color blue lstyle 4 text "\alpha=0.8"
- line color red lstyle 2 text "\alpha=0.5"
- line color black lstyle 0 text "\alpha=0.3"
- line color cyan lstyle 3 text "\alpha=0"
- end key
-
- amove xg(3) yg(3.14)
- set hei 0.35
- text fast
- amove xg(1.5) yg(5.5)
- set hei 0.35
- text slow
-
- amove 0 0
- begin graph
- nobox
- size 25 14
- xaxis min 83 max 210
- yaxis min -3 max 3
- data resp.dat
- d1 line smooth lstyle 1 color red
- !d2 line color blue
- !d3 line color green
- let d4 = d3*1.2
- d4 line smooth lstyle 2 color blue
- xtitle "time/a.u." hei 0.6
- ytitle "amplitude/a.u." hei .6
- title "driven Van der Pol oscillator - numerical results" hei .7
- end graph
- begin key
- hei 0.4
- position tr
- line lstyle 1 color red text "response"
- line lstyle 2 color blue text "driver"
- end key
- set hei 0.7
- amove 0.9 yg(2.6)
- text b
- set hei 1.3
- amove 1 yg(2.75)
- marker circle
-
- set hei 0.7
- amove 0.85 yg(9.1)
- text a
- set hei 1.3
- amove 1 yg(9.25)
- marker circle
-
- amove 9 yg(2.3)
- set hei 0.35
- text fast
- amove 6 yg(2.3)
- set hei 0.35
- text slow
-
- amove 5.5 2.5
- text out of phase
- amove 9.4 2.5
- text in phase
-
- end scale
-