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]);

Отфильтруйте аудио Используя параметрический эквалайзер 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)

Отфильтруйте аудио Используя параметрический эквалайзер 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 должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: 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.

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

designParamEQ

Синтаксис

Описание

Примеры

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Ссылки

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

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

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

Темы

Введенный в R2016a

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

Поддержка

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

designParamEQ

Синтаксис

Описание

Примеры

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Ссылки

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

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

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

Темы

Введенный в R2016a

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

Поддержка

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

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

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

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

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

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

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

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

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