sobota, 13 marca 2021 o 17:56:11 UTC+1 Piotr Gałka napisał(a):
ja nie wiem czy on moje posty widzi, niech mu ktoś zacytuje, bo takie podałem:
a tu w pythonie jak liczyć:
#-*- coding: utf-8 -*- import time import math
def WyliczKatRzeczywisty(BudSzerGeograficzna,SolarKat,SolarAzymut,nrdnia,godzina,minuta): godzina=godzina+float(minuta)/60 deklinacja=23.45*math.sin((math.pi/180)*360*(nrdnia+284)/365) katgodzinowy=15*(godzina-12) cosfi=math.sin((math.pi/180)*deklinacja)*(math.sin((math.pi/180)*BudSzerGeograficzna)*math.cos((math.pi/180)*SolarKat)-math.cos((math.pi/180)*BudSzerGeograficzna)*math.sin((math.pi/180)*SolarKat)*math.cos((math.pi/180)*SolarAzymut))+math.cos((math.pi/180)*deklinacja)*(math.cos((math.pi/180)*BudSzerGeograficzna)*math.cos((math.pi/180)*SolarKat)*math.cos((math.pi/180)*katgodzinowy)+math.sin((math.pi/180)*BudSzerGeograficzna)*math.sin((math.pi/180)*SolarKat)*math.cos((math.pi/180)*SolarAzymut)*math.cos((math.pi/180)*katgodzinowy)+math.sin((math.pi/180)*SolarKat)*math.sin((math.pi/180)*SolarAzymut)*math.sin((math.pi/180)*katgodzinowy)) kat=(180/math.pi)*math.acos(cosfi) return (katgodzinowy,90-kat) #funkcja zwraca azymut w stopniach, 0 dla poludnia, ujemne kier wschod, dodatnie zachod
BudSzerGeograficzna=51.1209 SolarKat=11 SolarAzymut=-65 ksmax=0
for nrdnia in range (1,365): for godzina in range (3,22): #for minuta in range (0,60): minuta=0 ks=WyliczKatRzeczywisty(SolarKat,SolarAzymut,nrdnia,godzina,minuta)[1] if ks<0: ks=0 else: pass print nrdnia,godzina,minuta,"kąt solarny:","%.2f" %ks if ks>ksmax: ksmax=ks print "największy znaleziony kąt:"nrdnia,godzina,minuta,ksmax
b.