| Zobacz poprzedni temat :: Zobacz następny temat |
| Autor |
Wiadomość |
Gość
|
 Wysłany: Nie 13:17, 10 Gru 2006 Temat postu: AstroSfer |
|
|
Dobra więc tak napisze tutaj wszysktie zadania domowe które mam z astronomii sferycznej - zastrzegam sobie że będa niepokoleii, jednak przy tych co mam date to dopisze ją w nawiasie. Rozwiązania zadań proponuje wklejać tutaj pisząc albo nr zadania który ja podam albo najpierw całą treść a poniżej rozwiązanie ja tak zrobie tych zadań co mam odpo wpisze tutaj:
No i jeszcze jedno przepraszam że dopiero dziś to umieszam - ale komp mi zdechł drugie to sorka za brak liter greckich ale nie mam jeszcze anie oo ani worda...
No więc zaczynam:
Zad 1
W Pułkowie fi=59st46' obserwowano słońce w prawdziw południe, zmieszono odległość zenitalną słńca Z0=39st33', czas gwiazdowy tej obserwacji wynosił S=3h57m41s.
Oblicz alfe(rekrastencje)słońca i delte(deklinacje)słońca
ODPOWIEDZI: als=3h52m52s dels=20st13'
Różnica czasów5m16s
Zad2
Oblicz refrakcję R zgodnie z modelem płaskiej atmosfery dla danych:
fi=52st15; Z0=45st<-wartość obsewowana
t=-2.5stC P=758.5mmHg
Zad3(nie mam dokładnej treści)
Podanie odległości między Irkuckiem a Poznaniem dla
Irkuck: fi=52st33'N i lambda=104st24'E
Poznań: fi=52st24'N i lambda=16st55'E
W domu mieliśmy wyliczyć tą odległość w km wychodzi 6374.140km
zad4
Rozwiąż układ 3 równań
6.2x+3.7y-5.3z=8.4
1,5x-6,3y+4,7z=8,4
3,6x-3,6y-3,6z=4,7
Zad5(16.XI.2006)
Korzystając z wyników z poprzedniego zadania oblicz max h nad horyzontem geodezyjnym(obserwatora) sztucznego satelity ziemi obserwowanego w tym miejscu. Satelita porusza się po orbicie kołowej geocentrycznej o r=8798km nachylonej do płaszczyzny równika ziemskiego pod katem 18st36'
Dane z poprzedniego zadania:ro:=6365km FI':=51.81
zad6(12.10.06)
a moment czasu 1985 01 07 14.42.00 UT planowana jest obserwacja obiektu 3D273 alfa=12h28.3m delta=2st8' teleskop ma montaż horyzontalny o współrzędnych lambda:=-83st31'fl:=40st15'
Wyznacz nastawę tego teleskopu
W domu trza było wyliczyć coś takiego
[x] [x1]
[y]=r(-180)q(90-fi)[y1]
[z](A,h) [z1](H,delta)
odp:h:=26,0564 A:=291,0779(st->?)
zad7
Mając dany trujkąt sferyczny oblicz a , b , c
Dane: A:=126st55'18'' B:=29st59'14'' C:=33st40'50''
Wzory z których korzystamy cosB=-cosAcosC+sinAsinCcosb
cosC=-cosAcosB+sinAsinBsinc
cosA=-cosBcosC+sinBsinCsina
odpowiedzi: a:=64st20'23'' b:=34st18'00'' c:=34st42'05''
zad8(23XI06)
Udowodnij że kąt wertykału v dany jest formułą
tgv=(e^2sin2FI)/2(1-e^2sin^2Fi) dodatkowe info e:=c/a c:=sqrt(a^2-b^2)
To chyba byłoby na tyle - jeśli ktoś widzi błedy to proszę pisać:
|
|
| Powrót do góry » |
|
 |
|
|
 |
Gość
|
| Wysłany: Nie 13:23, 10 Gru 2006 Temat postu: |
|
|
ODP DLA ZAD 4:
rownanie:={6.2*x+3.7*y-5.3*z=8.4, 1.5*x-6.3*y+4.7*z=8.4,
3.6*x-3.6*y-3.6*z=4.7};
niewiadome:={x,y,z};
solve(rownanie,niewiadome):
|
|
| Powrót do góry » |
|
|
Gość
|
| Wysłany: Nie 13:26, 10 Gru 2006 Temat postu: |
|
|
ZAD 6(mam nadzieje że dobre wklejam  ) :
print("Nastaw teleskopu");
rad:=PI/180:
rad1:=PI/12:
l:=(12+(28.3/60)):
d:=(2+(8/60)):
la:=(-83-(31/60)):
fi:=(40+(15/60)):
0ut:=(7+6/60+1/3600):
GS:=1.0027379094:
1ut:=14+42/60:
CGM0:=0ut+la;
CGM1:=CGM0+1ut*GS;
|
|
| Powrót do góry » |
|
|
Gość
|
| Wysłany: Nie 13:26, 10 Gru 2006 Temat postu: |
|
|
ZAd7
print("trojkat sferyczny");
RAD:=PI/180:
A:=(126+55/60+18/3600)*RAD:
B:=(29+59/60+14/3600)*RAD:
C:=(33+40/60+50/3600)*RAD:
ca:=(cos(A)+cos(B)*cos(C))/(sin(B)*sin(C)):
cb:=(cos(B)+cos(A)*cos(C))/(sin(A)*sin(C)):
cc:=(cos(C)+cos(B)*cos(A))/(sin(B)*sin(A)):
a:=(arccos(ca)):
a1:=a/RAD:
float(%);
b:=(arccos(cb)):
b1:=b/RAD:
float(%);
c:=(arccos(cc)):
c1:=c/RAD:
float(%);
akuku:
|
|
| Powrót do góry » |
|
|
Gość
|
| Wysłany: Nie 13:27, 10 Gru 2006 Temat postu: |
|
|
zad2:
rad:=3.14/180:
T:=-2.5:
z0:=45*rad:
P:=758.5:
R:=rad*(60+4/60)*((P/760)/(1+T/273))*tan(z0);
Akuku:
|
|
| Powrót do góry » |
|
|
Gość
|
| Wysłany: Wto 22:56, 19 Gru 2006 Temat postu: |
|
|
No więc Tak Tutaj są wszystkie zadania - 2 które zrobiła nasza dobra koleżanka Asia(wszyscy niech jej ładnie podziękują!!) umieszczam je w jej kolejności bo tak były zadawane przy każdym jest mała zmianka czego się tyczy więc pewnie się połapiecie...tzn mam taką nadzieję..jak nie to piszcie do mnie jak będę to pomoge...
zad1 Równanie
rownanie:={ 6.2*x+3.7*y-5.3*z=8.4 , 1.5*x-6.3*y+4.7*z=8.4 ,
3.6*x-3.6*y-3.6*z=4.7 }:
niewiadome:={x,y,z}:
solve(rownanie,niewiadome);
zad2 Trojkąt sferyczny
rad:=3.14/180:
A:=rad*(126+55/60+18/3600):
B:=rad*(29+59/60+14/3600):
C:=rad*(33+40/60+50/3600):
cosa:=(cos(A)+cos(B)*cos(C))/(sin(B)*sin(C)):
cosb:=(cos(B)+cos(A)*cos(C))/(sin(A)*sin(C)):
cosc:=(cos(C)+cos(A)*cos(B))/(sin(A)*sin(B)):
a:=arccos(cosa)/rad:
print("a",a):
b:=arccos(cosb)/rad:
print("b",b):
c:=arccos(cosc)/rad:
print("c",c):
zad3-Nastawa Teleskopu
rad:=3.14/180:
radgodz:=3.14/12:
alfa:=radgodz*(12+28.3/60):
delta:=rad*(2+8/60):
lambda:=-rad*(83+31/60):
fi:=rad*(40+15/60):
T1:=radgodz*(14+42/60):
T2:=radgodz*(7+6/60+1/3600):
S1:=T2+lambda:
M:=1.0027379094:
S:=S1+T1*M:
H:=S-alfa:
x1:=cos(H)*cos(delta):
y1:=sin(H)*cos(delta):
z1:=sin(delta):
q:=matrix(3,3,[[1,0,0],[0,cos(90*rad-fi),sin(90*rad-fi)],[0,-sin(90*rad-fi),cos(90*rad-fi)]]):
r:=matrix(3,3,[[1,0,0],[0,cos(-180*rad),sin(-180*rad)],[0,-sin(-180*rad),cos(-180*rad)]]):
P:=matrix(3,1,[[x1],[y1],[z1]]):
X:=r*q*P:
hrad:=arcsin(X[3]):
h:=hrad/rad;
Arad:=arctan(X[2]/X[1]):
A:=Arad/rad;
zad4 Refrakcja model płaskiej atmosfery
rad:=3.14/180:
T:=-2.5:
z0:=45*rad:
P:=758.5:
R:=rad*(60.4/3600)*((P/760)/(1+T/273))*tan(z0);
Akuku:
zad5 To jest to Coś z Pułkowem ten kto je ma mógłby je umieścić...
zad6 Max wysokość satelity nad horyzontem
rad:=3.14/180:
fiprim:=51.81300355*rad:
i:=rad*(18+36/60):
r:=8798:
p:=6364.912801:
OS:=sqrt(p*p+r*r-2*p*r*cos(fiprim-i)):
print("OS", OS):
cosx:=(p*p+OS*OS-r*r)/(2*p*OS):
x:=arccos(cosx):
H1:=x-90*rad:
H:=H1/rad;
akuku:
zad 8 Odległośc r' dla poszczególnych maximów
DIGITS:=20:
rad:=3.14/180:
lambda:=270*rad:
b:=45*rad:
AU:=1/206265:
r:=1/(1/3600*rad):
print("lambda slońca = 0"):
rprim0:=sqrt(AU*AU+r*r);
float(rprim0):
print("lambda slońca = 180"):
rprim180:=sqrt(AU*AU+r*r);
float(rprim180):
print("lambda slońca = 270"):
rprim270:=sqrt(r*r+AU*AU-2*AU*r*cos(180*rad-b));
float(rprim270):
print("lambda slońca = 90"):
rprim90:=sqrt(r*r+AU*AU-2*r*AU*cos(b));
float(rprim90):
No to mam nadzieję że się cieszycie I PAMIĘTAJCIE PODZIĘKUJCIE ASI !!
|
|
| Powrót do góry » |
|
|
Ucja
Dołączył: 20 Gru 2005
Posty: 41
Przeczytał: 0 tematów
Skąd: Poznań
|
| Wysłany: Wto 23:48, 19 Gru 2006 Temat postu: |
|
|
to ja zupełnie nie na temat zamieszczę pełne rozwiązanie OSTATNIEGO zadania domowego
| Kod: |
alf:= 2 + (15 + 54.6/60)/60:
delt:= 89 + (11 + 39/60)/60:
rad:= PI/180:
dd := x -> ( 0.6406161 * x + 0.0000839 * x^2 + 0.0000050 * x^3) * rad:
zz := x -> ( 0.6406161 * x + 0.0003041 * x^2 + 0.0000051 * x^3) * rad:
ff := x -> ( 0.5567530 * x - 0.0001185 * x^2 - 0.0000116 * x^3) * rad:
p := x -> matrix(3,3, [ [ 1, 0, 0],
[ 0, cos(x), sin(x)],
[ 0, -sin(x), cos(x)]]):
q := x -> matrix(3,3, [ [ cos(x), 0, -sin(x)],
[ 0, 1, 0],
[ sin(x), 0, cos(x)]]):
r := x -> matrix(3,3, [ [ cos(x), sin(x), 0],
[ -sin(x), cos(x), 0],
[ 0, 0, 1]]):
Alf := alf * rad * 15:
Delt := delt * rad:
x:= 1984.5: //epoka poczatkowa
T1 := (x - 2000) /100:
z1:= zz(T1):
f1:= ff(T1):
d1:= dd(T1):
P:= matrix(3,3):
P:= r( -z1 ) * q( f1 ) * r( -d1 ):
s0 := matrix(3,1, [ cos (Delt) * cos (Alf) , cos(Delt) * sin(Alf),sin(Delt)]):
sJ2000 := matrix(3,1):
sJ2000 := ( P^(-1) ) * s0 :
x:= 2005.5: //epoka docelowa
T2 := (x - 2000) /100:
P := r( -1*zz( T2 ) ) * q( ff( T2 ) ) * r( -1*dd(T2) ):
sk := matrix(3,1):
sk := P * sJ2000:
float( (arctan(sk[2]/sk[1])) / (rad*15) );
float(arcsin(sk[3]) / rad );
JD1984 := 2445883.625:
JD1985 := 2446248.875:
JD2005 := 2453553.875:
T:= (JD2005 - 2451545)/36525: //JDepoki docelowej
l := 134 + (57 + 46.733/60)/60 + (1325*360 + 198 + (52 + 2.633/60)/60 )*T + 31.31/3600 * T^2 + 0.064/3600 * T^3:
lprim := 357 + (31 + 39.804/66)/60 + ( 99*360 + 359 + (3 + 1.224/60)/60 ) * T - 0.577/3600 * T^2 -
0.012/3600 * T^3:
eF := 93 + (16 + 18.877/60)/60 + (1342*360 + 82 + (1 + 3.137/60)/60 )* T - 13.257/3600 * T^2 + 0.011/3600 * T^3:
De := 297 + (51 + 1.307/60)/60 + (1236*360 + 307 + (6 + 41.328/60)/60 )*T - 6.891/3600 * T^2 + 0.019 * T^3:
Om := 125 + (2 + 40.28/60)/60 - (5*360 + (134 + 8 + 10.539/60)/60 )*T +
7.455/3600 * T^2 + 0.008 * T^3:
a := array(1..10, [ 0, 0, -2, 2, -2, 1, 0, 2, 0, 0]):
b := array(1..10, [ 0, 0, 0, 0, 0, -1, -2, 0, 0, 1]):
c := array(1..10, [ 0, 0, 2, -2, 2, 0, 2, -2, 2, 0]):
d := array(1..10, [ 0, 0, 0, 0, 0, -1, -2, 0, -2, 0]):
e := array(1..10, [ 1, 2, 1, 0, 2, 0, 1, 1, 2, 0]):
S := array(1..10, [ -171996-174.2*T,
2062 +0.2*T,
46 +0.2*T,
11 +0.2*T,
-3 +0.2*T,
-3 +0.2*T,
-2 +0.2*T,
1 +0.2*T,
-13187 -1.6*T,
1426 -3.4*T]):
C := array(1..10, [ 92025 +8.9*T,
-895 +0.5*T,
-24 +0.5*T,
0 +0.5*T,
1 +0.5*T,
0 +0.5*T,
1 +0.5*T,
0 +0.5*T,
5736 -3.1*T,
54 -0.1*T]):
dnut := 0:
deps := 0:
for i from 1 to 10 do
A:= (a[i]* l + b[i] * lprim + c[i] * eF + d[i] * De + e[i] * Om)*rad:
dnut:= dnut + S[i] * sin(A) * 0.0001 /3600:
deps:= deps + C[i] * cos(A) * 0.0001 /3600:
end_for:
eps0:= 23 + (26 + 21.448/60)/60 - 46.8150/3600 * T + 0.00059/3600 * T^2 +
0.001813 * T^3:
eps:= eps0 + deps:
N := matrix(3,3):
N := p(-eps * rad) * r(-dnut * rad ) * p(eps0 * rad):
skp := N * sk:
float( (arctan(skp[2]/skp[1])) / (rad*15) );
float(arcsin(skp[3]) / rad ):
|
wynik działania programu wygląda tak:
| Kod: |
2.63532419
89.28794485
2.623807563
89.2887507
|
pierwsze dwie liczby to rektascencja i deklinacja gwiazdy obliczone metodą "dokładną" kolejne dwie to rektascencja i deklinacja liczone z uwzględnieniem wpływów nutacyjnych
Wynik wygląda prawdopodobnie, więc z dokładnością do literówek w danych winno być dobrze
Post został pochwalony 0 razy
|
|
| Powrót do góry » |
|
|
Ucja
Dołączył: 20 Gru 2005
Posty: 41
Przeczytał: 0 tematów
Skąd: Poznań
|
| Wysłany: Pon 19:40, 22 Sty 2007 Temat postu: |
|
|
Czy ktoś spisał co jest do zrobienia na czwartek?
Post został pochwalony 0 razy
|
|
| Powrót do góry » |
|
|
Gość
|
| Wysłany: Pon 20:16, 22 Sty 2007 Temat postu: |
|
|
Zad z 4.01.2007r
Rok zwrotnikowy=365,2422
Rok gwiazdowy=365,2564
Rok anamolistyczny=365,2596
a) oszacuje tempo i znam precesji osi obrotu ziemi <- na zajęciach
b) Oszacuj tempo i znak ruchu perychelium orbity ziemi <- dom
|
|
| Powrót do góry » |
|
|
Gość
|
| Wysłany: Pon 20:22, 22 Sty 2007 Temat postu: |
|
|
zad z 7.12.2006
Gwiazda Polarna o współrzędnych alfa=2h 15m 54,46s delta=89st 11,39' na epokę E=1984,5
a)oblicz w sposób przybliżony jej współrzędne na epokę 1985 i 2085
b) oblicz w sposób ścisły - > dom
Potem jedno z tych co ty podałaś..
Jedno zamieszcze jak odzyskam baterie do aparatu..dajmy dziś wieczorem albo jutro jak wróce z uczelni..jest za długie a nie mam sił tego pisać...
a do 1 zaraz zamieszcze rozwiązanie..No dobra tzn jak nawiąże połączenie z hoth'em
|
|
| Powrót do góry » |
|
|
Gość
|
| Wysłany: Pon 20:45, 22 Sty 2007 Temat postu: |
|
|
A to jest jedno z zadań - jednak nie pamiętam do niego treści wiem że siedzieliśmy nad nim jakoś 11.01.2007
[link widoczny dla zalogowanych]
|
|
| Powrót do góry » |
|
|
sekhet
Astronom
Dołączył: 20 Gru 2005
Posty: 8
Przeczytał: 0 tematów
Skąd: Poznań
|
| Wysłany: Pon 23:57, 22 Sty 2007 Temat postu: |
|
|
joł pipol
Jutro wrzucę to co mam/////
pozdro
3oL
Post został pochwalony 0 razy
|
|
| Powrót do góry » |
|
|
|
|
