home *** CD-ROM | disk | FTP | other *** search
/ c't freeware shareware 1999 February / CT_SW9902.ISO / mac / software / wissen / daten / gnuplot.hqx / gnuplot.2.0b4 / Scripts / demo / controls.dem < prev    next >
Text File  |  1997-05-03  |  771b  |  32 lines
  1. #
  2. # $Id: controls.dem,v 1.2 1993/09/27 17:11:18 alex Exp $
  3. #
  4. # warning:  this demo is SLOW on PCs without math coprocessors!
  5. #
  6. # From _Automatic_Control_Systems_, fourth ed., figure 6-14
  7. # transient response of a second-order system to a unit step input function
  8. #
  9. damp(t) = exp(-s*wn*t)/sqrt(1.0-s*s)
  10. per(t) = sin(wn*sqrt(1.0-s**2)*t - atan(-sqrt(1.0-s**2)/s))
  11. c(t) = 1-damp(t)*per(t)
  12. #
  13. #    wn is natural undamped frequency
  14. #    s is damping factor
  15. #
  16. wn = 1.0
  17. set xrange [0:13]
  18. set samples 50
  19. set dummy t
  20. #
  21. # plot c(t) for several different damping factors s
  22. #
  23. plot s=.1,c(t),s=.3,c(t),s=.5,c(t),s=.7,c(t),s=.9,c(t),s=1.0,c(t),s=1.5,c(t),s=2.0,c(t)
  24. pause -1 "Hit return to continue"
  25.  
  26. # undo what we have done
  27. set xrange [-10:10]
  28. set autoscale xy
  29. set samples 160
  30. set dummy x
  31.  
  32.