Моделирование движения системы солнце-земля-луна по законам кеплера

@K95 · Регистрация: 13.05.2016

Студворк — интернет-сервис помощи студентам

Помогите пожалуйста не могу никак написать программу, необходимо смоделировать движение солнца-земли и луны по законам кеплера, написала программу, но луна не вращается вокруг земли, не знаю что делать, знаю в программе ошибка, но не могу понять где

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
program dvizhellips;
 
uses
  GraphABC;
 
var
  GM, vx, vy, r, t, tk, dt, rx, Sx, Sy, Hu, Hv, ax, ay,  x, y, xl, yl, R1, vxl, vyl, axl, ayl,xlae,ylae: real;
 
procedure Init;
begin
  GM := 4 * pi * pi;
  vx := 0;
  vy := 2 * pi;
  vxl := 0;
  vyl := pi/2; 
  t := 0;
  tk := 30;
  rx :=6;
  dt := 0.0001;
  Sx := Window.Width / rx;
  Sy := Sx;
  Hu := Window.Width / 2;
  Hv := Window.Height / 2;
  x := 1;
  y := 0;
  xl := 1.1;
  yl := 0;
  line(0, round(Hv), Window.Width, round(Hv));
  line(round(Hu), 0, round(Hu), Window.Height);
end;
 
procedure OutRes(x, y: real; c: color);
begin
  var u := round(x * Sx + Hu);
  var v := round(-y * Sy + Hv);
  PutPixel(u, v,  c);
end;
 
procedure Eiler;
begin
  while (t < tk ) do
  begin
    r := sqrt(x * x + y * y );
    writeln(r);
    ax := -(GM / (r * r * r)) * x;
    ay := -(GM / (r * r * r)) * y;
    vx := vx + ax * dt;
    vy := vy + ay * dt;
    x := x + vx * dt;
    y := y + vy * dt;
    
    R1 := sqrt(xl * xl + yl * yl);
    writeln(R1);
    axl := -((GM) / ((R1 + r) * (R1 + r))) * xl;
    ayl := -((GM) / ((R1 + r) * (R1 + r))) * yl;
    vxl := vxl + axl * dt;
    vyl := vyl + ayl * dt;
    xl := xl + vxl * dt;
    yl := yl + vyl * dt;
   { xlae:=(R1*r)/xl;
    ylae:=(R1*r)/xl;
    xlae:=xlae+x;
    ylae:=ylae+y;}
    t := t + dt;
    OutRes(x, y, clblue);
    OutRes(xl, yl, clgreen);
  end;
end;
 
begin
  Init;
  Eiler;
end.

@K95 · 08.06.2016, 11:52 **[ТС]**

Помогите пожалуйста не могу никак написать программу, необходимо смоделировать движение солнца-земли и луны по законам кеплера, написала программу, но луна не вращается вокруг земли, не знаю что делать, знаю в программе ошибка, но не могу понять где

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
program dvizhellips;
 
uses
  GraphABC;
 
var
  GM, vx, vy, r, t, tk, dt, rx, Sx, Sy, Hu, Hv, ax, ay,  x, y, xl, yl, R1, vxl, vyl, axl, ayl,xlae,ylae: real;
 
procedure Init;
begin
  GM := 4 * pi * pi;
  vx := 0;
  vy := 2 * pi;
  vxl := 0;
  vyl := pi/2; 
  t := 0;
  tk := 30;
  rx :=6;
  dt := 0.0001;
  Sx := Window.Width / rx;
  Sy := Sx;
  Hu := Window.Width / 2;
  Hv := Window.Height / 2;
  x := 1;
  y := 0;
  xl := 1.1;
  yl := 0;
  line(0, round(Hv), Window.Width, round(Hv));
  line(round(Hu), 0, round(Hu), Window.Height);
end;
 
procedure OutRes(x, y: real; c: color);
begin
  var u := round(x * Sx + Hu);
  var v := round(-y * Sy + Hv);
  PutPixel(u, v,  c);
end;
 
procedure Eiler;
begin
  while (t < tk ) do
  begin
    r := sqrt(x * x + y * y );
    writeln(r);
    ax := -(GM / (r * r * r)) * x;
    ay := -(GM / (r * r * r)) * y;
    vx := vx + ax * dt;
    vy := vy + ay * dt;
    x := x + vx * dt;
    y := y + vy * dt;
    
    R1 := sqrt(xl * xl + yl * yl);
    writeln(R1);
    axl := -((GM) / ((R1 + r) * (R1 + r))) * xl;
    ayl := -((GM) / ((R1 + r) * (R1 + r))) * yl;
    vxl := vxl + axl * dt;
    vyl := vyl + ayl * dt;
    xl := xl + vxl * dt;
    yl := yl + vyl * dt;
   { xlae:=(R1*r)/xl;
    ylae:=(R1*r)/xl;
    xlae:=xlae+x;
    ylae:=ylae+y;}
    t := t + dt;
    OutRes(x, y, clblue);
    OutRes(xl, yl, clgreen);
  end;
end;
 
begin
  Init;
  Eiler;
end.

@JuriiMW · 08.06.2016, 12:15

Вы думаете, что мы сейчас прямо так вспомним законы Кеплера и „без базара“ поймём что, почему и зачем у вас описано:

Сообщение от K95

GM, vx, vy, r, t, tk, dt, rx, Sx, Sy, Hu, Hv, ax, ay, x, y, xl, yl, R1, vxl, vyl, axl, ayl,xlae,ylae: real;

… Уже не говоря о способе реализации в стиле „a la Turbo Pascal — мы из восьмидесятых“… ;–(

@K95 · 10.06.2016, 18:52 **[ТС]**

GM-гравитационная постоянная и масса
vx, vy-скорость земли по x и y
r-радиус земли
ax, ay-ускорение земли по x и y
x, y-координаты земли
rx, Sx, Sy, Hu, Hv - для рисования системы координат
xl, yl-координаты луны
R1-радиус луны
vxl, vyl-скорость луны по x и y
axl, ayl-ускорение луны по x и y
xlae,ylae-перевод лунных в астрономические единицы

@JuriiMW · 14.06.2016, 04:45

А теперь теоретическую составляющую: формулы и константы…
… или, хотя бы, ссылку ;–)

Новые блоги и статьи Все статьи Все блоги /
Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Рецензия / Мнение/ Перевод https:/ / **********/ gallery/ thinkpad-x220-tablet-porn-gzoEAjs . . .	PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .
От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .

@K95 1 / 1 / 0 Регистрация: 13.05.2016 Сообщений: 27
	10.06.2016, 18:52 [ТС]
	GM-гравитационная постоянная и масса vx, vy-скорость земли по x и y r-радиус земли ax, ay-ускорение земли по x и y x, y-координаты земли rx, Sx, Sy, Hu, Hv - для рисования системы координат xl, yl-координаты луны R1-радиус луны vxl, vyl-скорость луны по x и y axl, ayl-ускорение луны по x и y xlae,ylae-перевод лунных в астрономические единицы 0

@JuriiMW 5095 / 2661 / 2355 Регистрация: 10.12.2014 Сообщений: 10,059
	14.06.2016, 04:45
	А теперь теоретическую составляющую: формулы и константы… … или, хотя бы, ссылку ;–) 0