designfilt

Спроектируйте цифровые фильтры

Описание

пример

d = designfilt(resp,Name,Value) проектирует digitalFilter объект, d, с ответом вводят resp. Задайте фильтр далее с помощью набора Name,Value пары. Позволенные наборы спецификации зависят от типа ответа, resp, и состойте из комбинаций следующего:

  • Ограничения частоты соответствуют частотам, на которых фильтр предоставляет желаемое поведение. Примеры включают 'PassbandFrequency' и 'CutoffFrequency'. (См. полный список в соответствии с Аргументами в виде пар имя-значение.) Необходимо всегда задавать ограничения частоты.

  • Ограничения величины описывают поведение фильтра в конкретных частотных диапазонах. Примеры включают 'PassbandRipple' и 'StopbandAttenuation'. (См. полный список в соответствии с Аргументами в виде пар имя-значение.) designfilt обеспечивает значения по умолчанию для ограничений величины, оставленных незаданными. В проектах произвольной величины необходимо всегда задавать векторы из желаемых амплитуд.

  • 'FilterOrder'. Некоторые методы разработки позволяют вам задать порядок. Другие производят проекты минимального порядка. Таким образом, они генерируют самые маленькие фильтры, которые удовлетворяют заданным ограничениям.

  • 'DesignMethod' алгоритм, используемый, чтобы спроектировать фильтр. Примеры включают метод наименьших квадратов с ограничениями ('cls') и работа с окнами Кайзера ('kaiserwin'). Для некоторых наборов спецификации существует несколько методов разработки, доступных, чтобы выбрать из. В других случаях можно использовать только один метод, чтобы выполнить желаемым техническим требованиям.

  • Проектные решения являются параметрами, характерными для данного метода разработки. Примеры включают 'Window' для 'window' метод и оптимизация 'Weights' для произвольной величины equiripple проекты. (См. полный список в соответствии с Аргументами в виде пар имя-значение.) designfilt обеспечивает значения по умолчанию для проектных решений, оставленных незаданными.

  • 'SampleRate' частота, на которой действует фильтр. designfilt имеет частоту дискретизации по умолчанию 2 Гц. Используя это значение эквивалентно работе с нормированными частотами.

Примечание

Если вы задаете неполный или противоречивый набор пар "имя-значение" в командной строке, designfilt предложения открыть Ассистент Создания фильтра. Ассистент помогает вам спроектировать фильтр и вставляет откорректированный код MATLAB® по командной строке.

Если вы вызываете designfilt из скрипта или функции с неправильным набором технических требований, designfilt выдает ошибку сообщение со ссылкой, чтобы открыть Ассистент Создания фильтра. Ассистент помогает вам спроектировать фильтр, комментирует дефектный код в функции или скрипт, и вставляет откорректированный код MATLAB по следующей строке.

  • Использование filter в форме dataOut = filter(d,dataIn) отфильтровать сигнал с digitalFilterD.

  • Используйте FVTool, чтобы визуализировать digitalFilterD.

  • Введите d.Coefficients получить коэффициенты digitalFilterD. Для БИХ-фильтров коэффициенты описываются как секции второго порядка.

  • Смотрите digitalFilter для списка функций фильтрации и анализа, доступных для использования с digitalFilter объекты.

designfilt(d) позволяет вам отредактировать существующий цифровой фильтр, d. Это открывает Ассистент Создания фильтра, заполненный с техническими требованиями фильтра, которые можно затем изменить. Это - единственный способ, которым можно отредактировать digitalFilter объект. Его свойства в противном случае только для чтения.

Примеры

свернуть все

Спроектируйте минимальный порядок КИХ-фильтр lowpass с нормированной частотой полосы пропускания 0.25π rad/s, частота полосы задерживания 0.35π rad/s, неравномерность в полосе пропускания 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.25π rad/s, частота полосы пропускания 0.35π rad/s, неравномерность в полосе пропускания 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. Задайте нормированную ширину перехода 0.25π рад/с. Отобразите в линейных модулях ответ величины фильтра. Используйте его, чтобы отфильтровать вектор с 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);

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

свернуть все

Отфильтруйте ответ и тип в виде вектора символов или строкового скаляра. Кликните по одному из возможных значений resp расширять таблицу позволенных наборов спецификации.

 'lowpassfir'

 'lowpassiir'

 'highpassfir'

 'highpassiir'

 'bandpassfir'

 'bandpassiir'

 'bandstopfir'

 'bandstopiir'

 'differentiatorfir'

 'hilbertfir'

 'arbmagfir'

Типы данных: char | string

Цифровой фильтр в виде digitalFilter объект сгенерирован designfilt. Используйте этот вход, чтобы изменить технические требования существующего digitalFilter.

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: 'FilterOrder',20,'CutoffFrequency',0.4 достаточен, чтобы задать КИХ-фильтр lowpass.

Не все комбинации Name,Value пары допустимы. Допустимые комбинации зависят от ответа фильтра, в котором вы нуждаетесь и от ограничений частоты и величины вашего проекта.

Порядок фильтра

свернуть все

Порядок фильтра в виде разделенной запятой пары, состоящей из 'FilterOrder' и положительный целочисленный скаляр.

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

Порядок числителя БИХ проектирует в виде разделенной запятой пары, состоящей из 'NumeratorOrder' и положительный целочисленный скаляр.

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

Порядок знаменателя БИХ проектирует в виде разделенной запятой пары, состоящей из 'DenominatorOrder' и положительный целочисленный скаляр.

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

Ограничения частоты

свернуть все

Частота полосы пропускания в виде разделенной запятой пары, состоящей из 'PassbandFrequency' и положительная скалярная величина. Значение частоты должно быть в области значений Найквиста.

'PassbandFrequency1' более низкая частота полосы пропускания для проекта bandstop или полосы пропускания.

'PassbandFrequency2' более высокая частота полосы пропускания для проекта bandstop или полосы пропускания.

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

Частота полосы задерживания в виде разделенной запятой пары, состоящей из 'StopbandFrequency' и положительная скалярная величина. Значение частоты должно быть в области значений Найквиста.

'StopbandFrequency1' более низкая частота полосы задерживания для проекта bandstop или полосы пропускания

'StopbandFrequency2' более высокая частота полосы задерживания для проекта bandstop или полосы пропускания.

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

Частота на 6 дБ в виде разделенной запятой пары, состоящей из 'CutoffFrequency' и положительная скалярная величина. Значение частоты должно быть в области значений Найквиста.

'CutoffFrequency1' более низкая частота на 6 дБ для проекта bandstop или полосы пропускания.

'CutoffFrequency2' более высокая частота на 6 дБ для проекта bandstop или полосы пропускания.

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

Частота на 3 дБ в виде разделенной запятой пары, состоящей из 'HalfPowerFrequency' и положительная скалярная величина. Значение частоты должно быть в области значений Найквиста.

'HalfPowerFrequency1' более низкая частота на 3 дБ для проекта bandstop или полосы пропускания.

'HalfPowerFrequency2' более высокая частота на 3 дБ для проекта bandstop или полосы пропускания.

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

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

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

Частоты ответа в виде разделенной запятой пары, состоящей из 'Frequencies' и вектор. Используйте эту переменную, чтобы перечислить частоты, на которых фильтр произвольного ответа величины желал амплитуд. Частоты должны монотонно увеличиваться и лежать в области значений Найквиста. Первый элемент вектора должен быть или 0 или –fs/2, где fs является частотой дискретизации, и ее последним элементом должен быть fs/2. Если вы не задаете частоту дискретизации, designfilt использует значение по умолчанию 2 Гц.

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

Количество полос в многополосном проекте в виде разделенной запятой пары, состоящей из 'NumBands' и положительный целочисленный скаляр, не больше, чем 10.

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

Многополосные частоты ответа в виде разделенных запятой пар, состоящих из 'BandFrequenciesi' и числовой вектор. 'BandFrequenciesi', куда i запускается от 1 до NumBands, вектор, содержащий частоты, на которых i th полоса многополосного проекта имеет требуемые значения, 'BandAmplitudesi'. NumBands может быть самое большее 10. Частоты должны лечь в области значений Найквиста и должны быть заданы в монотонно увеличивающемся порядке. Смежные диапазоны частот должны иметь ту же амплитуду на своем перекрестке.

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

Ограничения величины

свернуть все

Неравномерность в полосе пропускания в виде разделенной запятой пары, состоящей из 'PassbandRipple' и положительная скалярная величина описывается в децибелах.

'PassbandRipple1' неравномерность в полосе пропускания нижней полосы для проекта bandstop.

'PassbandRipple2' неравномерность в полосе пропускания более высокой полосы для проекта bandstop.

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

Затухание в полосе задерживания в виде разделенной запятой пары, состоящей из 'StopbandAttenuation' и положительная скалярная величина описывается в децибелах.

'StopbandAttenuation1' затухание в полосе задерживания нижней полосы для полосового проекта.

'StopbandAttenuation2' затухание в полосе задерживания более высокой полосы для полосового проекта.

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

Желаемые амплитуды ответа произвольного ответа величины фильтруют в виде разделенной запятой пары, состоящей из 'Amplitudes' и вектор. Опишите амплитуды в линейных модулях. Вектор должен иметь ту же длину как 'Frequencies'.

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

Многополосные амплитуды ответа в виде разделенных запятой пар, состоящих из 'BandAmplitudesi' и числовой вектор. 'BandAmplitudesi', куда i запускается от 1 до NumBands, вектор, содержащий желаемые амплитуды в i th полоса многополосного проекта. NumBands может быть самое большее 10. Опишите амплитуды в линейных модулях. 'BandAmplitudesi' должен иметь ту же длину как 'BandFrequenciesi'. Смежные диапазоны частот должны иметь ту же амплитуду на своем перекрестке.

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

Метод разработки

свернуть все

Метод разработки в виде разделенной запятой пары, состоящей из 'DesignMethod' и вектор символов или строковый скаляр. Выбор метода разработки зависит на съемочной площадке ограничений частоты и величины, которые вы задаете.

  • 'butter' проектирует БИХ-фильтр Баттерворта. Фильтры Баттерворта имеют сглаженную монотонную частотную характеристику, которая является максимально плоской в полосе пропускания. Они жертвуют крутизной спада за плоскость.

  • 'cheby1' проектирует Чебышевский БИХ-фильтр типа I. Чебышевский тип , который я фильтрую, имеет частотную характеристику, которая является equiripple в полосе пропускания и максимально плоский в полосе задерживания. Их неравномерность в полосе пропускания увеличивается с увеличивающейся крутизной спада.

  • 'cheby2' проектирует Чебышевский БИХ-фильтр типа II. Чебышевские фильтры типа II имеют частотную характеристику, которая является максимально плоской в полосе пропускания и equiripple в полосе задерживания.

  • 'cls' проектирует КИХ-фильтр с помощью метода наименьших квадратов с ограничениями. Метод минимизирует несоответствие между заданной произвольной кусочной линейной функцией и ответом величины фильтра. В то же время это позволяет вам установить ограничения на неравномерность в полосе пропускания и затухание в полосе задерживания.

  • 'ellip' проектирует эллиптический БИХ-фильтр. Эллиптические фильтры имеют частотную характеристику, которая является equiripple и в полосе пропускания и в полосе задерживания.

  • 'equiripple' проектирует equiripple КИХ-фильтр с помощью алгоритма Парков-McClellan. Фильтры Equiripple имеют частотную характеристику, которая минимизирует максимальную величину пульсации по всем полосам.

  • 'freqsamp' проектирует КИХ-фильтр произвольного ответа величины путем выборки частотной характеристики однородно и взятия обратного преобразования Фурье.

  • 'kaiserwin' проектирует КИХ-фильтр с помощью метода окна Кайзера. Метод обрезает импульсную характеристику идеального фильтра и использует окно Кайзера, чтобы ослабить получившиеся колебания усечения.

  • 'ls' проектирует КИХ-фильтр с помощью наименьших квадратов. Метод минимизирует несоответствие между заданной произвольной кусочной линейной функцией и ответом величины фильтра.

  • 'maxflat' проектирует максимально плоский КИХ-фильтр. Эти фильтры имеют сглаженную монотонную частотную характеристику, которая является максимально плоской в полосе пропускания.

  • 'window' использует приближение наименьших квадратов, чтобы вычислить коэффициенты фильтра и затем сглаживает импульсную характеристику с 'Window'.

Типы данных: char | string

Опции метода разработки

свернуть все

Окно в виде разделенной запятой пары, состоящей из 'Window' и вектор из длины N + 1, где N является порядком фильтра. 'Window' может также быть соединен с именем окна или указателем на функцию, который указывает, что функция раньше генерировала окно. Любая такая функция должна взять N + 1, как сначала введено. Дополнительные входные параметры могут быть переданы путем определения массива ячеек. По умолчанию, 'Window' пустой вектор для 'freqsamp' метод разработки и @hamming для 'window' метод разработки.

Для списка доступных окон смотрите Windows.

Пример: 'Window',hann(N+1) и 'Window',(1-cos(2*pi*(0:N)'/N))/2 оба задают окно Hann, чтобы использовать с фильтром порядка N.

Пример: 'Window','hamming' задает Окно Хэмминга необходимого порядка.

Пример: 'Window',@mywindow позволяет вам задать свою собственную функцию окна.

Пример: 'Window',{@kaiser,0.5} задает окно Кайзера необходимого порядка параметром формы 0.5.

Типы данных: double | char | string | function_handle | cell

Полоса, чтобы соответствовать точно в виде разделенной запятой пары, состоящей из 'MatchExactly' и любой 'stopband', 'passband', или 'both'. 'both' доступно только для эллиптического метода разработки, где это - значение по умолчанию. 'stopband' значение по умолчанию для 'butter' и 'cheby2' методы. 'passband' значение по умолчанию для 'cheby1'.

Типы данных: char | string

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

Пример: 'PassbandOffset',0 результаты в фильтре с модулем получают в полосе пропускания.

Пример: 'PassbandOffset',2 результаты в фильтре с усилением полосы пропускания 2 дБ или 1.259.

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

Масштабируйте полосу пропускания в виде разделенной запятой пары, состоящей из 'ScalePassband' и логический скаляр. Когда вы устанавливаете 'ScalePassband' к true, полоса пропускания масштабируется после работы с окнами, так, чтобы фильтр имел модульное усиление на нулевой частоте.

Пример: 'Window',{@kaiser,0.1},'ScalePassband',true справка задает фильтр, ответ величины которого на нулевой частоте - точно 0 дБ. Дело обстоит не так, когда вы задаете 'ScalePassband',false. Чтобы проверить, визуализируйте фильтр с fvtool и увеличение.

Типы данных: логический

Нулевая фаза в виде разделенной запятой пары, состоящей из 'ZeroPhase' и логический скаляр. Когда вы устанавливаете 'ZeroPhase' к true, нулевой фазовый отклик получившегося фильтра всегда положителен. Это позволяет вам выполнить спектральную факторизацию на результате и получить фильтр минимальной фазы из него.

Типы данных: логический

Вес оптимизации полосы пропускания в виде разделенной запятой пары, состоящей из 'PassbandWeight' и положительная скалярная величина.

'PassbandWeight1' вес оптимизации полосы пропускания нижней полосы для bandstop КИХ-проекта.

'PassbandWeight2' вес оптимизации полосы пропускания более высокой полосы для bandstop КИХ-проекта.

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

Вес оптимизации полосы задерживания в виде разделенной запятой пары, состоящей из 'StopbandWeight' и положительная скалярная величина.

'StopbandWeight1' вес оптимизации полосы задерживания нижней полосы для полосового КИХ-проекта.

'StopbandWeight2' вес оптимизации полосы задерживания более высокой полосы для полосового КИХ-проекта.

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

Веса оптимизации в виде разделенной запятой пары, состоящей из 'Weights' и положительная скалярная величина или вектор из той же длины как 'Amplitudes'.

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

Многополосные веса в виде разделенных запятой пар, состоящих из 'BandWeightsi' и набор положительных скалярных величин или векторов. 'BandWeightsi', куда i запускается от 1 до NumBands, скаляр или вектор, содержащий веса оптимизации i th полоса многополосного проекта. Если задано как вектор, 'BandWeightsi' должен иметь ту же длину как 'BandAmplitudesi'.

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

Частота дискретизации

свернуть все

Частота дискретизации в виде разделенной запятой пары, состоящей из 'SampleRate' и положительная скалярная величина описывается в герц. Чтобы работать с нормированными частотами, установите 'SampleRate' к 2, или просто не используют его.

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

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

свернуть все

Цифровой фильтр, возвращенный как digitalFilter объект.

Больше о

свернуть все

Ассистент создания фильтра

Если вы задаете неполный или противоречивый набор расчетных параметров, designfilt предложения открыть Ассистент Создания фильтра.

(В описании аргумента для resp существует полный список допустимых наборов спецификации для всех доступных типов ответа.)

Ассистент ведет себя по-другому, если вы вызываете designfilt в командной строке или в рамках скрипта или функции.

Ассистент создания фильтра в командной строке

Вам дают сигнал, произведенный на уровне 2 кГц. Вас просят спроектировать КИХ-фильтр lowpass, который подавляет частотные составляющие выше, чем 650 Гц. “Частота среза” походит на хорошего кандидата на параметр спецификации. В командной строке MATLAB вы вводите следующее.

Fsamp = 2e3;
Fctff = 650;
dee = designfilt('lowpassfir','CutoffFrequency',Fctff, ...
                 'SampleRate',Fsamp);

Что-то, кажется, неправильно, потому что это диалоговое окно появляется на вашем экране.

Вы нажимаете Yes и получаете новое диалоговое окно, которое предлагает генерировать код. Вы видите, что переменные, которые вы задали прежде, были вставлены, где ожидается.

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

dee = designfilt('lowpassfir', 'FilterOrder', 10, ...
                 'CutoffFrequency', Fctff, 'SampleRate', Fsamp);

Введение имени фильтра повторяет информацию от диалогового окна.

dee
dee = 
 digitalFilter with properties:

    Coefficients: [1x11 double]
   Specifications:
    FrequencyResponse: 'lowpass'
      ImpulseResponse: 'fir'
           SampleRate: 2000
          FilterOrder: 10
      CutoffFrequency: 650
         DesignMethod: 'window'
 Use fvtool to visualize filter
 Use filter function 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', Fctff, ...
                 'PassbandRipple', 1, 'StopbandAttenuation', 60, ...
                 'SampleRate', Fsamp);
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(N,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);
Order = 5;
sgout = hipassfilt(Order,sgin);

Когда вы запускаете скрипт в командной строке, вы получаете сообщение об ошибке.

Сообщение об ошибке дает вам выбор открытия ассистента, чтобы откорректировать код MATLAB. Нажмите Click here получить Ассистент Создания фильтра на вашем экране.

Вы видите проблему: Вы не задавали ограничение частоты. Вы также забыли устанавливать частоту дискретизации. После экспериментирования вы находите, что можно задать Frequency units как Hz, Passband frequency, равный 700 Гц и Input Fs, равняются 2 000 Гц. Design method изменяется от Butterworth к Chebyshev type I. Вы нажимаете OK и получаете следующее.

Ассистент правильно идентифицировал файл, где вы вызываете designfilt. Нажмите Yes, чтобы принять изменение. Функция имеет откорректированный код MATLAB.

function dataOut = hipassfilt(N,dataIn)
% hpFilter = designfilt('highpassiir','FilterOrder',N);
hpFilter = designfilt('highpassiir', 'FilterOrder', N, ...
                'PassbandFrequency', 700, 'PassbandRipple', 1, ...
                'SampleRate', 2000);
dataOut = filter(hpFilter,dataIn);
end

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

Ассистент создания фильтра настройки

Можно установить 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 в поле диалогового окна подтверждения.

Поиск и устранение проблем

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

  • Вам не предлагают ассистент, если вы используете оценку секции кода, или от Панели инструментов MATLAB или путем нажатия Ctrl+Enter. (См., Делят Ваш Файл на Секции кода для получения дополнительной информации.)

  • Вам не предлагают ассистент, если ваш код имеет множественные вызовы designfilt, по крайней мере один из тех вызовов является неправильным, и

    • Вы вставляете код по командной строке и выполняете его путем нажатия Enter.

    • Вы выбираете код в Редакторе и выполняете его путем нажатия F9.

  • Вам не предлагают ассистент, если вы запускаете designfilt использование анонимной функции. (См. Анонимные функции для получения дополнительной информации.), Например, этот вход предлагает ассистенту.

    d = designfilt('lowpassfir','CutoffFrequency',0.6)
    Этот вход не делает.
    myFilterDesigner = @designfilt;
    d = myFilterDesigner('lowpassfir','CutoffFrequency',0.6)

  • Вам не предлагают ассистент, если вы запускаете designfilt использование eval. Например, этот вход предлагает ассистенту.

    d = designfilt('lowpassfir','CutoffFrequency',0.6)
    Этот вход не делает.
    myFilterDesigner = ...
        sprintf('designfilt(''%s'',''CutoffFrequency'',%f)', ...
                            'lowpassfir',0.6);
    d = eval(myFilterDesigner)

Ассистент Создания фильтра требует, чтобы программное обеспечение Java® и рабочий стол MATLAB запустились. Это не поддерживается, если вы запускаете MATLAB с -nojvm, -nodisplay, или -nodesktop опции.

Введенный в R2014a
Для просмотра документации необходимо авторизоваться на сайте