Return to controls.dem CVS log

Up to [local] / OpenXM_contrib / gnuplot / demo

Annotation of OpenXM_contrib/gnuplot/demo/controls.dem, Revision 1.1.1.1

1.1       maekawa     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: