реклама на сайте
подробности

 
 
 
Reply to this topicStart new topic
> Как отобразить АЧХ fir1 матлаб
Acvarif
сообщение Apr 13 2018, 16:27
Сообщение #1


Знающий
****

Группа: Участник
Сообщений: 998
Регистрация: 27-08-08
Пользователь №: 39 850



Подсобите пожалуйста с пониманием работы функции fir1
Например
Код
b = fir1(50, 0.1);

Здесь http://194.81.104.27/~brian/DSP/FIRFilterM...gnPractical.pdf говорится что в b будут возвращены коэффициенты НЧ фильтра 50 порядка.
Но вот со вторым членом функции совсем не врубаюсь.
Если необходима частота среза 100 Гц то какое значение должен иметь второй член функции и как сделать plot АЧХ фильтра?
Go to the top of the page
 
+Quote Post
ataradov
сообщение Apr 14 2018, 05:13
Сообщение #2


Профессионал
*****

Группа: Участник
Сообщений: 1 014
Регистрация: 8-01-07
Из: San Jose, CA
Пользователь №: 24 202



https://www.mathworks.com/help/signal/ref/fir1.html

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


Там же есть примеры как строить.

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

Сообщение отредактировал ataradov - Apr 14 2018, 05:16
Go to the top of the page
 
+Quote Post
Lmx2315
сообщение Apr 14 2018, 09:54
Сообщение #3


отэц
*****

Группа: Свой
Сообщений: 1 699
Регистрация: 18-09-05
Из: Москва
Пользователь №: 8 684



Цитата(Acvarif @ Apr 13 2018, 19:27) *
Подсобите пожалуйста с пониманием работы функции fir1
Например
Код
b = fir1(50, 0.1);


Если необходима частота среза 100 Гц то какое значение должен иметь второй член функции и как сделать plot АЧХ фильтра?


как уже было сказано - частота среза зависит от тактовой частоты фир фильтра, в вашем случае 0.1*Ft.
Чтобы построить АЧХ фильтра , нужно через ваш фир фильтр с полученными выше коэффициентами пропустить дельтафункцию и взять от результата БПФ . В качестве дельта функции можно применить единичный всплеск длительностью один отсчёт с большой амплитудой .


--------------------
b4edbc0f854dda469460aa1aa a5ba2bd36cbe9d4bc8f92179f 8f3fec5d9da7f0
SHA-256
Go to the top of the page
 
+Quote Post
ViKo
сообщение Apr 14 2018, 10:12
Сообщение #4


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 476
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Цитата(Lmx2315 @ Apr 14 2018, 12:54) *
нужно через ваш фир фильтр с полученными выше коэффициентами пропустить дельтафункцию

Импульсная характеристика фильтра называется.
Go to the top of the page
 
+Quote Post
Acvarif
сообщение Apr 14 2018, 15:58
Сообщение #5


Знающий
****

Группа: Участник
Сообщений: 998
Регистрация: 27-08-08
Пользователь №: 39 850



Цитата(ataradov @ Apr 14 2018, 08:13) *
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');

Частота дискретизации понятно. Но при чем тут частота Найквиста?
Go to the top of the page
 
+Quote Post
Grizzzly
сообщение Apr 14 2018, 20:16
Сообщение #6


Знающий
****

Группа: Свой
Сообщений: 537
Регистрация: 22-02-13
Пользователь №: 75 748



Цитата(Acvarif @ Apr 14 2018, 18:58) *
Частота дискретизации понятно. Но при чем тут частота Найквиста?

Частота Найквиста есть половина частоты дискретизации. В зарубежных источниках в основном её принято использовать для нормировки.
Go to the top of the page
 
+Quote Post
Stanislav
сообщение Apr 21 2018, 00:52
Сообщение #7


Гуру
******

Группа: Свой
Сообщений: 4 278
Регистрация: 13-05-05
Из: Москва
Пользователь №: 4 987



Цитата(Acvarif @ Apr 14 2018, 18:58) *
......
Частота дискретизации понятно. Но при чем тут частота Найквиста?
Попробуйте поиграться сперва с fdatool.
Это крайне полезное занятие, для уяснения сущности фильтров.

В расчётах ЦФ, как и писали коллеги, часто используется нормированная частота (1 соответствует частоте дискретизации).
То есть, по большому счёту, частоты в привычном понимании там нет. sm.gif
Нормированное значение частоты 0,5 для фильтра с вещественным коэффициентами соответствует точке зеркальной симметрии АЧХ фильтра.
В отечественной литературе она называлась частотой свёртывания, а в зарубежной - частотой Найквиста.


--------------------
Там, где начинается свобода слова, заканчивается свобода мысли.©
Go to the top of the page
 
+Quote Post

Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 21st July 2018 - 03:36
Рейтинг@Mail.ru


Страница сгенерированна за 0.12075 секунд с 7
ELECTRONIX ©2004-2016