- •Математическоеикомпьютерноемоделированиепроцессовисистемвсредеmatlab/simulink.
- •Содержание
- •6.Упражнениядлясамостоятельнойработы.....................................................67
- •7.Иллюстрациикупражнениям..........................................................................71
- •8.Литература.........................................................................................................89 введение
- •1.Общиесведенияопрограммномкомплексеsimulink
- •1.1.Решающиеэлементы
- •1.2.Построениепростейшихмоделей
- •3.Аппроксимациясигналов
- •2.Моделиалгебраическихобъектов
- •3.1.Методнаименьшихквадратов
- •3.2.Модифицированныйметодравныхплощадей
- •1].Интервалразбиваетсянапятьподинтервалов[0;0,1],[0;0,2],[0;
- •3.3.Блочно-импульснаяаппроксимация
- •4.Моделидинамическихобъектов
- •5.Обратимостьмоделейизадачиоптимизации
- •5.1.Обратимостьрешающихэлементов
- •5.2.Виртуальныерешающиеэлементысистемы
- •5.3.НеобратимыйлинейныйпреобразовательвсистемеSimulink
- •5.5.Виртуальныйаналогоперационногоусилителя
- •5.6.Обратимыйлинейныйпреобразователь
- •5.7.Обратимаямодельсистемылинейныхалгебраическихуравнений
- •5.8.Обратимаямодельсистемылинейныхнеравенств
- •5.9.Обратимаямодельзадачилинейногопрограммирования
- •5.10.Модельтранспортнойзадачилинейногопрограммирования
- •6.Упражнениядлясамостоятельнойработы
- •6.1.Сформироватьивизуализироватьсигналызаданнойформы.
- •6.2.Сформироватьструктурныесхемы.
- •6.3.Построитьвиртуальныемодели.
- •7.Иллюстрациикупражнениям
- •8.Литература
- •1.АнуфриевИ.Е.,СмирновА.Б.,СмирноваЕ.Н.Matlab7.–сПб.:
- •4.А.С.304600ссср.Квазианалоговоемоделирующееустройство/
- •6.АхиезерН.И.Лекциипотеорииаппроксимации.–м.-л.:Гос.Изд-во
- •12.ДьяконовВ.П.,КругловВ.Математическиепакетырасширения
- •14.ДэбниДж.,ХартманТ.Simulink4.Секретымастерства.–м.:
- •15.ЗагидуллинР.Ш.LabViewвисследованияхиразработках.–м.:
2.Моделиалгебраическихобъектов
Ваналоговойвычислительнойтехникесуществуетнесколькоспособовмоделированияобъектов,описываемыхсистемамиалгебраическихуравнений[22–25].Одинизних,наиболеераспространенный,сводитсякрешениюсистемыобыкновенныхдифференциальныхуравнений,установившеесярешениекоторойдаетрешениеалгебраическойзадачи.
Пустьбезинерционныйобъектописываетсясистемойуравнений
a11x1+a12x2+...+a1nxn=b1
a21x1+a22x2+...+a2nxn=b2
L
an1x1+an2x2+...+annxn=bn
Дляпостроениямоделиданногообъектаприменимметодсведениякэквивалентнойсистемедифференциальныхуравнений.
Введемследующуюсистемудифференциальныхуравнений:
dx1+ax
+ax
+...+ax
−b=0
dt 111 122 1nn 1
dx2+ax
+ax
+...+ax
−b=0
dt 211 222 2nn 2
L
dxn
ax ax ax b
dt+n11+
n22+...+
nnn−
n=0
Кактольковсепроизводныезатухнут(dxi
dt
=0),получаемсигнал
решения{x1,x2,...,xn}.
Дляэквивалентностиэтихдвухсистемуравненийдолжнообеспечиватьсязатухающеерешениесистемыдифференциальныхуравнений.Достаточнымусловием,обеспечивающимзатухающеерешение,являетсяположительнаяопределенностьматрицыкоэффициентовлинейнойсистемыуравнений.Этовозможно,в
n
частности,приусловии,когдаaii≥∑aij,i≠j.
j=1
Пример2.1.
Найтирешениесистемылинейныхалгебраическихуравнений:
4x1+2x2=14.
2x1+5x2=−5
Перейдем к эквивалентной системе дифференциальныхуравнений:
dx1=14−4x
−2x
dt 1 2
dx2=−5−2x
−5x
dt 1 2
Структурнаясхемамоделиданнойсистемыприведенанарис.
2.1.Онапостроениасиспользованиемклассическихметодованалоговойвычислительнойтехники.Переходныйпроцессустановлениярешенияизображеннаэкраневиртуального
осциллографа(рис.2.2).
Рис.2.1.Структурнаясхемамоделисистемыдифференциальныхуравнений,
эквивалентнойсистемелинейныхалгебраическихуравненийвторогопорядка
Нарис.2.2видно,чтопосле
t=2
навыходахвиртуальных
интеграторовустанавливаютсясигналы,соответствующиерешениюсистемылинейныхалгебраическихуравнений:
x1=5,
x2=−3.
Расширениеклассаалгебраическихобъектов,матрицакоэффициентовкоторыхимеетпроизвольныйвид,возможнопутемиспользованияметодовквазианалоговогомоделирования[22,23].Примерытакихметодовбудутрассмотреныпозднее.
Рис.2.2.Переходныйпроцессустановлениярешениясистемылинейныхалгебраическихуравненийпутемсведениякэквивалентнойдифференциальнойсистемеуравнений
Пример3.2.
Найтирешениесистемы:
5x1+2x2+x3=2,5
3x1+7x2+2x3=−1,5.
4x1+0,5x2+5x3=11,5
Модельсистемыирезультатыееработыпредставленынарис.
2.3,2.4соответственно.
Рис.2.3.Модельсистемылинейныхалгебраическихуравнений3порядка.
Решениесистемыуравнений,отображенноенацифровыхрегистраторах(Display1–Display3)равно:
x1=0,5x2=−1.x3=2
Рис.2.4.Переходныйпроцессустановлениярешениясистемыуравнений.
3.Аппроксимациясигналов
Вэтойразделерассмотримширокоизвестныйаппроксимационныйметоднаименьшихквадратовдляфункцииоднойпеременной,атакжеальтернативныйметод,которыйпозволяетснизитьвычислительныезатратыиупроститьреализациюнаПЭВМ–модифицированныйметодравныхплощадей.Вчастности,вSimulink-структураханализадинамическихсистемудобнееиспользоватьименнометодравныхплощадейблагодаряегопреимуществам.
Пустьнепрерывныйсигнал
x(t)
заданнанекотороминтервале
измененияаргумента
v
[0,T].Выбираетсянекотораясистемафункций
S={s(t)}n
,заданныхнатомжеинтервалеизмененияаргумента.Из
k k=0
теорииаппроксимацииизвестно[6,18],чтосигналможетбытьпредставленнаопределенноминтервалеизмененияаргументаобобщеннымполиномом
n
xa(t)=∑X(k)sk(t).
k=0
Систему коэффициентов
(3.1)
X(k)обычно называют
аппроксимирующимполиномиальнымспектром(АПС)сигналаx(t).
