home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 22 gnu / 22-gnu.zip / GNUPDEMO.ZIP / electron.demo < prev    next >
Text File  |  1992-03-25  |  2KB  |  102 lines
  1. #
  2. # $Id: electron.demo,v 3.26 92/03/24 22:32:20 woo Exp Locker: woo $
  3. #
  4. # Electronics demo
  5. #
  6. # Bipolar Transistor (NPN) Mutual Characteristic
  7. Ie(Vbe)=Ies*exp(Vbe/kT_q)
  8. Ic(Vbe)=alpha*Ie(Vbe)+Ico
  9. alpha = 0.99
  10. Ies = 4e-14
  11. Ico = 1e-09
  12. kT_q = 0.025
  13. set dummy Vbe
  14. set grid
  15. set offsets
  16. set nolog
  17. set nopolar
  18. set samples 160
  19. set title "Mutual Characteristic of a Transistor"
  20. set xlabel "Vbe (base emmitter voltage)"
  21. set xrange [0 : 0.75]
  22. set ylabel "Ic (collector current)"
  23. set yrange [0 : 0.005]
  24. set key .2,.0045
  25. set format y "%.4f"
  26. plot Ic(Vbe)
  27. set format "%g"
  28.  
  29. pause -1 "Hit return to continue"
  30.  
  31. # Junction Field Effect Transistor (JFET) Mutual Characteristic
  32. # drain current above pinch off
  33. Ida(Vd)=Ido*(1-Vg/Vp)**2
  34. # drain current below pinch off
  35. Idb(Vd)=Ido*(2*Vd*(Vg-Vp)-Vd*Vd)/(Vp*Vp)
  36. # drain current
  37. Id(Vd)= (Vd>Vg-Vp) ? Ida(Vd) : Idb(Vd)
  38. # drain current at zero gate voltage
  39. Ido = 2.5
  40. # pinch off voltage
  41. Vp = -1.25
  42. # gate voltage
  43. Vg = 0
  44. set dummy Vd
  45. set nogrid
  46. set nokey
  47. set offsets 0, 1, 0, 0
  48. set title "JFET Mutual Characteristic"
  49. set xlabel "Drain voltage Vd (V)"
  50. set xrange [0 : 4]
  51. set ylabel "Drain current Id (mA)"
  52. set yrange [0 : 5]
  53. set label "-0.5 Vp" at 4.1,0.625
  54. set label "-0.25 Vp" at 4.1,1.4
  55. set label "0" at 4.1,2.5
  56. set label "Vg = 0.5 Vp" at 4.1,3.9
  57. plot Vg=0.5*Vp,Id(Vd),Vg=0.25*Vp,Id(Vd),Vg=0,Id(Vd),Vg=-0.25*Vp,Id(Vd)
  58. set nolabel
  59.  
  60. pause -1 "Hit return to continue"
  61.  
  62. # amplitude frequency response
  63. A(jw) = ({0,1}*jw/({0,1}*jw+p1)) * (1/(1+{0,1}*jw/p2))
  64. p1 = 10
  65. p2 = 10000
  66. set dummy jw
  67. set grid
  68. set key
  69. set logscale xy
  70. set offsets 0, 0, 0, 0
  71. set title "Amplitude Frequency Response"
  72. set xlabel "jw (radians)"
  73. set xrange [1.1 : 90000.0]
  74. set ylabel "magnitude of A(jw)"
  75. set autoscale  y
  76. plot abs(A(jw))
  77.  
  78. pause -1 "Hit return to continue"
  79.  
  80. # phase frequency response
  81. set nolog y
  82. set logscale x
  83. set title "Phase Frequency Response"
  84. set ylabel "Phase of A(jw) (degrees)"
  85. plot 180/pi*arg(A(jw))
  86.  
  87. pause -1 "Hit return to continue"
  88.  
  89. # undo what we've done
  90. set dummy x
  91. set nogrid
  92. set offsets 0,0,0,0
  93. set title ""
  94. set ylabel ""
  95. set xlabel ""
  96. set xrange [-10:10]
  97. set autoscale xy
  98. set key
  99. set format xy "%g"
  100. set nolabel
  101. set nolog
  102.