Спроектируйте цифровые фильтры
designfilt больше не помогает в исправлении вызовов функции в рамках скрипта или функции. Для получения дополнительной информации см. Вопросы совместимости.
d = designfilt(resp,Name,Value)digitalFilter объект, d, с ответом вводят resp. Примеры resp 'lowpassfir' и 'bandstopiir'. Задайте фильтр далее с помощью набора Аргументов в виде пар имя-значение. Позволенные наборы спецификации зависят от resp и состойте из комбинаций их:
Frequency Constraints соответствует частотам, на которых фильтр предоставляет желаемое поведение. Примеры включают 'PassbandFrequency' и 'CutoffFrequency'. Необходимо всегда задавать ограничения частоты.
Magnitude Constraints описывает поведение фильтра в конкретных частотных диапазонах. Примеры включают 'PassbandRipple' и 'StopbandAttenuation'. designfilt обеспечивает значения по умолчанию для ограничений величины, оставленных незаданными. В проектах произвольной величины необходимо всегда задавать векторы из желаемых амплитуд.
Filter Order. Некоторые методы разработки позволяют вам задать порядок. Другие производят проекты минимального порядка. Таким образом, они генерируют самые маленькие фильтры, которые удовлетворяют заданным ограничениям.
Design Method является алгоритмом, используемым, чтобы спроектировать фильтр. Примеры включают метод наименьших квадратов с ограничениями ('cls') и работа с окнами Кайзера ('kaiserwin'). Для некоторых наборов спецификации существует несколько методов разработки, доступных, чтобы выбрать из. В других случаях можно использовать только один метод, чтобы выполнить желаемым техническим требованиям.
Design Method Options является параметрами, характерными для данного метода разработки. Примеры включают 'Window' для 'window' метод и оптимизация 'Weights' для произвольной величины equiripple проекты. designfilt обеспечивает значения по умолчанию для проектных решений, оставленных незаданными.
Sample Rate является частотой, на которой действует фильтр. designfilt имеет частоту дискретизации по умолчанию 2 Гц. Используя это значение эквивалентно работе с нормированными частотами.
Примечание
Если вы задаете неполный или противоречивый набор аргументов name-value в командной строке, designfilt предложения открыть Ассистент Создания фильтра. Ассистент помогает вам спроектировать фильтр и вставляет откорректированный MATLAB® код по командной строке.
Если вы вызываете designfilt из скрипта или функции с неправильным набором технических требований, designfilt выдает ошибку сообщение со ссылкой, чтобы открыть Ассистент Создания фильтра. Ассистент помогает вам спроектировать фильтр и вставляет откорректированный код MATLAB по командной строке. Спроектированный фильтр сохранен в рабочую область.
Использование filter в форме dataOut = filter(d,dataIn) отфильтровать сигнал с digitalFilterD. Для БИХ-фильтров, filter функционируйте использует прямую форму II реализаций.
Используйте FVTool, чтобы визуализировать digitalFilterD.
Введите d.Coefficients получить коэффициенты digitalFilterD. Для БИХ-фильтров коэффициенты описываются как секции второго порядка.
Смотрите digitalFilter для списка функций фильтрации и анализа, доступных для использования с digitalFilter объекты.
designfilt( позволяет вам отредактировать существующий цифровой фильтр, d)d. Это открывает Ассистент Создания фильтра, заполненный с техническими требованиями фильтра, которые можно затем изменить. Это - единственный способ, которым можно отредактировать digitalFilter объект. Его свойства в противном случае только для чтения.
Спроектируйте минимальный порядок КИХ-фильтр lowpass с нормированной частотой полосы пропускания рад/отсчет, частота полосы задерживания рад/отсчет, неравномерность в полосе пропускания 0,5 дБ и затухание в полосе задерживания 65 дБ. Используйте окно Кайзера, чтобы спроектировать фильтр. Визуализируйте его ответ величины. Используйте его, чтобы отфильтровать вектор из случайных данных.
lpFilt = designfilt('lowpassfir','PassbandFrequency',0.25, ... 'StopbandFrequency',0.35,'PassbandRipple',0.5, ... 'StopbandAttenuation',65,'DesignMethod','kaiserwin'); fvtool(lpFilt)
dataIn = rand(1000,1); dataOut = filter(lpFilt,dataIn);
Спроектируйте БИХ-фильтр lowpass с порядком 8, частота полосы пропускания 35 кГц и неравномерность в полосе пропускания 0,2 дБ. Задайте частоту дискретизации 200 кГц. Визуализируйте ответ величины фильтра.
lpFilt = designfilt('lowpassiir','FilterOrder',8, ... 'PassbandFrequency',35e3,'PassbandRipple',0.2, ... 'SampleRate',200e3); fvtool(lpFilt)
Используйте фильтр, который вы спроектировали, чтобы отфильтровать случайный сигнал с 1000 выборками.
dataIn = randn(1000,1); dataOut = filter(lpFilt,dataIn);
Выведите коэффициенты фильтра, описанные как секции второго порядка.
sos = lpFilt.Coefficients
sos = 4×6
0.2666 0.5333 0.2666 1.0000 -0.8346 0.9073
0.1943 0.3886 0.1943 1.0000 -0.9586 0.7403
0.1012 0.2023 0.1012 1.0000 -1.1912 0.5983
0.0318 0.0636 0.0318 1.0000 -1.3810 0.5090
Спроектируйте минимальный порядок highpass КИХ-фильтр с нормированной частотой полосы задерживания рад/отсчет, частота полосы пропускания рад/отсчет, неравномерность в полосе пропускания 0,5 дБ и затухание в полосе задерживания 65 дБ. Используйте окно Кайзера, чтобы спроектировать фильтр. Визуализируйте его ответ величины. Используйте его, чтобы отфильтровать 1 000 выборок случайных данных.
hpFilt = designfilt('highpassfir','StopbandFrequency',0.25, ... 'PassbandFrequency',0.35,'PassbandRipple',0.5, ... 'StopbandAttenuation',65,'DesignMethod','kaiserwin'); fvtool(hpFilt)
dataIn = randn(1000,1); dataOut = filter(hpFilt,dataIn);
Спроектируйте highpass БИХ-фильтр с порядком 8, частота полосы пропускания 75 кГц и неравномерность в полосе пропускания 0,2 дБ. Задайте частоту дискретизации 200 кГц. Визуализируйте ответ величины фильтра. Примените фильтр к вектору с 1000 выборками из случайных данных.
hpFilt = designfilt('highpassiir','FilterOrder',8, ... 'PassbandFrequency',75e3,'PassbandRipple',0.2, ... 'SampleRate',200e3); fvtool(hpFilt)
dataIn = randn(1000,1); dataOut = filter(hpFilt,dataIn);
Спроектируйте КИХ-фильтр полосы пропускания 20-го порядка с более низкой частотой среза, 500 Гц и более высокой частотой среза 560 Гц. Частота дискретизации составляет 1 500 Гц. Визуализируйте ответ величины фильтра. Используйте его, чтобы отфильтровать случайный сигнал, содержащий 1 000 выборок.
bpFilt = designfilt('bandpassfir','FilterOrder',20, ... 'CutoffFrequency1',500,'CutoffFrequency2',560, ... 'SampleRate',1500); fvtool(bpFilt)
dataIn = randn(1000,1); dataOut = filter(bpFilt,dataIn);
Выведите коэффициенты фильтра.
b = bpFilt.Coefficients
b = 1×21
-0.0113 0.0067 0.0125 -0.0445 0.0504 0.0101 -0.1070 0.1407 -0.0464 -0.1127 0.1913 -0.1127 -0.0464 0.1407 -0.1070 0.0101 0.0504 -0.0445 0.0125 0.0067 -0.0113
Спроектируйте БИХ-фильтр полосы пропускания 20-го порядка с более низкой частотой на 3 дБ, 500 Гц и более высокой частотой на 3 дБ 560 Гц. Частота дискретизации составляет 1 500 Гц. Визуализируйте частотную характеристику фильтра. Используйте его, чтобы отфильтровать случайный сигнал с 1000 выборками.
bpFilt = designfilt('bandpassiir','FilterOrder',20, ... 'HalfPowerFrequency1',500,'HalfPowerFrequency2',560, ... 'SampleRate',1500); fvtool(bpFilt)
dataIn = randn(1000,1); dataOut = filter(bpFilt,dataIn);
Спроектируйте 20-й порядок bandstop КИХ-фильтр с более низкой частотой среза, 500 Гц и более высокой частотой среза 560 Гц. Частота дискретизации составляет 1 500 Гц. Визуализируйте ответ величины фильтра. Используйте его, чтобы отфильтровать 1 000 выборок случайных данных.
bsFilt = designfilt('bandstopfir','FilterOrder',20, ... 'CutoffFrequency1',500,'CutoffFrequency2',560, ... 'SampleRate',1500); fvtool(bsFilt)
dataIn = randn(1000,1); dataOut = filter(bsFilt,dataIn);
Спроектируйте 20-й порядок bandstop БИХ-фильтр с более низкой частотой на 3 дБ, 500 Гц и более высокой частотой на 3 дБ 560 Гц. Частота дискретизации составляет 1 500 Гц. Визуализируйте ответ величины фильтра. Используйте его, чтобы отфильтровать 1 000 выборок случайных данных.
bsFilt = designfilt('bandstopiir','FilterOrder',20, ... 'HalfPowerFrequency1',500,'HalfPowerFrequency2',560, ... 'SampleRate',1500); fvtool(bsFilt)
dataIn = randn(1000,1); dataOut = filter(bsFilt,dataIn);
Спроектируйте фильтр дифференциатора полной полосы порядка 7. Отобразите его нулевой фазовый отклик. Используйте его, чтобы отфильтровать вектор с 1000 выборками из случайных данных.
dFilt = designfilt('differentiatorfir','FilterOrder',7); fvtool(dFilt,'MagnitudeDisplay','Zero-phase')
dataIn = randn(1000,1); dataOut = filter(dFilt,dataIn);
Спроектируйте Гильбертов трансформатор порядка 18. Задайте нормированную ширину перехода рад/отсчет. Отобразите в линейных модулях ответ величины фильтра. Используйте его, чтобы отфильтровать вектор с 1000 выборками из случайных данных.
hFilt = designfilt('hilbertfir','FilterOrder',18,'TransitionWidth',0.25); fvtool(hFilt,'MagnitudeDisplay','magnitude')
dataIn = randn(1000,1); dataOut = filter(hFilt,dataIn);
Вам дают сигнал, произведенный на уровне 1 кГц. Спроектируйте фильтр, который останавливает частоты между 100 Гц и 350 Гц и частоты, больше, чем 400 Гц. Задайте порядка фильтра 60. Визуализируйте частотную характеристику фильтра. Используйте его, чтобы отфильтровать случайный сигнал с 1000 выборками.
mbFilt = designfilt('arbmagfir','FilterOrder',60, ... 'Frequencies',0:50:500,'Amplitudes',[1 1 1 0 0 0 0 1 1 0 0], ... 'SampleRate',1000); fvtool(mbFilt)
dataIn = randn(1000,1); dataOut = filter(mbFilt,dataIn);
Если вы задаете неполный или противоречивый набор расчетных параметров, designfilt предложения открыть Ассистент Создания фильтра.
(В описании аргумента для resp существует полный список допустимых наборов спецификации для всех доступных типов ответа.)
Ассистент ведет себя по-другому, если вы вызываете designfilt в командной строке или в рамках скрипта или функции.
Вам дают сигнал, произведенный на уровне 2 кГц. Вас просят спроектировать КИХ-фильтр lowpass, который подавляет частотные составляющие выше, чем 650 Гц. “Частота среза” походит на хорошего кандидата на параметр спецификации. Вы вводите этот код в командной строке MATLAB.
Fsamp = 2e3; Fctff = 650; dee = designfilt('lowpassfir','CutoffFrequency',Fctff, ... 'SampleRate',Fsamp);
Что-то, кажется, неправильно, потому что это диалоговое окно появляется на вашем экране.
Вы нажимаете Yes и получаете новое диалоговое окно, которое предлагает генерировать код. Вы видите, что переменные, которые вы задали прежде, были вставлены, где ожидается.
После исследования некоторых опций предложил, вы решаете протестировать откорректированный фильтр. Вы нажимаете OK и получаете этот код по командной строке.
designfilt('lowpassfir','FilterOrder', 10, ... 'CutoffFrequency',Fctff,'SampleRate',2000);
Введение имени фильтра повторяет информацию от диалогового окна.
dee
dee =
digitalFilter with properties:
Coefficients: [-0.0036 0.0127 -0.0066 -0.0881 0.2595 ...
0.6521 0.2595 -0.0881 -0.0066 0.0127 -0.0036]
Specifications:
FrequencyResponse: 'lowpass'
ImpulseResponse: 'fir'
SampleRate: 2000
CutoffFrequency: 650
FilterOrder: 10
DesignMethod: 'window'
Use fvtool to visualize filter
Use designfilt to edit filter
Use filter to filter dataВы вызываете FVTool и получаете график deeчастотная характеристика.
fvtool(dee)
Сокращение не выглядит особенно резким. Ответ выше 40 дБ для большинства частот. Вы помните, что у ассистента была опция, чтобы настроить “ограничение величины”, вызвал “затухание в полосе задерживания”. Откройте ассистент путем вызова designfilt с фильтром называют как вход.
designfilt(dee)
Нажмите Magnitude constraints выпадающее меню и выбирает Passband ripple and stopband attenuation. Вы видите, что метод разработки изменился от Window к FIR constrained least-squares. Значение по умолчанию для затухания составляет 60 дБ, который выше, чем 40. Нажмите OK и визуализируйте получившийся фильтр.
dee = designfilt('lowpassfir','FilterOrder',10, ... 'CutoffFrequency',650,'PassbandRipple',1, ... 'StopbandAttenuation',60,'SampleRate',2000); fvtool(dee)
Сокращение все еще не выглядит резким. Затухание - действительно 60 дБ, но для частот выше 900 Гц.
Снова вызовите designfilt с вашим фильтром, как введено.
designfilt(dee)
Ассистент вновь появляется.
Чтобы сузить различие между принятыми и отклоненными частотами, увеличьте порядок фильтра или измените Frequency constraints от Cutoff (6dB) frequency к Passband and stopband frequencies. Если вы изменяете порядка фильтра от 10 до 50, вы получаете более резкий фильтр.
dee = designfilt('lowpassfir','FilterOrder',50, ... 'CutoffFrequency',650,'PassbandRipple',1, ... 'StopbandAttenuation',60,'SampleRate',2000); fvtool(dee)
Немного экспериментирования показывает, что можно получить подобный фильтр путем установки полосы пропускания и частот полосы задерживания соответственно к 600 Гц и 700 Гц.
dee = designfilt('lowpassfir','PassbandFrequency',600, ... 'StopbandFrequency',700,'PassbandRipple',1, ... 'StopbandAttenuation',60,'SampleRate',2000); fvtool(dee)
Вам дают сигнал, произведенный на уровне 2 кГц. Вас просят спроектировать фильтр highpass, который останавливает частоты ниже 700 Гц. Вы не заботитесь о фазе сигнала, и необходимо работать с фильтром младшего разряда. Таким образом БИХ-фильтр кажется соответствующим. Вы не уверены, какой порядок фильтра является лучшим, таким образом, вы пишете функцию, которая принимает порядок как вход. Откройте редактор MATLAB и создайте файл.
function dataOut = hipassfilt(Order,dataIn) hpFilter = designfilt('highpassiir','FilterOrder',N); dataOut = filter(hpFilter,dataIn); end
Чтобы протестировать вашу функцию, создайте сигнал, состоявший из двух синусоид с частотами 500 и 800 Гц, и сгенерируйте выборки в течение 0,1 с. Фильтр 5-го порядка кажется разумным как исходное предположение. Создайте скрипт под названием driveHPfilt.m.
% script driveHPfilt.m
Fsamp = 2e3;
Fsm = 500;
Fbg = 800;
t = 0:1/Fsamp:0.1;
sgin = sin(2*pi*Fsm*t)+sin(2*pi*Fbg*t);
N = 5;
sgout = hipassfilt(N,sgin);Когда вы запускаете скрипт в командной строке, вы получаете сообщение об ошибке.
Сообщение об ошибке дает вам выбор открытия ассистента, чтобы откорректировать код MATLAB. Нажмите Click here получить Ассистент Создания фильтра на вашем экране.
Вы видите проблему: Вы не задавали ограничение частоты. Вы также забыли устанавливать частоту дискретизации. После экспериментирования вы находите, что можно задать Frequency units как Hz, Passband frequency, равный 700 Гц и Input Fs, равняются 2 000 Гц. Design method изменяется от Butterworth к Chebyshev type I. Вы нажимаете OK и получаете это на командной строке.
hp = designfilt('highpassiir','FilterOrder',N, ... 'PassbandFrequency',700,'PassbandRipple',1, ... 'SampleRate',2000);
Новый digitalFilter объект hp сохранено в рабочую область. В зависимости от ваших конструктивных ограничений можно изменить набор спецификации.
Можно установить designfilt никогда не предлагать Ассистенту Создания фильтра. Это действие устанавливает настройку MATLAB, которая может быть сброшена с setpref:
Используйте setpref('dontshowmeagain','filterDesignAssistant',false) предлагаться ассистент каждый раз. С этой командой можно получить ассистент, снова отключив его.
Используйте setpref('dontshowmeagain','filterDesignAssistant',true) отключить ассистент постоянно. Можно также нажать Do not show this message again в начальном диалоговом окне.
Можно установить designfilt всегда корректировать дефектные технические требования без выяснения. Это действие устанавливает настройку MATLAB, которая может быть сброшена при помощи setpref:
Используйте setpref('dontshowmeagain','filterDesignAssistantCodeCorrection',false) иметь designfilt откорректируйте свой код MATLAB, не прося подтверждение. Можно также нажать Always accept в поле диалогового окна подтверждения.
Используйте setpref('dontshowmeagain','filterDesignAssistantCodeCorrection',true) гарантировать тот designfilt корректирует ваш код MATLAB только, когда вы подтверждаете, что хотите изменения. С этой командой можно отменить эффект наличия, нажал Always accept в поле диалогового окна подтверждения.
digitalFilter | double | fftfilt | filt2block | filter | filtfilt | filtord | firtype | freqz | FVTool | grpdelay | impz | impzlength | info | isallpass | isdouble | isfir | islinphase | ismaxphase | isminphase | issingle | isstable | phasedelay | phasez | single | ss | stepz | tf | zerophase | zpk | zplane