designParamEQ

Спроектируйте параметрический эквалайзер

Синтаксис

[B,A] =
designParamEQ(N,gain,centerFreq,bandwidth)

[B,A] =
designParamEQ(N,gain,centerFreq,bandwidth,mode)

[B,A] =
designParamEQ(___,Name,Value)

Описание

пример

[B,A] = designParamEQ(N,gain,centerFreq,bandwidth) проектирует Nпараметрический эквалайзер th-порядка с заданным усилением, сосредоточьте частоту и полосу пропускания. B и A матрицы числителя и коэффициентов знаменателя, со столбцами, соответствующими каскадным фильтрам секции второго порядка (SOS).

пример

[B,A] = designParamEQ(N,gain,centerFreq,bandwidth,mode) задает, реализован ли параметрический эквалайзер с секциями второго порядка или разделами четвертого порядка (FOS).

[B,A] = designParamEQ(___,Name,Value) задает опции с помощью одного или нескольких Name,Value парные аргументы.

Примеры

свернуть все

Спроектируйте параметрический эквалайзер 2D полосы

Скрипт Open Live Script

Укажите, что порядок фильтра, пиковое усиление в дБ, нормировал центральные частоты и нормировал полосу пропускания полос вашего параметрического эквалайзера.

N = [2, ...
     4];
 
получите = [6, ...
      -4]; 
  
centerFreq = [0.25, ...
              0.75]; 
          
полоса пропускания = [0.12, ...
              0.1];

Сгенерируйте коэффициенты фильтра с помощью заданных параметров.

[B,A] = designParamEQ(N,gain,centerFreq,bandwidth,"Orientation","row");

Визуализируйте свое создание фильтра.

fvtool([B,A]);

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

Отфильтруйте аудио Используя параметрический эквалайзер SOS

Скрипт Open Live Script

Спроектируйте параметрический эквалайзер секций второго порядка (SOS) с помощью designParamEQ и отфильтруйте аудиопоток.

Создайте считывающее устройство звукового файла и Системные объекты средства записи аудио устройства. Используйте частоту дискретизации читателя как частота дискретизации средства записи.

frameSize = 256;

fileReader = dsp.AudioFileReader("RockGuitar-16-44p1-stereo-72secs.wav","SamplesPerFrame",frameSize);

sampleRate = fileReader.SampleRate;

deviceWriter = audioDeviceWriter("SampleRate",sampleRate);

Проигрывайте звуковой сигнал через свое устройство.

count = 0;
while count < 2500
    audio = fileReader();
    deviceWriter(audio);
    count = count + 1;
end
reset(fileReader)

Спроектируйте параметрический эквалайзер SOS, подходящий для использования с dsp.BiquadFilter.

N = [4,4];
gain = [-25,35];
centerFreq = [0.01,0.5];
bandwidth = [0.35,0.5];
[B,A] = designParamEQ(N,gain,centerFreq,bandwidth);

Визуализируйте свое создание фильтра. Вызовите designParamEQ с теми же техническими требованиями проекта. Задайте выходную ориентацию как "row" так, чтобы это подошло для использования с fvtool.

[Bvisualize,Avisualize] = designParamEQ(N,gain,centerFreq,bandwidth,"Orientation","row");
fvtool([Bvisualize,Avisualize], ...
    "Fs",fileReader.SampleRate, ...
    "FrequencyScale","Log");

Создайте фильтр biquad.

myFilter = dsp.BiquadFilter( ...
    "SOSMatrixSource","Input port", ...
    "ScaleValuesInputPort",false);

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

scope = dsp.SpectrumAnalyzer( ...
    "SampleRate",sampleRate, ...
    "PlotAsTwoSidedSpectrum",false, ...
    "FrequencyScale","Log", ...
    "FrequencyResolutionMethod","WindowLength", ...
    "WindowLength",frameSize, ...
    "Title","Original and Equalized Signals", ...
    "ShowLegend",true, ...
    "ChannelNames",{'Original Signal','Equalized Signal'});

Проигрывайте отфильтрованный звуковой сигнал и визуализируйте исходные и отфильтрованные спектры.

count = 0;
while count < 2500
    originalSignal = fileReader();
    equalizedSignal = myFilter(originalSignal,B,A);
    scope([originalSignal(:,1),equalizedSignal(:,1)]);
    deviceWriter(equalizedSignal);
    count = count + 1;
end

Как лучшая практика, выпустите свои объекты, однажды сделанные.

release(deviceWriter)
release(fileReader)
release(scope)

Figure Spectrum Analyzer contains an axes object and other objects of type uiflowcontainer, uimenu, uitoolbar. The axes object with title Original and Equalized Signals contains 2 objects of type line. These objects represent Original Signal, Equalized Signal.

Отфильтруйте аудио Используя параметрический эквалайзер FOS

Скрипт Open Live Script

Спроектируйте параметрический эквалайзер разделов четвертого порядка (FOS) с помощью designParamEQ и отфильтруйте аудиопоток.

frameSize = 256;

fileReader = dsp.AudioFileReader( ...
    "RockGuitar-16-44p1-stereo-72secs.wav", ...
    "SamplesPerFrame",frameSize);

sampleRate = fileReader.SampleRate;

deviceWriter = audioDeviceWriter( ...
    "SampleRate",sampleRate);

Проигрывайте звуковой сигнал через свое устройство.

count = 0;
while count < 2500
    x = fileReader();
    deviceWriter(x);
    count = count + 1;
end
reset(fileReader)

Спроектируйте коэффициенты параметрического эквалайзера FOS.

N = [2,4];
gain = [5,10];
centerFreq = [0.025,0.65];
bandwidth = [0.025,0.35];
mode = "fos";

[B,A] = designParamEQ(N,gain,centerFreq,bandwidth,mode,"Orientation","row");

Создайте БИХ-фильтры FOS.

myFilter = dsp.FourthOrderSectionFilter(B,A);

Визуализируйте частотную характеристику своего параметрического эквалайзера.

fvtool(myFilter)

scope = dsp.SpectrumAnalyzer( ...
    "SampleRate",sampleRate, ...
    "PlotAsTwoSidedSpectrum",false, ...
    "FrequencyScale","Log", ...
    "FrequencyResolutionMethod","WindowLength", ...
    "WindowLength",frameSize, ...
    "Title","Original and Equalized Signals", ...
    "ShowLegend",true, ...
    "ChannelNames",{'Original Signal','Equalized Signal'});

Проигрывайте отфильтрованный звуковой сигнал и визуализируйте исходные и отфильтрованные спектры.

count = 0;
while count < 2500
    x = fileReader();
    y = myFilter(x);
    
    scope([x(:,1),y(:,1)]);
    
    deviceWriter(y);
    
    count = count + 1;
end

Как лучшая практика, выпустите свои объекты, однажды сделанные.

release(fileReader)
release(deviceWriter)
release(scope)

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

свернуть все

`N` — Порядок фильтра
скаляр | вектор-строка

Порядок фильтра в виде скалярного или вектора-строки та же длина как centerFreq. Элементами вектора должны быть даже целые числа.

`gain` — Пиковое усиление (дБ)
скаляр | вектор-строка

Пиковое усиление в дБ в виде скалярного или вектора-строки та же длина как centerFreq. Элементы вектора должны быть с действительным знаком.

`centerFreq` — Нормированная центральная частота полос эквалайзера
скаляр | вектор-строка

Нормированная центральная частота полос эквалайзера в виде скалярного или вектора-строки из вещественных значений в области значений от 0 до 1, где 1 соответствует частоте Найквиста (π рад/отсчет). Если centerFreq задан как вектор-строка, отдельные эквалайзеры спроектированы для каждого элемента centerFreq.

`bandwidth` — Нормированная полоса пропускания
скаляр | вектор-строка

Нормированная полоса пропускания в виде скалярного или вектора-строки та же длина как centerFreq. Элементы вектора указаны как вещественные значения в области значений от 0 до 1, где 1 соответствует частоте Найквиста (π рад/отсчет).

Нормированная полоса пропускания измеряется на уровне дБ усиления/2. Если усиление установлено в -Inf (отметьте фильтр), нормированная полоса пропускания измеряется в точке затухания на 3 дБ: $10 \times \log_{10} (0.5)$ .

Чтобы преобразовать полосу пропускания октавы в нормированную полосу пропускания, вычислите связанный Q - фактор как

$Q = \frac{\sqrt{2^{(o c t a v e b a n d w i d t h)}}}{2^{(o c t a v e b a n d w i d t h)} - 1} .$

Затем преобразуйте в полосу пропускания

$b a n d w i d t h = \frac{c e n t e r F r e q}{Q} .$

`mode` — Режим проектирования
`'sos'` (значение по умолчанию) | `'fos'`

Режим проектирования в виде 'sos' или 'fos'.

'sos' – Реализует ваш эквалайзер как расположенные каскадом фильтры второго порядка.
'fos' – Реализует ваш эквалайзер как расположенные каскадом фильтры четвертого порядка. Поскольку разделы четвертого порядка не требуют расчета корней, они обычно более в вычислительном отношении эффективны.

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

Аргументы name-value

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

Пример: 'Orientation',"row"

`Orientation` — Ориентация возвращенных коэффициентов фильтра
`"column"` (значение по умолчанию) | `"row"`

Ориентация возвращенных коэффициентов фильтра в виде разделенной запятой пары, состоящей из 'Orientation' и "column" или "row":

Установите 'Orientation' к "row" для функциональной совместимости с FVTool, dsp.DynamicFilterVisualizer, и dsp.FourthOrderSectionFilter.
Установите 'Orientation' к "column" для функциональной совместимости с dsp.BiquadFilter.

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

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

свернуть все

`B` — Коэффициенты фильтра числителя
матрица

Коэффициенты фильтра числителя, возвращенные как матрица. Размер и интерпретация B зависит от Orientation и mode:

Если 'Orientation' установлен в "column" и mode установлен в "sos", затем B возвращен как L-by-3 матрица. Каждый столбец соответствует коэффициентам числителя ваших каскадных секций второго порядка.
Если 'Orientation' установлен в "column" и mode установлен в "fos", затем B возвращен как L-by-5 матрица. Каждый столбец соответствует коэффициентам числителя ваших каскадных разделов четвертого порядка.
Если 'Orientation' установлен в "row" и mode установлен в "sos", затем B возвращен как 3 L матрицей. Каждая строка соответствует коэффициентам числителя ваших каскадных секций второго порядка.
Если 'Orientation' установлен в "row" и mode установлен в "fos", затем B возвращен как 5 L матрицей. Каждая строка соответствует коэффициентам числителя ваших каскадных разделов четвертого порядка.

`A` — Коэффициенты фильтра знаменателя
матрица

Коэффициенты фильтра знаменателя, возвращенные как матрица. Размер и интерпретация A зависит от Orientation и mode:

Если 'Orientation' установлен в "column" и mode установлен в "sos", затем A возвращен как L-by-2 матрица. Каждый столбец соответствует коэффициентам знаменателя ваших каскадных секций второго порядка. A не включает ведущие коэффициенты единицы.
Если 'Orientation' установлен в "column" и mode установлен в "fos", затем A возвращен как L-by-4 матрица. Каждый столбец соответствует коэффициентам знаменателя ваших каскадных разделов четвертого порядка. A не включает ведущие коэффициенты единицы.
Если 'Orientation' установлен в "row" и mode установлен в "sos", затем A возвращен как 3 L матрицей. Каждая строка соответствует коэффициентам знаменателя ваших каскадных секций второго порядка.
Если 'Orientation' установлен в "row" и mode установлен в "fos", затем A возвращен как 5 L матрицей. Каждая строка соответствует коэффициентам знаменателя ваших каскадных разделов четвертого порядка.

Ссылки

[1] Orfanidis, Софокл Дж. "Старший Цифровой Проект Параметрического эквалайзера". Журнал Общества звукоинженеров. Издание 53, ноябрь 2005, стр 1026–1046.

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

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

Темы

Введенный в R2016a

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

designParamEQ

Синтаксис

Описание

Примеры

Спроектируйте параметрический эквалайзер 2D полосы

Отфильтруйте аудио Используя параметрический эквалайзер SOS

Отфильтруйте аудио Используя параметрический эквалайзер FOS

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

`N` — Порядок фильтра
скаляр | вектор-строка

`gain` — Пиковое усиление (дБ)
скаляр | вектор-строка

`centerFreq` — Нормированная центральная частота полос эквалайзера
скаляр | вектор-строка

`bandwidth` — Нормированная полоса пропускания
скаляр | вектор-строка

`mode` — Режим проектирования
`'sos'` (значение по умолчанию) | `'fos'`

Аргументы name-value

`Orientation` — Ориентация возвращенных коэффициентов фильтра
`"column"` (значение по умолчанию) | `"row"`

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

`B` — Коэффициенты фильтра числителя
матрица

`A` — Коэффициенты фильтра знаменателя
матрица

Ссылки

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

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

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

Темы

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

Поддержка

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

designParamEQ

Синтаксис

Описание

Примеры

Спроектируйте параметрический эквалайзер 2D полосы

Отфильтруйте аудио Используя параметрический эквалайзер SOS

Отфильтруйте аудио Используя параметрический эквалайзер FOS

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

N — Порядок фильтра скаляр | вектор-строка

gain — Пиковое усиление (дБ) скаляр | вектор-строка

centerFreq — Нормированная центральная частота полос эквалайзера скаляр | вектор-строка

bandwidth — Нормированная полоса пропускания скаляр | вектор-строка

mode — Режим проектирования 'sos' (значение по умолчанию) | 'fos'

Аргументы name-value

Orientation — Ориентация возвращенных коэффициентов фильтра "column" (значение по умолчанию) | "row"

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

B — Коэффициенты фильтра числителя матрица

A — Коэффициенты фильтра знаменателя матрица

Ссылки

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

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

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

Темы

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

Поддержка

`N` — Порядок фильтра
скаляр | вектор-строка

`gain` — Пиковое усиление (дБ)
скаляр | вектор-строка

`centerFreq` — Нормированная центральная частота полос эквалайзера
скаляр | вектор-строка

`bandwidth` — Нормированная полоса пропускания
скаляр | вектор-строка

`mode` — Режим проектирования
`'sos'` (значение по умолчанию) | `'fos'`

`Orientation` — Ориентация возвращенных коэффициентов фильтра
`"column"` (значение по умолчанию) | `"row"`

`B` — Коэффициенты фильтра числителя
матрица

`A` — Коэффициенты фильтра знаменателя
матрица

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