Interessante Theorie, aber für die Erde liefert sie, sofern ich mich jetzt nicht selbst verrechnet habe, einen um [vm_earth/Differenz] = 366.31983 fach zu niedrigen Wert:
Code:
syms kg m sek positive
ha_earth = 1.4960e+011*m;
r_earth = 6363000*m;
v_earth = 29780*m/sek;
vm_earth = 463.9973*m/sek;
vw_earth = vm_earth/r_earth;
hw_earth = v_earth/ha_earth;
Code:
v_innen = (ha_earth-r_earth/2)*hw_earth;
v_aussen = (ha_earth+r_earth/2)*hw_earth;
Code:
Differenz = vpa(v_aussen-v_innen)
Differenz = 1.2666 m/sek
Hier nochmal mit allen Planeten, Matrix zeilenweise von Merkur bis Pluto:
Code:
ha_solar =
1.0e+012 m *
0.0579 0.1082 0.1496
0.2279 0.7787 1.4268
2.8724 4.4950 5.9063
r_solar =
1.0 m *
2439700 6051800 6363000
3392845 70701500 59250000
25460500 24693000 1151000
v_solar =
1.0 m/sek *
47870 35030 29780
24130 13070 9690
6810 5430 4720
vw_solar =
1.0e-003/sek *
0.0012 -0.0003 0.0729
0.0709 0.1759 0.1619
-0.1012 0.1093 -0.0114
hw_solar = v_solar./ha_solar;
Code:
v_innen = (ha_solar-r_solar./2).*hw_solar;
v_aussen = (ha_solar+r_solar./2).*hw_solar;
Faktor = vpa((v_aussen-v_innen)./vm_solar)
Code:
Faktor =
[ 0.666854, -1.08231, 0.00273]
[ 0.001494, 0.00009, 0.00004]
[ -0.00002, 0.00001, -0.00007]
würde ich statt der Winkelgeschwindigkeit der Bahn [hw_] die Winkelgeschwindigkeit der Planeten [vw_] einsetzen, würde die Rechnung übrigens stimmen und beim Faktor
überall 1 herauskommen...
Variablen.txt