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π rad/s, неравномерность в полосе пропускания 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)

Спроектируйте тот же фильтр с помощью 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)

Спроектируйте lowpass equiripple фильтр с нормированной частотой полосы пропускания 0.45π rad/s, нормированная частота полосы задерживания 0.55π rad/s, неравномерность в полосе пропускания 1 дБ и затухание в полосе задерживания 60 дБ. Отобразите фазовый отклик фильтра.

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

Спроектируйте эллиптический БИХ-фильтр lowpass с нормированной частотой полосы пропускания 0.4π rad/s, нормированная частота полосы задерживания 0.5π rad/s, неравномерность в полосе пропускания 1 дБ и затухание в полосе задерживания 60 дБ. Отобразите фазовый отклик фильтра.

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

Входные параметры

свернуть все

Коэффициенты передаточной функции в виде векторов. Опишите передаточную функцию в терминах 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ω++(N)ej(N1)ω.

Пример: b = [1 3 3 1]/6 и a = [3 0 1 0]/3 задайте третий порядок Фильтр Баттерворта с нормированной частотой на 3 дБ 0.5π рад/выборка.

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

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

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

Коэффициенты секции второго порядка в виде матрицы. sos K-by-6 матрица, где количество разделов, K, должно быть больше или быть равно 2. Если количество разделов меньше 2, функция обрабатывает вход как вектор числителя. Каждая строка sos соответствует коэффициентам второго порядка (biquad) фильтр. i th строка 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 должен иметь по крайней мере два элемента, потому что в противном случае функция интерпретирует его как nW = π соответствует частоте Найквиста.

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

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

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

свернуть все

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

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

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

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

Смотрите также

| | | | |

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