Exponenta Banner
exponenta event banner

фазы

Фазовый отклик цифрового фильтра

свернуть все на странице

Синтаксис

[phi, w] = фазез (b, a, n)
[phi, w] = фазез (sos, n)
[phi, w] = фазез (d, n)
[phi, w] = фазез (___, n, 'целое')
[phi, f] = фазез (___, n, fs)
phi = фазез (___, f, fs)
phi = фазез (___, w)
фазез (___)

Описание

[phi,w] = phasez(b,a,n) возвращает значение n-точный вектор фазовой характеристики phi и соответствующий вектор угловой частоты w для цифрового фильтра с коэффициентами передаточной функции, сохраненными в b и a.

[phi,w] = phasez(sos,n) возвращает значение n-точная фазовая характеристика, соответствующая матрице сечений второго порядка sos.

[phi,w] = phasez(d,n) возвращает значение n-точный фазовый отклик для цифрового фильтра d.

[phi,w] = phasez(___,n,'whole') возвращает отклик фазы в n выборочные точки вокруг всей единичной окружности. Этот синтаксис может включать любую комбинацию входных аргументов из предыдущих синтаксисов.

[phi,f] = phasez(___,n,fs) возвращает частотный вектор.

phi = phasez(___,f,fs) возвращает вектор фазового отклика phi оценивается на физических частотах, поставляемых в f. Этот синтаксис может включать любую комбинацию входных аргументов из предыдущих синтаксисов.

phi = phasez(___,w) возвращает отклик развернутой фазы в радианах на частотах, указанных в w.

пример

phasez(___) без выходных аргументов отображает фазовый отклик фильтра.

Примеры

свернуть все

Открыть сценарий в реальном времени

Использовать designfilt для конструирования КИХ-фильтра порядка 54 нормированная частота отсечки 0,3δ рад/с, пульсация полосы пропускания 0,7 дБ и затухание полосы останова 42 дБ. Используйте метод ограничения наименьших квадратов. Отображение фазовой характеристики фильтра.

Nf = 54;
Fc = 0.3;
Ap = 0.7;
As = 42;

d = designfilt('lowpassfir','CutoffFrequency',Fc,'FilterOrder',Nf, ...
               'PassbandRipple',Ap,'StopbandAttenuation',As, ...
               'DesignMethod','cls');
phasez(d)

Figure Filter Visualization Tool - Phase Response contains an axes and other objects of type uitoolbar, uimenu. The axes with title Phase Response contains an object of type line.

Спроектируйте один и тот же фильтр с помощью fircls1. Помните, что fircls1 использует линейные единицы измерения для измерения пульсации и затухания.

pAp = 10^(Ap/40);
Apl = (pAp-1)/(pAp+1);

pAs = 10^(As/20);
Asl = 1/pAs;

b = fircls1(Nf,Fc,Apl,Asl);
phasez(b)

Figure contains an axes. The axes with title Phase Response contains an object of type line.

Открыть сценарий в реальном времени

Сконструировать низкочастотный фильтр с нормированной частотой полосы пропускания 0,45δ рад/с, нормализованной частотой полосы останова 0,55λ рад/с, пульсацией полосы пропускания 1 дБ и ослаблением полосы останова 60 дБ. Отображение фазовой характеристики фильтра.

d = designfilt('lowpassfir', ...
               'PassbandFrequency',0.45,'StopbandFrequency',0.55, ...
               'PassbandRipple',1,'StopbandAttenuation',60, ...
               'DesignMethod','equiripple');
phasez(d)

Figure Filter Visualization Tool - Phase Response contains an axes and other objects of type uitoolbar, uimenu. The axes with title Phase Response contains an object of type line.

Открыть сценарий в реальном времени

Сконструировать эллиптический БИХ-фильтр нижних частот с нормированной частотой полосы пропускания 0,4δ рад/с, нормализованной частотой полосы останова 0,5δ рад/с, пульсацией полосы пропускания 1 дБ и ослаблением полосы останова 60 дБ. Отображение фазовой характеристики фильтра.

d = designfilt('lowpassiir', ...
               'PassbandFrequency',0.4,'StopbandFrequency',0.5, ...
               'PassbandRipple',1,'StopbandAttenuation',60, ...
               'DesignMethod','ellip');
phasez(d)

Figure Filter Visualization Tool - Phase Response contains an axes and other objects of type uitoolbar, uimenu. The axes with title Phase Response contains an object of type line.

Входные аргументы

свернуть все

Коэффициенты передаточной функции, указанные как векторы. Выразить передаточную функцию с точки зрения b и a как

H (ejω) =B (ejω) (ejω) =b (1) +b (2) e−jω + b (3) e−j2ω + ⋯ + b (M) e−j (M−1) ωa (1) +a (2) e−jω + (3) e−j2ω + ⋯ + (N) e−j (N−1) ω.

Пример: b = [1 3 3 1]/6 и a = [3 0 1 0]/3 укажите фильтр Баттерворта третьего порядка с нормированной частотой 3 дБ 0,5δ рад/выборка.

Типы данных: double | single
Поддержка комплексного номера: Да

Число точек вычисления, указанное в виде положительного целого скаляра не менее 2. Когда n отсутствует, по умолчанию - 512. Для получения наилучших результатов установите n до значения, большего, чем порядок фильтра.

Коэффициенты сечения второго порядка, заданные как матрица. sos является K-by-6 матрицей, где число сечений, K, должно быть больше или равно 2. Если число сечений меньше 2, функция обрабатывает входные данные как числительный вектор. Каждая строка sos соответствует коэффициентам фильтра второго порядка (биквад). I-й ряд sos соответствует [bi(1) bi(2) bi(3) ai(1) ai(2) ai(3)].

Пример: s = [2 4 2 6 0 2;3 3 0 6 0 0] задает фильтр Butterworth третьего порядка с нормализованной частотой 3 дБ 0,5δ рад/выборка.

Типы данных: double | single
Поддержка комплексного номера: Да

Цифровой фильтр, указанный как digitalFilter объект. Использовать designfilt формируют цифровой фильтр на основе характеристик частотного отклика.

Пример: d = designfilt('lowpassiir','FilterOrder',3,'HalfPowerFrequency',0.5) задает фильтр Butterworth третьего порядка с нормализованной частотой 3 дБ 0,5δ рад/выборка.

Частота выборки, заданная как положительный скаляр. Если единица времени - секунды, fs выражается в герцах.

Типы данных: double

Угловые частоты, определенные как вектор и выраженные в рад/выборке. w должен иметь по крайней мере два элемента, поскольку в противном случае функция интерпретирует его как n. w = λ соответствует частоте Найквиста.

Частоты, заданные как вектор. f должен иметь по крайней мере два элемента, поскольку в противном случае функция интерпретирует его как n. Если единица времени - секунды, f выражается в герцах.

Типы данных: double

Выходные аргументы

свернуть все

Фазовый отклик, возвращенный как вектор. При указании n, то phi имеет длину n. Если не указать n, или укажите n как пустой вектор, то phi имеет длину 512.

Если вход в phasez - единичная точность, функция вычисляет фазовый отклик, используя арифметику с единичной точностью. Продукция phi - единичная точность.

Угловые частоты, возвращаемые в виде вектора. w имеет значения в диапазоне от 0 до δ. При указании 'whole' во входных данных значения в w диапазон от 0 до 2δ. При указании n, w имеет длину n. Если не указать n, или укажите n в качестве пустого вектора, w имеет длину 512.

Частоты, возвращаемые в виде вектора, выраженного в герцах. f имеет значения в диапазоне от 0 до fs/ 2 Гц. При указании'whole' во входных данных значения в f диапазон от 0 до fs Гц. При указании n, f имеет длину n. Если не указать n, или укажите n как пустой вектор, f имеет длину 512.

См. также

| | | | |

Представлен до R2006a

Документация по инструментам обработки сигналов

Поддержка