Star Test simulator, (addestriamoci a fare lo star test)

Rugositàààà... adesso è tutto così semplice....

E' lambda/30 rms.

Nicola: il Suiter!

Un altro trucco dello star test è farlo muovendo uniformemente il fuoco per osservare come la luce collassa nel fuoco e riemerge. In questo esempio, incidentalmente, il momento in cui si passa per il fuoco è vicino a un minimo della turbolenza (vedi sotto).
L'astigmatismo che era difficile da vedere ora è evidente. La stella collassa da una parte e riemerge dall'altra.



E questo è uno star test sul mio ex 40 cm. http://autocostruttori.blogspot.it/2008/06...-telescope.html
Si nota un po' di astigmatismo, ma poco, (dovuto agli appoggi laterali che poi ho regolato).
Da notare che la turbolenza NON è seeing atmosferico.


______


PS mi raccomando non dite niente ai rifrattoristi perché queste cose nel loro mondo non esistono. Le stelle sono sempre così: grandi e tonde (ma tonde che più tonde non si può... per la grandezza si rimedia usando ingrandimenti a due cifre):-)
N.B. non c'è niente da mettere a fuoco: è già a fuoco.

Un po' alla volta. x Mauro: uso "diffract.r" che sono alcune classi in R. www.visi.com/~mkoehler/optics_tools/

Lo script che fa l'ultima animazione è questo:

size <- 101
t <- turb.ft(size*11,size*11)
t <- t/sd(as.vector (t))
t2 <- turb.ft(size*11,size*11)
t2 <- t2/sd(as.vector (t2))
p <- pupil(size)
p$obstruct(0.20)
step <- 4
step2 <- 2
nsteps <- floor((size*10)/step)
def=7
star <- star.test(p,def,size=100,pad=4,draw=FALSE)
zmax <- 2.5*max(star)
gamma <- 2.2
rmsi <- vector(mode="numeric",length=nsteps)
library("Cairo")

p$add.astig(0.07,1)
p$add.roughness(0.03,4)
p$add.micro(0.02)

for (i in 1:nsteps) {
istr=format(100000+i)
Cairo(paste("turbulence",istr,".png",sep=""), type="png", width=400, height=400)
tt <- 0.5*t[((i-1)*step+1) : ((i-1)*step+size),(size*1+1):(size*2)]+0.5*t2[((i-1)*step2+1) : ((i-1)*step2+size),(size*1+1):(size*2)]
rms <- sd(as.vector (tt))*2
rmsi[i] <- sd(as.vector(p$wf()),na.rm=TRUE)
p$set.turb2(tt,rms)
star <- star.test(p,def*(1-2*i/nsteps),size=100,pad=4,draw=FALSE)
zmax = 0.9*zmax + 0.1*max(star)
image ((star/zmax)^(1/gamma), zlim=c(0,1), col=gray256, asp=1, bty='n', axes=FALSE)
dev.off()
}

plot(rmsi,type="l")
mean(rmsi)

Poi bisogna comporre le immagini .png in un filmato (lo faccio con Mathematica).

Serve anche questo (il modello di turbolenza di Kolmogorov/Fried)


turb.ft <- function (size, friedratio=1, pad=2) {
n <- nextn (size*pad)
b <- matrix (rnorm (n^2)+1i*rnorm(n^2), n, n)
rho <- make.rho.fdom (n)
b <- b * sqrt (0.023*rho^(-11/3)/friedratio^(5/3))
b[rho == 0] <- 0
resize (Re (ift (b)), size, size)/ (2*pi)
}



Intendi la stella del rifrattore vero

Altre due paperette in palio. Un difetto è facile, l'altro difficile.



PS il movimento del fuocheggiatore non torna indietro. Va da extra ad intrafocale.

Bisogna riuscire a leggere attraverso la turbolenza. Eccolo senza turbolenza. E' scollimato (lambda/14 rms) e c'è un po' di rugosità di tipo leggero e ondulazioni a granfe scale (che di fatto non si vedono nemmeno, pensavo si vedessero meglio).

Questo è un difetto grave, e però non è facile da indovinare: