MATLAB.Exponenta
MATLAB и Simulink на русском
Технологии разработки и отладки
		сложных технических систем

MATLAB\MATLAB

В.Г.Потемкин "Справочник по MATLAB"
Анализ и обработка данных

  В оглавление книги \ К следующему разделу \ К предыдущему разделу

Свертка и фильтрация

UNWRAP
Корректировка фазовых углов

Синтаксис:

            Q = unwrap(P)
            Q = unwrap(P, cutoff)

Описание:

Функция Q = unwrap(P) корректирует фазовые углы элементов одномерного массива P при переходе через значение p, дополняя их значениями ±2p для того, чтобы убрать разрывы функции; если P - двумерный массив, то соответствующая функция применяется к столбцам.

Функция Q = unwrap(P, cutoff) позволяет пользователю изменить значение cutoff критического угла; по умолчанию cutoff = p .

Пример:

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

             image843.gif (407 bytes) .

Вычислим частотную характеристику этого фильтра в диапазоне w = 0.1:0.01:10, используя функцию freqresp пакета Control System Toolbox [1].

             w = 0.1:0.01:10;
             g = freqresp(num, den, sqrt(-1) * w);

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

             ph1 = (180./pi) * (atan2(imag(g), real(g)));
             ph2 = (180./pi) * unwrap(atan2(imag(g), real(g)));

w ph1 ph2
0.6000 126.7350 126.7350
0.7000 141.9270 141.9270
0.8000 -158.7123 201.2877
0.9000 -135.6888 224.3112
1.0000 135.0000 -225.0000
1.1000 -139.3435 220.6565
1.2000 -145.0993 214.9007
1.3000 -151.1794 208.8206
1.4000 -157.1667 202.8333
1.5000 -162.8839 197.1161
1.6000 -168.2575 191.7425
1.7000 -173.2642 186.7358
1.8000 -177.9068 182.0932
1.9000 177.7986 177.7986
2.0000 173.8298 173.8298

Из таблицы следует, что при входе и выходе из диапазона частот 0.8-2.0 Гц фазовые углы ph1 при достижении критического угла p терпят разрывы, которые устраняются функцией unwrap.

Построим фазовые частотные характеристики ph1(w) и ph2(w)-360, которые подтверждают сделанные выводы.

         semilogx(w, ph1), hold on, grid
         semilogx(w, ph2-360, 'b')

image844.gif (2910 bytes)

Сопутствующие функции: ANGLE, ABS, Control System Toolbox [1].

Ссылки:

1. Signal Processing Toolbox User’s Guide. Natick: The MathWorks, Inc., 1990

  В оглавление книги \ К следующему разделу \ К предыдущему разделу

 


Поиск по сайту:

Система Orphus

Яндекс.Метрика