Черных И.В. Simulink Среда для создания инженерных приложений (2003) (1152002), страница 42
Текст из файла (страница 42)
Результатом работы Б-функции в этих случаях % является пустой массив. % саве ( 1. 2, 4, Я ) зув=[); % Неизвестное значение переменной Е1ад % % аснегызве еггог(['Опьапй1ей Е1ад = ',пив2зсг(Е1ад))); % епй % Конец функции Сыпевсмо % 6. 5-фунхцнн н нх примсненнс =% % в1пев = випвххевз в1хев.ицвзсопсвсасев = О: % Число непрерывных переменных % состояния. в»*в» ьз»»»»» в1яев.плдп01всвсасев = О; в1хев.ицвбцсрцев = -1з в1яев.ншптпрцев — -)з в1яев.01гуеес(снгоцц)з = 1з в1яев.ишпвавр1еТипев = 1) вув = в1вв1яев(в1пев) вег [1) хО = [) Рис. бЗ ев = [-1 О) вххев = вхввххевз И. В.
Черных. ЯМЫ.ВЧК: среда создания нюкенпрных прняовеннй % вс)11п1С1а11яе51хев % зпс)11п161а11яе51гев % Функция инициализации % Расчет начальных условий, тактов дискретности, % размеров векторов сигналов, массива тактов и аадерзпек % Еццсс1оп [вув,хо,вег,ев) = вб11пхейа1хаевхяев() % % Число дискретных переменных % состояния. % Число выходных переменных. % Задается динамически.
% Число входных переменных % Задается динамически. % Прямая передача входного % сигнала на выход. % Размер вектора тактов % квантования. % Параметр зарезервирован для % будущего использования. % Задание вектора начальных % значений переменных состояния. % Такт квантования и вадерзпка. % Такт наследуется от % прелществующего блока. % Конец процедуры вс)11п1%1а11пе51пев % % пх[10церцсв % Функция вычисления вектора выходных сигналов % Ецпсс1оп вув = вз)10церцев(с,х,ц) % вув -" ц " 2з % Выходной сигнал блока есть вхолной сигнал, % умноженный на коэффициент 2. % % Конец процедуры вб10церцев Пример: Схема модели с Б-функцией йтез(во приведена на рис. 6.3.
На осцнллозрамме видно, что амплитуда выходного сигнала в два раза больше входного. Модель непрерывной системы Модель непрерывной системы, описываемой уравнениями в пространстве состояний, представлена файлом свйзпс.пх Соответствующая Б-функцня моделирует непрерывную систему с двумя вхпдами, двумя выходамн и двумя персменнымн состояния Параметры модели (значения матриц А, В, С, [1) задаются в теле Б-функцин и передаются в процедуры ответного вызова через их заголовки в вцле дополнительных аргументов. Ниже приведен текст этой Б-функции.
Ецпссхоп [вув,хо,вег,ев) = свЕцпс(г,х,ц,Е1ад) % 05РОНС Пример 5-функции. 0 помощью уравнений пространства Ъ состояния моделируется непрерывная систеыаз И. В. Черных. Б!МЛ ВЧК: среда создания инженерных врняоженнй Ъ Ъ х' = Ах + Вц У = Сх + Пц Ъ Ъ Ъ Ъ шаблон для создания Б-функции - файл вЕцпсир1.и. Ъ Значения матриц передаются в са11Ьаск-методы через их заголовки в качестве дополнительных параметров СоругхдЬС 1990-2001 'ХЪзе МаСЬИогйв, 1пс. фйеч1зфоп: 1.8 $ Авторский перевал комментариев: Черных И.В. Ъ Ъ Ъ Ъ Ъ Задание Ъ .А= [-О. 09 1 матриц: -0.01 0)з Ъ Матрица системы.
Ъ В=[ 1 0 Ъ С=[ 0 1 -7 -2] Ъ Матрица входа. 2 — 5]т Ъ Матрица выхода. Ъ П=[-3 1 0 О] Ъ Матрица обхода Ъ Инициализация Ъ Ъ вЂ” Ъ саве О, [зув, хО, всг, св] =ид11п1с[а11хеБ1хев (А, В, с, и): Ъ Вычисление производных Ъ саве 1, зув=и61Пе гфча С яме в ( С, х, ц, А. В, С, П); Ъ внфгсЬ Е1ад, Ъ В зависимости от значения переменной Е1ад Ъ происходит вызов того или иного метода: 6 Б ункцнн н нхпрнмененне Ъ==-- Ъ Вычисление вектора выходных сигналов Ъ Ъ саве 3 вуз=хх)10ЦСПНСз (С. х, ц. А, В, С, П); Ъ Неиспользуемые значения переменной Е1ад Ъ вЂ” — — -"--" — — — --= — - — -" — — -"-=== — ===В Ъ Ъ Неизвестное значение переменной Е1ад Ъ Ъя яЪ о сЬеги1зе еггог(['ппьапд1еб е1ад = ',пцзя2всг(Е1ад)])з епб Ъ )(анен свЕцпс В иб11пфсфа11хедххев Ъ Функция инициализации Ъ Расчет начальных условий, тактов дискретности, Ъ размеров векторов сигналов, массива тактов и задержек Ъ ецпсс1оп [вув, хО, зсг, св) =иб11пфсфа11хейххев (А, В, с, и) Ъ вфхев = в1из1хез; в1хев.нашсапсБсасев 2з Ъ Число непрерывных переменных Ъ состояния.
0; Ъ Число дискретных переменных Ъ состояния. вЕхев.НииПЕвсБСагев Ъ Ъ В примере не используются методы для завершения работы Ъ Б-функции и нет дискретных переменных состояния, позтому Ъ значения переменной Е1ад 2, 4, 9 не используются. Ъ Результатом Б-функции в зтом случае является пустая Ъ матрица. Ъ саве [ 2. А, 9 в) эд4 ЭВБ 6. 3-функции н нх применение в1аев в1зев выев выев НитпоиСриев Нищ1приев Р1греебг)«гоид)« Ииа«Яаа«р1еттщев в1щвыев(вязев); зеков(2,1): [)' аув ко [О О); вув = А*к + В*и« $ Рис б.я выходных сигналов Л,В,С,Р) РЕ В. Чарных ЯМ()).ПЧК; среда создания инженерных прил«аканнй 2; $ Число выходных переменных. 2< $ Число входным переменных.
$ Прямая передача. 1« $ Размер вектора тактов В дискретности. В Задание вектора начальных условий $ Параметр зарезервирован дпя $ будущега использования. $ Массив,определяющий такты $ дискретности и задержки. В Конец процедуры щб11п1ста11зевязев В В=-= В вю)1Регтнагянев В Вы числение значений производных вектора $ состояния непрерывной части систеьа« $ — — В $ ЕипсС1ап вув=п«З1Регянае1нев(е,х,и,й,В,С,Р) В В Конец процедуры щб1Регзнаехнев В В щб10иериев В Вычисление значений вектора Еипссзоп вув=щб10исрисв(с,х,и В вув = С*к + 0*из В $ Конец процедуры щб10иериев Пример: Схема модели с Б-функцией «зуцпс приведена на рис. 6.4. Медаль диенретиой системы Модель дискретной системы, описываемой уравнениями в пространстве состояний, представлена М-файлом «[зйн)с.ю.
Эта Б-функция, папи«щипая иа языке МКП.АВ, позволяет смоделировать дискретную систему с двумя входами, двумя выходами н двумя переменными состояния, Параметры модели [значения матриц А, В, С и Р) задаются в теле Б-функции и передаются процедурам ответного вызова в качестве входных аргументов через их заголовки. ЕипссЕоп [вув,хо,всг,св) = бвЕипс<с,к,и,Е1ад) В РЗРПНС Пример З-функции. В Модель дискретной систеыы в пространстве состояний: В к(п+1) = Хк(п) + Ви(п) В у(п) = Ск(п) + Ри(п) В В Значения матриц передаются в процедуры ответного вызова $ в вице входных аргументов М.
В. Черных. ЯМЕ)(.(МК: среда создания нвженерных нраложений 6, 5-функцнн н нх применение % % Шаблон для создания Б-функции - файл вЕцььсир1.п. % % СоРуг1дпе 1990-2001 ТЬе Мае)ьМогйв, Епс. % ЬйечЕвфопь 1.8 8 % Авторский перевод комментариевь Черных И. В. % % Задание матриц". А = [0.9135 0.1594 -0.7971 0.5947)) % Матрица системы. % В я (0.05189 0 0.4782 0)ь % С= [01 1 0)) % Р я [0.01 0 О -0.02)) % виЕасЬ Е1ад, % Матрица входа. % Матрица выхода. % Матрица прямой передачи.
% В зависимости от значения переменной Е1ад % происходит обрашение к тай или иной % процедуре ответного вызова: =Я=== =як= я †% —--- % Инициализация % %== яяяян я=я я % савв О, [вув,хО,всг,св) я вхИХп1%1а11зедфзев(А,В,С,Р) ) % Вычисление вектора состояний дискретной части системы % саве 2, вув пьй10Рйаее(С.Х,ц,А,В,С,Р)ь % Вычисление вектора % непрерывной час~и с %--- саве 3, вув пй1бцерцев( % выходььых сигналов истеыы Е,х,ц,А,С,Р)) — — — — — -я я % %- % Неиспользуемые значения переменной Е1ад % % % В примере не используются процедуры завершения, % нет непрерызнмх переменных состояния.
% Значения переменной Е1ад ), 4, 9 не используются. % Результатом 8-функции з этих случаях является пустой массив. % саве ( 1 4 9 ) аул=[1 % Неизвестное значение переменной Е1ад % ос)ьегы1ве еггог( ['цп)ьапй1ей Е1ад = ',ппш2зсг(Е1ад) ) ) ь епй % % Конец процедуры йвЕцпс % я % шй11п)%1а11ге81зев % % Функция инициализации % % Расчет начальнмх условий, тактов дискретности, % % размеров векторов сигналов, массива тактов и задержек Сз) = шй11пЕЕЕа11зе81хев(А,В,С,Р) вЕзе(Р,1) ) вЕзе(Р,2) з1хев.Нцшбцерцев везен.мше1прцсв % Ецпссьоп [вув хО зег % везен = вьшвезев в1зев.НпжСопЕЗСасев вЕзез.МциРЕзсБЕасев = зЕзе(А, 11) % Число непрерывных % переменных. % Число дискретных % переменных.
% Число выходных % переменных. % Числа входных % переменных. 6. Б функции н нх и ененис взкез.Жгрееделгоцпп = 1; в1хез.ицвзмхр1етьиев = 1 1 Ъ Прямая передача сигнала. Ъ Размер вектора Ъ тактов дискретности. Ъ ак)1црбаее Ъ Вычисление вектора состояния Ъ дискретной части систеым Ъ Ъ Конец процедуры пк31црс)аее Ъ Рнс. 6.5 391 И. В.
Черных. Б(М()о(ЧК: с да создания инженерных ~ няожсннй Ъ зув = выхзьзез(в[зев); Ъ хО = зеков(ззгез.нциПхзсзеаеез. 1)з Ъ Бодание вектора начальных условий ло переменным состояния. всг - [): Ъ Параметр зарезервирован. Сз [0.2 0)1 Ъ Массив, задающий, Ъ такт дискретности и задержку. Ъ Ъ Конец процедуры вд11п1с[а1[хезькев Ъ хцпссзоп вув = пк)1прс(асе(с,х.ц.А,В,О,П) вуз = А*к+В*цз Ъ Ъ ак(10церцез Ъ Вычисление вектора выходных сигналов Ъ 1олсезоп зув = ис(10церцев(г,х.ц,А,П,В) ЗУВ СХХ+Р*Ц( Ъ Ъ Конец процедуры мс(10церцсв Пример: Схема модели с 3-функцией ([зйщс приведена иа рнс.
6.5. Модель гибращной системы Модель гибридной системы, содержащей как непрерывные, так н дискретные компоненты, представлена файлом щ[хеб(п.щ. Рассматриваемая $-функция моделирует систему, состоящую из последовательно включенных интегратора ([Ь) и блока задержки () Й). Особенность'Я-функции для гибридной системы атом, что обрмцение к процедурам ответного вызова при моделировании дискретных компонентов выполняется в те же моменты времени, что н для непрерывных компонентов.