Общие сведения о численном интегрировании обыкновенного дифференциального уравнения
Многие задачи
сводятся к решению одного или нескольких
обыкновенных дифференциальных уравнений
первого или более высокого порядка.
Причем дифференциальное уравнение
высокого порядка почти всегда (если
может выделена старшая производная),
может быть сведено к системе дифференциальных
уравнений первого порядка. Решение
такой системы уравнений можно рассматривать
как решение одного уравнения
.
Формальное отличие состоит в том, что вместо скалярных величин для системы уравнений необходимо вводить векторы. В нашем случае вектор y имеет следующую форму:
.
Рассмотрим простейшее
обыкновенное дифференциальное уравнение
первого порядка
.
Необходимо найти
,
т.е.решение этого уравнения.
Полученное уравнение
описывает множество решений. Для
выделения конкретного решения необходимы
дополнительные условия. Задача значительно
упрощается, если начальные условия
задаются. Постановка задачи при данных
условиях есть задача Коши. При определенных
условиях решение этой задачи существует,
и оно является единственным. К сожалению
большинство таких задач точно не
решаются. Их можно решить приближенно
аналитически или численно.
В настоящее время разработано большое количество различных численных методов интегрирования обыкновенных дифференциальных уравнений первого прядка. Выбор конкретной формулы всегда является компромиссом между весьма противоречивыми желаниями. Все существующие методы можно классифицировать, анализируя различные стороны этих методов. При этом необходимо учитывать:
точность;
устойчивость вычислительной схемы;
трудоемкость (затрачиваемое время);
необходимую память для хранения информации;
какая нужна информация перед каждым шагом (никакой, значения функций или производных в m уже вычисленных точках) и т. д.
При выполнении численного интегрирования в проекте используется три метода: метод Эйлера, модифицированный метод Эйлера и метод Рунге-Кутта.
Метод Эйлера
Сущность метода
Эйлера заключается в том что производная
на каждом шаге интегрирования принимается
постоянной, равной значению
в начале каждого шага. Таким образом
действительная интегральная кривая,
проходящая через заданную точку,
заменяется ломаной, состоящей из
прямолинейных отрезков (ломаная Эйлера):
Здесь и далее h — шаг интегрирования.
Метод Эйлера обладает
малой точностью и дает сравнительно
удовлетворительные результаты лишь
при очень малых шагах интегрирования.
Он является абсолютно точным при
,
т.е. для линейных зависимостей. Метод
является одношаговым, т.е.позволяет в
любой момент изменять шаг.
Для заданного момента
времени значения Н, параметры атмосферы
, а также величины
меняются не значительно, поэтому их
изменениями можно пренебречь.
Таблица 4.3. – Метод Эйлера
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
t,s |
V,m/s |
X,m |
H,m |
Xc,m |
Rc,m |
β° |
q° |
e° |
13,5 |
474 |
3093 |
145 |
7600 |
4507,2246 |
-0,5720489 |
-1 |
0,4279511 |
14 |
487,9609 |
3329,964 |
140,8638 |
7400 |
4070,2412 |
-0,575239 |
-1,017635114 |
0,4423961 |
14,5 |
501,4803 |
3573,906 |
136,5306 |
|
|
|
-1,034958803 |
|
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
cos(q) |
sin(q) |
V*cos(q) |
V*sin(q) |
m,kg |
P(H) |
a(H) |
g(H) |
mg |
0,999848 |
-0,017452 |
473,9278 |
-8,272441 |
31,025 |
99600 |
339,737 |
9,806 |
304,23115 |
0,999842 |
-0,01776 |
487,884 |
-8,666259 |
30,6 |
99640 |
339,753 |
9,806 |
300,0636 |
19 |
20 |
21 |
22 |
23 |
24 |
25 |
26 |
27 |
M |
qS |
Cxtr |
Cxд |
Cxво |
Cx0 |
Cyalp |
Alp1 |
Cx1 |
1,395197 |
10721,5 |
0,029 |
0 |
0,048325 |
0,077325 |
3,28 |
0,457817034 |
0,000927 |
1,436223 |
11365,86 |
0,029 |
0 |
0,047334 |
0,076334 |
3,314 |
0,4228063 |
0,0008976 |
28 |
29 |
30 |
31 |
32 |
33 |
34 |
35 |
36 |
Cxa |
Xa |
Nxa |
Vdot |
DCY/AR |
Cya |
Ya |
Nya |
TetaDot |
0,078252 |
838,9786 |
2,829977 |
27,92189 |
0,006704 |
0,0262602 |
281,5486 |
0,970091863 |
-0,03527 |
0,077232 |
877,8037 |
2,739919 |
27,03879 |
0,006059 |
0,0244981 |
278,44217 |
0,969750979 |
-0,0346 |
Таблица 4.7. – Определение угла атаки при t=13,5 с
Teta(dot) |
-0,000402 |
|
|
Nya |
0,9701 |
|
|
А |
295,1042 |
|
|
АВ |
0,027525 |
|
|
Ср |
0,15856 |
Dcy(Alp0) |
DCY/AR1 |
Cya(alp) |
3,28 |
0 |
0,006134 |
Alp1 |
Alp2 |
|
|
0,008005 |
0,00799 |
Радианы |
|
0,458634 |
0,457817 |
Градусы |
|
Таблица 4.8. – Определение угла атаки при t=14 с
Teta(dot) |
-0,000389 |
-0,001 |
|
Nya |
0,9500861 |
0,97 |
|
А |
291,06169 |
Dcy(Alp0) |
0 |
АВ |
0,0256084 |
|
|
Ср |
0,1495707 |
|
DCY/AR1 |
Cya(alp) |
3,314 |
|
0,006704 |
Alp1 |
Alp2 |
|
|
0,007394 |
0,0073794 |
Радианы |
|
0,423625 |
0,4228063 |
градусы |
|
Применение метода Эйлера заключается в:
