Контрольная работа
По основам алгоритмирования и программирования
на тему:
«Исследование алгоритмов фильтрации и управления»
Содержание
Введение
1 Исследование алгоритма фильтрации
2 Исследование алгоритма управления
Приложение
Введение
Существует два метода проектирования систем автоматического управления (САУ) : экспериментальный и аналитический.
Экспериментальный метод заключается в создании различных систем без привлечения математических расчетов на основании опыта разработчика.
Использования аналитического метода при проектировании позволяет на
математических моделях просмотреть различные схемные решения систем
автоматического управления, проверить эффективность использования раз-
личных принципов и законов управления.
1 Исследование алгоритма фильтрации
Моделирование замкнутой системы управления.
Данная схема используется для слежения, программного управления и стабилизации. В такой системе регулятор в процессе управления учитывает как задание, так и реальное состояние объекта, а, кроме того, косвенно учитывает и возмущение.
شَيкِиоيальيаے ٌُهىа заىкيٍَой رہذ ٌ َпًавлهيиهى по оٍклоيهيи.
e(t) – оٍклоيهيиه (оّиلка ٌлهوهيиے, ًهمَлиًоваيиے) َпًавлےهىой вهличиيы y(t) оٍ заنаيиے ُз(t). خٌيовيыه лهىهيٍы ٌُهىы: оلْهкٍ َпًавлهيиے, коيًٍَ مлавيой оلًаٍيой ٌвےзи.
رًавيиващهه ًٌٍَойٌٍво (ٌَىىаٍоً) ًٌавيиваهٍ заنащَ и َпًавлےهىыه вهличиيы и вычиٌлےهٍ оٍклоيهيиه, оّиلкَ ًهمَлиًоваيиے
e(t) = ُз(t) - y (t).
ذهمَлےٍоً – выًаلаٍываهٍ ٍакоه َпًавлےщهه возنهйٌٍвиه u(t) يа оلْهкٍ َпًавлهيиے, коٍоًоه ٌвоنиٍ оّиلкَ к يَл или نопٌٍَиىоىَ ىиيиىَىَ. آ иنهалه, коمنа e = 0, ُз(t) = y (t)
ذаلоٍа ٌиٌٍهىы в ٌٍаٍикه
آ ٌٍаٍикه пهًهُоنيыه пًоِهٌٌы َوه закоيчилиٌь. آٌه ٌиميалы ًَавيовهّهيы и ٌلалаيٌиًоваيы. دًи ٍهкَщهى зيачهيии заنаيиے ُз(t) = const оٍклоيهيиه e (оّиلка ًهمَлиًоваيиے) поٌٍоےييа и هه вهличиيа ٍакова, чٍо ًهمَлےٍоً выًаلаٍываهٍ َпًавлےщهه возنهйٌٍвиه u, оلهٌпهчиващهه зيачهيиه َпًавлےهىой вهличиيы, نащهه поٌлه вычиٍаيиے из заنаيиے ٍоٍ ٌиميал оّиلки, коٍоًый и هٌٍь. دًичиي نлے изىهيهيиے ٌиميалов يهٍ и вهличиيа ٌиميала оّиلки ىала. دоٍоىَ َпًавлےهىаے вهличиيа пًиلлизиٍهльيо (или ٍочيо) ًавيа или пًопоًِиоيальيа заنаيи.
ذаلоٍа ٌиٌٍهىы в نиيаىикه
رлهوهيиه
دٌٍَь ٌهйчаٌ ٌиٌٍهىа يаُоنиٌٍے в ٌٍَаيовивّهىٌے ًهوиىه и e = 0. آ какой-ٍо ىоىهيٍ вًهىهيи заنаيиه يачиيаهٍ ُз(t) ًаٌٍи. دоٌколькَ ًهمَлےٍоً и оلْهкٍ َпًавлهيиے иيهًِиоييыه, ٍо َпًавлےهىаے вهличиيа в пهًвыه ىоىهيٍы вًهىهيи يه ٌَпهваهٍ возًаٌٍи. دоٍоىَ оٍклоيهيиه ه возًаٌٍаهٍ. سчиٍываے ٍо возًаٌٍаيиه ًهمَлےٍоً َвهличиваهٍ َпًавлےщَ вهличиيَ, оلْهкٍ воٌпًиيиىаهٍ ٍо َвهличهيиه и ًهаمиًَهٍ يа يهمо َвهличهيиهى َпًавлےهىой вهличиيы. فٍо пًивоنиٍ к َىهيьّهيи оٍклоيهيиے. تоيًٍَ оًٍиِаٍهльيой оلًаٍيой ٌвےзи (خخر) заىкيَлٌے. ذهمَлےٍоً пًоنолوаهٍ изىهيےٍь َпًавлےщَ вهличиيَ نо ٍهُ поً, пока оٍклоيهيиه يه َىهيьّиٌٍے نо ىалой и نаوه يَлهвой вهличиيы. زهى ٌаىыى َпًавлےهىаے вهличиيа поنٍےمиваهٌٍے к заنаيи, ٌлهنиٍ за заنаيиهى. تоيًٍَ оٌَщهٌٍвлےهٍ خخر, чٍо пًи повыّهيии оّиلки ًهمَлиًоваيиے пًивоنиٍ, ىоوهٍ لыٍь поٌٍهпهييо во вًهىهيи, к هه َىهيьّهيи.
ذهمَлےٍоً иىههٍ иيфоًىаِи о ٌоٌٍоےيии оلْهкٍа, поٌколькَ оٍклоيهيиه опًهنهлےهٌٍے вычиٍаيиهى َпًавлےщهй вهличиيы из заنаيиے.
رٍаلилизаِиے (коىпهيٌаِиے возىَщهيиے)
بзىهيهيиه возىَщهيиے пًивоنиٍ к оٍклоيهيи َпًавлےهىой вهличиيы оٍ هه ًٍهلَهىоمо зيачهيиے. فٍо يهىهنлهييо ٌказываهٌٍے يа оٍклоيهيии (оّиلкه ًهمَлиًоваيиے). ذهمَлےٍоً يа оٌيоваيии ٍоمо изىهيهيиے оٍклоيهيиے ٍак изىهيےهٍ َпًавлےщَ вهличиيَ, чٍо оٍклоيهيиه َىهيьّаهٌٍے. تоيًٍَ оًٍиِаٍهльيой оلًаٍيой ٌвےзи (خخر) заىкيَлٌے. آлиےيиه изىهيهيиے возىَщهيиے коىпهيٌиًَهٌٍے, ٍ.ه. оٌَщهٌٍвлےهٌٍے ٌٍаلилизаِиے.
بٍак, пًиيِип ًаلоٍы رہذ ٌ َпًавлهيиهى по оٍклоيهيи оٌيовываهٌٍے يа ٌвойٌٍвه коيًٍَа, оلًазоваييоمо оًٍиِаٍهльيой оلًаٍيой ٌвےзь ٌ выُоنа оلْهкٍа, ًٌٍَهىлےٍь к يَл ٌиميал يа выُоنه ٌَىىаٍоًа.
جаٍهىаٍичهٌкиه ىоنهли ٌиٌٍهى, возنهйٌٍвий и ًهакِий
جаٍهىаٍичهٌкиه ىоنهли оلْهкٍов и ٌиٌٍهى авٍоىаٍичهٌкоمо َпًавлهيиے ًٌٍоےٌٍے ٌ ِهль иُ иٌٌлهنоваيиے и поًٌٍоهيиے يа همо оٌيовه ىоنهлهй опٍиىальيыُ, оٍвهчащиُ ًٍهلоваيиےى заказчика ٌиٌٍهى.
آ оلщهى виنه оلْهкٍы и ٌиٌٍهىы опиٌываٌٍے يهлиيهйيыى نиффهًهيِиальيыى ًَавيهيиهى, ٌвےзыващиى возنهйٌٍвиه u(t) يа ٌиٌٍهىَ и همо пًоизвоنيыه по вًهىهيи ٌ ًهакِиهй v(t) ٌиٌٍهىы и هه пًоизвоنيыىи.
اачаٌٍَ пًи َпًавлهيии يهлиيهйيыىи оلْهкٍаىи пًиًащهيиه َпًавлےщиُ возنهйٌٍвий и иُ пًоизвоنيыُ ىалы по ًٌавيهيи ٌо ٌٍаِиоيаًيыىи зيачهيиےىи, поٍоىَ ٌиٌٍهىа ىоوهٍ لыٍь лиيهаًизоваيа, ٍ.ه. يهлиيهйيоه نиффهًهيِиальيоه ًَавيهيиه ىоوهٍ لыٍь заىهيهيо лиيهйيыى оٍيоٌиٍهльيо ٍиُ пًиًащهيий:
منه пًиًащهيиے x, y – ىалы, نлے пًиًащهيий يачальيыه ٌَловиے يَлهвыه.
تоффиِиهيٍы ai и kj опًهنهлےٌٍے ًаلочهй ٍочкой, в коٍоًой оٌَщهٌٍвлهيа лиيهаًизаِиے.
آ ًَавيهيии извهٌٍيы коффиِиهيٍы, оيи опًهنهлےٌٍے ٌвойٌٍваىи ٌиٌٍهىы, возنهйٌٍвиه x(t) – фَيкِиے вًهىهيи, и يачальيыه ٌَловиے – оيи يَлهвыه. بٌкоىой вهличиيой ےвлےهٌٍے фَيкِиے вًهىهيи y(t) – ًهакِиے, оٍклик ٌиٌٍهىы يа возنهйٌٍвиه x(t).
رиميалы – ٍо ىаٍهىаٍичهٌкиه ىоنهли возنهйٌٍвий и ًهакِий ٌиٌٍهى и оلْهкٍов, ٍ.ه. ٍо ىоنهли физичهٌкиُ вهличиي. رиميалы пًهنٌٍавлےٍ ٌоلой фَيкِии вًهىهيи и ےвлےٌٍے يоٌиٍهлےىи иيфоًىаِии.
Воздействия и соответственно сигналы бывают:
- непрерывными и дискретными,
- детерминированными и случайными.
دًоٌٍهйّиه ٌиميалы ًٌавيиٍهльيо лهمко مهيهًиًٌٍَے и иٌпользٌٍَے в качهٌٍвه пًоليыُ пًи иٌٌлهنоваيии оلْهкٍов и ٌиٌٍهى َпًавлهيиے.
ر поىощь пًоٌٍهйّиُ ٌиميалов ىоويо пًهنٌٍавиٍь пًоизвольيый (ٌлоويый ٌиميал).
دًоٌٍهйّиه ٌиميалы:
гармонический (синусоидальный);
ступенчатый единичный сигнал (10(t) );
дельта - функция Дирака (δ(t));
степенная функция (tν).
أаًىоيичهٌкий ٌиميалأаًىоيичهٌкий ٌиميал изىهيےهٌٍے во вًهىهيи по ٌиيٌَоиنальيоىَ закоيَ:
X(t) = Xm *sin(ωt + φ)
ر поىощь ٍакоمо ٌиميала или يаلоًа (ٌَىىы) ٍакиُ ٌиميалов َنоليо ىоنهлиًоваٍь пهًиоنичهٌкиه возنهйٌٍвиے يа ٌиٌٍهىы, يапًиىهً виلًаِии, а ٍакوه ىоويо ىоنهлиًоваٍь ٌиميалы пًоизвольيоمо виنа. تًоىه ٍоمо, ٌиيٌَоиنальيый ٌиميал иٌпользَهٌٍے в качهٌٍвه пًоليоمо пًи иٌٌлهنоваيии ٌٍَаيовивّهمоٌے ًهوиىа ًаلоٍы رہذ.
حهзавиٌиىаے пهًهىهييаے - вًهىے t ىهيےهٌٍے оٍ ىиيٌَ نо плٌ لهٌкоيهчيоٌٍи.
أаًىоيичهٌкий ٌиميал ُаًакٍهًизَهٌٍے ًٍهىے паًаىهًٍаىи (чиٌлаىи):
- амплитудой Xm – максимальной величиной сигнала,
- круговой частотой ω [рад/сек] и
- начальной фазой φ [градусов].
حаًےنَ ٌ يазваييыىи оٌيовيыىи паًаىهًٍаىи ٌиيٌَоиنальيоمо ٌиميала иٌпользٌٍَے и نополيиٍهльيыه паًаىهًٍы:
- f [Гц] = ω/2π - циклическая частота;
- Т [сек] = 1/f – период.
ہлمоًиٍى фильًٍаِии.
ؤлے поٌлهنоваٍهльيой оلًаلоٍки ًٍهلَهٌٍے ٍهкَщаے оِهيка ٌиميала يа каونоى ٍакٍه ٌ َчهٍоى иيфоًىаِии, поٌٍَпащهй يа вُоن фильًٍа в пًоِهٌٌه يаلлنهيиے. آ يаٌٍоےщهه вًهىے ّиًокоه ًаٌпًоًٌٍаيهيиه полَчили аنапٍивيыه фильًٍы, в коٍоًыُ поٌٍَпащаے يоваے иيфоًىаِиے иٌпользَهٌٍے نлے يهпًهًывيой коًًهкٍиًовки ًаيهه ٌنهлаييой оِهيки ٌиميала (ٌопًовоونهيиه ِهли в ًаنиолокаِии, ٌиٌٍهىы авٍоىаٍичهٌкоمо ًهمَлиًоваيиے в َпًавлهيии и ٍ.ن).
Алгоритм фильтрации технологического параметра
Исходные данные
АХ=2 К=1,5
АЕ=0,15 Т=40
ТХ=50 τ=15
ТЕ=8 Т0=1
Задача:
Отфильтровать входящий сигнал
Выделить полезную составляющую сигнала
Исследовать алгоритм настройки пропорционально-интегрального модулятора для построения замкнутой системы управления
Для выделения полезного сигнала, перед использованием его в качестве сигнала обратной связи регулятора. Выделяют несколько видов фильтров:
- Фильтр текущего среднего
- Экспоненциальный фильтр
Расчет фильтра текущего среднего
Находим:
x(t)=Ax*sinώxt - полезный сигнал
Е(t)=AE*sinώEt - помеха
z(t)=x(t)+E(t) - входной сигнал
ώx= ώе=
y(t) - выходной сигнал
y2(t)=[z(t)+z(t-1)]
y3(t)=[z(t)+z(t-1)+z(t-2)]
y4(t)=[z(t)+z(t-1)+z(t-2)+z(t-3)]
Расчет экспоненциального фильтра
y(t)=б*y(t-1)+(1-б)*z(t-1)
б=exp(-)
Тб=5ч10 для расчета выбираем 5
Расчеты сводим в таблицу:
Алгоритм фильтрации | |||||||
Ax= | 2 | ||||||
Ae= | 0,15 | ||||||
Tx= | 50 | юx= | 0,13 | ||||
Te= | 8 | юe= | 0,79 | ||||
Tб= | 5 | б= | 0,82 | ||||
t | Z(t) | E(t) | X(t) | Y2(t) | Y3(t) | Y4(t) | Ya |
0 | 0,00 | 0,00 | 0,00 | 0,00 | 0,00 | 0,00 | 0,00 |
1 | 0,36 | 0,11 | 0,25 | 0,18 | 0,18 | 0,18 | 0,00 |
2 | 0,65 | 0,15 | 0,50 | 0,50 | 0,50 | 0,50 | 0,06 |
3 | 0,84 | 0,11 | 0,74 | 0,74 | 0,62 | 0,46 | 0,17 |
4 | 0,96 | 0,00 | 0,96 | 0,90 | 0,82 | 0,70 | 0,29 |
5 | 1,07 | -0,11 | 1,18 | 1,02 | 0,96 | 0,88 | 0,41 |
6 | 1,22 | -0,15 | 1,37 | 1,14 | 1,08 | 1,02 | 0,53 |
7 | 1,43 | -0,11 | 1,54 | 1,33 | 1,24 | 1,17 | 0,66 |
8 | 1,69 | 0,00 | 1,69 | 1,56 | 1,45 | 1,35 | 0,80 |
9 | 1,92 | 0,11 | 1,81 | 1,80 | 1,68 | 1,56 | 0,96 |
10 | 2,05 | 0,15 | 1,90 | 1,98 | 1,89 | 1,77 | 1,13 |
11 | 2,07 | 0,11 | 1,96 | 2,06 | 2,01 | 1,93 | 1,30 |
12 | 2,00 | 0,00 | 2,00 | 2,03 | 2,04 | 2,01 | 1,44 |
13 | 1,89 | -0,11 | 2,00 | 1,94 | 1,99 | 2,00 | 1,54 |
14 | 1,81 | -0,15 | 1,96 | 1,85 | 1,90 | 1,94 | 1,60 |
15 | 1,80 | -0,11 | 1,90 | 1,81 | 1,83 | 1,87 | 1,64 |
16 | 1,81 | 0,00 | 1,81 | 1,80 | 1,81 | 1,83 | 1,67 |
17 | 1,79 | 0,11 | 1,69 | 1,80 | 1,80 | 1,80 | 1,70 |
18 | 1,69 | 0,15 | 1,54 | 1,74 | 1,77 | 1,77 | 1,71 |
19 | 1,48 | 0,11 | 1,37 | 1,58 | 1,65 | 1,69 | 1,71 |
20 | 1,18 | 0,00 | 1,18 | 1,33 | 1,45 | 1,53 | 1,67 |
21 | 0,86 | -0,11 | 0,96 | 1,02 | 1,17 | 1,30 | 1,58 |
22 | 0,59 | -0,15 | 0,74 | 0,72 | 0,87 | 1,02 | 1,45 |
23 | 0,39 | -0,11 | 0,50 | 0,49 | 0,61 | 0,75 | 1,29 |
24 | 0,25 | 0,00 | 0,25 | 0,32 | 0,41 | 0,52 | 1,13 |
25 | 0,11 | 0,11 | 0,00 | 0,18 | 0,25 | 0,33 | 0,97 |
26 | -0,10 | 0,15 | -0,25 | 0,00 | 0,09 | 0,16 | 0,81 |
27 | -0,39 | 0,11 | -0,50 | -0,25 | -0,13 | -0,03 | 0,65 |
28 | -0,74 | 0,00 | -0,74 | -0,56 | -0,41 | -0,28 | 0,46 |
29 | -1,07 | -0,11 | -0,96 | -0,90 | -0,73 | -0,57 | 0,24 |
30 | -1,33 | -0,15 | -1,18 | -1,20 | -1,04 | -0,88 | 0,00 |
31 | -1,48 | -0,11 | -1,37 | -1,40 | -1,29 | -1,15 | -0,24 |
32 | -1,54 | 0,00 | -1,54 | -1,51 | -1,45 | -1,35 | -0,46 |
33 | -1,58 | 0,11 | -1,69 | -1,56 | -1,53 | -1,48 | -0,66 |
34 | -1,66 | 0,15 | -1,81 | -1,62 | -1,59 | -1,56 | -0,82 |
35 | -1,80 | 0,11 | -1,90 | -1,73 | -1,68 | -1,64 | -0,98 |
36 | -1,96 | 0,00 | -1,96 | -1,88 | -1,81 | -1,75 | -1,12 |
37 | -2,10 | -0,11 | -2,00 | -2,03 | -1,95 | -1,88 | -1,28 |
38 | -2,15 | -0,15 | -2,00 | -2,12 | -2,07 | -2,00 | -1,43 |
39 | -2,07 | -0,11 | -1,96 | -2,11 | -2,11 | -2,07 | -1,56 |
40 | -1,90 | 0,00 | -1,90 | -1,99 | -2,04 | -2,06 | -1,65 |
41 | -1,70 | 0,11 | -1,81 | -1,80 | -1,89 | -1,96 | -1,70 |
42 | -1,54 | 0,15 | -1,69 | -1,62 | -1,71 | -1,80 | -1,70 |
43 | -1,43 | 0,11 | -1,54 | -1,49 | -1,56 | -1,64 | -1,67 |
44 | -1,37 | 0,00 | -1,37 | -1,40 | -1,45 | -1,51 | -1,63 |
45 | -1,28 | -0,11 | -1,18 | -1,33 | -1,36 | -1,41 | -1,58 |
46 | -1,11 | -0,15 | -0,96 | -1,20 | -1,25 | -1,30 | -1,53 |
47 | -0,84 | -0,11 | -0,74 | -0,98 | -1,08 | -1,15 | -1,45 |
48 | -0,50 | 0,00 | -0,50 | -0,67 | -0,82 | -0,93 | -1,34 |
49 | -0,14 | 0,11 | -0,25 | -0,32 | -0,49 | -0,65 | -1,19 |
50 | 0,15 | 0,15 | 0,00 | 0,00 | -0,16 | -0,33 | -1,00 |
По полученным данным строим график
2 Исследование алгоритма управления
Находим:
x(t) – сигнал задания
у(t) = б*у(t-1)+k(1-б)*U(t-1-n) – выходной сигнал
где: б=exp(-)
n=
U(t)=U(t-1)+g0*E(t)+g*E(t-1) - сигнал на исполнительный механизм
Где: g0=kR g= -kR(1-)
E(t)=x(t)-y(t) -ошибка регулирования
Для расчета изначально принимаем:
kR=0.05 - коэффициент усиления регулятора
ТI= 999 - время интегрирования
Расчетные данные сводим в таблицу : Приложение А.
По полученным данным строим график:
Настраиваем регулятор:
Отключаем интегрирующую составляющую ТI=999
Постепенно увеличиваем коэффициент усиления регулятора, выводим систему в режим автоколебания.
Запоминаем значение коэффициента усиления регулятора соответствующее
режиму автоколебания
kRкр= 0,2486
kR= kRкр Определяем период автоколебания
ТА = 50 Вычисляем параметры:
kR= kRкр/2=0,1243
ТI=(0.3ч0.8)TA
TI=0.5*50=25
Строим график:
Литература
1. Рокотян С.С. и Шапиро И.М. - «Справочник по фильтрациям систем» 3-е издание, переработанное и дополненное, Москва, Энергоатомиздат, 1985 г; [1]
2.Солдаткина Л.А. – Автоматизированные информационные системы:учебное пособие для вузов. Москва, Энергия, 1978 г.;