
- •В.И. Королев, в.В. Сахаров, о.В. Шергина оценка параметров, моделирование динамических систем и электрических цепей в среде 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.
2.5. Модель твердого тела с двумя неподвижными точками
Классические модели статики твердых тел и механических конструкций при воздействии на них произвольной системы сил хорошо известны из курса теоретической механики. Их построение базируется на пяти аксиомах статики и системе правил, определяющих необходимые и достоверные условия равновесия сил в пространстве относительно выбранной системы координат.
В общем случае для равновесия системы сил, приложенных к твердому телу, необходимо и достаточно, чтобы сумма проекций всех этих сил на каждую из трех произвольно выбранных координатных осей равнялась нулю и чтобы сумма их моментов относительно каждой из этих осей также равнялась нулю.
Получив систему уравнений, отвечающую условиям равновесия, мы можем представить ее в векторно-матричной форме и решить в среде MatLAB.
Рассмотрим следующую задачу.
На горизонтальном
валу трансмиссии (рис. 2.11) находятся два
жестко соединенных с ним шкива. Радиусы
большого и малого шкивов, соответственно,
равны
и
.
Известны силы натяжения ветвей ремня,
охватывающего передний шкив,
и
.
Сила
вертикальна, а сила
,
как следует из рис. 2.11 б), образует с
вертикалью угол
.
Натяжения ветвей ремня, охватывающего
второй шкив, равны
и
,
причем силы
и
параллельны и образуют с вертикалью
угол
.
Необходимо
определить реакции подшипников
и
,
перпендикулярные к оси вращения вала,
и модули сил
и
,
при условии, что
.
Реакции подшипников и силы
,
,
,
уравновешиваются (вал вращается
равномерно). Расстояние между шкивами
и расстояния шкивов от подшипников
приведены на рис. 2.3 а). Весом вала и
шкивов пренебрегаем.
Решение задачи
начнем с выбора начала координат (точки
)
и направления осей
,
,
.
Силы
,
,
,
перпендикулярны оси
.
Cоставляющие по осям
и
реакций подшипников в точках
и
обозначим, соответственно, через
,
и
,
.
Таким образом, определению подлежат
шесть неизвестных:
,
,
,
,
,
.
Для их нахождения необходимо составить
шесть уравнений: пять уравнений равновесия
и шестое уравнение – соотношение
,
данное в условии задачи. Заметим, что
все силы перпендикулярны к оси
;
их проекции на эту ось равны нулю. Поэтому
уравнение проекций на ось
не составляется. Проекции сил на оси
и
представим в виде уравнений:
Рис. 2.11. Схема трансмиссии
(0.0)
(0.0)
Составим уравнения
моментов всех сил относительно каждой
из координатных осей. Моменты сил,
приложенных в точках
и
,
относительно оси
равны нулю, так как линии действия этих
сил пересекают ось. Рассмотрим рис. 2.3
б). Чтобы получить моменты сил
,
,
,
относительно оси
,
достаточно модуль каждой из этих сил
умножить на радиус соответствующего
шкива и взять их произведения с
соответствующим знаком:
(2.0)
Моменты сил
относительно оси
определим по формуле:
(2.0)
Моменты сил
относительно оси
:
(2.0)
Исходные данные: |
|
|
|
|
|
|
|
Введем вектор
искомых переменных:
После подстановки
исходных данных в формулы (2.492.53),
с учетом равенства
,
получим систему уравнений:
.
Модель трансмиссии в векторно-матричной форме имеет вид:
, (2.0)
где
,
.
Для решения системы (2.54) составим простую программу, приведенную ниже.
Файл
%File ‘sah21.m’. example from theoretical mechanics.
A= [1 1 0 0 1.5
0 0 1 1 –1.5*sqrt(3)
0 0 0 –1.5 1.5* sqrt(3)
0 0 0 0 –0.2
0 1.5 0 0 1.5]
pause,
B=[–100*sqrt(2); 340+100*sqrt(2); –85–25*sqrt(2); –56; –25*sqrt(2)]
pause
C=inv(A)
pause
X=C*B
Решение задачи, выполненное на компьютере в среде MatLAB, выведено на печать с экрана дисплея путем нажатия клавиши <Print Screen>.
Здесь вектор
,
матрица
,
в чем легко убедиться, если умножить
на
.
Рассмотренная задача может быть решена в режиме прямых вычислений. С этой целью, в порядке упражнения, выполните операции, предусмотренные приведенной программой.
п sah21
A =
1.0000 1.0000 0 0 1.5000
0 0 1.0000 1.0000 –2.5981
0 0 0 –1.5000 2.5981
0 0 0 0 –0.2000
0 1.5000 0 0 1.5000
B =
–141.4214
481.4214
–120.3553
–56.0000
–35.3553
C =
1.0000 0 0 2.5000 –0.6667
0 0 0 5.0000 0.6667
0 1.0000 0.6667 –4.3301 0
0 0 –0.6667 –8.6603 0
0 0 0 –5.0000 0
X =
–257.8511
–303.5702
643.6716
565.2111
280.0000
п