
- •Анализ устойчивости и расчет переходного процесса дискретной системы
- •230100 - Информатика и вычислительная техника 230101 - Вычислительные машины, комплексы, системы и сети 230102 - Автоматизированные системы обработки информации и управления
- •Оглавление
- •Введение
- •1. Описание системы
- •2. Порядок выполнения домашнего задания
- •3. Методические указания
- •Библиографический список
- •Анализ устойчивости и расчет переходного процесса дискретной системы
- •620002, Г. Екатеринбург, ул. Мира, 19
2. Порядок выполнения домашнего задания
1. Для нескорректированной системы найти дискретную передаточную функцию разомкнутой системы Wнс(z) и Wнс().
2. Определить критическое значение передаточного коэффициента непрерывной части нескорректированной системы Ккр.
3. Рассчитать весовую решетчатую функцию разомкнутой нескорректированной системы wнс[i] при К = 0,5Ккр. Построить график функции wнс(t).
4. Рассчитать и построить логарифмические псевдочастотные характеристики нескорректированной разомкнутой системы при К = 0,5Ккр. Подтвердить правильность расчетов.
5. Рассчитать дискретную передаточную функцию последовательного корректирующего звена Wk(z), обеспечивающего в замкнутой системе переходный процесс минимальной длительности при отсутствии установившейся ошибки. Составить уравнение алгоритма цифровой реализации корректирующего звена.
6. Рассчитать и построить переходную функцию h(t) замкнутой скорректированной системы.
3. Методические указания
Указания даются по каждому из пунктов задания, при этом сохраняется порядок их следования. Анализ и синтез цифровой системы рекомендуется проводить методами теории импульсных САУ.
1. Импульсным элементом в данной системе следует считать ЦАП, так как на его выходе ступенчато изменяющийся сигнал u(t) эквивалентен последовательности прямоугольных импульсов с длительностями, равными периоду квантования Т0.
Следовательно, ЦАП работает как фиксирующий импульсный элемент и осуществляет амплитудно-импульсную модуляцию первого рода с единичным коэффициентом преобразования. Такой импульсный элемент называют также экстраполятором нулевого порядка. Его описание приводится к расчетной схеме с идеальным импульсным элементом (ИИЭ) и формирующим звеном, который имеет передаточную функцию
Расчетная схема нескорректированной системы представлена на рис. 2. Здесь Wф(s) и W(s) - передаточные функции формирующего звена и непрерывной части соответственно.
На выходе ИИЭ указан сигнал *(t), который представляет последовательность идеальных импульсов, интенсивности (площади) и полярности которых определяются решетчатой функцией [i]:
*(t)
=
i(t
– iT0)
Рис. 2. Расчетная схема нескорректированной системы
В работе необходимо использовать передаточную функцию приведенной непрерывной части нескорректированной системы
.
Пользуясь правилами определения дискретной передаточной функции разомкнутой нескорректированной системы с учетом формирователя и таблицей z-преобразований, можно получить выражение для Wнс(z) в общем виде:
где
К
- коэффициент непрерывной части (пока
численно не определен),
.
Таким образом структурную схему нескорректированной системы (для моментов съема) можно представить в виде рис. 3.
Рис. 3. Структурная схема нескорректированной системы
2. Применяя билинейное преобразование, найдите выражение передаточной функции разомкнутой нескорректированной дискретной системы в зависимости от переменной ρ:
Wнс(ρ)=Rнс(ρ)/Qнс(ρ).
Тогда характеристическое уравнение замкнутой системы можно записать в виде
Qз(ρ)=Rнс(ρ)+Qнс(ρ).
Определите величину Ккр, используя алгебраический критерий Гурвица.
3. Значения дискрет весовой решетчатой функции wнс[i] могут быть найдены по передаточной функции разомкнутой дискретной системы Wнс(z), которую необходимо для этого представить в виде числового ряда по отрицательным степеням z:
Wнс(z)= wнс[0]+wнс[1]z -1+wнс[2]z -2+… .
Такое представление можно получить, например, в результате численного деления полинома Rнс(z) на Qнс(z).
4. Для расчета псевдочастотных характеристик используйте выражение Wнс(ρ)=Rнс(ρ)/Qнс(ρ), которое целесообразно привести к виду
Wнс(ρ)= K(τ1ρ+1)(1-ρ)/2ρ(τ2ρ+1).
После подстановки К=0,5Ккр и ρ=i v, где v - псевдочастота, рассчитайте псевдочастотные характеристики. Постройте эти характеристики.
Подтвердите правильность расчетов, доказав по критерию Найквиста, что характеристики соответствуют устойчивой замкнутой системе, имеющей двукратный запас устойчивости по коэффициенту усиления.
5. Требуемая передаточная функция корректирующего эвена W к(z) определяется формулой
W к(z)= Qнс(z)/(zn Rнс(1) - Rнс(z)),
где Rнс(z) и Qнс(z) - полиномы передаточной функции разомкнутой дискретной системы W к(z),
n - степень полинома Qнс(z),
Rнс(1) = Rнс(z)| z=1.
Найдите передаточную функцию дискретного корректирующего звена и, разделив числитель и знаменатель полученного выражения на (z – 1), представьте W к(z) в виде
W к(z)=Ko(z+b1)/K(z+b2).
Здесь величины Ко, b1, b2 определяются численно, а коэффициент усиления непрерывной части К можно принять равным Ко, в результате чего получаем удобную для реализации функцию
W к(z)=(z+b1)/(z+b2).
Структурная схема скорректированной системы приведена на
рис. 4.
Рис. 4. Структурная схема скорректированной системы
С
целью получения алгоритма цифровой
реализации звена с передаточной функцией
W
к(z)
запишем
определяемые ею соотношения для
изображений
и
:
,
.
Перейдя к оригиналам - решетчатым функциям, найдите искомый алгоритм реализации цифрового корректирующего звена в виде формулы для вычисления текущих значений скорректированного сигнала ψ[i].