Какие функции реализуют операторы freqs feedback

freqz

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

[ h , w ] = freqz( sos , n ) возвращает n — укажите комплексную частотную характеристику, соответствующую матрице секций второго порядка sos .

[ h , w ] = freqz( d , n ) возвращает n — укажите комплексную частотную характеристику для цифрового фильтра d .

[ h , w ] = freqz( ___ , n ,’whole’) возвращает частотную характеристику в n точки выборки вокруг целого модульного круга.

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

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

h = freqz( ___ , w ) возвращает вектор частотной характеристики h оцененный на нормированных частотах предоставляется в w .

h = freqz( ___ , f , fs ) возвращает вектор частотной характеристики h оцененный на физических частотах предоставляется в f .

freqz( ___ ) без выходных аргументов строит частотную характеристику фильтра.

Примеры

Частотная характеристика от передаточной функции

Вычислите и отобразите ответ величины БИХ третьего порядка фильтр lowpass, описанный следующей передаточной функцией:

H ( z ) = 0 . 0 5 6 3 4 ( 1 + z — 1 ) ( 1 — 1 . 0 1 6 6 z — 1 + z — 2 ) ( 1 — 0 . 6 8 3 z — 1 ) ( 1 — 1 . 4 4 6 1 z — 1 + 0 . 7 9 5 7 z — 2 ) .

Опишите числитель и знаменатель как полиномиальные свертки. Найдите частотную характеристику в 2 001 точке, охватывающей полный модульный круг.

b0 = 0.05634; b1 = [1 1]; b2 = [1 -1.0166 1]; a1 = [1 -0.683]; a2 = [1 -1.4461 0.7957]; b = b0*conv(b1,b2); a = conv(a1,a2); [h,w] = freqz(b,a,'whole',2001);

Постройте ответ величины, описанный в децибелах.

plot(w/pi,20*log10(abs(h))) ax = gca; ax.YLim = [-100 20]; ax.XTick = 0:.5:2; xlabel('Normalized Frequency (\times\pi rad/sample)') ylabel('Magnitude (dB)')

Частотная характеристика от секций второго порядка

H ( z ) = 0 . 0 5 6 3 4 ( 1 + z — 1 ) ( 1 — 1 . 0 1 6 6 z — 1 + z — 2 ) ( 1 — 0 . 6 8 3 z — 1 ) ( 1 — 1 . 4 4 6 1 z — 1 + 0 . 7 9 5 7 z — 2 ) .

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

b0 = 0.05634; b1 = [1 1]; b2 = [1 -1.0166 1]; a1 = [1 -0.683]; a2 = [1 -1.4461 0.7957]; sos1 = [b0*[b1 0] [a1 0]]; sos2 = [b2 a2]; [h,w] = freqz([sos1;sos2],'whole',2001);

Постройте ответ величины, описанный в децибелах.

plot(w/pi,20*log10(abs(h))) ax = gca; ax.YLim = [-100 20]; ax.XTick = 0:.5:2; xlabel('Normalized Frequency (\times\pi rad/sample)') ylabel('Magnitude (dB)')

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

Частотная характеристика КИХ-фильтра

Спроектируйте КИХ фильтр lowpass порядка 80 с помощью окна Кайзера с β = 8 . Задайте нормированную частоту среза 0 . 5 π рад/отсчет. Отобразите величину и фазовые отклики фильтра.

b = fir1(80,0.5,kaiser(81,8)); freqz(b,1)

Figure contains 2 axes objects. Axes object 1 contains an object of type line. Axes object 2 contains an object of type line.

Спроектируйте тот же фильтр с помощью designfilt . Отобразите его величину и фазовые отклики с помощью fvtool .

d = designfilt('lowpassfir','FilterOrder',80, . 'CutoffFrequency',0.5,'Window','kaiser',8>); freqz(d)

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

Частотная характеристика КИХ-полосового фильтра

Спроектируйте КИХ-полосовой фильтр с полосой пропускания между 0 . 3 5 π и 0 . 8 π рад/отсчет и 3 дБ пульсации. Первая полоса задерживания идет от 0 к 0 . 1 π рад/отсчет и имеет затухание 40 дБ. Вторая полоса задерживания идет от 0 . 9 π рад/отсчет к частоте Найквиста и имеет затухание 30 дБ. Вычислите частотную характеристику. Постройте его величину и в линейных модулях и в децибелах. Подсветите полосу пропускания.

sf1 = 0.1; pf1 = 0.35; pf2 = 0.8; sf2 = 0.9; pb = linspace(pf1,pf2,1e3)*pi; bp = designfilt('bandpassfir', . 'StopbandAttenuation1',40, 'StopbandFrequency1',sf1,. 'PassbandFrequency1',pf1,'PassbandRipple',3,'PassbandFrequency2',pf2, . 'StopbandFrequency2',sf2,'StopbandAttenuation2',30); [h,w] = freqz(bp,1024); hpb = freqz(bp,pb); subplot(2,1,1) plot(w/pi,abs(h),pb/pi,abs(hpb),'.-') axis([0 1 -1 2]) legend('Response','Passband','Location','South') ylabel('Magnitude') subplot(2,1,2) plot(w/pi,db(h),pb/pi,db(hpb),'.-') axis([0 1 -60 10]) xlabel('Normalized Frequency (\times\pi rad/sample)') ylabel('Magnitude (dB)')

Figure contains 2 axes objects. Axes object 1 contains 2 objects of type line. These objects represent Response, Passband. Axes object 2 contains 2 objects of type line.

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

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

Коэффициенты передаточной функции в виде векторов. Опишите передаточную функцию в терминах b и a как

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

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

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

n — Количество точек оценки
512 (значений по умолчанию) | положительный целочисленный скаляр

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

sos — Коэффициенты секции второго порядка
матрица

Коэффициенты секции второго порядка в виде матрицы. 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
Поддержка комплексного числа: Да

d — Цифровой фильтр
digitalFilter объект

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

Пример: d = designfilt(‘lowpassiir’,’FilterOrder’,3,’HalfPowerFrequency’,0.5) задает третий порядок Фильтр Баттерворта с нормированной частотой на 3 дБ 0.5π рад/отсчет.

fs — Частота дискретизации
положительная скалярная величина

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

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

w — Угловые частоты
вектор

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

f — Частоты
вектор

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

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

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

h — Частотная характеристика
вектор

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

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

w — Угловые частоты
вектор

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

f — Частоты
вектор

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

Алгоритмы

Частотная характеристика цифрового фильтра может быть интерпретирована как передаточная функция, оцененная в z = e jω [1].

freqz определяет передаточную функцию из (действительный или комплексный) полином числителя и полином знаменателя, который вы задаете, и возвращает комплексную частотную характеристику, H (e jω ), цифрового фильтра. Частотная характеристика оценена в точках выборки, определенных синтаксисом, который вы используете.

freqz обычно использует Алгоритм бпф, чтобы вычислить частотную характеристику каждый раз, когда вы не предоставляете вектор из частот как входной параметр. Это вычисляет частотную характеристику как отношение преобразованного числителя и коэффициентов знаменателя, дополненных нулями к желаемой длине.

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

Ссылки

[1] Оппенхейм, Алан V, Рональд В. Шафер и Джон Р. Бак. Обработка сигналов дискретного времени. 2-й Эд. Верхний Сэддл-Ривер, NJ: Prentice Hall, 1999.

Расширенные возможности

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.

Указания и ограничения по применению:

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

Какие функции реализуют операторы freqs feedback

[ h , w ] = freqz( b,a , n ) returns the n -point frequency response vector, h , and the corresponding angular frequency vector, w , for the digital filter with numerator and denominator polynomial coefficients stored in b and a , respectively.

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

[ h , w ] = freqz(___, n ,’whole’) returns the frequency response at n sample points around the entire unit circle.

[ h , w ] = freqz(___, n ,’whole’) Возвращает частотную характеристику n точек выборки вокруг всего единичного круга.

Frequency Response from Transfer Function

Compute and display the magnitude response of the third-order IIR lowpass filter described by the following transfer function:

Express the numerator and denominator as polynomial convolutions. Find the frequency response at 2001 points spanning the complete unit circle.

[ h , w ] = freqz( d , n ) returns the n -point complex frequency response for the digital filter, d .

Frequency Response of an FIR Bandpass Filter

Design an FIR bandpass filter with passband between and rad/sample and 3 dB of ripple. The first stopband goes from 0 to rad/sample and has an attenuation of 40 dB. The second stopband goes from rad/sample to the Nyquist frequency and has an attenuation of 30 dB. Compute the frequency response. Plot its magnitude in both linear units and decibels. Highlight the passband.
sf1 = 0.1;
pf1 = 0.35;
pf2 = 0.8;
sf2 = 0.9;
pb = linspace(pf1,pf2,1e3)*pi;

bp = designfilt(‘bandpassfir’, .
‘StopbandAttenuation1’,40, ‘StopbandFrequency1’,sf1.
‘PassbandFrequency1′,pf1,’PassbandRipple’,3,’PassbandFrequency2′,pf2, .
‘StopbandFrequency2′,sf2,’StopbandAttenuation2’,30);

[h,w] = freqz(bp,1024);
hpb = freqz(bp,pb);

subplot(2,1,1)
plot(w/pi,abs(h),pb/pi,abs(hpb),’.-‘)
axis([0 1 -1 2])
legend(‘Response’,’Passband’,’Location’,’South’)
ylabel(‘Magnitude’)

subplot(2,1,2)
plot(w/pi,db(h),pb/pi,db(hpb),’.-‘)
axis([0 1 -60 10])
xlabel(‘Normalized Frequency (\times\pi rad/sample)’)
ylabel(‘Magnitude (dB)’)

freqz(___) with no output arguments plots the frequency response of the filter.

Frequency Response of an FIR filter

Design an FIR lowpass filter of order 80 using a Kaiser window with . Specify a normalized cutoff frequency of rad/sample. Display the magnitude and phase responses of the filter.

b = fir1(80,0.5,kaiser(81,8));
freqz(b,1)

Design the same filter using designfilt. Display its magnitude and phase responses using fvtool.
d = designfilt(‘lowpassfir’,’FilterOrder’,80, .
‘CutoffFrequency’,0.5,’Window’, );
freqz(d)

Note: If the input to freqz is single precision, the frequency response is calculated using single-precision arithmetic. The output, h , is single precision.

[ h , w ] = freqz( sos , n ) returns the n -point complex frequency response corresponding to the second-order sections matrix, sos .

[ h , f ] = freqz(___, n , fs ) returns the frequency response vector, h , and the corresponding physical frequency vector, f , for the digital filter with numerator and denominator polynomial coefficients stored in b and a , respectively, given the sample rate, fs .

[ h , f ] = freqz(___, n ,’whole’, fs ) returns the frequency at n points ranging between 0 and fs .

h = freqz(___, w ) returns the frequency response vector, h , at the normalized frequencies supplied in w .

h = freqz(___, f , fs ) returns the frequency response vector, h , at the physical frequencies supplied in f .

Проектирование цифровых фильтров в MATLAB с использованием возможностей пакета Signal Processing Toolbox

Частотная характеристика дискретного фильтра рассчитывается с помощью функции freqz. В простейшем виде функция freqz используется следующим образом:

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

При расчете используются нормированные значения частот, измеряемые в радианах на отсчет (при такой нормировке частота дискретизации равна , а частота Найквиста — ). По умолчанию выбирается 512 частотных точек, равномерно распределенных в диапазоне .

При отсутствии выходных параметров функция freqz строит графики АЧХ (в децибелах) и ФЧХ (в градусах) фильтра.

Построение графика частотных характеристик фильтра с использованием функции freqz производится следующим образом:

Полученное в результате выполнения команд графическое окно с графиками АЧХ и ФЧХ цифрового фильтра имеет вид, приведенный на рисунке.

Рисунок — графическое окно с графиками АЧХ и ФЧХ цифрового фильтра, полученного с использованием функции freqz(b,a)

Если при вызове функции указаны выходные параметры, построение графика не производится:

В этом случае в векторе h возвращаются рассчитанные значения комплексного коэффициента передачи, а в векторе w – использованные при расчете значения нормированных частот.

Чтобы задать количество частотных точек для расчета, используется третий входной параметр n (частоты при этом по прежнему равномерно распределяются в диапазоне ):

Построение графика частотных характеристик фильтра в этом случае производится следующим образом:

Рисунок — графическое окно с графиками АЧХ и ФЧХ цифрового фильтра, полученного с использованием функции freqz(b,a,5)

Чтобы перейти от использования нормированных частот к использованию абсолютных значений частот в Гц, используется четвертый входной параметр Fs, равный частоте дискретизации (частоты при этом равномерно распределяются в диапазоне ):

Построение графика частотных характеристик фильтра в этом случае производится следующим образом:

Рисунок — графическое окно с графиками АЧХ и ФЧХ цифрового фильтра, полученного с использованием функции freqz(b,a,512, 2000)

Можно заставить функцию freqz использовать ряд частот, равномерно распределенных на интервале вплоть до частоты дискретизации. Это является необходимым при анализе фильтров с комплексными коэффициентами, частотные характеристики которых не являются симметричными относительно частоты Найквиста. Такое указание дается функции freqz с помощью четвертого входного параметра – строки ‘whole’:

freqz(b,a, n, ‘whole’);

Есть возможность задать вектор нормированных круговых частот w для расчета частотной характеристики:

3.1.1. Проектирование цифровых фильтров с бесконечной импульсной характеристикой

3.1.1.1. Использование аналоговых прототипов

Для проектирования цифровых фильтров с бесконечной импульсной характеристикой (БИХ, IIF) используются функции расчета фильтров, учитывающие аналоговые прототипы [3]:

[b,a]=butter(n,Wn,’ftype’) (фильтр Баттерворта)

[b,a]=cheby1(n,Rp,Wn,’ftype’) (фильтр Чебышева I)

[b,a]=cheby2(n,Rs,Wn,’ftype’) (фильтр Чебышева II)

[b,a]=ellip(n,Rp,Rs,Wn,’ftype’) (эллиптические фильтры).

Фильтры Баттерворта отличаются максимально гладкой и плоской АЧХ в полосе пропускания, однако спады и подъемы частотной характеристики этих фильтров менее крутые, чем у фильтров Чебышева и эллиптических фильтров, особенно при малых порядках [4]. Фильтры Чебышева I рода характеризуются равновеликими пульсациями АЧХ в полосе пропускания и монотонным изменением ослабления в полосе задержания. Фильтры Чебышева II рода характеризуются равновеликими пульсациями АЧХ в полосе задержания и монотонным изменением ослабления в полосе пропускания. Они не обеспечивают такого же крутого спада АЧХ, как фильтры I рода, однако в полосе пропускания их АЧХ меняются монотонно, т.е. отсутствуют пульсации. Эллиптические фильтры характеризуются более крутыми спадами и подъемами АЧХ, чем фильтры Баттерворта и Чебышева. Их АЧХ имеют равновеликие пульсации в полосе пропускания и в полосе задержания.

Параметр Wn задает частоту среза ФНЧ, если это скаляр. Следует помнить, что при обращении к указанным функциям при проектировании цифровых фильтров частота среза Wn задается в единицах нормированной частоты в диапазоне от 0 до 1, где 1 соответствует половине частоты дискретизации или частоте Найквиста. Частота среза Wn соответсвует значению амплитуды, равному . Если его задать как двухкомпонентный вектор Wn=[w1,w2], то функция будет возвращать характеристики полосового фильтра. При отсутствии параметра ’ftype’ получаются характеристики ФНЧ или полосового фильтра. Параметр ’ftype’ имеет два значения, что позволяет проектировать еще два важных типа фильтров: high (ФВЧ), stop (режекторный с заданием Wn=[w1,w2]).

Параметр Rp задает неравномерность АЧХ в децибелах в полосе пропускания. Параметр Rs задает минимальное затухание АЧХ в децибелах в полосе задержания.

3.1.1.2. Использование функций maxflat и yulewalk

В группу функций проектирования фильтров с БИХ и КИХ включены еще две функции [3]:

Функция maxflat задает проектирование обобщенного фильтра Баттерворта, у которого передаточная функция представлена полиномами числителя и знаменателя с разной степенью: nb и na соответственно. Частота среза этого фильтра Wn задается на уровне – 3 дБ. Она находится между 0 и 1, где 1 соответствует частоте Найквиста.

Расчёт цифрового фильтра с использованием функции maxflat осуществляется следующим образом:

b = 0.0642 0.3009 0.5436 0.4514 0.1407 -0.0156 -0.0097 0.0021

a = 1.0000 0.0000 0.4776

Полученное в результате выполнения команды графическое окно с графиками АЧХ синтезированного цифрового фильтра имеет вид, приведенный на рисунке.

Рисунок – графическое окно с графиками АЧХ цифрового фильтра ФНЧ, синтезированного с использованием функции [b,a]=maxflat(nb,na,Wn,’plots’)

Функция yulewalk служит для проектирования рекурсивного цифрового фильтра Юла-Уокера порядка n с заданной АЧХ, которая представляется векторами частот f и значений m. Частоты должны быть заданы в интервале , где 1 соответствует частоте Найквиста (частота, равная половине частоты дискретизации). Вектор частот должен иметь первый элемент 0 и последний 1, а в промежутке частоты указываются в возрастающем порядке. Если надо задать скачок АЧХ, то частота скачка указывается подряд дважды, а в векторе m указываются значения АЧХ до скачка и после скачка.

Рисунок – графическое окно с графиками АЧХ цифрового фильтра, синтезированного с использованием функции [b,a]=yulewalk(10,f,m)

Использование функции yulewalk иллюстрируется следующим примером:

В последних строках фрагмента строятся идеальная и реальная АЧХ фильтра, что облегчает их сравнение: рисунок.

3.1.1.3. Выбор минимального порядка фильтров с БИХ

В пакет Signal Processing Toolbox включены 4 функции [3], которые обеспечивают вычисление минимального порядка и граничной частоты Wn фильтров с БИХ по следующим заданным параметрам: Wp – частота пропускания (среза), Ws – частота задержания, Rp – максимально допустимые пульсации в области пропускания в децибелах и Rs – минимальное затухание в полосе задержания в децибелах. К этим функциям относятся:

Для ВЧ фильтров должно соблюдаться условие: Wp > Ws. Для полосовых и режекторных фильтров Wp и Ws должны задаваться двухэлементными векторами, элементы которых задают граничные частоты полос. При этом первой указывается более низкая частота. В этом случае параметр Wn также будет двухэлементным вектором.

Величина Wp нормирована и должна находиться в диапазоне , где 1 соответсвует частоте Найквиста или половине частоты дискретизации [2]. Величина Ws находится в том же диапазоне, что и Wp. Кроме минимального порядка также вычисляется частота среза Wn, которая далее может быть использована в качестве исходной величины для расчета параметров фильтров.

Расчет ФНЧ Баттерворта, обеспечивающего затухание не более 3 дБ в диапазоне от 0 до 100 Гц и не менее 15 дБ в диапазоне от 150 Гц до частоты Найквиста (частота дискретизации 1000 Гц) выполняется следующим образом >> Wp=100/500;

Полученное в результате выполнения команды графическое окно с графиками АЧХ и ФЧХ синтезированного цифрового фильтра имеет вид, приведенный на рисунке.

Рисунок – графическое окно с графиками АЧХ и ФЧХ цифрового фильтра, синтезированного с использованием функции [b,a]=butter(n,Wn)

3.1.1.4. Методы преобразования аналоговых фильтров в цифровые

Для преобразования аналоговых фильтров в цифровые используются функции[3]:

[numd,dend]=bilinear(num,den,fs) (метод билинейного преобразования)

[numd,dend]=bilinear(num,den,fs,Fp) (метод билинейного преобразования)

[bz,az]=impinvar(b,a,fs,tol]) (метод инвариантности имп. характеристики)

При методе билинейного преобразования передаточная характеристика аналогового прототипа преобразуется в передаточную функцию цифрового фильтра в соответствии с выражением:

где fs – частота дискретизации.

Параметр Fp является дополнительным и задает предискажения. Он определяет частоту в Гц, для которой значения АЧХ до и после преобразования совпадают. В этом случае билинейное преобразование видоизменяется следующим образом:

При методе инвариантного импульсного преобразования осуществляется преобразование аналогового фильтра в цифровой фильтр, имеющий ту же импульсную характеристику, что и у аналогового фильтра с точностью до множителя 1/f_s.

3.1.2. Проектирование цифровых фильтров с конечной импульсной характеристикой

3.1.2.1. Метод Ремеза

Для проектирования цифровых фильтров с КИХ, имеющих равные пульсации как в полосе пропускания, так и в полосе задержания, по методу ремеза служит функция [3]:

Эта функция вычисляет вектор b, содержащий (n+1) коэффициентов фильтра с КИХ, АЧХ которого наилучшим образом представляет некоторую функцию fresp. Вектор f должен содержать пары частот, граничных для областей пропускания и задержания. Частоты должны быть в диапазоне от -1 до +1, причем +1 соответствует частоте Найквиста (половине частоты дискретизации). Частоты должны идти в строго нарастающем порядке. Параметр fresp может принимать различные значения: lowpass (ФНЧ), highpass (ФВЧ), bandpass (полосовой), bandstop (задерживающий), differentiator (дифференцирующий), hilfilt (Гильберта).

Расчёт фильтра ФНЧ 20-го порядка по методу Ремеза осуществляется следующим образом:

>> b=cremez(20,[-1 -0.4 -0.3 0.6 0.7 1],’lowpass’);

Полученное в результате выполнения второй команды графическое окно с графиками АЧХ и ФЧХ синтезированного цифрового фильтра ФНЧ имеет вид, приведенный на рисунке. Для вывода значений коэффициентов фильтра необходимо убрать знак ; после первой команды, либо дополнительно ввести команду b.

Рисунок – графическое окно с графиками АЧХ и ФЧХ цифрового фильтра ФНЧ, синтезированного с использованием функции b=cremez(20,[-1 -0.4 -0.3 0.6 0.7 1],’lowpass’)

3.1.2.2. Метод взвешивания

Для проектирования фильтров с КИХ методом взвешивания служит функция [3]:

В простейшей форме b=fir1(n,Wn) функция создает вектор из n+1 коэффициентов ФНЧ с КИХ и линейной АЧХ для заданной частоты среза Wn (от 0 до 1). Групповое время задержки таких фильтров равно n/2. При проектировании используется весовое окно Хемминга. Для проектирования полосового фильтра надо задать Wn как вектор из двух частот [w1,w2], причем w1 < w2. Параметр ‘ftype’ задает тип фильтра: high (ФВЧ), stop (заграждающий). Вектор window размером n+1 позволяет задавать отсчеты окна.

Рисунок – графическое окно с графиками АЧХ и ФЧХ цифрового полосового фильтра, синтезированного с использованием функции b=fir1(32,[0.3 0.7])

Использование функции fir1 для проектирования цифрового полосового фильтра 32-ого порядка с граничными частотами 0.3 и 0.7 осуществляется следующим образом:

Полученное в результате выполнения второй команды графическое окно с графиками АЧХ и ФЧХ синтезированного цифрового фильтра имеет вид, приведенный на рисунке. Из рисунка видно, что ФЧХ фильтра в полосе пропускания действительно носит строго линейный характер.

3.1.2.3. Метод взвешивания для фильтра с произвольной АЧХ

Для проектирования фильтров с КИХ и произвольной частотной характеристикой методом взвешивания с использованием окон служит функция [3]:

Эту функцию не рекомендуется применять для построения обычных фильтров: полосовых, режекторных, ФНЧ и ФВЧ.

Функция b=fir2(n,f,m) возвращает вектор с (n+1) коэффициентами ФНЧ порядка n. АЧХ фильтра задается векторами f (значения частот от 0 до 1) и m (отсчеты АЧХ для заданных в векторе f частот). Для ступенчатых АЧХ задаются два подряд отсчета частоты для ступеньки и два значения АЧХ. Параметр npt (по умолчанию 512) задает число отсчетов сетки для интерполяции, параметр lap (по умолчанию 25) задает число точек в области скачков АЧХ, а вектор window задает тип окна (по умолчанию окно Хемминга).

Использование функции fir2 для проектирования цифрового заграждающего фильтра 40-ого порядка с прямоугольной АЧХ и скачками на частотами 0.3 и 0.6 осуществляется следующим образом:

На рисунке показаны АЧХ идеального фильтра и фильтра, построенного с помощью функции fir2.

Рисунок – графическое окно с графиками АЧХ цифрового полосового фильтра, синтезированного с использованием функции b=fir2(40,f,m)

3.1.2.4. Метод наименьших квадратов

Метод наименьших квадратов в пакете Signal Processing Toolbox реализуют три функции:

Для проектирования фильтров с КИХ и произвольной АЧХ служит функция:

Функция задает проектирование фильтра порядка n, заданного (n+1) коэффициентами вектора b с АЧХ, которая задается своими отсчетами частоты f и значениями amp. Векторы up и lo имеют то же число элементов, что и вектор amp, и задают верхние и нижние границы ЧХ для каждой из полос фильтра.

Для проектирования фильтров ФНЧ и ФВЧ с КИХ методами наименьших квадратов функция:

В простейшей форме b=fircls1(n,W0,dp,ds) функция создает вектор коэффициентов ФНЧ порядка n с частотой среза W0 между 0 и 1, где 1 соответствует половине частоты дискретизации, максимально допустимой функцией ошибки АЧХ в полосе пропускания dp и максимальной функцией ошибки в полосе задержания ds. Параметр ‘high’ задает проектирование ФВЧ. Задается дополнительный параметр wt, выше которого или ниже которого выполняются требования к АЧХ проектируемого фильтра. Если wt в полосе пропускания, то , где — функция ошибки. Если wt вне полосы, то пропускания . Векторы up и lo имеют то же число элементов, что и вектор amp, и задают верхние и нижние границы ЧХ для каждой из полос фильтра.

Использование функции fircls1 для проектирования цифрового ФВЧ 60-ого порядка имеет следующий вид:

Полученная при вычислении АЧХ приведена на рисунке. Из рисунка наглядно видно, что АЧХ укладывается в заданные нормы отклонений dp и ds.

Рисунок – графическое окно с графиками АЧХ цифрового ФВЧ, синтезированного с использованием функции b=fircls1(n,w0,dp,ds,’high’,’plots’)

Функция b=fircls(n,f,a,w, ‘ftype’) также использует метод наименьших квадратов для построения фильтра с КИХ и линейной ФЧХ для заданной АЧХ самой разнообразной формы. АЧХ задается векторами частот f и значений нужной АЧХ a. В векторе w можно задать вес для аппроксимации АЧХ, причем число элементов вектора w равно половине числа элементов вектора f. Параметр ‘ftype’ позволяет задать тип фильтра: hilbert (фильтр Гильберта типов III и IV с линейной ФЧХ) и differentiator (фильтр с АЧХ, линейно нарастающей с частотой).

Использование функции fircls для проектирования цифрового фильтра Гильберта 25-ого порядка с трапецеидальной АЧХ, отсчеты которой заданы в векторах частот f и значений a имеет следующий вид:

Полученная при вычислении АЧХ приведена на рисунке.

Рисунок – графическое окно с графиками АЧХ фильтра Гильберта, синтезированного с использованием функции b=firls(25,F,A,’hilbert’)

freqs

h = freqs( b , a , w ) возвращает комплексную частотную характеристику аналогового фильтра, заданного векторами коэффициентов b и a , оцененный на угловых частотах w .

[ h , wout ] = freqs( b , a , n ) использование n частота указывает, чтобы вычислить h и возвращает соответствующие угловые частоты в wout .

freqs( ___ ) без выходных аргументов строит величину и фазовые отклики как функции угловой частоты в окне текущей фигуры. Можно использовать этот синтаксис с любым из предыдущих входных синтаксисов.

Примеры

Частотная характеристика от передаточной функции

Найдите и изобразите частотную характеристику в виде графика передаточной функции

H ( s ) = 0 . 2 s 2 + 0 . 3 s + 1 s 2 + 0 . 4 s + 1 .

Figure contains 2 axes objects. Axes object 1 contains an object of type line. Axes object 2 contains an object of type line.

Можно также сгенерировать графики путем вызова freqs без выходных аргументов.

Figure contains 2 axes objects. Axes object 1 contains an object of type line. Axes object 2 contains an object of type line.

Частотная характеристика аналога Lowpass фильтр Бесселя

Спроектируйте аналог 5-го порядка lowpass фильтр Бесселя с приблизительно постоянной групповой задержкой до 1 0 4 рад/с. Постройте частотную характеристику фильтра с помощью freqs .

Figure contains 2 axes objects. Axes object 1 contains an object of type line. Axes object 2 contains an object of type line.

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

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

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

Пример: [b,a] = butter(5,50,’s’) задает пятый порядок Фильтр Баттерворта с частотой среза 50 рад/секунда.

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

w — Угловые частоты
положительный вектор действительных чисел

Угловые частоты в виде положительного вектора действительных чисел описываются в рад/секунда.

Читать:

Как в 1с повысить оклад сотруднику

Пример: 2*pi*logspace(6,9) задает 50 логарифмически расположенных с интервалами угловых частот от 1 МГц ( 2π × 10 6 рад/секунда) и 1 ГГц ( 2π × 10 9 рад/секунда).

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

n — Количество точек оценки
200 (значений по умолчанию) | положительный целочисленный скаляр

Количество оценки указывает в виде положительного целочисленного скаляра.

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

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

h — Частотная характеристика
вектор

Частотная характеристика, возвращенная как вектор.

wout — Угловые частоты
вектор

Угловые частоты, в который h вычислен, возвращен как вектор.

Алгоритмы

freqs возвращает комплексную частотную характеристику аналогового фильтра, заданного b и a . Функция оценивает отношение полиномов Преобразования Лапласа

H ( s ) = B ( s ) A ( s ) = b ( 1 ) s n + b ( 2 ) s n − 1 + ⋯ + b ( n + 1 ) a ( 1 ) s m + a ( 2 ) s m − 1 + ⋯ + a ( m + 1 )

вдоль мнимой оси на частоте указывает s = jω :

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

Открытый пример

У вас есть модифицированная версия этого примера. Вы хотите открыть этот пример со своими редактированиями?

Документация Signal Processing Toolbox

Поддержка

1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.

2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.

3. Сохраняйте структуру оригинального текста — например, не разбивайте одно предложение на два.

4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.

5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.

Похожие публикации:

Как установить офис на планшет андроид бесплатно
Как установить принтер hp laserjet 1005 на windows 10
Как установить скоростной usb хост контроллер на windows xp
Как установить фубар 2000 на русском языке

Какие функции реализуют операторы freqs feedback

Функции matlab для расчета линейных систем (цепей).

Пакет Signal Processing содержит функции позволяющие выполнять расчет линейных систем. В частности они позволяют выполнить расчет комплексного коэффициента передачи, амплитудно-частотную и фазо-частотную характеристики системы, а также осуществить преобразование способов описания линейных систем.

Для расчета комплексного коэффициента передачи в функцию передачи необходимо подставить мнимый аргумент׃

S = Jω.

Соответствующие расчеты выполняются с помощью функции freqs. Простейший вид синтаксиса: freqs (b,a);

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

Для расчета характеристик по умолчанию выбирается 200 частот, логарифмически равномерно распределенных в диапазоне от 0,1 до 10.

При отсутствии выходных параметров функция freqs строит графики АЧХ и ФЧХ.

АЧХ выводится в логарифмическом масштабе (без пересчета в децибелы), ФЧХ – в градусах.

Если вместо графика нужно получить вектор рассчитанных значений комплексного коэффициента передачи, нужно присвоить результат, возвращаемый функцией freqs, какой либо переменной:

h = freqs (b,a);

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

[h, w] = freqs (b,a);

Можно построить вектор значений частотной характеристики АЧХ в линейном , а не в логарифмическом масштабе:

» [h, w] = freqs (b, a)

» plot (w, abs (h))

»grid on

Частотный диапазон можно задать принудительно также с помощью третьего входного параметра, которым является вектор круговых частот w׃

h = freqs (b, a, w);

Функция unwrap позволяет избавиться лишних разрывов ФЧХ. Она ищет в переданном ей векторе скачки между соседними элементами, превышающими заданную пороговую величину (по умолчанию π) и сдвигает соответствующие фрагменты вектора на ±2π нужное число раз

» [h, w] = freqs (b, a); % — ФЧХ

» ph_i = angle (h)

» ph_i = unwrap (ph_i) % — устранение скачков

» plot (w, ph_i * 180/ p_i) % отображается в градусах

» grig on

Функции преобразования способов описания систем.

В пакете Signal Processing системы MATLAB имеется ряд функций, предназначенных для преобразования описаний из одной формы в другую.

Имена функций имеют вид XX2YY, где XX – обозначение исходной формы описания, а YY – обозначение целевой формы описания системы.

Формы описания цепей в именах функций׃

— tf – (transfer function) – коэффициенты полинома числителя и знаменателя функции передачи H (s);

— zp – нули и полюсы (zeros and poles);

— ss – описание в пространстве состояний (state-space).

Необходимость в преобразовании описаний часто возникает из-за того, что функции, например фильтров, дают результат расчета в одной форме, а функция построения, например частотной характеристики, в другой форме.

Для входных и выходных параметров используются следующие обозначения:

а) функция передачи – в – вектор-строке коэффициентов (в порядке убывания степеней) числителя функции передачи;

а – вектор-строке коэффициентов (в порядке степеней) знаменателя функции передачи;

б) нули и полюсы׃

z – вектор нулей (столбец);

p – вектор полюсов (столбец);

K — коэффициент усиления (скаляр).

в) Пространство состояний׃

A— квадратная матрица связи вектора состояний и его производной;

B — вектор-столбец связи входных сигналов и производный вектора состояний;

C — вектор-строка связи выходных сигналов и вектора состояний;

D — скалярный коэффициент связи выходных и входных сигналов.

Дадим пояснения функций преобразования.

Функция tf2zp преобразует наборы коэффициентов полиномов числителя, знаменателя функции передачи в векторы нулей и полюсов, рассчитывая также значения общего коэффициента усиления

[z, p, K] = tf2zp (b, a);

Коэффициент усиления рассчитывается из соотношения

Функция zp2tf является обратной по отношению к функции tf2zp׃ она осуществляет преобразование коэффициента усиления, векторов нулей и полюсов функций передачи в коэффициенты полиномов её числителя и знаменателя׃

[b, a] = zp2tf (z, p, K);

Функция tf2ss преобразует наборы коэффициентов полиномов числителя и знаменателя функции передачи в параметры представления системы в пространстве состояний׃

[A, B, C, D] = tf2ss (b, a);

Функция ss2tf является обратной по отношению к функции tf2ss׃ она преобразует параметры пространства состояний в коэффициенты полиномов функции передачи системы׃

[b, a] = ss2tf (A, B, C, D);

Функция zp2ss преобразует нули, полюсы и коэффициенты усиления в ее параметры пространства состояний׃

[A, B, C, D] = zp2ss;

Функция ss2zp является обратной по отношению к предыдущей, т.е. она преобразует параметры пространства состояний в нули, полюсы и коэффициенты усиления׃

Русские Блоги

[MATLAB] Введение в функцию freqz (частотная характеристика цифрового фильтра)

freqz

Frequency response of digital filter

Syntax

Description

[ h , w ] = freqz(___, n ,’whole’) returns the frequency response at n sample points around the entire unit circle.

[ h , w ] = freqz(___, n ,’whole’) Возвращает частотную характеристику n точек выборки вокруг всего единичного круга.

Frequency Response from Transfer Function

Compute and display the magnitude response of the third-order IIR lowpass filter described by the following transfer function:

Express the numerator and denominator as polynomial convolutions. Find the frequency response at 2001 points spanning the complete unit circle.

[ h , w ] = freqz( d , n ) returns the n -point complex frequency response for the digital filter, d .

Frequency Response of an FIR Bandpass Filter

[h,w] = freqz(bp,1024);
hpb = freqz(bp,pb);

subplot(2,1,1)
plot(w/pi,abs(h),pb/pi,abs(hpb),’.-‘)
axis([0 1 -1 2])
legend(‘Response’,’Passband’,’Location’,’South’)
ylabel(‘Magnitude’)

subplot(2,1,2)
plot(w/pi,db(h),pb/pi,db(hpb),’.-‘)
axis([0 1 -60 10])
xlabel(‘Normalized Frequency (\times\pi rad/sample)’)
ylabel(‘Magnitude (dB)’)

freqz(___) with no output arguments plots the frequency response of the filter.

Frequency Response of an FIR filter

Design an FIR lowpass filter of order 80 using a Kaiser window with . Specify a normalized cutoff frequency of rad/sample. Display the magnitude and phase responses of the filter.

b = fir1(80,0.5,kaiser(81,8));
freqz(b,1)

Note: If the input to freqz is single precision, the frequency response is calculated using single-precision arithmetic. The output, h , is single precision.

[ h , w ] = freqz( sos , n ) returns the n -point complex frequency response corresponding to the second-order sections matrix, sos .

[ h , f ] = freqz(___, n ,’whole’, fs ) returns the frequency at n points ranging between 0 and fs .

h = freqz(___, w ) returns the frequency response vector, h , at the normalized frequencies supplied in w .

h = freqz(___, f , fs ) returns the frequency response vector, h , at the physical frequencies supplied in f .

Какие функции реализуют операторы freqs feedback

Frequency response of analog filters

Description

freqs returns the complex frequency response H(j) (Laplace transform) of an analog filter

given the numerator and denominator coefficients in vectors b and a .

h = freqs(b,a,w) returns the complex frequency response of the analog filter specified by coefficient vectors b and a . freqs evaluates the frequency response along the imaginary axis in the complex plane at the angular frequencies in rad/sec specified in real vector w , which must contain more than one frequency.

[h,w] = freqs(b,a) automatically picks a set of 200 frequency points w on which to compute the frequency response h .

[h,w] = freqs(b,a,f) picks f number of frequencies on which to compute the frequency response h .

freqs with no output arguments plots the magnitude and phase response versus frequency in the current figure window.

freqs works only for real input systems and positive frequencies.

Find and graph the frequency response of the transfer function given by:

You can also create the plot with:

To convert to hertz, degrees, and decibels, use:

freqs evaluates the polynomials at each frequency point, then divides the numerator response by the denominator response:

[MATLAB] Введение в функцию freqz (частотная характеристика цифрового фильтра)

freqz

Frequency response of digital filter

Syntax

Description

[ h , w ] = freqz(___, n ,’whole’) returns the frequency response at n sample points around the entire unit circle.

Frequency Response from Transfer Function

Compute and display the magnitude response of the third-order IIR lowpass filter described by the following transfer function:

Express the numerator and denominator as polynomial convolutions. Find the frequency response at 2001 points spanning the complete unit circle.

[ h , w ] = freqz( d , n ) returns the n -point complex frequency response for the digital filter, d .

Frequency Response of an FIR Bandpass Filter

[h,w] = freqz(bp,1024);
hpb = freqz(bp,pb);

subplot(2,1,2)
plot(w/pi,db(h),pb/pi,db(hpb),’.-‘)
axis([0 1 -60 10])
xlabel(‘Normalized Frequency (\times\pi rad/sample)’)
ylabel(‘Magnitude (dB)’)

freqz(___) with no output arguments plots the frequency response of the filter.

Frequency Response of an FIR filter

Design an FIR lowpass filter of order 80 using a Kaiser window with . Specify a normalized cutoff frequency of rad/sample. Display the magnitude and phase responses of the filter.

b = fir1(80,0.5,kaiser(81,8));
freqz(b,1)

Note: If the input to freqz is single precision, the frequency response is calculated using single-precision arithmetic. The output, h , is single precision.

[ h , w ] = freqz( sos , n ) returns the n -point complex frequency response corresponding to the second-order sections matrix, sos .

[ h , f ] = freqz(___, n ,’whole’, fs ) returns the frequency at n points ranging between 0 and fs .

h = freqz(___, w ) returns the frequency response vector, h , at the normalized frequencies supplied in w .

h = freqz(___, f , fs ) returns the frequency response vector, h , at the physical frequencies supplied in f .

How to use `feedback` function in Matlab?

Matlab’s feedback function is used to obtain the closed loop transfer function of a system. Example:

enter image description here

sys = feedback(sys1,sys2) returns a model object sys for the negative feedback interconnection of model objects sys1 , sys2 . To compute the closed-loop system with positive feedback, use sign = +1, for negative feedback we use -1.

My question arises when we have a system of the following type:

enter image description here

According to these docs, we can use feedback to create the negative feedback loop with G and C.

This is a source of confusion, shouldn’t the above be: sys = feedback(G*C,1,-1)? These are not the same.

However, looking at these docs, for a unit loop gain k, you can compute the closed-loop transfer function T using:

enter image description here

Why are we using 1 and not -1 ? This is still negative feedback and not positive feedback.

Какие функции реализуют операторы freqs feedback

Frequency response of digital filter

Syntax

Description

[ h , w ] = freqz(___, n ,’whole’) returns the frequency response at n sample points around the entire unit circle.

Frequency Response from Transfer Function

Compute and display the magnitude response of the third-order IIR lowpass filter described by the following transfer function:

Express the numerator and denominator as polynomial convolutions. Find the frequency response at 2001 points spanning the complete unit circle.

[ h , w ] = freqz( d , n ) returns the n -point complex frequency response for the digital filter, d .

Frequency Response of an FIR Bandpass Filter

[h,w] = freqz(bp,1024);
hpb = freqz(bp,pb);

subplot(2,1,2)
plot(w/pi,db(h),pb/pi,db(hpb),’.-‘)
axis([0 1 -60 10])
xlabel(‘Normalized Frequency (\times\pi rad/sample)’)
ylabel(‘Magnitude (dB)’)

freqz(___) with no output arguments plots the frequency response of the filter.

Frequency Response of an FIR filter

Design an FIR lowpass filter of order 80 using a Kaiser window with . Specify a normalized cutoff frequency of rad/sample. Display the magnitude and phase responses of the filter.

b = fir1(80,0.5,kaiser(81,8));
freqz(b,1)

Note: If the input to freqz is single precision, the frequency response is calculated using single-precision arithmetic. The output, h , is single precision.

[ h , w ] = freqz( sos , n ) returns the n -point complex frequency response corresponding to the second-order sections matrix, sos .

[ h , f ] = freqz(___, n ,’whole’, fs ) returns the frequency at n points ranging between 0 and fs .

h = freqz(___, w ) returns the frequency response vector, h , at the normalized frequencies supplied in w .

h = freqz(___, f , fs ) returns the frequency response vector, h , at the physical frequencies supplied in f .

Проектирование цифровых фильтров в MATLAB с использованием возможностей пакета Signal Processing Toolbox

b = 0.0642 0.3009 0.5436 0.4514 0.1407 -0.0156 -0.0097 0.0021

a = 1.0000 0.0000 0.4776

Похожие публикации:

Как отключить автозапуск торрента
Как сделать меню в тильде
Как смотреть нхл в россии после санкций
Как снять задний диван bmw x3

How do you use feedback frequency and timing to motivate and empower others?

Feedback frequency and timing affect how people process and act on the information you provide. Too much feedback can overwhelm or demoralize them, while too little feedback can leave them confused or uncertain. The right balance depends on several factors, such as the nature of the task, the level of skill, the stage of development, and the personality of the receiver. You want to give feedback that is timely, relevant, specific, and actionable, but also respectful, supportive, and encouraging.