Студопедия — Если ввести аппроксимацию производной
Студопедия Главная Случайная страница Обратная связь

Разделы: Автомобили Астрономия Биология География Дом и сад Другие языки Другое Информатика История Культура Литература Логика Математика Медицина Металлургия Механика Образование Охрана труда Педагогика Политика Право Психология Религия Риторика Социология Спорт Строительство Технология Туризм Физика Философия Финансы Химия Черчение Экология Экономика Электроника

Если ввести аппроксимацию производной






, (2.4.15) при разбиении времени t на малые отрезки D t, а также использовать дискретизацию времени t = k D t по целочисленным отсчетам k = 0, 1, 2, …, то уравнение (2.4.4) можно записать в рекурсивной форме

x (k +1) @ (I + D t A) x (k) + D t B u (k). (2.4.16)

Точность решения уравнения (2.4.16) повышается с уменьшением D t.

Метод дискретизации уравнения состояния оказывается чрезвычайно полезным при вычислении временных характеристик нелинейных систем. В этом случае уравнение состояния имеет общий вид

, (2.4.17)

где f есть функция вектора состояния x и вектора входа u.

Используя (2.4.15) легко получить

x (k +1) = x (k) + D t f [ x (k), u (k), k ]. (2.4.18)

2.5. Идентификация и устойчивость дискретных моделей

линейных систем.

Предположим, что наблюдается некоторый дискретный сигнал yn. Рассмотрим его дискретную модель

yn = a10 yn-1 + a20 yn-2 + a30 yn-3 + + hn,

n = 0, 1, 2, …, N –1. (2.5.1)

В данной модели предполагается, что входной сигнал системы

xn , z =

имеет некоторую задержку z и описывается набором L импульсов dn с неизвестными амплитудами D xn, также подлежащими идентификации.

Будем считать, что помеха hn имеет нормальный гауссовый характер и некоррелирована (см. Модуль 1).

Будем искать параметры a10, a20, a30 и D xn дискретно-разностной динамической модели (2.51) методом наименьших квадратов, минимизируя энергию помехи

СКО = = [ yn a1 yn-1 a2 yn-2a3 yn-3

]2Þ min. (2.5.2)

Тогда, из условия ¶ОШ/¶ aij = ¶ОШ/¶D xk = 0, получим следующую систему 3-х линейных уравнений для оценивания параметров динамической модели

Z 11 a 1 + Z 12 a 2 + Z 13 a 3= R 1, (2.5.3)

Z 21 a 1 + Z 22 a 2+ Z 23 a 3= R 2,

Z 31 a 1 + Z 32 a 2 + Z 33 a 3 = R 3, где

Z 11 = = (2.5.4)

= (y0 y0 + y1 y1 +…+ yN-1 yN-1) – (yz-1 yz-1 + yz-2 yz-2 +…+ yz-L-1 yz-L-1),

Z12 = =

= (y1 y0 + y2 y1 +…+ yN-1 yN-2)– (yz-1 yz-2 + yz-2 yz-3 +…+ yz-L-1 yz-L-2),

Z13 = =

= (y2 y0 + y3 y1 +…+ yN-1 yN-3)– (yz-1 yz-3 + yz-2 yz-4 +…+ yz-L-1 yz-L-3),

Z21 = = Z12,

Z22 = =

= (y0 y0 + y1 y1 +…+ yN-2 yN-2) ‑ (yz-2 yz-2 + yz-3 yz-3 +…+ yz-L-2 yz-L-2),

Z23 = =

= (y1 y0 + y2 y1 +…+ yN-2 yN-3) – (yz-2 yz-3 + yz-3 yz-4 +…+ yz-L-2 yz-L-3),

Z31 = = Z13,

Z32 = = Z23,

Z33 = =

= (y0 y0 + y1 y1 +…+ yN-3 yN-3) – (yz-3 yz-32 + yz-4 yz-4 +…+ yz-L-3 yz-L-3).

Правые части системы уравнений (2.5.4) определяются как

R1 = = (2.5.5)

= (y1 y0 + y2 y1 +…+ yN yN-1)–(yz yz-1 + yz-1 yz-2 + yz-2 yz-3 +…+ yz-L yz-L-1),

R2 = =

= (y2 y0 + y3 y1 +…+ yN yN-2) –(yz yz-2 + yz-1 yz-3 + yz-2 yz-4 +…+ yz-L yz-L-2),

R3 = =

= (y3 y0 + y4 y1 +…+ yN yN-3)–(yz yz-3 + yz-1 yz-4 + yz-2 yz-5 +…+ yz-L yz-L-3).

Из (2.5.3) следует, что Z21 = Z12, Z31 = Z13, Z32 = Z23. Поэтому детерминант системы и ее алгебраические дополнения соответственно равны

Detz = Z11 Z22 Z33 Z11 Z232 Z22 Z132Z33 Z122 + 2 Z12 Z13 Z23 ,

Deta1 = R1 (Z22 Z33 Z232) + (2.5.6)

+ R2 (Z13 Z23Z12 Z33) + R3 (Z12 Z23Z13 Z22),

Deta2 = R1 (Z13 Z23 Z12 Z33) +

+ R2 (Z11 Z33Z132) + R3 (Z12 Z13Z11 Z23),

Deta3 = R1 (Z12 Z33 Z13 Z22) +

+ R2 (Z12 Z13Z11 Z23) + R3 (Z11 Z22Z122).

Параметры модели оцениваются как

a1 = Deta1 / Detz = (2.5.7)

= ;

a2 = Deta2 / Detz =

= ;

a3 = Deta3 / Detz =

= .

 

Коэффициенты D xn, вычисляются по оцененным параметрам модели как

D xn = yna10 yn-1a20 yn-2a30 yn-3, (2.5.8)

n = 0, 1, 2, …, N -1.

Спрашивается – как произвести оценку устойчивости полученной модели (2.5.1).

Возможно два подхода.

Теоретическая устойчивость дискретно-разностной модели (2.5.1) проверяется на основании известных в теории автоматического регулирования условий Рауса-Гурвица

a10 + a20 + a30 < 1,

a10a20 + a30 > –1, (2.5.9)

a302 + 3a102 + 3a20 (a30a10) – 4a30 a10 +

+ 9a30 + 8a209a10 > – 8

или более «мягких» условий

a10 + a20 + a30 < 1, (2.5.10)

a10a20 + a30 > –1,

3(a30 + 1) + a20 > a1,

3(a30 –1) < a10 + a20.

Практическая устойчивость дискретно-разностной модели (2.5.1) проверяется по условию

ОШмод / Энергия сигнала = [ ynzn ] 2 / yn 2 < D , (2.5.11)

где D – некоторый порог устойчивости (D £ 1); а zn – модель, построенная в соответствии с (2.5.1) для оцененных параметров.

2.6. Экспертные корреляционно-экстремальные

системы управления для спектральной

компьютерной квалиметрии.

Корреляционное сравнение данных.

При контроле качества технологических процессовили самих изделий необходимо проводить сравнение наблюдаемых сигналов (данных) X = [ X 1, X 2, …, XN ] c некоторой конечной совокупностью { S m } = {[ X 1 m, X 2 m, …, XNm ]} опорных сигналов (m = = 1, 2, …, M). На основании результатов сравнения вырабатывается управляющее воздействие u (X, { S m }). В таких системах основная тяжесть задачи падает на синтез методов, алгоритмов и систем сравнения наблюдаемых и эталонных сигналов.

Для примера, на рис. 2.6.1 приведена система контроля качества жидких сред по оптическому светорассеянию. Луч лазера Л освещает кювету Кс исследуемой жидкостью. Рассеянный свет попадает на детектор Д, формирующий спектр X рассеяния (угловой, цветовой, частотный, оптический, инфракрасный, комбинационный – в зависимости от конкретной задачи и исследуемой жидкости). Данный спектр поступает в устройство сравнения УС (аналоговое, цифровое, цифроаналоговое), куда поступают эталонные спектры S m, хранящиеся в базе данных БД. Само УС и БД могут быть выполнены как на основе цифровой техники, так и аналоговой голографической техники, проводить сравнение последовательно или параллельно. УС вырабатывает сигналы сравнения mm = m (X,{ S m }), на основе которых блок управления БУ формирует управление u ({ mm }).

 

БД

 


Рис. 2.6.1

Как видно из рисунка, управление, вырабатываемое в контуре обратной связи, осуществляет специализированный поиск в БД. Например, в простейшем случае ищется максимальное значение mm.

Предположим, что наблюдаемые спектральные данные X описываются аддитивной моделью в виде

X = S + H, (2.6.1)

где S – вектор-столбец опорного сигнала, а H – вектор-столбец аддитивных нормально распределенных нестационарных помех с некоррелированными компонентами, имеющими нулевые математические ожидания á H ñ = 0 и дисперсии sn2 = á Hn2 ñ.

Оказывается, в зависимости от конкретного вида модели (2.6.1) возможны различные методы оптимального сравнения неизвестных и опорных данных. При этом, оптимальность понимается в смысле экстремума критерия различительной способности при сравнении.

Рассмотрим следующие случаи применения оптимальных корреляторов – алгоритмов, реализующих сравнение данных.

Амплитуда сигнала и мощности помех á H2n ñ известны.

Тогда УС формирует корреляционную меру сходства

m (X, S ú F) = S T F T FX = å Nn (F) n 2(S) n (Х) n, (2.6.2) (7)

(F) n = 1/ sn = 1 / á H2n ñ1/2. (8)

Амплитуда сигнала известна, а мощности помех á H2n ñ

неизвестны.

Тогда УС формирует корреляционную меру сходства

m (X, S ú F) = S T F T FX / [(X S) T F T F (X S) + 1], (2.6.3)

(F) n = 1/ sn *, где sn *2 оценки мощности помех.

Амплитуда сигнала и мощности помех неизвестны.

Тогда УС формирует корреляционную меру сходства

m (X, S ú F) = 1/ [1 (S T F T FX)2 /(S T F T FS) (X T F T FX)], (2.6.4) (10)

(F) n = 1 / sn *, где sn *2 оценки мощности помех.

Методы спектральной компьютерной квалиметрии.

Практически все стандарты по сертификации качества продукции различных отраслей пищевых производств включают мероприятия по оценке органолептических свойств (признаков или характеристик) пищевых сред – сырья и готовых продуктов, которые гораздо больше, чем физико-химические свойства и пищевая ценность, влияют на выбор потребителей и, в конечном счете, формируют их спрос. Органолептическими свойствами пищевых продуктов являются внешний вид, текстура, вкус и запах, аромат. К сожалению, осуществить грамотную оценку органолептических свойств пищевых продуктов, тем более оперативную и в условиях массового производства, могут далеко не все, даже подготовленные эксперты. Необходимы некоторые «эталонные» качества, основанные на опыте и природных способностях. В то же время, создать подобное количество датчиков различных органолептических свойств, какое имеется у органов осязания и обоняния человека практически невозможно. В этих условиях и была предложена информационная технология спектральной компьютерной квалиметрии (СКК).

Идея СКК основана на разработке специальных или же использовании уже известных инструментальных средств (методов, датчиков и приборов) оперативного контроля разнообразных физико-химических, оптических, реологических и биологических характеристик пищевых сред совместно с автоматизированными (компьютерными) экспертными системами для связывания «спектров» (совокупностей) получаемых инструментальных характеристик с соответствующими органолептическими признаками и другими (потребительскими) характеристиками, оцениваемыми или инструментально измеряемыми опытными экспертами. В дальнейшем, с накоплением данных и знаний, такие экспертные системы должны самостоятельно давать оперативные оценки органолептических и потребительских свойств исследуемых пищевых сред, как отклик на входные запросы в виде соответствующих инструментальных характеристик данных сред. Реализация же СКК требует привлечения современных «синтетических» знаний, а также разработки специальных методов и средств.

Идея СКК наглядно показана на рисунке 2.6.1.

 


Рис. 2.6.1

СКК особенно эффективна в области управления качеством продукции, основанной на органических средах (растительные и минеральные масла, нефтепродукты), а также вин и ликероводочной продукции, что связано, в первую очередь, с возможностью практического осуществления оперативного спектрального анализа данных жидких сред в широком диапазоне электромагнитных волн.

Для оперативного контроля разнообразных сред представляют интерес следующие методы получения спектральных данных.

1. Метод теплового удара мгновенное изменение локальной (поверхностной) температуры среды (например, излучением полупроводникового лазера) и измерение температурного отклика среды.

2. Метод акустического удара мгновенное возбуждение в среде акустических (ультразвуковых) шумов и измерение их спектральных характеристик.

3. Метод электровозбуждения возбуждение электротока в среде и измерение спектра мощности его флуктуаций.

4. Метод инфракрасной и СВЧ (субмиллиметровой, миллиметровой и сантиметровой) спектроскопии измерение спектральных свойств среды в диапазонах длин волн от (0,75¸80) мкм до (0,30¸30,00) смсоответственно.

5. Метод оптической спектроскопии измерение спектральных свойств среды в диапазонах длин волн от (0,20¸ 1000) мкм.

6. Методы жидкостной и газовой хроматографии – разделение вещества на фракции в соответствие их молекулярного состава.

В любом из перечисленных методов формируют спектр S (w) отклика исследуемых сред на внешнее воздействие (см. стр. 44 45). Данные спектры несут косвенную информацию о микро- и макроскопических, физико-химических и структурных свойствах сред. При этом количество M = W / D W потенциальных свойств определяется максимальной шириной W полосы отклика (0 £ w £ W) и разрешающей способностью D W метода.

В результате исследуемая среда с различными потребительскими свойствами 1 и 2 потенциально может откликаться близкими спектральными распределениями S 1(w) и S 2(w), носящими нечетко выраженные различия, как показано на рис. 2. 6.3. Приведенные спектральные распределения измерены с 10% точностью (± H). Поэтому необходимы специальные высокоразрешающие методы численной обработки для их различения.

 

Рис. 2.6.3

К сожалению, даже не сплошные, а линейчатые спектры пищевых сред носят нечеткий характер, что обусловлено их многокомпонентным составом. Так, на рис. 2.6.4 приведен спектр аминокислотного состава растительного масла, полученный хроматографическим методом, при котором последовательно во времени, путем диффузии, выделяются аминокислотные фракции в соответствии с их молекулярной массой аминокислоты с легкими молекулярными массами выделяются первыми.

Из рисунка видно, что хроматографические спектры при наличии всех аминокислот могут различаются только своими огибающими.

 
 

Причина нечеткости спектральных данных пищевых сред заключается в том, что относительно узкие резонансные линии их многочисленных компонентов накладываются друг на друга, что приводит к «смазыванию» спектральных распределений.

Рис. 2.6.4

Это легко пояснить, опираясь на теорию спектрального анализа, опирающегося на рассмотренную выше теорию частотных характеристик систем (см. стр. 44 45). Спектральное распределение функции отклика среды для «частиц» (атома, молекулы, молекулярных соединений) одного «к»-го сорта, соответствующее зависимости мощности поглощения от частоты w возбуждения, вблизи резонанса wк, описывается функцией Лоренца

Sк (w) @ S0 wк Dwк /[(wкw)2 + Dwк 2], (2.6.5)

где Dwк ширина поглощения частицы «к»-го сорта.

Отсюда следует, что спектральное распределение функции отклика среды, содержащей различные частицы (к = 1, 2, …, К) имеет вид

Så (w) = åк Sк (w) @ åк S0 wк Dwк / [(wкw)2 + Dwк 2]. (2.6.6)

Простые вещества (газы, жидкости), состоящие из небольшого количества разнотипных невзаимодействующих «частиц», имеют лорентцевские спектры показанные на рисунке 2.6.5.

 


Рис. 2.6.5

Пищевые среды, как и другие многокомпонентные вещества, находящиеся в конденсированном состоянии, состоят из огромного количества различных «частиц», лорентцевские спектры которых значительно уширены из-за связей, удерживающих данные частицы в едином комплексе.

В качестве наглядного примера на рис. 2.6.6. приведен спектр поглощения среды состоящей из частиц трех разных типов.

Видно, что из-за наложения лорентцевских спектральных линий и большой ширины каждой из них результирующий спектр приобретает «смазанный» или нечеткий характер. Вполне естественно, что спектральное разрешение уменьшается с уменьшением полосы частот используемых воздействий. С падением разрешения падает и способность различения близких спектральных распределений, особенно в нерезонансной области спектра поглощения.

 

 


Рис. 2.6.6

Основная проблема спектральной компьютерной квалиметрии – уверенное различение «близких» спектральных распределений. Для решения этой проблемы возможно использовать рассмотренные выше корреляторы.

Другой проблемой является формирование базы данных «опорных» спектров { S m = [ S 1 m, S 2 m, …, SNm ]} (m = 1, 2, …, M) и «эталонных» векторов { Y m = [ Y 1 m, Y 2 m, , YNm ]} описания свойств пищевых продуктов, а также – построение экспертных правил (знаний), на основании которых будет синтезироваться вектор отклика Y на неизвестный входной спектр X.

Построение экспертных правил возможно осуществить на основе моделей регрессионных зависимостей, минимизирующих СКО соответствующих выходов модели

Y k (X) = (å MmÁ;[ X, S m ]) - 1å Mm Y km Á;[ X, S m ], (2.6.6)

где Á;[ X, S m ] – ядро преобразования, определяемое на основании мер сходства X и S m, а также – статистических предположений о распределении «эталонных» образов { S m, Y m }.

Вопросы для самопроверки к главе 2.

1. Каково назначение ЦАП и АЦП в цифровых системах

управления? Изобразите функциональную схему

одноконтурной цифровой системы управления.

2. Что такое квантователь или ключ? Какое

функциональное преобразование он осуществляет?

3. Как работает ЦАП? Изобразите функциональную схему

ключа и фиксатора (экстраполятора нулевого уровня).

Чему равна передаточная функция экстраполятора?

4. Как выражается z -преобразование сигнала?

5. Как построить z -преобразование передаточной функции

разомкнутой дискретной системы?

6. Как определить отклик дискретной системы на

единичный импульсный входной сигнал, на единичный

ступенчатый входной сигнал?

7. Как определить устойчивость линейной непрерывной

системы с обратной связью по ее непрерывной

передаточной функции П (s)?

8. Как определить устойчивость линейной непрерывной

системы с обратной связью по ее дискретной

передаточной функции П (z)?

9. В чем принципиальная разница по устойчивости между

дискретной и непрерывной системами второго порядка?

10. Как получить дискретную аппроксимацию цифрового

регулятора?

11. Как получить уравнение непрерывной динамической

системы в переменных ее состояния?

12. Как описывается в общем случае состояние линейной

динамической системы в переменных ее состояния?

Приведите векторное описание.

13. Приведите общий вид решения уравнения состояния.

14. Как определяется фундаментальная или переходная

матрица состояния?

15. Выпишите вид линейного и нелинейного уравнений

состояний в дискретно-разностной форме.

16. Объясните смысл и СКО метод идентификации

дискретных моделей линейных систем.

17. Как провести оценку устойчивости дискретно-

разностной модели линейной системы.

18. Приведите функциональную схему системы контроля

качества жидких сред по оптическому светорассеянию.

19. Приведите примеры корреляционных мер сходства.

Чем они отличаются?

20. В чем заключаются методы спектральной

компьютерной квалиметрии?

21. Какие методы получения спектральных данных Вы

знаете?

22. Что лежит в основе теории спектрального анализа?

23. Почему спектры многокомпонентных веществ,

находящихся в конденсированном состоянии носят

нечеткий размытый характер?

24. Какие две проблемы спектральной компьютерной

квалиметрии Вы знаете?

25. Объясните смысл построения экспертных правил в

компьютерной квалиметрии.







Дата добавления: 2015-09-19; просмотров: 379. Нарушение авторских прав; Мы поможем в написании вашей работы!



Кардиналистский и ординалистский подходы Кардиналистский (количественный подход) к анализу полезности основан на представлении о возможности измерения различных благ в условных единицах полезности...

Обзор компонентов Multisim Компоненты – это основа любой схемы, это все элементы, из которых она состоит. Multisim оперирует с двумя категориями...

Композиция из абстрактных геометрических фигур Данная композиция состоит из линий, штриховки, абстрактных геометрических форм...

Важнейшие способы обработки и анализа рядов динамики Не во всех случаях эмпирические данные рядов динамики позволяют определить тенденцию изменения явления во времени...

Растягивание костей и хрящей. Данные способы применимы в случае закрытых зон роста. Врачи-хирурги выяснили...

ФАКТОРЫ, ВЛИЯЮЩИЕ НА ИЗНОС ДЕТАЛЕЙ, И МЕТОДЫ СНИЖЕНИИ СКОРОСТИ ИЗНАШИВАНИЯ Кроме названных причин разрушений и износов, знание которых можно использовать в системе технического обслуживания и ремонта машин для повышения их долговечности, немаловажное значение имеют знания о причинах разрушения деталей в результате старения...

Различие эмпиризма и рационализма Родоначальником эмпиризма стал английский философ Ф. Бэкон. Основной тезис эмпиризма гласит: в разуме нет ничего такого...

Метод архитекторов Этот метод является наиболее часто используемым и может применяться в трех модификациях: способ с двумя точками схода, способ с одной точкой схода, способ вертикальной плоскости и опущенного плана...

Примеры задач для самостоятельного решения. 1.Спрос и предложение на обеды в студенческой столовой описываются уравнениями: QD = 2400 – 100P; QS = 1000 + 250P   1.Спрос и предложение на обеды в студенческой столовой описываются уравнениями: QD = 2400 – 100P; QS = 1000 + 250P...

Дизартрии у детей Выделение клинических форм дизартрии у детей является в большой степени условным, так как у них крайне редко бывают локальные поражения мозга, с которыми связаны четко определенные синдромы двигательных нарушений...

Studopedia.info - Студопедия - 2014-2024 год . (0.014 сек.) русская версия | украинская версия