pi = 3.14
ws = 23.45 * pi / 180 '弧度，黃赤交角
A = 5.15 * pi / 180 '白黃交角，弧度值
y = Val(Text1.Text) '計算日年份
D = Val(Text2.Text) '從年初至該日期的總天數
l = (y - 2000) / 4 + 1 'J2000系統
l = Int(l) 'J2000系統，從2000年首至計算年年首的閏年數
t = ((y - 2000) * 365 + D + l - 0.5) / 36525 '從2000年1月12時起算的儒略世紀
s = (218.31643 + 481267.88128 * t + 0.00161 * t ^ 2 + 0.000005 * t ^ 3) / 360
x = Int(s)
s = (s - x) * 360 * pi / 180 '月球平均經度，弧度值
h = (280.46607 + 36000.7698 * t + 0.0003 * t ^ 2) / 360
X1 = Int(h)
h = (h - X1) * 360 * pi / 180 '從春分點起算的太陽平均經度
P = (83.35345 + 4069.01388 * t - 0.01031 * t ^ 2 - 0.00001 * t ^ 3) / 360
X2 = Int(P)
P = (P - X2) * 360 * pi / 180 '月球近地點平均經度
ps = (282.93835 + 1.71946 * t + 0.00046 * t ^ 2 + 0.000003 * t ^ 3) / 360
x3 = Int(ps)
ps = (ps - x3) * 360 * pi / 180 '太陽近地點的平均經度
N = (125.04452 - 1934.13626 * t + 0.00207 * t ^ 2 + 0.000002 * t ^ 3) / 360
x4 = Int(N)
N = (N - x4) * 360 * pi / 180 '弧度值,月球升交點的平均經度
COSI = Cos(ws) * Cos(A) - Sin(ws) * Sin(A) * Cos(N)
sinI = Sqr(1 - (COSI) ^ 2) '白赤交角
tanI = sinI / COSI
I = Atn(tanI) '結果為弧度值
sinv = Sin(N) * Sin(A) / Sin(I)
cosv = Sqr(1 - (sinv) ^ 2)
v = Atn(sinv / cosv)
'結果為弧度值，γA
Sinz1 = Sin(N) * Sin(ws) / Sin(I) '有些問題,AΩ
Cosz1 = Cos(N) * Cos(v) + Sin(N) * Sin(v) * Cos(ws)
z1 = Atn(Sinz1 / Cosz1) '結果為弧度值AΩ
z2 = z1 + s - N
'弧度值，月球在白道上從A點起算的平均經度，σ
e = 0.0549
'月球軌道偏心率
m = 0.074804
'太陽與月球運行速度的比值
u = z2 + 2 * e * Sin(s - P) + (5 / 4) * e ^ 2 * Sin(2 * s - 2 * P) + m * e * (15 / 4 + 263 / 16 * m) * Sin(s - 2 * h + P) + m ^ 2 * (11 / 8 + 59 / 12 * m + 75 / 16 * (e ^ 2 / m)) * Sin(2 * s - 2 * h) + 17 / 8 * m ^ 2 * e * Sin(3 * s - 2 * h - P) + 77 / 16 * m ^ 2 * e * Sin(2 * s - h + ps) '弧度值,從輔助春分點起算的月球在白道上的真經度
sinc = Sin(l) * Sin(I)
'δ月球赤緯
cosc = Sqr(1 - (sinc) ^ 2)
c = Atn(sinc / cosc)
'弧度值
Label5.Caption = c * 180 / pi
z3 = Val(Text3.Text)
'度數φ，地理緯度
b = 0 To 23
'太陰時
總公式：
N北方向，e東西方向
fn = 1 / 2 * Sin(2 * z3 * pi / 180) * (3 * (sinc) ^ 2 - 1) - Cos(2 * z3 * pi / 180) * Sin(2 * c) * Cos(15 * b * pi / 180) - 1 / 2 * Sin(2 * z3 * pi / 180) * (cosc) ^ 2 * Cos(30 * b * pi / 180)
fe = -Sin(z3 * pi / 180) * Sin(2 * c) * Sin(15 * b * pi / 180) + Cos(z3 * pi / 180) * (cosc) ^ 2 * Sin(30 * b * pi / 180)
這是代碼 幫我寫個程式 要有圖示和曲線 填充繪制text和按鈕
uj5u.com熱心網友回復：
圖表和曲線 - - 大神幫幫忙轉載請註明出處，本文鏈接：https://www.uj5u.com/gongcheng/283813.html
標籤：VB基礎類
上一篇：DirectShow的MJPEG Decompressor Filter如何修改輸出引腳的MediaType？
下一篇：Duilib繪制一個圖