- •«Исследование динамических систем»
- •1 Идентификация обыкновенного линейного дифференциального уравнения 1го порядка
- •1.1 Постановка задачи
- •1.2 Описание используемых методов
- •1.2.1 Аппроксимация на смежных отрезках
- •1.2.1 Аппроксимация на скользящих интервалах
- •1.3 . Результаты решения задачи аппроксимации
- •1.3.1 Проведение идентификации в среде Excel
- •1.3.2 Проведение идентификации в среде Delphi
- •1.3.3 Проверка и сравнение результатов идентификации
- •1.4 Заключение
- •2 Исследование динамики системы
- •2.1 Постановка задачи
- •2.2 Запись конечно-разностных аналогов дифференциальных уравнений
- •2.3 Решение в среде Excel
- •2.4 Решение в среде Delphi
- •2.5 Заключение
- •Литература
- •Приложение 1
- •Приложение 2
2.4 Решение в среде Delphi
Блок-схема алгоритма численного моделирования динамических процессов в системе автоматического регулирования изображена на рисунке 2.3. Код программы находятся в приложении2.
Исходные данные: | |||||
B= |
1 |
Z0= |
0 |
|
|
С= |
4 |
a0= |
1 |
|
|
k1= |
4 |
a1= |
1 |
|
|
k2= |
-1,5 |
a2= |
1 |
|
|
k3= |
2 |
a3= |
1 |
|
|
k4= |
1,5 |
a4= |
-0,1 |
|
|
k5= |
1 |
b0= |
2 |
|
|
dt= |
0,1 |
b1= |
0 |
|
|
t0= |
0 |
b2= |
0 |
|
|
Y0= |
0 |
|
|
|
|
|
|
|
|
|
|
Коэффициенты, получены в результате расчета: | |||||
D1= |
0,90909 |
D5= |
0,41707 |
|
|
D2= |
0,36364 |
D6= |
-0,3659 |
|
|
D3= |
-0,1364 |
D7= |
-0,4171 |
|
|
D4= |
0,97561 |
D8= |
0,36585 |
|
|
Рисунок 2.3 – Блок-схема алгоритма численного моделирования динамических процессов
Таблица 2.2 – Решение системы уравнений (2.4) в среде Delphi
|
|
|
|
|
|
I |
τ |
X |
G |
Y |
Z |
0 |
0,0 |
1,841 |
2 |
0,000 |
0,000 |
1 |
0,1 |
1,882 |
2 |
0,661 |
0,173 |
2 |
0,2 |
1,914 |
2 |
1,249 |
0,346 |
3 |
0,3 |
1,935 |
2 |
1,769 |
0,516 |
4 |
0,4 |
1,947 |
2 |
2,223 |
0,681 |
5 |
0,5 |
1,949 |
2 |
2,616 |
0,839 |
6 |
0,6 |
1,941 |
2 |
2,949 |
0,989 |
7 |
0,7 |
1,925 |
2 |
3,227 |
1,130 |
8 |
0,8 |
1,899 |
2 |
3,452 |
1,259 |
9 |
0,9 |
1,865 |
2 |
3,629 |
1,376 |
10 |
1,0 |
1,823 |
2 |
3,760 |
1,481 |
11 |
1,1 |
1,773 |
2 |
3,848 |
1,572 |
12 |
1,2 |
1,717 |
2 |
3,898 |
1,649 |
13 |
1,3 |
1,655 |
2 |
3,912 |
1,712 |
14 |
1,4 |
1,587 |
2 |
3,894 |
1,760 |
15 |
1,5 |
1,515 |
2 |
3,846 |
1,794 |
16 |
1,6 |
1,439 |
2 |
3,772 |
1,814 |
17 |
1,7 |
1,361 |
2 |
3,676 |
1,821 |
18 |
1,8 |
1,280 |
2 |
3,560 |
1,814 |
19 |
1,9 |
1,198 |
2 |
3,427 |
1,794 |
20 |
2,0 |
1,116 |
2 |
3,281 |
1,762 |
21 |
2,1 |
1,034 |
2 |
3,124 |
1,720 |
22 |
2,2 |
0,953 |
2 |
2,959 |
1,667 |
23 |
2,3 |
0,875 |
2 |
2,790 |
1,604 |
24 |
2,4 |
0,799 |
2 |
2,617 |
1,534 |
|
|
|
|
|
|
|
|
|
|
|
|
I |
τ |
X |
G |
Y |
Z |
25 |
2,5 |
0,727 |
2 |
2,445 |
1,456 |
26 |
2,6 |
0,659 |
2 |
2,275 |
1,373 |
27 |
2,7 |
0,596 |
2 |
2,110 |
1,284 |
28 |
2,8 |
0,538 |
2 |
1,951 |
1,192 |
29 |
2,9 |
0,485 |
2 |
1,800 |
1,098 |
30 |
3,0 |
0,439 |
2 |
1,660 |
1,002 |
31 |
3,1 |
0,400 |
2 |
1,531 |
0,906 |
32 |
3,2 |
0,367 |
2 |
1,414 |
0,812 |
33 |
3,3 |
0,341 |
2 |
1,312 |
0,719 |
34 |
3,4 |
0,323 |
2 |
1,224 |
0,630 |
35 |
3,5 |
0,311 |
2 |
1,152 |
0,545 |
36 |
3,6 |
0,307 |
2 |
1,095 |
0,464 |
37 |
3,7 |
0,309 |
2 |
1,055 |
0,390 |
38 |
3,8 |
0,319 |
2 |
1,031 |
0,322 |
39 |
3,9 |
0,335 |
2 |
1,023 |
0,261 |
40 |
4,0 |
0,357 |
2 |
1,032 |
0,208 |
41 |
4,1 |
0,386 |
2 |
1,056 |
0,164 |
42 |
4,2 |
0,420 |
2 |
1,095 |
0,128 |
43 |
4,3 |
0,459 |
2 |
1,149 |
0,101 |
44 |
4,4 |
0,502 |
2 |
1,216 |
0,082 |
45 |
4,5 |
0,550 |
2 |
1,295 |
0,073 |
46 |
4,6 |
0,601 |
2 |
1,386 |
0,073 |
47 |
4,7 |
0,656 |
2 |
1,487 |
0,082 |
48 |
4,8 |
0,713 |
2 |
1,598 |
0,100 |
49 |
4,9 |
0,771 |
2 |
1,715 |
0,126 |
Рисунок 2.3 – результат решения в среде Delphi