Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.naic.edu/~emuller/cor_rms.pro Дата изменения: Mon Aug 9 16:29:44 2004 Дата индексирования: Mon Oct 1 21:48:02 2012 Кодировка: Поисковые слова: http astrokuban.info astrokuban

pro cor_rms,dat,x1=x1,x2=x2
;program to measure some fundamental stats for the data input

N1=float(dat.b1.h.cor.lagsbcout)
if not keyword_set(x1) then x1=0


if n_elements(dat.b1.d(0,*)) eq 2 then begin
b1_2pol='T'
x2=n_elements(dat.b1.d(*,0))-1
rms1a=sqrt(total((dat.b1.d(x1:x2,0))^2)/N1)
rms1b=sqrt(total((dat.b1.d(x1:x2,1))^2)/N1)
max1a=max((dat.b1.d(x1:x2,0)))
max1b=max((dat.b1.d(x1:x2,1)))

endif else begin
x2=n_elements(dat.b1.d(*))-1
rms1a=sqrt(total((dat.b1.d(x1:x2))^2)/N1)
max1a=max((dat.b1.d(x1:x2)))
endelse

N2=float(dat.b2.h.cor.lagsbcout)
if n_elements(dat.b2.d(0,*)) eq 2 then begin
b2_2pol='T'
rms2a=sqrt(total((dat.b2.d(x1:x2,0))^2)/N2)
rms2b=sqrt(total((dat.b2.d(x1:x2,1))^2)/N2)
max2a=max((dat.b2.d(x1:x2,0)))
max2b=max((dat.b2.d(x1:x2,1)))
endif else begin
rms2a=sqrt(total((dat.b2.d(x1:x2))^2)/N1)
max2a=max((dat.b2.d(x1:x2)))
endelse

N3=float(dat.b2.h.cor.lagsbcout)
if n_elements(dat.b3.d(0,*)) eq 2 then begin
b3_2pol='T'
rms3a=sqrt(total((dat.b3.d(x1:x2,0))^2)/N3)
rms3b=sqrt(total((dat.b3.d(x1:x2,1))^2)/N3)
max3a=max((dat.b3.d(x1:x2,0)))
max3b=max((dat.b3.d(x1:x2,1)))
endif else begin
rms3a=sqrt(total((dat.b3.d(x1:x2))^2)/N3)
max3a=max((dat.b3.d(x1:x2)))
endelse

N4=float(dat.b2.h.cor.lagsbcout)
if n_elements(dat.b4.d(0,*)) eq 2 then begin
b4_2pol='T'
rms4a=sqrt(total((dat.b4.d(x1:x2,0))^2)/N4)
rms4b=sqrt(total((dat.b4.d(x1:x2,1))^2)/N4)
max4a=max((dat.b4.d(x1:x2,0)))
max4b=max((dat.b4.d(x1:x2,1)))
endif else begin
rms4a=sqrt(total((dat.b4.d(x1:x2))^2)/N4)
max4a=max((dat.b4.d(x1:x2)))
endelse


print,strcompress("RMS pol 1 board 1: "+string(rms1a))
if b1_2pol eq'T' then print,strcompress("RMS pol 2 board 1: "+string(rms1b))
print,strcompress("RMS pol 1 board 2: "+string(rms2a))
if b1_2pol eq 'T' then print,strcompress("RMS pol 2 board 2: "+string(rms2b))
print,strcompress("RMS pol 1 board 3: "+string(rms3a))
if b1_2pol eq 'T' then print,strcompress("RMS pol 2 board 3: "+string(rms3b))
print,strcompress("RMS pol 1 board 4: "+string(rms4a))
if b1_2pol eq 'T' then print,strcompress("RMS pol 2 board 4: "+string(rms4b))
print,"-------------------"
print,strcompress("Max pol 1 board 1: "+string(max1a))
if b1_2pol eq 'T' then print,strcompress("Max pol 2 board 1: "+string(max1b))
print,strcompress("Max pol 1 board 2: "+string(max2a))
if b1_2pol eq 'T' then print,strcompress("Max pol 2 board 2: "+string(max2b))
print,strcompress("Max pol 1 board 3: "+string(max3a))
if b1_2pol eq 'T'then print,strcompress("Max pol 2 board 3: "+string(max3b))
print,strcompress("Max pol 1 board 4: "+string(max4a))
if b1_2pol eq 'T' then print,strcompress("Max pol 2 board 4: "+string(max4b))

end