Электронный научный журнал
Современные проблемы науки и образования
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,791

АНАЛИЗ ЭФФЕКТИВНОСТИ ПРИМЕНЕНИЯ ДЕЛЬТА-ПРЕОБРАЗОВАНИЯ ПРИ МОДЕЛИРОВАНИИ ЗВЕНЬЕВ ВТОРОГО ПОРЯДКА

Андреев В.С. 1 Бутусов Д.Н. 1 Каримов Т.И. 1 Липкин С.М. 2 Сотнин М.И. 3
1 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Санкт-Петербургский государственный электротехнический университет «ЛЭТИ» им. В.И. Ульянова (Ленина)»
2 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Южно-Российский государственный технический университет (Новочеркасский политехнический институт)»
3 Государственное образовательное учреждение высшего профессионального образования «Новгородский государственный университет им. Ярослава Мудрого»
В статье производится оценка эффективности применения дельта-преобразования при переходе от непрерывного описания систем дифференциальных уравнений к дискретной форме на примере моделирования различных типовых звеньев второго порядка. Дается краткое описание дельта- оператора и описываются способы эквивалентного перехода от передаточных функций в z- и s-форме к передаточным функциям в дельта-форме. Методом совместного интегрирования реализована обобщенная модель динамической системы n-го порядка, описываемая дельта-оператором. Для оценки дельта-систем во временной и частотной области использовалась среда модельного проектирования NI LabVIEW, где был реализован ряд моделей передаточных функций в различных формах представления. Результаты моделирования представлены графически. Показано, что применение дельта- преобразования и совместного интегрирования при вычислениях с фиксированной точкой позволяет повысить степень эквивалентности цифровой и непрерывной систем.
погрешность интегрирования.
z-преобразование
дельта-преобразование
звено второго порядка
моделирование
1. Бутусов Д. Н., Жуков К. Г. Коррекция погрешности решения дифференциальных уравнений второго порядка // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета. – 2010. – № 6 (113). – С. 121-126.
2. Бутусов Д. Н., Жуков К. Г. Реализация цифровых фильтров методом совместного интегрирования // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета. – СПб., 2009. – № 6. – С. 26-34.
3. Жуков К. Г., Бутусов Д. Н. Коррекция погрешности решения уравнения гармонического осциллятора методом последовательного интегрирования // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета. – 2010. – № 6 (113). – С. 135-144.
4. Ту Ю. Современная теория управления. «Машиностроение», 1971.
5. Newman M. J., Holmes D. G. Delta operator digital filter for high performance inverter applications // IEEE Transactions on Power Electronics. – 2003. – Vol.18. – No.1. – Part 2.

Введение

Дельта-преобразование является, наряду с более известным z-преобразованием, одним из методов математического описания дискретных систем. Разработанное практически одновременно с z-преобразованием, оно долгое время оставалось мало востребованным, так как невысокое быстродействие вычислительных машин прошлого, применяемых при реализации систем автоматического управления и при цифровой обработке сигналов, не позволяло проявляться главному недостатку z-преобразования, а именно: при стремлении периода дискретизации к нулю корни и полюса системы в z-области стремятся к единице [4]. При ограниченной точности машинного представления чисел различные корни стремятся «слиться» друг с другом и с единицей; динамическая характеристика такой цифровой системы не повторяет характеристику непрерывной и может существенно отличаться от нее [5]. Этим эффектом можно пренебречь, если использовать представление чисел с плавающей точкой. Однако на аппаратном уровне его поддерживает довольно ограниченное число контроллеров. Использование же чисел с плавающей точкой на программном уровне приводит к существенному увеличению времени выполнения программы. В результате разработка цифровой системы управления с малым периодом дискретизации становится труднореализуемой задачей.

В то же время дельта-преобразование не имеет вышеописанного недостатка. При уменьшении периода дискретизации динамическая характеристика такой цифровой системы стремится к характеристике непрерывной. Еще больше повысить эффективность дельта- преобразования можно путем реализации программы не каноническими формами программирования, а так называемым совместным интегрированием, описанным в работе [2]. При этом становится возможным использование 16-битного и даже 8-битного контроллера там, где в случае z-преобразования приемлемую точность давал бы только 32-битный.

Преобразование непрерывных и z-систем в дельта-системы

Основная идея дельта-преобразования – использование метода Эйлера для вычисления производной:

, (1)

где – период дискретизации.

Как аналог оператора введем оператор , который имеет следующий смысл:

, (2)

Выведем связь с z-преобразованием. Как известно, его суть состоит в переходе от операции дифференцирования к операции сдвига во времени. Переменная при этом обозначает сдвиг во времени на 1 шаг вперед:

(3)

Исходя из (2) и (3), можно записать связь между переменными и :

(4)

Одним из наиболее удобных способов аппроксимации непрерывной системы, сохраняющей ее устойчивость, является билинейное преобразование, или преобразование Тастина. С учетом выражения (3), можно записать формулу перехода по методу Тастина из s-области в дельта-область:

(5)

В данном случае – период дискретизации при дифференцировании методом Эйлера, а – период дискретизации интегрирования, причем они необязательно равны.

С точки зрения инженерной практики, удобно сперва перевести непрерывную систему в z-систему, а затем преобразовать ее коэффициенты таким образом, чтобы заменить оператор на. Допустим, у нас имеется дискретная передаточная функция второго порядка:

(6)

Из формулы (4) известна связь и , откуда следует связь и :

(7)

Используя выражения (4) или (7), мы можем преобразовать передаточную функцию дискретной системы к виду:

(8)

Связь коэффициентов функций в и форме представлена в таблице 1.

Таблица 1. Связь коэффициентов дельта и для передаточных функций 2-го порядка

Функции более высокого порядка, чем 2-й, на практике лучше реализовать как последовательное соединение функций 1-го и 2-го порядка для уменьшения эффектов округления и переполнения разрядной сетки [1].

Совместное интегрирование

В работе [2] показано, что совместное интегрирование имеет более высокую точность при реализации численных методов интегрирования. Рассмотрим этот метод подробно.

Пусть имеется передаточная функция -го порядка, которая была преобразована в систему дифференциальных уравнений 1-го порядка:

(9)

Приведем ее к виду для совместного интегрирования:

(10)

Теперь перейдем от непрерывной системы к дельта-системе. Для этого оператор интегрирования заменим на . Учитывая, что , можно записать:

(11)

Подставляя выражение для из системы (10), получим окончательно:

(12)

Проделав аналогичные выкладки для всех уравнений системы (10), нетрудно записать общий вид матриц пространства состояний совместного интегрирования для дельта-системы:

,,, (13)

Сравнение z-преобразования и дельта-преобразования на простейшем звене

Возьмем в качестве примера простейшее колебательное звено второго порядка:

(14)

Пусть период дискретизации и . Тогда в z-области система будет выглядеть так:

(15)

Отсюда, используя формулы из таблицы 1, получим дельта-форму:

(16)

Точные значения коэффициентов переведем в значения с фиксированной точкой при длине машинного слова 16 бит (при этом 14 бит отводятся на запись единицы, то есть ).

1. Реакция на ступенчатое воздействие

Для анализа дельта-оператора были построены графики z-системы в форме прямого программирования, дельта-системы в форме совместного интегрирования и непрерывной системы. График ошибки моделирования цифровых систем представлен на рисунке 1.

 

Рисунок 1. Погрешность моделирования переходной характеристики цифровых систем относительно непрерывной системы

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

2. Реакция на гармоническое воздействие

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

Пусть период дискретизации и параметр; частота входного синусоидального сигнала выбрана (период равен 100 периодам дискретизации). Моделирование проводится при длине машинного слова 16 бит.

На рисунке 2 приведены графики погрешности моделирования z- и дельта-систем относительно выходного сигнала непрерывной системы.

 

Рисунок 2. Погрешность реакции цифровых систем относительно реакции непрерывной

Для дельта-системы используется совместное интегрирование как более точный метод. Ошибка в этом случае меньше на порядок, как и в случае реакции на единичный скачок.

Сравнения частотных характеристик z- и дельта-систем

Одной из часто используемых систем второго порядка является режекторный фильтр. Он задается передаточной функцией:

(17)

где - вырезаемая частота (рад/с), а - ширина полосы заграждения. На практике режекторный фильтр применяются чаще всего для отсечения помехи частотой 50±2 Гц. Соответственно, возьмем и . Тогда функция (34) примет вид:

(18)

Используя преобразование Тастина, найдем функцию в z-форме:

(19)

Произведем дельта-преобразование, опираясь на выражение (19):

(20)

Период дискретизации и параметр . Моделирование вновь проводится при длине машинного слова 16 бит.

 

Рисунок 3. АЧХ непрерывного и цифровых режекторных фильтров

Анализ графика АЧХ показывает, что z-система ведет себя неудовлетворительно:

а) подавляемая ею частота составляет 48 Гц, а не 50 Гц;

б) на частоте 51 Гц (входит в область подавления) ее коэффициент усиления больше единицы;

в) АЧХ неровная, имеет побочные пики (например, на частоте 93 Гц).

Таким образом, при 16-битном машинном представлении режекторный фильтр, построенный на основе z-преобразования, имеет неудовлетворительную АЧХ. В то же самое время АЧХ фильтра, построенного на основе дельта-преобразования, очень близка к АЧХ непрерывного фильтра. Более того, в работе [5] показано, что, используя дельта-преобразование, возможна реализация подобного режекторного фильтра даже на 8-битном контроллере.

Заключение

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

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

Преимущества дельта-преобразования отчетливо проявляются при малых периодах дискретизации. Так, с помощью дельта-преобразования был реализован цифровой режекторный фильтр, который показал значительно лучшие характеристики, чем фильтр на основе z-преобразования.

Работа выполнена при поддержке Министерства образования и науки РФ в рамках Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России на 2009–2013 годы» (Государственный контракт № 14.B37.21.2021 от 11 ноября 2012 г.).

Рецензенты:

Анисимов Владимир Иванович, д.т.н., профессор кафедры cистем автоматизированного проектирования. Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Санкт-Петербургский государственный электротехнический университет «ЛЭТИ» им. В. И. Ульянова (Ленина)», г. Санкт-Петербург.

Сольницев Ремир Иосифович, д.т.н., профессор кафедры cистем автоматизированного проектирования. Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Санкт-Петербургский государственный электротехнический университет «ЛЭТИ» им. В. И. Ульянова (Ленина)», г. Санкт-Петербург.


Библиографическая ссылка

Андреев В.С., Бутусов Д.Н., Каримов Т.И., Липкин С.М., Сотнин М.И. АНАЛИЗ ЭФФЕКТИВНОСТИ ПРИМЕНЕНИЯ ДЕЛЬТА-ПРЕОБРАЗОВАНИЯ ПРИ МОДЕЛИРОВАНИИ ЗВЕНЬЕВ ВТОРОГО ПОРЯДКА // Современные проблемы науки и образования. – 2013. – № 1.;
URL: http://science-education.ru/ru/article/view?id=8128 (дата обращения: 09.12.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074