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.2

1.1       maekawa     1: #
1.1.1.2 ! maekawa     2: # $Id: controls.dem,v 1.1.1.1.2.1 1999/10/11 13:25:20 lhecking Exp $
1.1       maekawa     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
1.1.1.2 ! maekawa    20: set key box
1.1       maekawa    21: #
                     22: # plot c(t) for several different damping factors s
                     23: #
                     24: 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)
                     25: pause -1 "Hit return to continue"
1.1.1.2 ! maekawa    26: reset
1.1       maekawa    27: