Реклама Google — средство выживания форумов :)
REM /характеристич. скорости для межпланетных перелетов LPO - LPO/ R#(1) = .3871: R#(2) = .72333: R#(3) = 1: R#(4) = 1.52369: R#(5) = 5.2028 R#(6) = 9.53884: R#(7) = 19.18185: R#(8) = 30.058: R#(9) = 39.52 V#(1) = 47.87: V#(2) = 35.02: V#(3) = 29.79: V#(4) = 24.13: V#(5) = 13.06 V#(6) = 9.65: V#(7) = 6.8: V#(8) = 5.43: V#(9) = 4.7 Mp#(1) = 3.3E+11: Mp#(2) = 4.87E+12: Mp#(3) = 5.973E+12: Mp#(4) = 6.42E+11 Mp#(5) = 1.9E+15: Mp#(6) = 5.68E+14: Mp#(7) = 8.7E+13 Mp#(8) = 1.03E+14: Mp#(9) = 4.7E+11 Re#(1) = 4878: Re#(2) = 12104: Re#(3) = 12756: Re#(4) = 6785: Re#(5) = 142796 Re#(6) = 120000: Re#(7) = 52400: Re#(8) = 48600: Re#(9) = 4900 Rc#(1) = 111620: Rc#(2) = 615800: Rc#(3) = 924550: Rc#(4) = 577650 Rc#(5) = 48176000: Rc#(6) = 54529000: Rc#(7) = 51741000: Rc#(8) = 86746000 Rc#(9) = 33923000: GE# = 66.742# CLS : i$(1) = "Афелий или радиус начальной орбиты астероида (a.е.)__ " i$(2) = "Диаметр астероида (км) " i$(3) = "Высота круговой LPO у планеты назначения, км " i$(4) = "Планета отправления ... Земля - 3, Марс - 4, ... __ " i$(5) = "Планета назначения Меркурий - 1, Венера - 2, ... __ " in: u = u + 1: PRINT u MOD 10; i$(u MOD 10); : INPUT i#(u MOD 10) IF u < 5 THEN GOTO in nn% = (ABS(i#(5) - 1) MOD 9) + 1: COLOR 14 exn% = (ABS(i#(4) - 1) MOD 9) + 1 Rp# = R#(nn%): Vor# = V#(nn%) Vcaf# = V#(3) * SQR(1 / R#(exn%)) Vaff# = Vcaf# * SQR(2 - 2 * R#(exn%) / (R#(exn%) + Rp#)) Tg# = SQR(((R#(exn%) + Rp#) / 2) ^ 3) / 2 Vex# = Vaff# - Vcaf#: Vcd# = Vex# n% = exn%: HLPO = 200: GOSUB dd dVc# = Vimp# - Vc#: dVp# = (SQR(2) - 1) * Vc#: dVg# = dVcd# PRINT "Скорость выхода из С.Д. планеты №"; exn%; ", км/c: "; USING "+###.###"; Vex#; PRINT " _(Vsol="; USING "###.###"; Vaff#; : PRINT ")" PRINT "Скорость старта с круговой (h="; CINT(Hat# / 2); "км):_____ "; USING "###.###"; Vimp#; PRINT " _(dVcirc="; USING "###.###"; dVc#; : PRINT ")" PRINT "Приращение с круговой до параболической:__ "; USING "###.###"; dVp# PRINT "dV до гиперболической (для h стартовой LPO): "; USING "###.###"; dVg# PRINT "*": PRINT "Время перелета по гомановской, лет, дней:"; FIX(Tg#); "a"; PRINT CINT(365.24 * (Tg# - FIX(Tg#))); "d": PRINT "*" kvp# = SQR(2 - 2 * Rp# / (R#(exn%) + Rp#)) Vcd# = Vor# * (kvp# - 1#) Vin# = Vor# * kvp# PRINT "Скорость входа в С.Д. планеты №"; nn%; ", км/c ___"; USING "+###.###"; Vcd#; PRINT " _(Vsol="; USING "###.###"; Vin#; : PRINT ")" n% = nn%: HLPO = i#(3): GOSUB dd Vimp0# = Vimp#: dVcd0# = dVcd# PRINT "Скорость входа в атмосферу планеты __ Vgip "; USING "###.###"; Vimp#; Vcd# = 0: GOSUB dd dVgip# = Vimp0# - Vimp# 'Vper(вход)-Vper(орб а=СД) dVelp# = Vimp# - Vc# PRINT "Скорость в перицентре переходной орбиты __ "; USING "###.###"; Vimp#; PRINT " _(dVgip="; USING "###.###"; dVgip#; : PRINT ")" PRINT "XC перехода на круговую орбиту _______ dVelp "; USING "###.###"; dVelp# IVqc# = dVp# + dVgip#: IVqe# = dVg# + dVelp# PRINT "Характеристическая скорость для ЖРД ________ "; USING "###.###"; IVqc# PRINT "Характеристическая скорость для ЭРД ________ "; USING "###.###"; IVqe#; PRINT " _(всего="; USING "###.###"; IVqc# + IVqe#; : PRINT ")" COLOR 7 INPUT " Выход 0, изменить входной параметр №(1-5) "; u u = (u MOD 6) + 9 IF u = 9 THEN SYSTEM ELSE GOTO in END dd: Hat# = HLPO * 2 Rec# = 500 * (Re#(n%) + Hat#) '200 -высота LPO Vc# = SQR(GE# * Mp#(n%) / Rec#) / 1000 'круговая w# = SQR((Vcd# / Vc#) ^ 2 - (Re#(n%) + Hat#) / Rc#(n%) + 2) Vimp# = w# * Vc#: dVt# = Vimp# - Vc# dVcd# = Vimp# - SQR(2) * Vc# RETURN