acvarif 0 13 апреля, 2018 Опубликовано 13 апреля, 2018 · Жалоба Подсобите пожалуйста с пониманием работы функции fir1 Например b = fir1(50, 0.1); Здесь http://194.81.104.27/~brian/DSP/FIRFilterM...gnPractical.pdf говорится что в b будут возвращены коэффициенты НЧ фильтра 50 порядка. Но вот со вторым членом функции совсем не врубаюсь. Если необходима частота среза 100 Гц то какое значение должен иметь второй член функции и как сделать plot АЧХ фильтра? Цитата Поделиться сообщением Ссылка на сообщение Поделиться на другие сайты Поделиться
ataradov 0 14 апреля, 2018 Опубликовано 14 апреля, 2018 (изменено) · Жалоба https://www.mathworks.com/help/signal/ref/fir1.html If Wn is a scalar, then fir1 designs a lowpass or highpass filter with cutoff frequency Wn. The cutoff frequency is the frequency at which the normalized gain of the filter is –6 dB. Там же есть примеры как строить. Для цифровых фильтров конкретная частота зависит от частоты дискретизации, поэтому значение задается в радиан/отсчет в нормализованном виде. Изменено 14 апреля, 2018 пользователем ataradov Цитата Поделиться сообщением Ссылка на сообщение Поделиться на другие сайты Поделиться
Lmx2315 2 14 апреля, 2018 Опубликовано 14 апреля, 2018 · Жалоба Подсобите пожалуйста с пониманием работы функции fir1 Например b = fir1(50, 0.1); Если необходима частота среза 100 Гц то какое значение должен иметь второй член функции и как сделать plot АЧХ фильтра? как уже было сказано - частота среза зависит от тактовой частоты фир фильтра, в вашем случае 0.1*Ft. Чтобы построить АЧХ фильтра , нужно через ваш фир фильтр с полученными выше коэффициентами пропустить дельтафункцию и взять от результата БПФ . В качестве дельта функции можно применить единичный всплеск длительностью один отсчёт с большой амплитудой . Цитата Поделиться сообщением Ссылка на сообщение Поделиться на другие сайты Поделиться
ViKo 1 14 апреля, 2018 Опубликовано 14 апреля, 2018 · Жалоба нужно через ваш фир фильтр с полученными выше коэффициентами пропустить дельтафункцию Импульсная характеристика фильтра называется. Цитата Поделиться сообщением Ссылка на сообщение Поделиться на другие сайты Поделиться
acvarif 0 14 апреля, 2018 Опубликовано 14 апреля, 2018 · Жалоба https://www.mathworks.com/help/signal/ref/fir1.html Там же есть примеры как строить. Для цифровых фильтров конкретная частота зависит от частоты дискретизации, поэтому значение задается в радиан/отсчет в нормализованном виде. Нашел. Да. cutoff frequency Wn определяется частотой дискретизации и частотой Найквиста filtord = 48; % Filter Order Fs = 144000; % Sampling Frequency Fn = Fs/2; % Nyquist Frequency Fc = 100; % Cutoff Frequency b = fir1(filtord, Fc/Fn); % Calculate Filter Coefficients [h,f]=freqz(b,1); plot(f*fs/(2*pi),20*log10(abs(h))) xlabel('frequency/Hz');ylabel('gain/dB');title('The gain response of lowpass filter'); Частота дискретизации понятно. Но при чем тут частота Найквиста? Цитата Поделиться сообщением Ссылка на сообщение Поделиться на другие сайты Поделиться
Grizzly 0 14 апреля, 2018 Опубликовано 14 апреля, 2018 · Жалоба Частота дискретизации понятно. Но при чем тут частота Найквиста? Частота Найквиста есть половина частоты дискретизации. В зарубежных источниках в основном её принято использовать для нормировки. Цитата Поделиться сообщением Ссылка на сообщение Поделиться на другие сайты Поделиться
Stanislav 0 21 апреля, 2018 Опубликовано 21 апреля, 2018 · Жалоба ...... Частота дискретизации понятно. Но при чем тут частота Найквиста? Попробуйте поиграться сперва с fdatool. Это крайне полезное занятие, для уяснения сущности фильтров. В расчётах ЦФ, как и писали коллеги, часто используется нормированная частота (1 соответствует частоте дискретизации). То есть, по большому счёту, частоты в привычном понимании там нет. :) Нормированное значение частоты 0,5 для фильтра с вещественным коэффициентами соответствует точке зеркальной симметрии АЧХ фильтра. В отечественной литературе она называлась частотой свёртывания, а в зарубежной - частотой Найквиста. Цитата Поделиться сообщением Ссылка на сообщение Поделиться на другие сайты Поделиться