hold_on_tf.pro

来自「IDL语言编写的用于天文自适应光学仿真的软件CAOS V6.0的第一部分。」· PRO 代码 · 共 39 行

; $Id: hold_on_tf.pro,v 1.2 2002/03/14 11:49:11 riccardi Exp $

;+
; 
; NAME
;
;   HOLD_ON_TF
;
;
;   tf = hold_on_tf(s, T)
;
; Return the (Laplace) Tranfer Function associated to the Hold-on
; (DAC) effect for a digital signal with period T.
; s vector of complex frequency [rad/s] (s = sigma + i*omega)
; T period [s]
;
;            1 - exp(-s*T)      sinh(s*T/2)
;   TF(s) = ---------------- = ----------- * T*exp(-s*T/2)
;                  s               s*T/2
;
;-

function hold_on_tf, s, T

iu = complex(0.0, 1.0)
idx = where(s ne 0.0, count)

sT2 = s*T/2.0

tf = exp(-sT2) * T

if count ne 0 then begin
    tf[idx] = (0.5*(exp(sT2)-exp(-sT2))[idx]/sT2[idx])*tf[idx]
endif

return, tf
end