Черных И.В. Simulink Среда для создания инженерных приложений (2003) (1152002), страница 43
Текст из файла (страница 43)
Поэтому пользователь, при написании Я-функцни, должен предусмотреть проверку правильности обращения к процедурам ответного вызова„моделирующих дискретную часть системы. а-функция ([ябзпс: хцпсг1оп [зуз,хО,вег, Св) изхе(эи(е,х,ц, 61ад) Ъ Н?ХЕЕИ Пример Б-Функции. Ъ Б-Функция моделирует систему, которая состоит Ъ из последовательности интегратора (1!в) Ъ и блока задержки (1/х). 6. Я-функции и нх применение епб % % Окончание в1хебв яы1ссп Е1ад % Инициализация % )(шп()1всБгагев = 1( я%зев % Вычисление производных % в%зев в1зез в%зев в1хез НивбиЕ.риси = 1) Ншп?присв = 1) 01греейсйгоидИ = 0; Е)шпзавр1еТ1вея = 2: саве 1 вуз=тпй10ег1нас1нев(с,х,и) % в1вя1зез(в%зев) пегов(2,1); зуз хО саве 2, (1: [О 0; брегхой боЕЕяег)) саве 3 % — - — — — -% % Завершение моделирования % саве 9 зуз = []; % З9З З92 И. В. Черных.
Б! М(Л2(ЧК: среда соэдпних инженерных приложений % Шаблон для создания Я-функции — файл яЕипсвр1.в. % % СоРугЕБЬГ 1990-2001 Т)пе МаГЫМогХя, Тпс. % Яйен1в1оп: 1.27 $ % Авторский перевод комментариез: Черных И. В. % % Такт дискретности и задержка для дискретной части систеыы: йрегфой = 1п йоЕЕвес = 0: % % В зависимости от значения переменной Е1ад % происходит обращение к тои или иной % процедуре ответного вызова: саве 0 !вуз,хО,вгг.гв]=вб11п1%1а11хеБ1зея(йрег1ой,йоЕЕвег)„. % % Вычисление вектора состояний дискретной части системы вуя=вй1црйасе[с,х,и,брег%ой,йоЕЕвес); % % Вычисление вектора выходных сигналов вуз=вй10иприсв(с,х,и,йоЕЕвег,йрегхой); % % Неизвестное значение переменной Е1ад % оспегы1ве еггог(('иппапй1ей Е1ад = ',пив2всг(Е1ад)])п % вй11пхеЕа]1зеЯ1зез % % Функция инициализации % % Расчет начальных условий, тактов дискретности, % % размеров векторов сигналов, массива тактов и задержек % % — — -- — — — — — — — -- — — — — — — — — — — ===================% ЕипсСЕоп (вув,хб,вгг,св]=вй11п1г1а11зея1зея(брег%ой,йоЕЕвег) % в1зея = я1вв1хев; я1хез.йивСолгзгапев 1; % Число непрерывных пеРеменных % состояния.
% Число дискретных переменных % состояния. % Число выкодных переменных % Число входных переменных % Прямая передача сигнала. % Размер вектора % тактов дискретности. % Задание вектора начальных % значений переменных состояния. % Параметр зарезервирован. % Такт дискретности для % непрерывной части системы. % Такт дискретности для % дискретной части системы.
% Конец процедуры вй11п1%1а11зеБ1зея % % пк510ег1нап1нез % Вычисление производных вектора состояния % непрерывной части систеиы %— Р1 В. Черных. 3[М[Л.! К[К: среда создания инженерных приложений б. д-фу))каин и их применение гцпссхоп зуз=шд10ег1чаС1чев(г,х,ц) В зув цс В % Конец процедуры шд10егунагхнев В В— В шд1[)рдаге В В Вычисление вектора состояния В $ дискретной части системы В В яя чч в я я я я я яяя я яВ В ГцПСС1ОП Зуаянд1[[рдасв(С,Х,П)драхзсд,дОГГВЕС) В В Вычисление переменных состояния дискретной части В выполняется в моменты времени, определяемые В с погрешностью, которая иа превышает 1е-15.
$ Если условие, заданное оператором 1г истинно, тс выхслной В сигнал блока изменяется. В противном случае выходной $ сигнал остается равным значению на прелыдущеы шаге. В 11 апв[гоцпд[[С - доставя)/дрег1од) . (à — дотбвег)/дрегьод) с 1е-15 вуз х [1); е1ве зуз = [)[ епд $ Конец процедуры )ад10рдасе В В— В вк)10церцев В В Вычисление вектора выходных сигналов $ В" — — — — —- я===яяяяяяяя = === яя====а=яя я я ЫВ В гцпсе1оп зуз=шд10цгрцгз(е,хсц.доггзее,дрег1од) В $ Вычисление выходного сигнала системы выполняется в моменты В времени, определяемые с погрешностью, которая не превышает $ 1е-15.
В Если условие, заданное оператором 1г истинно, то выходной В сигнал блока изменяется. В противном случае выходной В сигнал остается равным значению на предыдущем шаге. В 1[ а)хс [гоцлп[ [с. дос)яес) /дрег1од) - [с до[свес ) /ориг[оси с 1е 15 аув = х(2) ) е1зе вув = []) епд В $ Конец процедуры шд10церисз Ирннср[ Схема модели с а-функцией пнхед[н приведена нв рис. 6.Ь. сяясяя Б е мьел ясяя Рис б.б Медаль даснретаей системы е аеремеиаь[м тамтам нааатевааая Модель дискретной системы с переменным тактом квантования представлена в файле тяйпсзп.
Рассматриваемая Я-функция моделирует устройство задержки сигнала. Время задержки определяется величиной входного сигнала. Для вычисления момента времени срабатывания блока используется процедура ответного вызова пийбсПппеОИЧех[УагНВ, вьпываемая прн значении переменой баб = 4. Поскольку время срабатывания блока зависит от входного сигнала„то параметр О[гресхйЬгопйй должен иметь значение равное ). Ниже приведен текст этой Я-функции. а-функция Йьйвнг б. ункцнн н нк ~ мененне % Завершение моделирования % %-- — — — — -" е=========== ======= =====% саве 9. вув=зк]1тегв1пасе(с,х,ц)г вув = [)~ %= -е= саве 2, % 8 ецпсс1оп (вуз,хО, эсг В эфгев вфввфгев; % в[ген.ццвСопедсагев ВУВеВО10рдаСЕ(В,Х,ц]; В Св] =в611п1%1а11 геЯ[гез саве Э, вуз=вс)10цгрцсв(е,х,ц); В в[зев.Нцв01всзсагез вфгев.ышкбцгрцгэ = 1) в%зев.пцвХпрцгв 2( в1гев.ОЕгуеес)гцгоцдп )( 396 И.
В. Черных. Я(М(](.НЧК: среде создания ннженернык приложений ецпссЕоп (вуэ,хО.эсг,сз] = чвЕипс(с,х,ц,е1ад) В УЯРОМС Пример Я-функции. В Я-функция моделирует устройство залержки сигнала % с переменным тактом квантования % Время задержки определяется величиной сигнала, % поступающего на второй вход". % % бс = и(2) В у(С+ОС) = ц(С] В % Сопутствуюе(ие функции: ЯРОНТМРЬ, СЯРОВ(С. ПЯРОНС.
В В Соруг1д)ьг 1990-2001 Тпе Иагпыог)св. Хпс. В Яйеч1эфоп: 1.9 $ % Авторский перевод комментариев; Черных И. В. В зм1ссЬ Е1ад, % В эависиыости от значения переменной Е1ад % происходит обращение к той или иной % процедуре ответного вызова: % Инициализация % В-- — - — - — — -— — — — -- — В саве О, [вув,хО,всг,св]=шб11пфс1а11геЯЕгеэ % % Вычисление вектора состояний дискретной части системы % % Вычисление вектора выходных сигналов В Вычисление времени срабатывания для следующей % расчетной точки дискретной части системы саве Я, эуэ=лх)1ОеСТЕвеОЕНехеуагН1С(г,х,ц)р % Неиспользуеыые значения переменной Е1ад % % — — — — — - —— =====я==я= = == ====я====ее=я== =я=я=ее==В % % В примере нет непрерывных переменных состояния, поэтому % значение переменной Е1ад = 1 не используется. Результатом % Я-функции в этом случае является пустой массив, саве 1, % Неизвестное значение переменной Е1ад В В- — -еВ огьегм1ве егтог(('Опцапс)1ес) Е1ад ',пцв2вгт(ЕХад)]) % епс] % % Конец функции чвЕцпс В %=-- — —— =я====я== = ============е=яе=е= = я==я===я==янезе == = ==В В вс]11п1с1а11гед1гев В % Функция инициализации % % Расчет начальных условий, тактов дискретнфсти, % % размеров векторов сигналов, массива тактов и задержек % -=В % Число непрерывных переменных % состояния.
% Число дискретных переменных % состояния. % Число выкодных переменных. % Число входных переменных. % Прямая перелача сигналад. % Используется в процедуре В вй16егТЕвеОЕНехгуатНЕВ. 6. 3-функции н нх применение вув = вкввкяев(вкхев) хо = [О] вег []; Св = [-2 О] Ъ в61текв1па Ъ Завершение моделирования Ъ Ъ— Рпс. б.7 39[) ;:е[[(я()е(ь)(с;: дп создания нпжене ных приложений кпсев = 1< Ъ размер вектора тактов Ъ дискретности. Ъ Задание вектора начальных Ъ значений.
Ъ Параметр зарезервирован. Ъ Переменный такт дискретности. Ъ При вызове Я-функции со значением Е1ад = 4 выполняется Ъ расчет слецукллей точки срабатывания. Ъ Ъ Конец процелуры вс)11п1сйа11яеэ1хев Ъ Ъ як)10рс)асе Ъ Ъ Вычисление вектора состояния дискретной части системы Ъ Ъ 0ипсе1оп вув=вс)10рс)аее(е,х,и) Ъ вув = и(1) ) Ъ Ъ конец процедуры вс)10рс)асе Ъ Ъ вс)10иериев Ъ Вычисление вектора выходных сигналов Ъ сипсе1оп вув=вс)10иерикя<е,х,и) Ъ вув = х(1) с Ъ Ъ Конец процедуры вс)10иериев Ъ Ъ пю)10еСТ1веОйнехЕЧагн1Е Ъ Вмчисление следующей расчетной точки Ъ для дискретной части системы. Ъ Функция используется только в случае моделирования Ъ дискретной части с переменным тактом дискретности.
Ъ В этом случае массив Ев должен быть задан как [-2 О]. Ъ сипсе1оп вув=ипд10еСТ1веОЖехСЧагн1С(е.х,и) Ъ вуэ = С + и(2)с Ъ Ъ конец процедуры пИ10естйвеОЕБехсчагн1с Ъ Еипсе1оп яув=пк)1Тегв1паее(е,х, и Ъ вув = []" Ъ Ъ Конец процедуры вб1Тесвйпасе Пример: Схема модели с Б-функцией чяйдзс показана на рис. 6.7. В примере время задержки сигнала увеличивается с 0.0005 с до 0.00[5 с.
б. 8-функции н нх црнмеаенне 1 — — -С Ф ю-г.11, ~я нт 1 (6.5) е) ~' = — (СН . Ф-ю'- М~-], Н 7 (6.6) М =См .Ф.1, (6.7) 4 и=е+к 1+Ь вЂ”, И г (6.1) (6.8) 4 о> М-М,=2 — ° г (6.2) +ие Мг +зй- (6.3) М =Си Ф.(, с=С„Ф.го, Выходные переменные: (6.4) + уе го з у~ . Переменные состояния: Ихо Ю См.Ф 1 ха ! + и О + О з(г (. А (. (6.9) 4Х! См Ф 1 — хо+О хз+О.ио — —.и!. 4г .7 (6.10) уо = Си 'Ф.хо+О х,, у, =0 х,+х,. (6.11) (6.12) Ых — =А.х+В и, сй (6.13) у =Сх. (6.14) 401 И.
В. Черных. яМШ НЧК: среда созданию ннжеиерных приложений 6.5. Ненрерывнав модель электродвнгвтелв ПОСТОЯННОГО ТОКО В этом разделе на примере двигателя постоянного тока независимого возбуждения (ДПТ НВ) рассмотрим создание 8-функции для модели реалы ного объекта. Математическое описание ДПТ НВ Двигатель постоянного тока независимого возбуждения (рнс. 6.8), описывается следующей системой дифференциальных и алгебраических уравнений в абсолютных единицах: где и — напряжение на якорной обмотке двигателя, е — эпектродвижущая сила (ЭДС) якоря, 1 — ток якоря, Ф вЂ” поток, создаваемый обмоткой возбуждения, М вЂ” электромагнитный момент двигателя, Мс — момент сопротивления движению, оз — скорость вращения вала двигателя, й — активное сопротивление якорной цепи, з.