.1.2 Перевірка адекватності дінамічної МОДЕЛІ
Перевірку адекватності МОДЕЛІ здійснюю такоже помощью програмного середовища Matlab.
t=[0 20 40 60 80 100 120 140 160 200 440 480 560]; _x=[3.9 3.92 4 4.1 4.25 4.4 4.5 4.6 4.7 4.8 4.9 4.9 4.9]; _n=3.9; dT=4.9-3.9;=37.522;i=1:length(t)(i)=(1-exp(-t(i)./T).*(1+t(i)./T+0.5.*(t(i)./T).^2))*dT+T_n;(i)=abs(Q_x(i)-h(i))/dT*100_pohubka =Max (delta)
Результатом Виконання програми є: _pohubka=1.9800-max відносна похібка,%
Відносна похібка НЕ ??перевіщує 3%, отже Знайду модель адекватна заданій експеріментальній крівій Розгон.
5.1.1 Знаходження дінамічної МОДЕЛІ допоміжного об'єкта регулювання
.1.1.1 Вибір структурованих МОДЕЛІ та розрахунок ее параметрів
будую експериментальну криву Розгон, отриманий при стрібкоподібній зміні регулюючої Дії.
t=[0 4 39 74 109 143 179 214 285 425 500]; _x=[30 31 33 34.5 36 38 40 41 43 44 44.0001]. * 10000; (t, T_x, « b »); grid;
рис.5.4 Крива Розгон, отримай при стрібкоподібній зміні регулюючої Дії у=20% ходу РВ.
а) Оберемо функцію передачі у вігляді:
де, T - стала годині, n - кількість аперіодічніх ланок.
б) Алгоритм знаходження параметрів математичної моделі:
. Для зручності розрахунків експериментальну криву Розгон обєкта регулювання нормо діленням ее значень на максимальну зміну віхідної величини.
В результаті отримайте нормовану перехідну функцію (рис.5.5).
рис.5.5 Нормована експериментальна перехідна функція ОР.
2. З нормованої експериментальної перехідної Функції знаходжусь Значення годині і. Розраховуйтесь відношення и перевіряю Виконання умови 0.430.68.
=126/270.8=0.46. Ця Умова віконується и найближче
відповідне відношенню Значення n=2.
. Для Отримання Значення n=2 з табл.1 визначавши відношення
и з якіх знаходжусь сталлю годині T.
=1.68=3.89=126/1.68=75;=270.8/3.89=69.6144;=(T1 + T2) / 2=72.3072 c
Отже функція передачі ЗР буде наступна:
в) Для Отримання значень n, T знаходжусь перехідну функцію:
Коефіцієнт передачі k =,
де - стрібкоподібна зміна вхідної розмірів
експериментальна и розраховану кріві Розгон будую помощью програмного середовища Matlab.
t=[0 4 39 74 109 143 179 214 285 425 500]; _x=[30 30.2 31.5 34.1 36 38 40 41 43 44 44]; _n=30; dT=44-30;=72.3072;=[ 0:0.01:500];=(1-exp (-t1. / T). * (1 + t1. / T)) * dT + T_n;
plot (t1, h, «b», t, T_x, «r *"); grid («t, c»); ylabel («м3/рік»);
ріс.5.6 порівняння кривих Розгон:
" *" - експериментальна крива Розгон;
" -- Знайду аналітично.
5.1.1.2 Перевірка адекватності дінамічної МОДЕЛІ
Перевірку адекватності МОДЕЛІ здійснюю такоже помощью програмного середовища Matlab.
t=[0 4 39 74 109 143 179 214 285 425 500]; _x=[30 30.2 31.5 34.1 36 38 40 41 43 44 44];_n=30;dT=44-30;=72.3072;i=1:length(t)(i)=(1-exp(-t(i)./T).*(1+t(i)./T))*dT+T_n;(i)=abs(T_x(i)-h(i))/dT*100;_pohubka=max(delta)
Результатом Виконання програми є: _pohubka=2.4529 - max відно...
