Программа для расчета дистанции

@Mindhunter · Регистрация: 13.12.2012

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

Здравствуйте. Проблема вот в чем, мы вводим данный в градусах (направление) и морских милях (дстанция) и узлах (скорость). А в программе эти все данные считаются как радианы, метры и м/с - по системе Си. Подаскажите пожалуйста, как можно прописать код, чтоб перед расчетом сама перевела все параметры в систему Си. Я попробовал дописать, для перевода, но походу не получилось. И не могу найти ошибку в коде, может кто заметит, считает не правильно.

C#
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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;
 
namespace Laba4
{
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
        }
 
        private void button1_Click(object sender, EventArgs e)
        {
            try
            {
                double E = Convert.ToDouble(textE.Text);
                double T2 = Convert.ToDouble(textT2.Text);
                double Vh = Convert.ToDouble(textVh.Text);
                double Kh = Convert.ToDouble(textKh.Text);
                double M = Convert.ToDouble(textM.Text);
                double D1 = Convert.ToDouble(textD1.Text);
                double P1 = Convert.ToDouble(textP1.Text);
                double D2 = Convert.ToDouble(textD2.Text);
                double P2 = Convert.ToDouble(textP2.Text);  
                double PI = 4.0 * Math.Atan(1);
                double PI1 = PI/180;
                double a = 1852.0;
                double b = 0.5144;
                 
                //1. Определение координат первой точки
                double x1, y1; 
                
                x1 = (D1*a) * (Math.Sin(P1*PI1));
                y1 = (D1 * a) * Math.Cos(P1 * PI1);
                labelX1.Text = "= "+ x1.ToString();
                labelY1.Text = "=" + y1.ToString();
 
                //2. Определение координат второй точки
                double x2, y2;
                x2 = (D2*a) * Math.Sin(P2 * PI1);
                y2 = (D2*a) * Math.Cos(P2 * PI1);
                
                labelX2.Text = "= "+ x2.ToString();
                labelY2.Text = "=" + y2.ToString();
 
                //3. Определение delta t
                double deltaT = T2/60;
                labelDeltaT.Text = "= "+ deltaT.ToString();
                
                //4. Составляющая вестора скорости Vн
                double deltaXvh = ((Vh*b) * deltaT * Math.Sin(Kh*PI1));
                double deltaYvh = ((Vh*b) * deltaT * Math.Cos(Kh*PI1));
                
                labelDeltaXvh.Text = "= "+ deltaXvh.ToString();
                labelDeltaYvh.Text = "= "+ deltaYvh.ToString();
 
                //5. Точка начала вектора Vн
                double Xh = x1 - deltaXvh;
                double Yh = y1 - deltaYvh;
                
                labelXh.Text = "= "+ Xh.ToString();
                labelYh.Text = "= "+ Yh.ToString();
 
                //6. Длина вектора
                double Vh_2 = Math.Sqrt((Math.Pow(deltaXvh, 2)) + (Math.Pow(deltaYvh, 2)));
                labelVh.Text = "= "+ Vh_2.ToString();
 
                //7. Составляющая вектора скорости цели Vц
                double deltaXc = x2 - Xh;
                double deltaYc = y2 - Yh;
                double Vc = Math.Sqrt((Math.Pow(deltaXc, 2)) + (Math.Pow(deltaYc, 2)));
                
                labelDeltaXc.Text = "= "+ deltaXc.ToString();
                labelDeltaYc.Text = "= "+ deltaYc.ToString();
                labelVc.Text = "= " + Vc.ToString(); 
 
                //8. Определение дополнительных параметров
                Double Y2h = y2 - Yh;
                Double X2h = x2 - Xh;
                double delta = Vh / Vc;
                
                labelY2h.Text = "= "+ Y2h.ToString();
                labelX2h.Text = "= "+ X2h.ToString();
                labelDelta.Text = "= " + delta.ToString();
 
 
                //9. Координаты точек касания Xk1 и Xk2 к годографу
                Double Xk1 = Xh + (Math.Pow(delta, 2))*X2h + delta*(Math.Sqrt((Math.Pow(delta, 2)) + (Math.Pow(X2h, 2)) + (Math.Pow(Y2h, 2)) - (Math.Pow(Vh, 2))));
                Double Xk2 = Xh + (Math.Pow(delta, 2)) * X2h - delta * (Math.Sqrt((Math.Pow(delta, 2)) + (Math.Pow(X2h, 2)) + (Math.Pow(Y2h, 2)) - (Math.Pow(Vh, 2))));
                
                labelXk1.Text = "= "+ Xk1.ToString();
                labelXk2.Text = "= "+ Xk2.ToString();
 
 
                //10. Вторая координата точки касания
                double Y11 = Yh + (Y2h / 2) + (Math.Sqrt(Math.Pow((Y2h / 2), 2) + (Xk1 - Xh) * (x2 - Xk1)));
                double Y12 = Yh + (Y2h / 2) - (Math.Sqrt(Math.Pow((Y2h / 2), 2) + (Xk1 - Xh) * (x2 - Xk1)));
                double Y21 = Yh + (Y2h / 2) + (Math.Sqrt(Math.Pow((Y2h / 2), 2) + (Xk2 - Xh) * (x2 - Xk2)));
                double Y22 = Yh + (Y2h / 2) - (Math.Sqrt(Math.Pow((Y2h / 2), 2) + (Xk2 - Xh) * (x2 - Xk2)));
                
                labelY11.Text = "= "+ Y11.ToString(); 
                labelY12.Text = "= "+ Y12.ToString();
                labelY21.Text = "= "+ Y21.ToString();
                labelY22.Text = "= "+ Y22.ToString();
 
                //11. Проверка полученного значения путем выполнения условия
                double Yk1=0;
                double Yk2=0;
                double E1 = Math.Abs(Math.Pow((Y11 - Yh), 2) + Math.Pow((Xk1 - Xh), 2) - Math.Pow(Vh, 2));
                if (E1 <= E)
                {
                    bool goodE1 = true;
                    labelXk1Y11.Text = "= " + goodE1.ToString();
                    Yk1 = Y11;
                }
                else 
                {
                    bool goodE1 = false;
                    labelXk1Y11.Text = "= " + goodE1.ToString();
                }
                /////////////////////////////////////////////
                
                double E2 = Math.Abs(Math.Pow((Y12 - Yh), 2) + Math.Pow((Xk1 - Xh), 2) - Math.Pow(Vh, 2));
                if (E2 <= E)
                {
                    bool goodE2 = true;
                    labelXk1Y12.Text = "= " + goodE2.ToString();
                    Yk1 = Y12;
                }
                else
                {
                    bool goodE2 = false;
                    labelXk1Y12.Text = "= " + goodE2.ToString();
                }
                
                /////////////////////////////////////////////
                double E3 = Math.Abs(Math.Pow((Y21 - Yh), 2) + Math.Pow((Xk2 - Xh), 2) - Math.Pow(Vh, 2));
                if (E3 <= E)
                {
                    bool goodE3 = true;
                    labelXk2Y21.Text = "= " + goodE3.ToString();
                    Yk2 = Y21;
                }
                else
                {
                    bool goodE3 = false;
                    labelXk2Y21.Text = "= " + goodE3.ToString();
                }
 
                /////////////////////////////////////////////
                double E4 = Math.Abs(Math.Pow((Y22 - Yh), 2) + Math.Pow((Xk2 - Xh), 2) - Math.Pow(Vh, 2));
                if (E4 <= E)
                {
                    bool goodE4 = true;
                    labelXk2Y22.Text = "= " + goodE4.ToString();
                    Yk2 = Y22;
                }
                else
                {
                    bool goodE4 = false;
                    labelXk2Y22.Text = "= " + goodE4.ToString();
                }
 
                //12. Вычисление курсов
                double deltaYkh1 = Yk1 - Yh;
                double deltaXkh1 = Xk1 - Xh;
 
                labelDeltaYkh.Text = "= " + deltaYkh1.ToString();
                labelDeltaXkh.Text = "= " + deltaXkh1.ToString();
 
                if (deltaXkh1 >= 0)
                {
                    deltaXkh1 = PI/2;
                }
                else
                {
                    deltaXkh1 = 3*2/PI;
                }
 
                double Kn1 = (Math.Atan((-deltaYkh1 / deltaXkh1)) + deltaXkh1);
                double K1 = Kn1/PI1;
 
                labelKh1.Text = "= " + K1.ToString();
 
                double deltaYkh2 = Yk2 - Yh; 
                double deltaXkh2 = Xk2 - Xh;
 
                labelDeltaYkh2.Text = "= " + deltaYkh2.ToString();
                labelDeltaXkh2.Text = "= " + deltaXkh2.ToString();
 
                if (deltaXkh2 >= 0)
                {
                    deltaXkh2 = PI/2;
                }
                else
                {
                    deltaXkh2 = 3*2-PI;
                }
                double Kn2 = (Math.Atan(-(deltaYkh2 / deltaXkh2)) + deltaXkh2);
                double K2 = (Kn2 / PI1);
 
                labelKh2.Text = "= " + K2.ToString();
 
                //13. Вычисляем координаты упреждающей точки
 
                double Xym = x2 + (M * (x2 - x1));
                double Yym = y2 + (M *(y2 - y1));
                labelXym.Text = "= " + Xym.ToString();
                labelYym.Text = "= " + Yym.ToString();
 
                //14. Вычисляем дистанцию кратчайщего сближения
 
                double Dne1 = ((- 1 /Vh) * (Yym * (Yk1 - Yh) + (Xym * (Xk1 - Xh))));
                double Dn1 = Dne1 / a;
                double Dne2 = ((- 1 /Vh) * (Yym * (Yk2 - Yh) + (Xym * (Xk2 - Xh))));
                double Dn2 = Dne2 / a;
                
                labelDn1.Text = "= " + Dn1.ToString();
                labelDn2.Text = "= " + Dn2.ToString();
 
 
                //15. Выбор максимальной дистанции
 
                if (Dn1 > Dn2)
                {
                   labelResultDn.Text = "= " + Dn1.ToString();
                   labelResultK.Text = "= " + K1.ToString();
                                   }
                else
                {
                    labelResultDn.Text = "= " + Dn2.ToString();
                    labelResultK.Text = "= " + K2.ToString();
                   
                }
                
            }
            catch (Exception ex)
            {
                MessageBox.Show(ex.Message);
            }
        }
 
        private void Form1_Load(object sender, EventArgs e)
        {
 
        }
    }
}

Алгоритм которым руководствуюсь:

Известные значение, которые вводит пользователь (оператор):
Т2 = 3 (но чтоб можно было и самому менять, рекомендуется) – время наблюдения.
Vн – скорость
Кн – курс
М – коэффициент (рекомендуется равен 5, но можно самому менять)
D1 – дистанция
П1 – направление на объект
D2 – дистанция
П2 – направление на объект
E – эпселент = 0,01

1. Определение координат первой точки
X1 = D1*sinП1
Y 1 = D1*cos П1
2. Определение координат второй точки
X2 = D2*sinП2
Н2 = D2*cosП2
3. Определение ∆t
∆t = (T2-T1)/60 = T2/60
4. Составляющая вестора скорости Vн
∆Xvн = Vн*∆t*sinKн
∆Yvн = Vн*∆t*cosKн
5. Точка начала вестора Vн
Xн = X1 - ∆Xvн
Yн = Y1 - ∆Yvн
6. Длина ветора
Vн = корень из ∆Xvн^2 + ∆Yvн^2
7. Составляющая вектора скорости цели Vц
∆Xц = X2-Xн
∆Yц = Y2 – Yн
Vц = корень из ∆Xц^2 + ∆Yц^2
8. Дополнительные параметры
Y2н = ∆Yц
X2н = ∆Xц
∂ = Vн/Vц
9. Координаты точек касания
Xк1 = Xн+(∂^2)*X2н + ∂*√(∂^2)+(X2y^2)+(Y2н^2)-(Vн^2)
Xк2 = Xн+(∂^2)*X2н -∂*√(∂^2)+(X2y^2)+(Y2н^2)-(Vн^2)
10. Вторая координата точек касания
Y11 = Yh + (Y2h/2) + √(Y2h/2)^2 + (Xk1 – Xh)*(X2 – Xk1);
Y12 = Yh + (Y2h/2) - √(Y2h/2)^2 + (Xk1 – Xh)*(X2 – Xk1);
Y21 = Yh + (Y2h/2) + √(Y2h/2)^2 + (Xk2 – Xh)*(X2 – Xk2);
Y22 = Yh + (Y2h/2) - √(Y2h/2)^2 + (Xk2 – Xh)*(X2 – Xk2);
11. Проверка полученного значения путем выполнения условия
(Xk1; Y11), (Xk1; Y12), (Xk2; Y21), (Xk2; Y22) – по условие отбрасываем значения. Оставляем только подходящие по условию.
│((Y11 – Yh)^2) + ((Xk1 – Xh)^2) – (Vh^2)│≤ E;
│((Y12 – Yh)^2) + ((Xk1 – Xh)^2) – (Vh^2)│≤ E;
│((Y21 – Yh)^2) + ((Xk2 – Xh)^2) – (Vh^2)│≤ E;
│((Y22 – Yh)^2) + ((Xk2 – Xh)^2) – (Vh^2)│≤ E;
12. Вычисление К
∆Ykh1 = Yk1 – Yh;
∆Xkh1 = Xk1 – Xh;
Kh1 = -arctg (∆Ykh1/∆Xkh1) +
+90 если ∆Xkh1≥0. +270 если ∆Xkh1<0;

∆Ykh2 = Yk2 – Yh;
∆Xkh2 = Xk2 – Xh;
Kh1 = -arctg (∆Ykh2/∆Xkh2) +
+90 если ∆Xkh1≥0. +270 если ∆Xkh1<0;

13. Вычисляем координаты упреждающей точки
Xym = X2 + (M*(X2-X1));
Ннь = Y2 + (M*(Y2 – Y1));
14. Вычисляем дистанцию
Dn1 = (-1/Vh) * (Yym * (Yk1 – Yh) + Xym * (Xk1-Xh));
Dn2 = (-1/Vh) * (Yym * (Yk2 – Yh) + Xym * (Xk2-Xh));
Если Dn1>Dn2, то выбираем Dn1. Если Dn2>Dn1, то выбираем Dn2.

15. В зависимости от выбранного D выбираем К.
Если Dn1, то K1. Если Dn2 то К2.
16. Вывод на экран полученное знаение Dn и К.

Данный по которым пытаюсь проверить:
T2 = 3 (минуты); Vh = 12(узлы); Kh = 50(градусы), П1 = 80,9(градусы); D1 = 9,5(морские мили); П2 = 80; D2 = 6,8. Ответы должны получится: Dn1 = 1,5 K1 = 329; Dn2 = 3,7 K2 = 232.
Спасибо большое))

@Sergey_B · 25.01.2013, 10:40

Может проще вводить/выводить данные в нужных единицах измерения, а расчеты все оставить так как есть (в системе СИ). Просто пусть пользователь на форме вводит то, что вам нужно (градусы, морских мили, узлы (скорость)), потом вы эти начальные данные переводите в единицы СИ, проводите расчеты, и при выводе на форму полученных значений делаете обратный перевод величин.

@Mindhunter · 25.01.2013, 11:33 **[ТС]**

Sergey_B, Я пытался это реализовать, что в коде прописывать нужные коэффициенты, а че то не получается. Но я еще кое что попробую. Спасибо

@Sergey_B · 25.01.2013, 13:06

Вот попытался сделать, но у меня что-то тоже ответ не сошелся. Может поможет...

@Mindhunter · 25.01.2013, 13:12 **[ТС]**

Sergey_B, Спасибо большое. Там в формуле опечатка. Я попробую подредактировать. А опечатка в 14 пункте.
Там Vh это другое значение, Vh = Math.sqr ((Yk - Yh)^2 + (Xk - Xh)^2)
Все равно спасибо большое. буду пробовать))

@Mindhunter · 26.01.2013, 11:46 **[ТС]**

Вот примерно так было бы удобнее для сравнения с личными расчетами.

@Sergey_B · 26.01.2013, 14:54

Не работает после 10 шага, может ошибка в формулах, или единицах измерения.

@IamRain · 26.01.2013, 16:02

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

X1 = D1*sinП1
Y 1 = D1*cos П1

По-моему надо наоборот: чтобы найти проекцию на x - x-координату, нужно гипотенузу умножать на косинус угла
, x = D1* cos(П1). и для y также.

Добавлено через 30 секунд
Y1.*

Добавлено через 50 секунд
то бишь Y1 - уже на синус. наоборот, вообщем.

@Mindhunter · 26.01.2013, 16:12 **[ТС]**

IamRain, Не, это полюбасу правильно. Доктор технических наук решал)))
И было проверенно на примере, графические построения были сделаны и на калькуляторе счетали.

@IamRain · 26.01.2013, 16:24

Надо просто смотреть геометрическое представление задачи. Я вслепую написал.

Добавлено через 44 секунды
и не факт что ваш доктор прав.

Добавлено через 7 минут
доктора тоже люди

@Mindhunter · 26.01.2013, 16:33 **[ТС]**

IamRain, Это козырный доктор))
Тут правильно все, и переправиряли, как графически так и аналитически, все сходиться.

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование igorrr37 13.02.2026 / * Логарифм записывается как: (x-2)log(x^2+2) - означает логарифм (x^2+2) по основанию (x-2). Унарный минус обозначается как ! */ #include <iostream> #include <stack> #include <cctype>. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .

@Sergey_B 174 / 174 / 57 Регистрация: 08.01.2013 Сообщений: 335
	25.01.2013, 10:40
	Может проще вводить/выводить данные в нужных единицах измерения, а расчеты все оставить так как есть (в системе СИ). Просто пусть пользователь на форме вводит то, что вам нужно (градусы, морских мили, узлы (скорость)), потом вы эти начальные данные переводите в единицы СИ, проводите расчеты, и при выводе на форму полученных значений делаете обратный перевод величин. 1

@Mindhunter 0 / 0 / 0 Регистрация: 13.12.2012 Сообщений: 25
	25.01.2013, 11:33 [ТС]
	Sergey_B, Я пытался это реализовать, что в коде прописывать нужные коэффициенты, а че то не получается. Но я еще кое что попробую. Спасибо 0

@Mindhunter 0 / 0 / 0 Регистрация: 13.12.2012 Сообщений: 25
	25.01.2013, 13:12 [ТС]
	Sergey_B, Спасибо большое. Там в формуле опечатка. Я попробую подредактировать. А опечатка в 14 пункте. Там Vh это другое значение, Vh = Math.sqr ((Yk - Yh)^2 + (Xk - Xh)^2) Все равно спасибо большое. буду пробовать)) 0

@Mindhunter 0 / 0 / 0 Регистрация: 13.12.2012 Сообщений: 25
	26.01.2013, 16:12 [ТС]
	IamRain, Не, это полюбасу правильно. Доктор технических наук решал))) И было проверенно на примере, графические построения были сделаны и на калькуляторе счетали. 0

@IamRain 4694 / 2702 / 734 Регистрация: 02.08.2011 Сообщений: 7,228
	26.01.2013, 16:24
	Надо просто смотреть геометрическое представление задачи. Я вслепую написал. Добавлено через 44 секунды и не факт что ваш доктор прав. Добавлено через 7 минут доктора тоже люди 0

@Mindhunter 0 / 0 / 0 Регистрация: 13.12.2012 Сообщений: 25
	26.01.2013, 16:33 [ТС]
	IamRain, Это козырный доктор)) Тут правильно все, и переправиряли, как графически так и аналитически, все сходиться. 0