Sphère céleste

Miroir en Caelostat ou Sidérostat

Latitude du lieu

Position du Soleil

Angle horaire

Déclinaison

Rotation miroir

Rot. Est-Ouest

Rot. Nord-Sud

Rotation diurne

Normale au miroir

PhM - Observatoire de Lyon
février 2015

Expérimenter et vérifier que :
- quelle que soit la direction du Soleil, on peut toujours ramener 
   le rayon réfléchi dans la direction horizontale Sud
- pour compenser la rotation diurne (curseur H), il suffit de faire 
   tourner le miroir d'Est en Ouest avec le curseur α.
- cet angle est la moitié de la rotation H.

β = 54 (slider from 0 to 180, step 0.1)
a = 0.2
Mir0 = Polygon((a, (-a)), (a, a), 4)
α_M = -30.6 (slider from -90 to 90, step 0.1)
Sph_C = Sphere((0, 0, 0), 1)
c1 = Circle((0, 0), 1)
φ = 45.6 (slider from -90 to 90, step 0.1)
axePole = Line((0, 0, 0), Vector[((-cos(φ°)), 0, sin(φ°))])
vPole = Vector(((-cos(φ°)), 0, sin(φ°)))
Mir1 = Rotate(Mir0, β°, yAxis)
vEqu = Vector((sin(φ°), 0, cos(φ°)))
vMir = UnitVector(Mir0)
Mir2 = Rotate(Mir1, ((-α_M))°, axePole)
H = 15.1 (slider from 0 to 360, step 0.1)
δ = 12 (slider from -23 to 23, step 0.1)
vInc = Rotate(Vector[(1; 0°; δ° + (90 - φ)°)], ((-H))°, axePole)
vMir2 = UnitVector(Mir2)
vRefl = Rotate(vInc, 180*°, Line[(0, 0, 0), vMir2])
rInc = Ray((0, 0, 0), vInc)
rRefl = Ray((0, 0, 0), vRefl)
b = true
S = Translate((0, 0, 0), vInc)
cdec0 = Rotate(Circle[(0, 0), cos(δ°)], (φ + 90)°, yAxis)
cdec1 = Translate(Rotate[Circle[(0, 0, 0), cos(δ°), xOyPlane], (φ + 90)°, yAxis], Vector[(vPole * sin(δ°))])
S' = Intersect(OrthogonalLine[S, axePole, space], axePole)
sSS' = Segment(S, S')
cmer = Circle((1, 0, 0), (0, 0, 1), (-1, 0, 0))
sEO = Segment((0, -1, 0), (0, 1, 0))
sSud = Segment((0, 0, 0), (1, 0, 0))
P = Rotate((0, 0, 1), ((-(90 - φ)))°, yAxis)
P' = Rotate((0, 0, -1), ((-(90 - φ)))°, yAxis)
cequat = Rotate(Circle[(0, 0), 1], (φ + 90)°, yAxis)
O = (0, -1, 0)
E = (0, 1, 0)
Sud = (1, 0, 0)
c = true
d = true
N = (-1, 0, 0)