Методы частотного анализа

Частотный анализ линейных и линеаризованных систем методом расчета матрицы Якоби системы.

Постановка задачи

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

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

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

Алгоритмы частотного анализа

При использовании метода переменных состояния расчет частотных характеристик сводится к выполнению следующих действий:

1. Сформировать математическую модель в виде совместной системы дифференциальных, алгебраических и разностных уравнений.
2. Линеаризовать полученные уравнения. При этом линеаризация будет происходить в окрестности заданной рабочей точки для дифференцирования в численном или аналитическом виде.
3. Выполнить замену операторов дифференцирования, сдвига по времени и запаздывания на соответствующие комплексные коэффициенты усиления. В результате будет сформирована система линейных алгебраических уравнений с комплексными коэффициентами.
4. Решив полученную систему, найти значения действительной Re(ω) и мнимой Im(ω) частей комплексной частотной характеристики.
5. По значениям Re(ω), Im(ω) определить значения частотных характеристик.

Действия 3…5 выполняются для каждого из значений частоты, на которых рассчитываются частотные характеристики.

Более подробное описание каждого этапа расчета представлено ниже.

Тогда уравнения динамической системы запишутся в виде:

где DIS(T_dis) - оператор векторного покомпонентного дискретного запаздывания; T_dis - вектор соответствующих периодов дискретизации; DEL(T_del), EXT(T_вxt) - операторы покомпонентного запаздывания и экстраполяции; T_del, T_ext - векторы соответствующих постоянных времени.

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

Линеаризация полученных уравнений сводится к определению матриц частных производных:

где f = [f₁, f₂, f₃, f₄, f₅]^T, x = [xdif, x_alg, x_dis, x_del, x_вxt]^T. Расчет этих матриц производиться численно, заменяя дифференциалы конечными приращениями. Алгоритм расчета матрицы A и вектора-строки C представлен ниже с использованием оператора языка Паскаль:

{ x0 - заданная рабочая точка }
f0 := f(x0, u);
g0 := g(x0, u);
for i := 1 to N do
begin
{ relinc - относительное приращение, absinc - абсолютное приращение }
dxi := relinc*abs(x[i]) + absink;
xsave := x0[i];
x0[i] := xsave + dxi;
f1 := f (x0, u);
g1 := g (x0, u);
for j := 1 to N do A[j, i] := ( f1[j] - f0[j] ) / dxi;
C[j] := (g1[j] - g0[j] ) / dxi;
x0[i] := xsave
end;

Определение матриц B и D происходит аналогичным способом, путем задания приращения переменной u.

Далее следует произвести замену оператора дифференцирования на коэффициент jω, а операторов дискретного и непрерывного запаздывания - на exp(-jωT), учитывая при этом, что передаточная функция экстраполятора имеет вид (*). В результате получим систему алгебраических уравнений вида:

где E = diag[E_dif, E_alg, E_dis, E_del, E_вxt] - комплексная диагональная матрица, которая состоит из пяти диагональных матриц, ненулевые элементы которых равны:

Принимая u = 1 и решая систему (*), требуется найти значения вещественной и мнимой частотных характеристик на заданной частоте ω: y = Re(ω) + jIm(ω).

Решение алгебраической системы с комплексными коэффициентами сводится к решению системы с вещественными коэффициентами, которая имеет вид:

где E_I, E_R - действительная и мнимая части матрицы E. Тогда значения вещественной и мнимой характеристик определяются по формулам:

Значения амплитудной A(ω), логарифмической амплитудной L(ω) и фазовой j(ω) частотных характеристик вычисляются в виде:

где arg - аргумент комплексного числа, который вычисляется по формуле:

где k - целое число, задаваемое так, чтобы исключить разрывы фазовой характеристики.

Реализация в SimInTech

Примеры применения данной методики анализа систем приведены в примерах, расположенных в каталоге: \SimInTech\Demo\Графики и исследования\Частотный анализ\

Однако данный метод ограничен возможными некорректными расчетами характеристик нелинейных систем, поскольку нелинейные элементы при малых приращениях не имеют возможности менять своего выхода (например, реле и прочие разрывные нелинейности) при недостаточной амплитуде приращений. Поэтому в разных траекториях системы метод выдает принципиально разные результаты (т.е., например в один момент времени АЧХ будет нулевой, а в другой – ненулевой). Данный метод рекомендуется применять для частотного анализа линейных систем и для получения мгновенных частотных характеристик нелинейных систем. Кроме того, при этом методе затруднительно учесть влияние амплитуды входного сигнала на АЧХ системы (для нелинейных систем иногда принципиально).

Для использования данного метода построения частотных характеристик в библиотеке присутствует блок Построение частотных характеристик, который находится в закладке «Анализ и оптимизация».

Для того, чтобы построить частотную характеристику при помощи данного блока надо:

1. Создать модель системы.
2. Поставить на схему блок «Построение частотных характеристик»
3. Первый вход этого блока подключить на выход блока, который установлен в качестве источника входного воздействия системы.
4. Второй вход блока подключить на выход системы.
5. Установить диапазон расчета частотных характеристик и их типы в свойствах блока (для этого надо выделить блок, вызвать контекстное меню нажатием на блок правой кнопкой мыши и в контекстном меню выбрать пункт «Свойства»).
6. Нажать кнопку «Пуск» для расчета частотной характеристики системы и вывода ее на график.

Схема подключения блока «Построение частотных характеристик» приведена ниже:

Частотный анализ систем методом гармонического анализа при подаче на вход системы синусоидального испытательного сигнала.

Алгоритм частотного анализа

На вход системы подается гармонический сигнал u(t) = a ⋅ sin(ωt). Установившийся сигнал на выходе системы возможно представить в виде суммы гармоник:

Тогда амплитуды и фазы гармоник рассчитываются по формулам:

Значения частотной характеристики находим по 1-й гармонике:

Интервал оценки АЧХ выбирается из расчета t_k = kT, а значения характеристик оцениваются на каждом интервале [t₁, t₂],...[t_k, t_k+1] до тех пор, пока разность оценок между соседними интервалами не будет удовлетворять заданной точности. Чтобы уменьшить время переходного процесса, на 1-м периоде [0, T] можно принять u(t) = at/T.

Реализация в SimInTech

Для построения частотной характеристики системы методом гармонического анализа в библиотеке блоков есть блок Построение частотных характеристик, который находится в закладке «Анализ и опитимизация». Блок представляет собой субмодель, реализующую алгоритм прямого расчета частотной характеристики произвольной системы при подаче на вход синусоидального сигнала с дискретно изменяющейся частотой и заданной амплитудой колебаний.

Пример использования данного блока находится в файле: \SimInTech\Demo\Графики и исследования\Частотный анализ\Гармонический анализ\Динамическая система\Анализ динамической системы.prt

Ниже показана схема включения блока для исследования частотных характеристик системы:

Свойства блока:

• Количество интервалов АЧХ (nfpoints) – количество участков с дискретно меняющейся частотой
• Минимальная круговая частота, 1/с (fmin) – начальная круговая частота
• Максимальная круговая частота, 1/с (fmax) – конечная круговая частота
• Количество периодов стабилизации (nstab) – множитель показывающий сколько периодов синусоидального сигнала на входе системы будет занимать измерение АЧХ на одном частоте.
• Амплитуда сигнала (Ampl) – амплитуда синусоидального сигнала на входе системы.

Для получения частотной характеристики необходимо запустить систему на моделирование. По окончании расчета частотных характеристик блок сам сформирует команду на остановку моделирования. Конечное время расчета в окне параметров расчета проекта необходимо задать большим числом, поскольку время для расчета полной АЧХ рассчитывается автоматически. Ниже приведен результат расчета АЧХ и ФЧХ колебательного звена данным методом при параметрах звена:

• Коэффициенты усиления k = [1].
• Постоянные времени T = 0.1.
• Коэффициент демпфирования b = 0.1.

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

Спектральный анализ выходного сигнала системы методом подачи на вход системы шумоподобного сигнала или ступенчатого сигнала.

Алгоритм спектрального анализа

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

Реализация в SimInTech

Для оценки спектральных характеристик систем в SimInTech существует блок Спектральная плотность, который находится в закладке «Статистика». Пример подключения данного блока приведен на рисунке ниже и реализован в демо-примерах каталога \SimInTech\Demo\Обработка сигналов\Спектральная плотность\.

Пример подключения данного блока для анализа выходного спектра системы при случайном входном воздействии приведен ниже:

На вход исследуемой системы подается шумоподобный сигнал в диапазоне значений от −0.5 до 0.5 с математическим ожиданием 0 и размахом 1. На выходе системы подключен блок расчёта спектральной характеристики, который оценивает спектр системы по выборке в 512 точек (настраивается в свойствах блока) с шагом 0.1 сек, с усреднением полученных спектров по всем сериям измерений.

Ниже приведен результат расчета спектральной характеристики системы (колебательное звено, параметры которого описаны выше, входной сигнал - шум).

Пример модели оценки АЧХ системы путем анализа спектра отклика на ступенчатое воздействие приведен ниже:

Результат расчета спектральной плотности отклика системы при ступенчатом воздействии:

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

Заключение

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