; Marcos van Dam ; calculate the tip-tilt mirror dynamics ; ; The dynamics were originally calculated by Microgate using the frequency response provided by LLNL for the Keck II mirror: ; mirror=8.52e10 / (8.52e10 + 5.8080e7 * s + 1.32e4 * s^2 + s^3) ; ; By trial and error using the CLMP and closed loop diagnostics, I came up with the form of the transfer function below FUNCTION TTMIRRORDYNAMICS,s,telescope=telescope IF N_ELEMENTS(telescope) EQ 0 THEN telescope='Keck II' CASE telescope OF 'Keck I': BEGIN ; Keck I coefficients a=4.e10 b=7.e7 c=0. END ELSE: BEGIN ; Keck II coefficients a=2.2e10 b=5.8e7 c=1.32e4 END ENDCASE mirror=(a / (a + b * s + c * s^2 + s^3)) RETURN,mirror END