filter_prod.pro

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

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

pro filter_prod, c_num_p, c_den_p, freq, input_tf, $
output_tf=output_tf, LOG=log, HI2LOW=hi2low

if keyword_set(hi2low) then begin
	c_num=rotate(c_num_p,2) 
	c_den=rotate(c_den_p,2) 
endif else begin
	c_num=c_num_p
	c_den=c_den_p
endelse

filter_tf = poly(2*!pi*complex(0,1)*freq, c_num)/ $
			poly(2*!pi*complex(0,1)*freq, c_den)

pmulti_save = !p.multi
!p.multi = [0,1,2]

output_tf=filter_tf*input_tf

plot_oo, freq, abs(output_tf), xtitle='Frequency [Hz]', ytitle='bitC/bitP'
plot_oi, freq, imaginary(alog(output_tf))*180/!pi, xtitle='Frequency [Hz]',$
								ytitle='Phase [deg]'
!p.multi=pmulti_save
!p.multi(0)=0
return
end