Проектируйте параметрический эквалайзер
[ проектирует B,A] =
designParamEQ(N,gain,centerFreq,bandwidth)Nпараметрический эквалайзер I порядка с заданным коэффициентом усиления, центральной частотой и шириной полосы пропускания. B и A являются матрицами коэффициентов числителя и знаменателя с столбцами, соответствующими каскадным фильтрам секции второго порядка (SOS).
[ задает опции с использованием одного или нескольких B,A] =
designParamEQ(___,Name,Value)Name,Value аргументы в виде пар.
Задайте порядок фильтра, пиковый коэффициент усиления в дБ, нормализованные центральные частоты и нормированную полосу пропускания полос вашего параметрического эквалайзера.
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]);
Спроектируйте параметрический эквалайзер секций второго порядка (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) с помощью 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)
Создайте анализатор спектра, чтобы визуализировать исходный аудиосигнал, и аудиосигнал прошел через ваш параметрический эквалайзер.
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)
[1] Orfanidis, Sophocles J. «High-Order Digital Parametric Equalizer Design». Журнал Общества аудиотехники. Том 53, ноябрь 2005, стр. 1026-1046.
