
- •В.И. Королев, в.В. Сахаров, о.В. Шергина оценка параметров, моделирование динамических систем и электрических цепей в среде 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)