- •В.И. Королев, в.В. Сахаров, о.В. Шергина оценка параметров, моделирование динамических систем и электрических цепей в среде matlab
- •Санкт-Петербург
- •Рецензенты:
- •Isbn 5-88964-073-X © Королев в. И.,
- •Содержание
- •1. Идентификация параметров моделей систем на основе квадратичных методов оценивания………..……………………………………………………………........................... 11
- •Модели установившихся режимов в электрических цепях и
- •4. Модели пространства состояний в электрических цепях и
- •5. Системы и цепи под воздействием параметрических возмущений: численные и
- •6. Моделирование динамических систем и электрических цепей средствами
- •8. Моделирование технологических процессов в системах с помощью
- •Введение
- •Идентификация параметров моделей систем на основе квадратичных методов оценивания
- •Линейные модели систем Рассмотрим систему линейных уравнений
- •Рассмотрим линейное уравнение
- •Для существования точного решения уравнения (1.1) должно выполняться условие
- •1.2. Переопределенные системы уравнений. Метод наименьших квадратов
- •Введем векторы
- •1.3. Оценивание параметров периодических сигналов по экспериментальным данным
- •Таким образом, требуется оценить амплитуды и фазы
- •Разделив период t на 24 равных интервала, мы получим . Так как , то каждому будет соответствовать угол (в градусах).
- •1.4. Рекуррентный метод оценивания параметров моделей
- •1.5.Алгоритм оценивания параметров моделей в условиях ограничений
- •Минимальное значение (1.43) соответствует условию
- •1.6. Нормы оценивания параметров в теории инверсных систем
- •1.7. Оценивание параметров по критерию минимума среднего абсолютного значения остатков
- •Уравнение измерителя
- •1.8. Нелинейное оценивание параметров моделей в режиме прямых вычислений
- •1.9. Алгоритм ортогонализации в оценке параметров динамических систем
- •Рассмотрим нелинейную дискретную модель системы
- •2. Модели установившихся режимов в электрических цепях и системах
- •2.1. Пример модели цепи постоянного тока
- •Пример модели цепи переменного тока
- •2.3.Несимметричные режимы в трехфазных электрических цепях: метод симметричных составляющих
- •2.4. Модель многоузлового разветвления русла в стационарном режиме. Аналогии с электрической цепью
- •2.5. Модель твердого тела с двумя неподвижными точками
- •2.6. Многопродуктовая модель экономики: межотраслевой баланс
- •2.7. Модель оценки элементов матрицы преобразования
- •3. Динамические звенья и системы под воздействием гармонических сигналов
- •3.1.Построение частотных характеристик простых динамических звеньев
- •3.2.Резонанс в динамических системах
- •3.4. Резонанс в электрических цепях
- •3.4. Режим биений
- •4. Модели пространства состояний в электрических цепях и системах
- •4.1. Понятие состояния
- •4.2. Уравнения состояния для электрических цепей
- •4.3. Приведение к резистивной форме моделей электрических цепей с одним накопителем энергии
- •4.4. Определение коэффициентов дифференциальных уравнений моделей электрических цепей с двумя накопителями энергии
- •4.5. Электрическая цепь с четырьмя накопителями энергии
- •4.6. Примеры составления уравнений состояния rlc -цепей в матричной форме
- •Исходная цепь
- •Краткие сведения о решателях дифференциальных уравнений
- •5.2. Об аналитическом решении дифференциальных уравнений
- •5.3.Матричная форма решения уравнений состояния динамических систем
- •Моделирование динамических систем с двумя накопителями энергии при параметрических возмущениях Для определенности рассмотрим модель динамической системы с двумя накопителями энергии:
- •Моделирование –цепей. Аналитические решения
- •Сопротивления ,
- •Переходный процесс в -цепи под воздействием синусоидального входного сигнала
- •Системы, находящиеся под воздействием периодических сигналов прямоугольной формы
- •Динамические системы с двумя накопителями энергии
- •Корни комплексно-сопряженные.
- •Кратные корни.
- •Вещественные неравные корни
- •Динамические системы с тремя накопителями энергии
- •6. Моделирование динамических систем и электрических цепей средствами символьной математики
- •6.1. Символьные выражения и алгебра
- •6.2. Алгебраические и трансцендентные уравнения. Расчет цепи постоянного тока
- •Символьной алгебры
- •6.3. Символьное дифференцирование и интегрирование
- •6.4. Решение дифференциальных уравнений в символьной форме
- •6.5. Моделирование переходных процессов в электрических цепях средствами пакета символьной математики
- •6.6. Операторный метод расчета переходных процессов с использованием пакета символьной математики
- •7. Модели детерминированного хаоса и их применение
- •7.1. Дискретные динамические модели первого порядка со сложной динамикой
- •7.2. Программное обеспечение и моделирование нелинейных дискретных систем. Хаос
- •7.3. О применении нелинейных моделей систем для получения псевдослучайных хаотических последовательностей
- •Моделирование технологических процессов в системах с помощью нейронных сетей
- •8.1.Общие положения
- •8.2. Последовательность операций при создании нейронной сети в среде matlab ( Neural Networks Toolbox)
- •3. Оценка погрешности нейронной модели. Для оценки используется функция моделирования, где в скобках, согласно синтаксису, приводятся сеть и входной сигнал:
- •Моделирование уровней воды в водной коммуникации на основе нейронных сетей
- •В водной коммуникации
- •8.4. Применение нейронной сети для определения химического состава песчано – гравийной смеси
- •Моделирование технологического процесса оценки и прогноза рыночных факторов, воздействующих на работу предприятия и бизнес в классе нейронных сетей
- •Библиографический список
- •Оценка параметров, моделирование динамических систем и электрических цепей в среде MatLab
- •165300, Г. Котлас, ул. Невского, 20.
4.3. Приведение к резистивной форме моделей электрических цепей с одним накопителем энергии
Рассмотрим электрическую цепь, изображенную на (рис. 4.2).
В момент коммутации (при замыкании ключа «k»), согласно первому закону коммутации, ток в цепи мгновенно изменяться не может. Поэтому в момент t=0+ ток в цепи , и напряжение источника уравновешивается напряжением на индуктивности.
Согласно правилам, изложенным в предыдущем параграфе, заменим индуктивность источником тока.
а) б)
в) г)
Рис. 4.2. RL – цепь и ее резистивная форма
Обозначим этот ток (рис. 4.2 б). Цепь содержит один накопитель энергии (индуктивность) и один источник ЭДС (внутреннее сопротивление источника ЭДС равно нулю). Поэтому, согласно (4.10), динамическая модель цепи может быть представлена уравнением [23]:
(4.13)
Заметим, что эта модель справедлива также для любой разветвленной цепи, содержащей только одну индуктивность и один источник электроэнергии (например, источник тока). Однако выражения для коэффициентов а и b будут отличаться, поскольку они определяются топологией конкретной цепи.
Цепь, изображенная на (рис. 4.2 б), приемлема только для момента t = 0+ при условии, что изменяется скачком. Она не содержит инерционных элементов. Поэтому (1.3-1) можно заменить алгебраическим уравнением:
(4.14)
Постоянный коэффициент а в уравнении (4.14) определяет приращение напряжения на индуктивности при изменении (скачком) тока , а коэффициент b - приращение напряжения на индуктивности при изменении ЭДС источника на величину.
Используя принцип суперпозиции (наложения), мы можем искусственно вводить лишь приращение одной переменной, а вторую при этом «замораживать». Например,
, (4.15)
. (4.16)
Для удобства определения неизвестных коэффициентов приращения переменных можно принять равными единице. Если увеличение тока на одну единицу (например, на 1 Ампер) будет приводить к уменьшению напряжения на индуктивности на величину , то коэффициент а должен приниматься со знаком «минус», а если к увеличению - со знаком «плюс». Аналогично определяется знак коэффициента b.
Согласно (4.15), «заморожено». Поэтому в резистивной форме цепи (рис.4.2 в) источник ЭДС отсутствует, однако сохранено его внутреннее сопротивление, равное нулю. Если в момент t = 0+ с помощью источника тока генерировать скачком ток, это вызовет падение напряжения на сопротивлении R, равное . С учетом принятых на схеме знаков, напряжение на источнике тока уменьшится на величину . Следовательно, коэффициент а по модулю равен , а с учетом знака приращения напряжения, он должен быть принят со знаком «-», то есть а = - R.
Теперь остановимся на уравнении (4.16). В момент времени t = 0+ ему соответствует резистивная цепь, приведенная на рис.4.2 г). Источник тока из цепи исключен, а его внутреннее сопротивление, равное бесконечности, оставлено.
Заметим, что при изменении скачком (на 1 вольт) напряжения , оно будет полностью приложено к источнику тока (ток в цепи отсутствует). С учетом знаков сигналов, напряжение на источнике тока возрастет до значения . Следовательно, коэффициент b = 1.
Дифференциальное уравнение (4.13), в отличие от алгебраического, позволяет описать процесс в любой момент времени 0+ t < , в том числе - в момент t = 0+. Поэтому полученные коэффициенты могут быть перенесены в (4.13). Тогда будем иметь:
(4.17)
Запись (4.17) в форме Коши:
(4.18)
показывает, что матрицы А и В модели (4.4) вырождаются в элементы
Рассмотрим цепь, содержащую одну индуктивность и источник тока (рис.4.3 а). Ее резистивная схема представлена на (рис. 4.3 б).
Модель цепи (рис.4.2) сохраняется, с учетом того, что источником электроэнергии является источник тока
(4.19)
а) б)
Рис.4.3. RL – цепь с источником тока и ее резистивная форма.
Для резистивной цепи (рис. 4.3 б) справедливо алгебраическое уравнение
(4.20)
Генерируя при , мы получим: и аналогично будем иметь .
Таким образом, динамическая модель цепи, определяемая (4.19), представляется дифференциальным уравнением
или в форме Коши:
(4.21)
На рис. (4.4 а) приведена RC-цепь, подключаемая к источнику постоянной ЭДС. Возвращаясь к (4.10), мы видим, что при наличии одного накопителя энергии (емкости) и одного источника, независимо от числа и схемы соединений активных сопротивлений, модель цепи должна иметь вид:
(4.22)
а) б)
в) г)
Рис. 4.4. RC – цепь и ее резистивная форма.
Для получения значений а и b необходимо, согласно правилам, определить резистивную форму цепи. С этой целью выберем направление тока IC , генерируемого источником, и заменим емкость источником ЭДС. Направление ЭДС этого источника должно быть встречным току IС. Использование источника ЭДС вместо емкости в момент замыкания ключа t =0+ обосновывается тем, что напряжение на емкости, согласно второму закону коммутации, не может изменяться скачком. В момент замыкания ключа k в цепи устанавливается максимальный ток (при условии, если ).
Резистивная форма RС-цепи приведена на рис. (4.4 б). Левая часть (4.22) имеет размерность тока. Поэтому, переходя от дифференциального уравнения к алгебраическому, для резистивной цепи (рис. 4.4 б) мы получим:
(4.23)
Для определения коэффициентов а и b будем поочередно варьировать одну из двух переменных, оставляя другую без изменения. Предположим, что
.
Коэффициент а определяет, во сколько раз изменится ток через емкость, если напряжение на емкости изменится на одну единицу (например, 1 Вольт). Знак коэффициента будем принимать отрицательным, если под действием напряжения будет генерироваться ток, направленный встречно току IС .
При совпадении направлений генерируемого тока через емкость и тока IC соответствующий коэффициент будем записывать со знаком «+».
Согласно обозначениям, используемым на рис. (4.4 в), приращение тока по модулю равно:
Так как направление тока является встречным относительно IС , знак коэффициента должен быть «-», то есть .
Рассмотрим цепь, представленную на рис. (4.4 г), где источник ЭДС исключен, но в цепи содержится его внутреннее сопротивление, равное нулю. Под действием в цепи возникает ток .
Поскольку этот ток совпадает по направлению с током IС , знак коэффициента b, устанавливающего соотношение между и , должен быть «+», то есть:
.
Коэффициент .
Динамическая модель цепи (1.3-10) принимает вид:
,
и для переменной состояния UС (t) мы получаем стандартную форму записи:
(4.24)