phasedelay

Задержка фазы цифрового фильтра

Описание

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

пример

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

[phi,w] = phasedelay(d,n) возвращает n- укажите ответ задержки фазы цифрового фильтра d.

пример

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

[phi,f] = phasedelay(___,n,fs) возвращает ответ задержки фазы и соответствующий n- укажите вектору частоты f для цифрового фильтра, спроектированного, чтобы отфильтровать сигналы, произвел на уровне fs.

[phi,f] = phasedelay(___,n,'whole',fs) возвращает вектор частоты f в n точки, располагающиеся между 0 и fs.

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

phi = phasedelay(___,f,fs) возвращает ответ задержки фазы, оцененный на частотах, заданных в f.

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

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

пример

phasedelay(___) строит ответ задержки фазы по сравнению с частотой.

Примеры

свернуть все

Используйте метод наименьших квадратов с ограничениями, чтобы спроектировать КИХ-фильтр lowpass порядка 54 и нормированной частоты среза 0.3. Задайте неравномерность в полосе пропускания и затухание в полосе задерживания как 0,02 и 0.08, соответственно, описанный в линейных модулях. Вычислите и постройте ответ задержки фазы фильтра.

Ap = 0.02;
As = 0.008;

b = fircls1(54,0.3,Ap,As);
phasedelay(b)

Figure contains an axes object. The axes object contains an object of type line.

Повторите пример с помощью designfilt. Следует иметь в виду, что эта функция описывает пульсации в децибелах.

Apd = 40*log10((1+Ap)/(1-Ap));
Asd = -20*log10(As);

d = designfilt('lowpassfir','FilterOrder',54,'CutoffFrequency',0.3, ...
               'PassbandRipple',Apd,'StopbandAttenuation',Asd);
phasedelay(d)

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

Спроектируйте третий порядок Фильтр Баттерворта lowpass с частотой среза 200 Гц. Частота дискретизации составляет 1 000 Гц.

fc = 200;
fs = 1000;

[z,p,k] = butter(3,fc/(fs/2),'low');

Используйте zp2sos функционируйте, чтобы преобразовать нули, полюса и усиление к секциям второго порядка. Вычислите ответ задержки фазы фильтра и определите номер точек оценки к 1024. Отобразите результат.

sos = zp2sos(z,p,k);
phasedelay(sos,1024)

Figure contains an axes object. The axes object contains an object of type line.

Спроектируйте эллиптический фильтр порядка 10 и нормированной частоты полосы пропускания 0.4. Задайте неравномерность в полосе пропускания 0,5 дБ и затухание в полосе задерживания 20 дБ. Отобразите ответ задержки фазы фильтра по полному модульному кругу.

[b,a] = ellip(10,0.5,20,0.4); 
phasedelay(b,a,512,'whole')

Figure contains an axes object. The axes object contains an object of type line.

Повторите пример с помощью designfilt.

d = designfilt('lowpassiir','DesignMethod','ellip','FilterOrder',10, ...
               'PassbandFrequency',0.4, ...
               'PassbandRipple',0.5,'StopbandAttenuation',20);
phasedelay(d,512,'whole')

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

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

свернуть все

Коэффициенты передаточной функции в виде векторов.

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

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

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

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

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

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

Угловые частоты, по поводу которых функция оценивает ответ задержки фазы в виде вектора и описала в рад/отсчете. Частоты обычно между 0 и π. w должен содержать по крайней мере два элемента.

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

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

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

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

свернуть все

Ответ задержки фазы, возвращенный как вектор из длины n. Ответ задержки фазы оценен в n равномерно распределенные точки вокруг верхней половины модульного круга.

Примечание

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

Угловые частоты в рад/отсчете, возвращенном как вектор. Если вы задаете nW имеет длину n. Если вы не задаете n или вы задаете n как пустой вектор, затем w имеет длину 512.

Частоты в герц, возвращенном как вектор. Если вы задаете nF имеет длину n. Если вы не задаете n или вы задаете n как пустой вектор, затем f имеет длину 512.

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

Алгоритмы

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

H(ejω)=B(ejω)A(ejω)=b(1)+b(2)ejω++b(m+1)ejmωa(1)+a(2)ejω++a(n+1)ejnω

учитывая числитель и коэффициенты знаменателя во входных параметрах b и a.

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

| | | | |

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