phasez

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

Описание

[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.

Проектируйте lowpass equiripple фильтр с нормализованной частотой полосы пропускания 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.

Проектируйте эллиптический lowpass БИХ с нормализованной частотой полосы пропускания 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ω)A(ejω)=b (1)+b (2)ejω+b (3)ej2ω++b (M)ej(M1)ωа (1)+а (2)ejω+а (3)ej2ω++одинej(N1)ω.

Пример: 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 соответствует коэффициентам фильтра второго порядка (biquad). 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] задает фильтр Баттерворта третьего порядка с нормализованной частотой 3 дБ 0.5, рад/выборка.

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

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

Пример: d = designfilt('lowpassiir','FilterOrder',3,'HalfPowerFrequency',0.5) задает фильтр Баттерворта третьего порядка с нормализованной частотой 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