- •В.И. Королев, в.В. Сахаров, о.В. Шергина оценка параметров, моделирование динамических систем и электрических цепей в среде 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.
6.4. Решение дифференциальных уравнений в символьной форме
Уравнение, содержащее производные неизвестной функции или неизвестных функций, называется дифференциальным уравнением. Вместо производных в уравнение могут входить дифференциалы.
В процессе моделирования динамических систем и электрических цепей часто приходится решать уравнения двух видов: обыкновенные дифференциальные уравнения и дифференциальные уравнения в частных производных.
Остановимся на рассмотрении решения в символьной форме обыкновенных дифференциальных уравнений, для которых характерна зависимость неизвестной функции от одного аргумента. Решение выполняется с помощью функции dsolve, которая возвращает результат в символьном виде. Согласно синтаксису функции, запись первой производной осуществляется посредством символа D, для второй производной используется D2, для третьей - D3 и т.д. В этой связи при объявлении символьных переменных не следует использовать букву D, если далее предполагается решать уравнения с помощью dsolve.
При выводе результата в решении должны содержаться постоянные интегрирования, которые представлены коэффициентами , и т.д. Их число равно порядку дифференциального уравнения. Порядком дифференциального уравнения называют порядок наивысшей производной, входящей в это уравнение.
В рассмотренных ниже примерах приводятся решения дифференциальных уравнений, выполненные с соблюдением синтаксиса функции-решателя dsolve. В первых двух примерах приведены решения линейных уравнений первого порядка без правой части и с правой частью. В третьем примере найдено общее решение линейного уравнения второго порядка с постоянными коэффициентами без правой части при вещественных корнях. В четвертом примере корни характеристического уравнения комплексно-сопряженные. Поэтому решение состоит из экспоненциальной и периодической составляющих. Разложение функции на простые множители выполнено в четвертом примере с помощью оператора факторизации factor, а решение представлено в формате, приближенном к математическому (использована функция pretty). Пятый пример содержит решение дифференциального уравнения второго порядка без правой части при мнимых корнях характеристического уравнения. Поэтому общее решение состоит из двух периодических составляющих, умноженных на постоянные интегрирования. В шестом и седьмом примерах представлен синтаксис решения линейного уравнения второго порядка с правой частью. В первом случае правая часть уравнения представляет собой экспоненту, а во втором – полином, приводящий к общему решению, состоящему из суммы экспоненциальной и полиномиальной составляющих. Наконец, в восьмом примере представлено общее решение дифференциального уравнения пятого порядка. Его характеристическое уравнение содержит однократный действительный корень и пару двукратных сопряженных мнимых корней. В решении в символьной форме представлены все составляющие, определяемые составом корней: экспоненциальная составляющая, периодические и множитель , появляющийся вследствие кратности мнимых корней.
Очевидно, аналогично можно получить решения для линейных уравнений любого порядка, в том числе – при наличии правых частей.
Примеры решения дифференциальных уравнений.
Пример1. Уравнение без правой части.
syms x t y
dsolve('Dy-t/(1+t^2)*y=0')
ans =
C1*(1+t^2)^(1/2)
Пример 2 . Уравнение с правой частью.
dsolve('Dy-t/(1+t^2)*y=t')
ans =
1+t^2+C1*(1+t^2)^(1/2)
Пример 3 . Уравнение второго порядка.
dsolve('D2y=-1/4*Dy+3/8*y')
ans =
C1*exp(-3/4*t)+C2*exp(1/2*t)
Пример 4 . Комплексно-сопряженные корни
характеристического уравнения.
f1=dsolve('D2y+Dy+y=0');
f2=factor(f1);
pretty(f2)
1/2 1/2
exp(- 1/2 t) (sin(1/2 3 t) C1 + C2 cos(1/2 3 t))
Пример 5. Мнимые корни характеристического
уравнения.
dsolve('D2y+y=0')
ans =
C1*cos(t)+C2*sin(t)
Пример 6 . Уравнение с правой частью.
dsolve('D2y-1/2*Dy-1/2*y=3*exp(1/2*t)')
ans =
-6*exp(1/2*t)+C1*exp(t)+C2*exp(-1/2*t)
Пример 7. Уравнение с правой частью, заданной полиномом.
dsolve('D2y-3*Dy=t^2+3*t')
ans =
-11/27*t-11/18*t^2-1/9*t^3+C1+C2*exp(3*t)
Пример 8 .Уравнение пятого порядка.
dsolve('D5y+D4y+2*D3y+2*D2y+Dy+y=0')
ans =
C1*cos(t)+C2*sin(t)+C3*exp(-t)+C4*cos(t)*t+C5*sin(t)*t
Для решения системы уравнений с помощью функции dsolve требуется использовать соответствующий синтаксис:
.
Функция возвращает решение в общем виде. Оно содержит соответствующее число постоянных интегрирования. Однако, если известны граничные условия, то их следует использовать для исключения постоянных интегрирования и получения из множества общих одного частного решения. Граничные условия записываются по следующей схеме:
,
Приведем несколько примеров.
Предположим, что необходимо в символьных выражениях решить систему дифференциальных уравнений
,
.
Собственные значения матрицы состояния eig(A) являются комплексно-сопряженными числами (), и решение имеет вид:
,
.
Аналогичный результат получим в пакете символьной математики (после объявления символьных величин):
>> [x,y]=dsolve('Dx=3*x+5*y','Dy=-5*x+3*y')
x =
exp(3*t)*(cos(5*t)*C1+sin(5*t)*C2)
y =
-exp(3*t)*(sin(5*t)*C1-cos(5*t)*C2)
>>
Предположим, что необходимо определить и при следующих граничных условиях , . В этом случае
,
.
Интегрирование дифференциальных уравнений с использованием dsolve при заданных начальных условиях позволяет получить тот же результат, что можно видеть из расположенного ниже текста.
>> [x,y]=dsolve('Dx=3*x+5*y','Dy=-5*x+3*y','x(0)=2','y(0)=0')
x =
2*exp(3*t)*cos(5*t)
y =
-2*exp(3*t)*sin(5*t)
Одним из замечательных свойств пакета символьной математики является возможность решения дифференциальных уравнений при расщепленных граничных условиях. Иначе говоря, решение может быть получено не только при заданных начальных условиях (в момент ), но и при других граничных условиях. Например, переменная может быть задана в момент , а в момент . В приведенном ниже фрагменте расчеты выполнены для двух граничных точек, расположенных на тех же траекториях и в моменты времени и : , .
>> [x,y]=dsolve('Dx=3*x+5*y','Dy=-5*x+3*y','x(0.5)=-7.1810','y(0)=0')
x =
-7181/1000*exp(3*t)*cos(5*t)/cos(5/2)/exp(3/2)
y =
7181/1000*exp(3*t)*sin(5*t)/cos(5/2)/exp(3/2)
>> -7181/1000/cos(5/2)/exp(3/2)
ans =
2.0000
>>
Граничные условия могут задаваться в символьной форме. Предположим, что в момент переменные состояния системы и . Тогда система уравнений имеет решение
,
,
что также подтверждается приведенным ниже фрагментом символьных вычислений:
>> [x,y]=dsolve('Dx=3*x+5*y','Dy=-5*x+3*y','x(0)=A','y(0)=B')
x =
exp(3*t)*(cos(5*t)*A+sin(5*t)*B)
y =
-exp(3*t)*(sin(5*t)*A-cos(5*t)*B)
>>
Построение графиков символьных функций осуществляется с помощью оператора ezplot, если функции не содержат таких неопределенных коэффициентов, например, как постоянные интегрирования , и т.п. По умолчанию, синтаксис ezplot(t) возвращает график заданной символьной функции в интервале [-2pi, 2pi] изменения независимой переменной . Диапазон изменения независимой переменной может быть задан пользователем. В этом случае применяется запись
,
позволяющая построить графики при изменении от минимального значения до максимума с нанесением координатной сетки. Примеры построения графиков приводятся в следующем параграфе.
В заключение отметим, что в символьной форме может решаться определенный класс нелинейных уравнений (в основном, первого порядка, если существует решение в замкнутом виде). В случае отсутствия решения, на дисплее появляется текст, информирующий пользователя о невыполнимости предписанных операций.