Получение импульсной передаточной функции с помощью Mathcad Текст научной статьи по специальности «Математика»

УДК 004.021

ПОЛУЧЕНИЕ ИМПУЛЬСНОЙ ПЕРЕДАТОЧНОЙ ФУНКЦИИ С ПОМОЩЬЮ MATHCAD

В.А. Фаерман, Е.М. Яковлева

Томский политехнический университет E-mail: TheUnit@sibmail.com

Предложено применение математического пакета Mathcad 15 к получению передаточных функций (в операторной z-форме) типовых одноконтурных импульсных систем, имеющих произвольный типовой импульсный элемент на выходе сравнивающее-суммирующего устройства, с известной передаточной функцией непрерывной части.

Ключевые слова:

Импульсные системы, z-преобразование, импульсные элементы, Mathcad, импульсные передаточные функции, формирующие элементы.

Key words: Pulse systems, z-transfer, pulse elements, Mathcad, pulse transfer functions, forming units.

Введение

Возможности современной вычислительной техники позволяют выполнять численными методами детальный анализ систем автоматического управления, сопоставимый с результатами физического макетирования. Для анализа или синтеза, необходима математическая модель адекватно описывающая систему. Современное программное обеспечение позволяет существенно упростить процесс моделирования различных систем автоматического управления. Линейной импульсной системой называется такая система автоматического управления, которая кроме звеньев, описываемых линейными дифференциальными уравнениями, содержит импульсный элемент, преобразующий непрерывное входное воздействие в последовательность импульсов. Одноконтурную импульсную систему автоматического управления можно представить как взаимодействующие друг с другом импульсную и непрерывную (НЧ) части систем автоматического управления (САУ). В непрерывную часть обычно входит объект управления, а также усилительное и исполнительное устройства. Импульсная часть (ИЧ), как правило, является управляющим устройством и объединяет функциональные элементы, участвующие в импульсном преобразовании сигнала. В процессе преобразования непрерывного сигнала в дискретный импульсный элемент выполняет две операции - квантование по времени и импульсную модуляцию. Первая заключается в том, что сигнал на выходе импульсного элемента появляется в определенные дискретные, обычно равноотстоящие моменты времени t I = Т ■ /. В результате импульсной модуляции изменяется какой-либо параметр выходного импульса (амплитуда, ширина и т. д.) в зависимости от входного сигнала. В зависимости от изменяемого параметра импульса различают разнообразные виды модуляции: широтноимпульсная модуляция (ШИМ), амплитудно-импульсная модуляция (АИМ) и другие. При амплитудно-импульсной модуляции изменяемым параметром служит амплитуда импульсов, которая обычно пропорциональна значению непрерывного сигнала g(Ti) в моменты времени Т\. При широтно-импульсной модуляции модулируемым параметром является ширина

Фаерман Владимир Андреевич, студент Института кибернетики ТПУ.

E-mail: TheUnit@sibmail.com Область научных интересов: математическое моделирование систем автоматического

управления, автоматизированное проектирование

электронных устройств. Яковлева Елена Максимов-на, канд. техн. наук, доцент кафедры автоматики и компьютерных систем Института кибернетики ТПУ.

E-mail: yakovlevaem@tpu.ru Область научных интересов: автоматизация проектиро-вания печатных плат электронных устройств в САПР DipTrace и Catia с 3-х мерным изображением.

(длительность) импульсов, амплитуда при этом остается постоянной. Системы, в которых присутствует ШИМ относятся к классу нелинейных импульсных систем [1].

Основным отличием импульсных систем от непрерывных является прохождение сигнала только в определенные, дискретные моменты времени, поэтому непрерывная часть реагирует лишь на дискретные значения непрерывного сигнала в моменты квантования пТ. Поэтому непрерывную функцию х(У), описывающую непрерывный сигнал, можно заменить соответствующей решетчатой функцией х(пТ), где п - целое. Если интервал квантования Т задан, то по функции х(0 решетчатая функция х(пТ) определяется однозначно. Соотношение между решетчатой функцией х(п) и ее разностью Д кх(п) определяет уравнение в конечных разностях или разностное уравнение. Если это соотношение линейно, то разностное уравнение называется линейным. Оно может быть представлено как

акДк(п) + ак_ 1Дк_ 1 (п) + —Ъ а0х(п) = /(п) (1)

Ъкх(п + к) + Ък_ 1х(п + к — 1) + —Ъ Ъ0х(п) = /(п), (2)

где Дп) - известная решетчатая функция, х(п) - искомая решетчатая функция, представляющая собой решение разностного уравнения. Данное разностное уравнение, содержащее х(п) и х(п+к), называется разностным уравнением £-ого порядка. В том случае, когда _Дп)=0, уравнения (1) и (2) называются однородными [2].

Классические методы решения разностных уравнений во многом аналогичны классическим методам решения дифференциальных уравнений. В частности, общее решение неоднородного линейного разностного уравнения может быть представлено в виде суммы общего решения однородного разностного уравнения (свободной составляющей) и частного решения неоднородного разностного уравнения (вынужденной составляющей).

х (п) = хсв (п)+Хп^р 00 . (3)

Свободная составляющая решения (3) определяется из характеристического уравнения и имеет вид

к

Хсв № = ^С1 21П ,

¿ = 1

где - произвольные постоянные, определяемые начальными условиями. Вынужденная составляющая определяется видом правой части уравнения и может быть найдена, для многих типовых функций, методом подбора. Решением разностного уравнения х(п) является решетчатая функция существующая в моменты времени Тп, что соответствует наличию ключа на выходе системы. В реальных системах ключ отсутствует, поэтому выходной сигнал является непрерывной функцией, которая может быть восстановлена из х(п). Несмотря на то, что решение уравнений состояния (как аналитическое, так и численное) не представляет сложности, широкого применения на практике оно не находит. Это объясняется трудностями, связанными с получением разностных уравнений из дифференциальных уравнений непрерывной части. В ряде случаев разностные уравнения непрерывных звеньев получают из их импульсных передаточных функций [1, 3].

Для исследования линейных импульсных систем, на практике обычно применяются передаточные функции, являющиеся отношением операторных изображений выходного и входного сигналов системы. Так как в импульсных системах сигналы существуют в моменты времени Т, то есть являются дискретными, то передаточные функции называются импульсными и определяются как [2]:

X в ых (Ю

вх (Ю

Для решетчатых функций может быть введено дискретное преобразование Лапласа, определяемое формулой [1]:

е (р) = 1ит^_11" . (4)

В случае, если (4) сходится при Яе(р) < со то оригиналу ^}) соответствует некоторое изображение F(p). Под 2-преобразованием понимают изображение функции, определяемое следующим образом

е 00=1 Т= о! (I )■? 1. (5)

Изображение существует и соответствует оригиналу, если ряд (5) сходится. Дискретные преобразования связаны соотношением е_^т = г. Формула 2-преобразования может быть записана также для непрерывной функции в виде

е (г)=Х 0= о! ( ^■г _1, 1^ = 1 Т. (6)

Согласно (5) и (6) оригиналом W(z) должна служить импульсная переходная функция м>(?)

№ (г) =X0= о™ (IТ ■ г _1. (7)

1. Формирующие элементы в импульсных системах

В импульсных системах, используются различные импульсные элементы, которые принято описывать как последовательное соединение идеального импульсного элемента (ключа) и формирующего элемента. Идеальный импульсный элемент под воздействием непрерывного входного сигнала g(t) формирует идеальные мгновенные импульсы х(Тп) вида -функций (рис. 1). Формирующий элемент преобразует эти импульсы в Рис. 1 Идеальный ключевой элемент сигналы х(0 нужной формы.

При описании импульсных систем используются различные формирующие элементы, имеющие отличные математические описания и характеристики, влияющие на импульсную передаточную функцию системы. Наиболее распространены экстраполяторы нулевого и первого порядков, а также интерполяторы первого порядка.

Экстраполятор нулевого порядка удерживает на выходе величину, равную в течение всего периода Т (рис. 2). Взяв реакцию элемента на отдельное замыкание ключа (рис. 3), получим прямоугольный импульс. Математически он может быть описан следующей функцией:

*п (О

_ [дп> ^ є (Тп, Тп+1) _ 0, 0ґИ.ЄГ]Л/І5Є

9п(.^-(.^п) 1(^71 Т')')

(8)

Рис. 2. Экстраполятор 0-ого порядка

Рис. 3. Реакция на один импульс

Используя выражения (8), (4-5) составим передаточную функцию формирующего элемента

щэ О)

1-е

-бТ

(9)

Совершив замену в (9) е^ = г получим

Щэ ( ^ ,г)

1-г-1 г-1 1

Таким образом, экстраполятор нулевого порядка можно представить как два блока, один из которых можно включить в непрерывную часть. Экстраполятор первого порядка может быть представлен через преобразователь нулевого порядка следующим образом

1

91 Щ —► і - 1 7 ►- Тя2 г

1 экстра погиятор х(1п)

г Б 0-ого порядка г

Рис. 4. Экстраполятор первого порядка

Аналогично экстраполятор второго порядка может быть представлен через экстраполятор первого порядка и так далее. Таким образом, в виде соединения непрерывного и импульсного блока, может быть представлен любой типовой формирующий элемент (интерполятор или экстраполятор любого порядка) [4].

/э (5, 2) = (5 ) ■ (Ю (10)

Представление (10) импульсного элемента упрощает получение импульсной передаточной функции системы с помощью МаШса^

2. Получение передаточной функции линейных импульсных систем в Mathcad

Для удобства нахождения импульсных передаточных функций замкнутой системы по задающему воздействию, по ошибке и др., удобно начать с передаточной функции разомкнутой системы. Пусть имеем одноконтурную систему (рис. 5), содержащую импульсный элемент, расположенный после сравнивающе-суммирующего устройства; а также непрерывную часть, включающую усиливающие устройства, объект управления в прямой цепи; корректирующие устройства, расположенные в цепи обратной связи.

Рис. 5. Одноконтурная импульсная система

Выходной сигнал непрерывной системы у(0, при нулевых начальных условиях, определяется как свертка входного х(0 сигнала (задающего воздействия) и импульсной передаточной функции (весовой функции) м>(?)

у ( 0 = /о( 1 -т) х (т) йт, (11)

где весовая функция м>(?) - реакция системы на функцию Дирака. Весовая функция системы м>(?) является оригиналом передаточной функции непрерывной системы Щя), то есть связана с ней обратным преобразованием Лапласа. Следовательно, для того, чтобы перейти от непрерывного преобразования Лапласа к дискретному преобразованию достаточно проделать следующие действия [3]:

1. Определить передаточную функцию разомкнутой системы, с учетом формирующего элемента (14):

Wp (s) = ИфГ (S) ■ WT (S) ■ wr (s) . (12)

2. С помощью обратного преобразования Лапласа (согласно (11)) можно найти импульсную переходную функцию разомкнутой системы

w ( 0 = L _ 1{WПрс ф } . (13)

3. Согласно (6), определить весовую последовательность системы (решетчатую функцию веса):

w(nT) = wft^t = Tn. (14)

4. Согласно (7), применить прямое дискретное преобразование Лапласа

W 0 рс (z) = Z{w (пГ) } . (15)

5. Согласно схеме (рис. 8) домножить выражение (19) на WJM ( z)

W- (z) = W 0 рс (z) ■ W|T (z) .

Получение передаточной функции разомкнутой линейной импульсной системы в ППП Mathcad, согласно действиям, перечисленным выше, выполняется следующим образом: первые четыре пункта реализуются участком программы, представленным на рис. 6, пятый пункт реализуется участком, представленным на рис. 7.

Рис. 6. s-z преобразование

Wpc(z) := W(z) - W*np(z) simplify —>

Рис. 7. Получение передаточной функции

где invlaplace (W(s)) - выполняется символьно и реализует (13);

substitute, t = Tn - выполняется символьно, реализует (14);

ztrans (w(Tn)) - выполняется символьно, реализует (15);

simplify (W(z)) - выполняется символьно, упрощает выражение-аргумент (в данном случае сокращает подобные множители).

Функции, указанные в рис. 6 и 7, могут быть найдены на символьной панели инструментов (вызывается через меню окна приложения следующим образом Вид - Панели Инструментов - Символьные) или в меню вставки функции (вызывается иконкой fx) в строке инструментов окна приложения) в разделе - Только аналитические (symbolic only).

Форма, в которой будет получено решение, неудобна для работы и нуждается в дальнейших преобразованиях (ввиду численного решения возникают коэффициенты, порядок которых может достигать нескольких десятков). Для того чтобы представить полученную функцию в удобном для дальнейшего исследования системы виде:

W* (z)= . (16)

Либо:

Ykh-7i

W- (z) =

(17)

достаточно найти коэффициенты а¿, / = 0, ж; Ъ¿, / = 0 , 0; и корни характеристического уравнения 2 ¿, / = 0 , ж Для того, чтобы представить передаточную функцию в виде (16) можно воспользоваться участком программы, представленным на рис. 8; а для (17) необходимо вопользоваться участком программы, представленным на рис. 9.

В := пишет (W(z)) coeffs А := denom(W(z)) coeffs ->

В := Round

А := Round

В

Arovv-s(A}—1 А

. 10

-7

¡-i

W(z) :=

Aron's (A) rows(B)—1 ,

V M

i = 0 rows (A)—1 , .

Z M

i = 0

Рис. 8. Получение передаточной функции в форме (16)

В := mmierCW(z)) coeffs -> А := denom(W(z)) coeffs В

Afflws(A) — 1

( А

А := Round

С := polyroots (В)

W(z) :=

Afows(A)-1 rows (В)— 1 .

V M

i= о

rows (С)— 1

П (z-q>

Рис. 9. Получение передаточной функции в форме(17)

где denom(W(z)) - выполняется символьно, возвращает знаменатель дробно-рационального выражения; numer(W(z)) - выполняется символьно, возвращает числитель дробно-

рационального выражения; coeffs (^(z)) - выполняется символьно, возвращает коэффициенты полинома; polyroots (A) - возвращает вектор корней полинома по его коэффициентам; round (x, е) - округляет x с точностью до Е. Формализация округления коэффициентов упрощает дальнейшее получение импульсной передаточной функции замкнутой системы.

Если необходимо сделать процесс преобразования более наглядным, можно предварительно использовать функцию parfrac (W(z)) (выполняется символьно, представляет дробно-рациональное выражение в виде суммы элементарных дробей). В дальнейшем предложенный алгоритм преобразования (кроме шага 5) может последовательно применяться к каждой из элементарных дробей.

Для того чтобы получить передаточную функцию замкнутой системы, можно воспользоваться выражением [3]

Шзам(z) = 1+W-pс (Z) , (18)

где Wp (z) _ импульсная передаточная функция прямой цепи, которая определяется по предложенному выше алгоритму, с учётом выражения:

Wp (s) = ИфТ (s) ■ WWjT (s) .

Для упрощения преобразования (18) как вручную, так и с помощью Mathcad, удобно представить Wnp (z) и 1 + W с (z) в виде (17), используя участок программы, представленный на рис. 10. Аналогичным образом может быть найдена импульсная передаточная функции

системы по ошибке регулирования по задающему воздействию. При этом вместо выражения (18) следует использовать выражение

We(z) =

1+WpC(z)' Wnp(z)

Wsaarfz) := ------------ simplify

1 + Wpc(z) '

Взам := numer(W3aM(z}) coeffs

Азам := denom(W3aM(z)) coeffs

Азам

Азам := Round

Взам := Round

А3 ™го,Лг с(д.іш) _ 1 Взам ^ As^rowsi Азам)-1

,Ю

- 131

,10

- 71

WsaM(z) :=

rowsiBsai)-1

X

i= о

Взам.

i"0

пш(Аза[)-1

s

i= 0

Азам

Рис. 10. Получение передаточной функции замкнутой системы

В качества примера найдем передаточную функцию разомкнутой импульсной системы автоматического регулирования (САР) частоты вращения двигателя постоянного тока, структурная схема которой приведена на рис. 11. Численные значения параметров системы приведены в таблице.

Рис. 11. Структурная схема САУ двигателем постоянного тока

Таблица. З начения па раметров САР

Параметр Кі К2 Кз К4 Ті Т2 Тз Т4 T

Значение 18,75 20 1,6 0,02 0,06 0,02 0,3 0,01 0,157

Согласно (12) найдем выражение для приведенной непрерывной части разомкнутой системы

Ш (з) = ___________КгЪЯгЪ______________

РСУ ' Б■ ОУГз-*2'+Г3-5+1) Т-б+1) Т -я+1) ■

Выполнив операции, приведенные на рис. 6, 7, получим следующий результат

Рис. 12. Получение передаточной функции

Проделав операции, приведенные на рис. 8 получим результат в форме (17) О (УЮ ■ 7 - ?. 0Я7 ■ 72 - 7 ■ ? - 7. 19Й х 1(Г 6

ТУрс(г) ----------

-2.848 х 10

-:

г - 0.042 ■ 7" - -0.643 - г + г + 4.325 х 10

12

Рис. 13. Передаточная функция в форме (17)

Данный результат был использован для получения импульсных передаточных функций замкнутой системы и по ошибке регулирования. Достоверность использования данного алгоритма получения импульсных передаточных функций может быть подтверждена моделированием в ППП Matlab.

Заключение

Сейчас с применением современного программного обеспечения возможно решение широкого спектра задач моделирования и анализа импульсных систем. В настоящей работе было продемонстрировано применение математического пакета Mathcad 15 к получению передаточной функции одноконтурной линейной импульсной системы автоматического управления по известной передаточной функции непрерывной части. Предложенный алгоритм может быть использован для упрощения решения задач синтеза линейных импульсных систем с заданными требованиями, а также для упрощения процесса получения разностного уравнения непрерывной части импульсной САУ. Последнее актуально при моделировании систем с широтно-импульсной модуляцией, требующем численного решения системы разностных уравнений.

СПИСОК ЛИТЕРАТУРЫ

1. Бесекерский В.А., Попов Е.П. Теория систем автоматического управления. - СПб.: Профессия, 2004. - 767 с.

2. Яковлева Е. М., Замятин С. В. Курсовое проектирование по теории автоматического управления. - Томск: Изд-во ТПУ, 2010. - 106 с.

3. Теория автоматического управления / под ред. А. А. Воронова. - М.: Высшая школа, 1986. - 367 с.

4. Чечурин Л.С. Курс лекций по ТАУ: дискретные динамические объекты. URL: http://chechurin.com/book_tau/contents/p1_lecture2.html/ (дата обращения 26.06.2011).

Поступила 30.06.2011 г.