Как оно?
pro intfunction, lenght, A, F
x=fltarr(854)
y=A(1)/(sqrt(2)*A(2))
x=(lenght-A(0))/(sqrt(2)*A(2))
f=A(3)*voigt(y,x)
end
set_plot, 'win'
; Voight.pro
device,dec=0 & loadct,27,/sil
dir='X:\Students|countor_analiz'
fileps='deviation.ps'
; create gaussian
4len=854
lenght=findgen(4len)
lensize=n_elements(lenght)
maximum=1000
center=150
fwhm_g=5
fwhm_l=3
fwhm_death=fltarr(100)
center_death=fltarr(100)
maximum_death=fltarr(100)
xplot=fltarr(100)
Sig_Mois=1
for m=0,99 do begin
Sig_Mois=Sig_Mois+2
xplot(m)=sig_Mois
fwhm_death(m)=0
center_death(m)=0
maximum_death(m)=0
for j=0,9 do begin
gauss1=fltarr(lensize-1)
gauss=gaussian(lenght,[maximum,center,fwhm_g])
for i=0,lensize-2 do gauss1(i)=gauss(i)
lors=30*fwhm_l
lorens=fltarr(lors)
lorens1=fltarr(lensize-1)
;create lorens
for i=0,lensize-2 do begin
lorens1(i)=maximum/(1+(2*(lenght(i)-center)/ (2*fwhm_l))^2)
endfor
for i=0,lors-1 do begin
lorens(i)=maximum/(1+(2*(lenght(i+center-15*fwhm_l)-center)/ (2*fwhm_l))^2)
endfor
;!p.multi=[0,1,3]
prof=convol(gauss,lorens,total(lorens))
for i=0,4len-1 do begin
ran=RANDomn(seed,1)/2
prof(i)=prof(i)+ran*maximum/Sig_mois
;print,ran
endfor
;dalee deconvolvim
;weights=1/prof
;weights=sqrt(abs(prof/maximum))
weights=make_array(xlen,val=1.)
A=[center, fwhm_l, fwhm_g, maximum]
A(0)=A(0)+RANDomn(SEED,1)*5
A(2)=A(2)+RANDomn(SEED,1)*2
A(3)=A(3)+RANDomn(SEED,1)*10
yfit = CURVEFIT(lenght, prof, weights, A, SIGMA,$
FUNCTION_NAME='intfunction',/NODERIVATIVE)
;plot, lenght,yfit,xr=[0,300], yr=[0,900]
;plot, lenght,prof,xr=[0,300], yr=[0,900]
;plot,lenght,lorens1,xr=[0,300]
;oplot,lenght,gauss,col=30
;print,A, sigma
fwhm_death(m)=(fwhm_death(m)+A(2)-fwhm_g)/fwhm_g
fwhm_g=fwhm_g+1
center_death(m)=(center_death(m)+A(0)-center)/center
maximum_death(m)=(maximum_death(m)+A(3)-maximum)/maximum
endfor
fwhm_death(m)=fwhm_death(m)/10
center_death(m)=center_death(m)/10
maximum_death(m)=maximum_death(m)/10
endfor
print,a
plot,xplot,abs(fwhm_death)/fwhm_g,xst=1
s=1000
set_plot,'PS'
device, file=fileps,xsize=18,ysize=28,xoffset=1,yoffset=0,/portrait
plot,xplot,maximum_death,pos=[2*s,2*s,17*s,6*s],yr=[-1,1],xst=1,xtitle='Signal/nois',ytitle='FWHM_deviation',/noerase,/device
plot,xplot,center_death,pos=[2*s,8*s,17*s,12*s],yr=[-1,1],xst=1,xtitle='Signal/nois',ytitle='FWHM_deviation',/noerase,/device
plot,xplot,fwhm_death,pos=[2*s,14*s,17*s,20*s],yr=[-1,1],xst=1,xtitle='Signal/nois',ytitle='center_deviation',/noerase,/device
xyouts,9*s, 3*s, 'deviations',alignment=0.5, charsize=1.5,/device
device,/close
set_plot, 'win'
!p.multi=[0,1,1]
end
![]()
The device consists of the Focal Reducer (F/2.9) and of the prime focus Adapter (with tip-tilt auto-guiding system and lamps for flat-field and wavelength calibrations). |
SCOPIO User manual (ver. 8 Nov 2004): PDF |
Design of the SCORPIO: The optical part contains of the collimator (F/2.2) and the camera (F/1.8). The total focal ratio it the prime focus of the telescope is (F/2.6). The optics of the reducer compensate for the aberrations of the main mirror of the telescope, all optical surfaces have antireflecting coating in a spectral range 3500-10000 A. SCORPIO has two filter's wheels for interferometric narrow-band and middle-band filters (Wheel 1 in the focal plane) and for broad-band filters (Wheel 2 between the field lens and the collimator). The dispersers (grisms or Queensgate IFPs ET-50) installed between the collimator and the camera where the exit pupil of the optical system is located. The diameter of the collimated beam is 35mm. SCORPIO are under full remote control from a PC. The Tecktronix CCD 1024x1024 px with was used for observations in 2000-2002 years. The CCD EEV 42-40 with size 2048x2048 px are used since 2003. The project just started at October,1999 in SAO RAS. The first light on 6m telescope was observed 21 September 2000. The similar device (without adapter) was constructed by us for 2.6m telescope in the Byurakan (Armenia). Scientific advisor of the project:V.L. Afanasiev |