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параметрический эквалайзер I порядка с заданным коэффициентом усиления, центральной частотой и шириной полосы пропускания. B и A являются матрицами коэффициентов числителя и знаменателя с столбцами, соответствующими каскадным фильтрам секции второго порядка (SOS).

пример

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

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

Примеры

свернуть все

Проектирование двухдиапазонного параметрического эквализатора

Открыть Live Script

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

N = [2, ...
     4];
 
коэффициент усиления = [6, ...
      -4]; 
  
centerFreq = [0.25, ...
              0.75]; 
          
bandwidth = [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 and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains an object of type line.

Фильтрация аудио с помощью параметрического эквалайзера SOS

Открыть 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");

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

Создайте фильтр 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 and other objects of type uiflowcontainer, uimenu, uitoolbar. The axes with title Original and Equalized Signals contains 2 objects of type line. These objects represent Original Signal, Equalized Signal.

Фильтрация аудио с помощью FOS Parametric Equalizer

Открыть 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");

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

myFilter = dsp.FourthOrderSectionFilter(B,A);

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

fvtool(myFilter)

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

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

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 -factor как

$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 должны находиться внутри кавычек. Можно задать несколько аргументов в виде пар имен и значений в любом порядке 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, Sophocles J. «High-Order Digital Parametric Equalizer Design». Журнал Общества аудиотехники. Том 53, ноябрь 2005, стр. 1026-1046.

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

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

См. также

designShelvingEQ | designVarSlopeFilter | dsp.BiquadFilter | multibandParametricEQ

Темы

Введенный в R2016a

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

designParamEQ

Синтаксис

Описание

Примеры

Проектирование двухдиапазонного параметрического эквализатора

Фильтрация аудио с помощью параметрического эквалайзера SOS

Фильтрация аудио с помощью FOS Parametric Equalizer

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

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

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

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

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

`mode` - Проект
`'sos'` (по умолчанию) | `'fos'`

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

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

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

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

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

Ссылки

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

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

См. также

Темы

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

Поддержка

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

designParamEQ

Синтаксис

Описание

Примеры

Проектирование двухдиапазонного параметрического эквализатора

Фильтрация аудио с помощью параметрического эквалайзера SOS

Фильтрация аудио с помощью FOS Parametric Equalizer

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

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

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

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

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

mode - Проект 'sos' (по умолчанию) | 'fos'

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

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

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

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

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

Ссылки

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

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

См. также

Темы

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

Поддержка

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

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

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

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

`mode` - Проект
`'sos'` (по умолчанию) | `'fos'`

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

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

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

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