- •«Исследование динамических систем»
- •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
1.3.3 Проверка и сравнение результатов идентификации
Проверка результатов идентификации осуществляется путем численного решения методом трапеции:
, (1.11)
где
где А - коэффициент, отражающий инерциальные свойства объекта,
K- коэффициент усиления,
Xi– воздействие на объект в текущий момент времени,
Xi-1– воздействие на объект в предыдущий момент времени,
Yi– реакция объекта в текущий момент времени,
Yi-1– реакция объекта в предыдущий момент времени,
время.
Исходные данные
AE= |
1,018187 |
KE= |
1,000887 |
AD= |
1,00296 |
KD= |
1,00337 |
t= |
0,15 |
где AE- коэффициент, отражающий инерциальные свойства объекта, полученный в средеExcel, kE - коэффициент усиления, полученный в средеExcel,AD- коэффициент, отражающий инерциальные свойства объекта, полученный в средеDelphi, kD- коэффициент усиления, полученный в средеDelphi.
Таблица 1.2 Проверка идентификации | |||||
i |
|
X |
Y |
YвостE |
YвостD |
1 |
0,00 |
0,079 |
0,502 |
0,5017 |
0,5017 |
2 |
0,15 |
0,106 |
0,448 |
0,4456 |
0,4448 |
3 |
0,30 |
0,138 |
0,387 |
0,4012 |
0,4000 |
4 |
0,45 |
0,138 |
0,354 |
0,3651 |
0,3636 |
5 |
0,60 |
0,189 |
0,332 |
0,3375 |
0,3359 |
6 |
0,75 |
0,197 |
0,301 |
0,3177 |
0,3161 |
7 |
0,90 |
0,231 |
0,313 |
0,3035 |
0,3020 |
8 |
1,05 |
0,280 |
0,289 |
0,2969 |
0,2956 |
9 |
1,20 |
0,304 |
0,308 |
0,2963 |
0,2952 |
10 |
1,35 |
0,340 |
0,312 |
0,2999 |
0,2991 |
11 |
1,50 |
0,375 |
0,299 |
0,3078 |
0,3074 |
12 |
1,65 |
0,399 |
0,308 |
0,3187 |
0,3187 |
13 |
1,80 |
0,441 |
0,326 |
0,3327 |
0,3330 |
14 |
1,95 |
0,475 |
0,362 |
0,3499 |
0,3505 |
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
X |
Y |
YвостE |
YвостD |
15 |
2,10 |
0,526 |
0,365 |
0,3706 |
0,3716 |
16 |
2,25 |
0,556 |
0,406 |
0,3940 |
0,3954 |
17 |
2,40 |
0,607 |
0,417 |
0,4198 |
0,4216 |
18 |
2,55 |
0,637 |
0,439 |
0,4476 |
0,4498 |
19 |
2,70 |
0,685 |
0,492 |
0,4770 |
0,4794 |
20 |
2,85 |
0,706 |
0,506 |
0,5070 |
0,5098 |
21 |
3,00 |
0,731 |
0,550 |
0,5361 |
0,5392 |
22 |
3,15 |
0,773 |
0,582 |
0,5658 |
0,5691 |
23 |
3,30 |
0,795 |
0,609 |
0,5959 |
0,5994 |
24 |
3,45 |
0,837 |
0,643 |
0,6261 |
0,6299 |
25 |
3,60 |
0,871 |
0,646 |
0,6575 |
0,6615 |
26 |
3,75 |
0,889 |
0,694 |
0,6881 |
0,6923 |
27 |
3,90 |
0,902 |
0,721 |
0,7167 |
0,7210 |
28 |
4,05 |
0,919 |
0,743 |
0,7433 |
0,7477 |
|
|
|
|
|
|
Продолжение таблицы 1.2 | |||||
i |
|
X |
Y |
YвостE |
YвостD |
29 |
4,20 |
0,949 |
0,782 |
0,7696 |
0,7741 |
30 |
4,35 |
0,953 |
0,809 |
0,7946 |
0,7991 |
31 |
4,50 |
0,976 |
0,836 |
0,8180 |
0,8226 |
32 |
4,65 |
0,997 |
0,840 |
0,8413 |
0,8459 |
33 |
4,80 |
0,984 |
0,855 |
0,8619 |
0,8665 |
34 |
4,95 |
1,012 |
0,880 |
0,8806 |
0,8852 |
35 |
5,10 |
0,985 |
0,889 |
0,8969 |
0,9014 |
36 |
5,25 |
0,999 |
0,910 |
0,9101 |
0,9145 |
37 |
5,40 |
0,989 |
0,936 |
0,9217 |
0,9260 |
38 |
5,55 |
0,987 |
0,938 |
0,9309 |
0,9350 |
39 |
5,70 |
0,973 |
0,956 |
0,9377 |
0,9418 |
40 |
5,85 |
0,962 |
0,935 |
0,9419 |
0,9457 |
41 |
6,00 |
0,951 |
0,949 |
0,9440 |
0,9477 |
42 |
6,15 |
0,942 |
0,946 |
0,9444 |
0,9479 |
43 |
6,30 |
0,909 |
0,958 |
0,9419 |
0,9452 |
44 |
6,45 |
0,889 |
0,929 |
0,9361 |
0,9392 |
45 |
6,60 |
0,882 |
0,918 |
0,9293 |
0,9321 |
46 |
6,75 |
0,859 |
0,912 |
0,9212 |
0,9239 |
47 |
6,90 |
0,816 |
0,909 |
0,9098 |
0,9122 |
48 |
7,05 |
0,799 |
0,905 |
0,8959 |
0,8980 |
49 |
7,20 |
0,744 |
0,885 |
0,8788 |
0,8807 |
50 |
7,35 |
0,720 |
0,856 |
0,8588 |
0,8604 |
51 |
7,50 |
0,687 |
0,830 |
0,8376 |
0,8389 |
52 |
7,65 |
0,646 |
0,825 |
0,8142 |
0,8153 |
53 |
7,80 |
0,615 |
0,779 |
0,7891 |
0,7899 |
54 |
7,95 |
0,579 |
0,776 |
0,7628 |
0,7633 |
55 |
8,10 |
0,539 |
0,750 |
0,7349 |
0,7351 |
56 |
8,25 |
0,505 |
0,717 |
0,7058 |
0,7057 |
57 |
8,40 |
0,460 |
0,667 |
0,6752 |
0,6748 |
58 |
8,55 |
0,432 |
0,645 |
0,6437 |
0,6432 |
59 |
8,70 |
0,392 |
0,602 |
0,6119 |
0,6112 |
60 |
8,85 |
0,363 |
0,587 |
0,5798 |
0,5788 |
61 |
9,00 |
0,332 |
0,545 |
0,5479 |
0,5467 |
62 |
9,15 |
0,293 |
0,528 |
0,5156 |
0,5143 |
63 |
9,30 |
0,265 |
0,491 |
0,4832 |
0,4817 |
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
X |
Y |
YвостE |
YвостD |
64 |
9,45 |
0,212 |
0,435 |
0,4497 |
0,4480 |
65 |
9,60 |
0,189 |
0,402 |
0,4156 |
0,4137 |
66 |
9,75 |
0,154 |
0,391 |
0,3821 |
0,3800 |
67 |
9,90 |
0,148 |
0,343 |
0,3504 |
0,3483 |
68 |
10,05 |
0,114 |
0,320 |
0,3203 |
0,3181 |
69 |
10,20 |
0,077 |
0,289 |
0,2895 |
0,2871 |
70 |
10,35 |
0,079 |
0,249 |
0,2604 |
0,2581 |
71 |
10,50 |
0,056 |
0,226 |
0,2340 |
0,2316 |
72 |
10,65 |
0,031 |
0,195 |
0,2079 |
0,2055 |
73 |
10,80 |
0,034 |
0,190 |
0,1838 |
0,1814 |
74 |
10,95 |
0,021 |
0,159 |
0,1623 |
0,1600 |
75 |
11,10 |
0,008 |
0,150 |
0,1421 |
0,1398 |
76 |
11,25 |
0,002 |
0,133 |
0,1230 |
0,1208 |
77 |
11,40 |
0,003 |
0,113 |
0,1063 |
0,1041 |
78 |
11,55 |
0,005 |
0,101 |
0,0922 |
0,0901 |
79 |
11,70 |
0,004 |
0,090 |
0,0802 |
0,0782 |
80 |
11,85 |
0,017 |
0,064 |
0,0706 |
0,0688 |
81 |
12,00 |
0,031 |
0,062 |
0,0642 |
0,0626 |
82 |
12,15 |
0,042 |
0,066 |
0,0604 |
0,0590 |
83 |
12,30 |
0,038 |
0,052 |
0,0577 |
0,0564 |
84 |
12,45 |
0,060 |
0,059 |
0,0566 |
0,0555 |
85 |
12,60 |
0,073 |
0,065 |
0,0580 |
0,0571 |
86 |
12,75 |
0,113 |
0,049 |
0,0628 |
0,0621 |
87 |
12,90 |
0,120 |
0,065 |
0,0702 |
0,0698 |
88 |
13,05 |
0,158 |
0,080 |
0,0796 |
0,0794 |
89 |
13,20 |
0,199 |
0,095 |
0,0932 |
0,0933 |
90 |
13,35 |
0,219 |
0,093 |
0,1091 |
0,1095 |
91 |
13,50 |
0,233 |
0,131 |
0,1252 |
0,1259 |
92 |
13,65 |
0,276 |
0,132 |
0,1430 |
0,1439 |
93 |
13,80 |
0,302 |
0,166 |
0,1630 |
0,1642 |
94 |
13,95 |
0,340 |
0,203 |
0,1848 |
0,1862 |
95 |
14,10 |
0,379 |
0,202 |
0,2088 |
0,2105 |
96 |
14,25 |
0,430 |
0,237 |
0,2357 |
0,2377 |
97 |
14,40 |
0,464 |
0,270 |
0,2648 |
0,2670 |
98 |
14,55 |
0,500 |
0,282 |
0,2946 |
0,2972 |
где - время, X - воздействие на объект, Y – экспериментальные данные реакция объекта, YвостE – проверка решения в среде Excel, YвостD - проверка решения в среде Delphi.
Рисунок 1.7 – сравнение результатов идентификации