Разработайте отлогий эквалайзер
[B,A] =
designShelvingEQ(gain,slope,Fc)[B,A] =
designShelvingEQ(gain,slope,Fc,type)Разработайте три БИХ-эквалайзера низкой полки второго порядка с помощью designShelvingEQ. Три отлогих эквалайзера используют три отдельных наклонных спецификации.
Задайте частоту дискретизации, пиковое усиление, наклонный коэффициент и нормированную частоту среза для трех отлогих эквалайзеров. Частота дискретизации находится в Гц. Пиковое усиление находится в дБ.
Fs = 44.1e3; gain = 5; slope1 = 0.5; slope2 = 0.75; slope3 = 1; Fc = 1000/(Fs/2);
Разработайте коэффициенты фильтра с помощью заданных параметров.
[B1,A1] = designShelvingEQ(gain,slope1,Fc); [B2,A2] = designShelvingEQ(gain,slope2,Fc); [B3,A3] = designShelvingEQ(gain,slope3,Fc);
Создайте матрицы фильтра, совместимые с fvtool.
SOS1 = [B1',[1,A1']]; SOS2 = [B2',[1,A2']]; SOS3 = [B3',[1,A3']];
Визуализируйте свой проект фильтра.
fvtool(... dsp.BiquadFilter('SOSMatrix',SOS1),... dsp.BiquadFilter('SOSMatrix',SOS2),... dsp.BiquadFilter('SOSMatrix',SOS3),... 'Fs',Fs,... 'FrequencyScale','Log'); legend('slope = 0.1',... 'slope = 0.5',... 'slope = 1');
Разработайте эквалайзер низкой полки, и затем используйте его, чтобы отфильтровать звуковой сигнал.
Создайте считывающее устройство звукового файла и объекты средства записи аудио устройства. Используйте частоту дискретизации читателя как частота дискретизации средства записи. Вызовите setup, чтобы уменьшать вычислительную загрузку инициализации в цикле аудиопотока.
frameSize = 256; fileReader = dsp.AudioFileReader(... 'RockGuitar-16-44p1-stereo-72secs.wav',... 'SamplesPerFrame',frameSize); sampleRate = fileReader.SampleRate; deviceWriter = audioDeviceWriter(... 'SampleRate',sampleRate); setup(fileReader); setup(deviceWriter,ones(frameSize,2));
Проигрывайте звуковой сигнал через свое устройство.
count = 0; while count < 2500 audio = step(fileReader); play(deviceWriter,audio); count = count+1; end reset(fileReader)
Разработайте эквалайзер низкой полки разделов второго порядка (SOS).
gain = 10; slope = 3; Fc = 0.025; [B,A] = designShelvingEQ(gain,slope,Fc);
Визуализируйте свой проект эквалайзера.
SOS = [B',[1,A']]; fvtool(dsp.BiquadFilter('SOSMatrix',SOS),... '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 Signal',... 'ShowLegend',true,... 'ChannelNames',{'Original Signal','Equalized Signal'});
Проигрывайте компенсируемый звуковой сигнал и визуализируйте исходные и компенсируемые спектры.
setup(scope,ones(frameSize,2)); count = 0; while count < 2500 originalSignal = fileReader(); equalizedSignal = myFilter(originalSignal,B,A); scope([originalSignal(:,1),equalizedSignal(:,1)]); deviceWriter(equalizedSignal); count = count+1; end release(fileReader) release(scope) release(deviceWriter)
Разработайте трех БИХ второго порядка высокие эквалайзеры полки с помощью designShelvingEQ. Три отлогих эквалайзера используют три отдельных спецификации усиления.
Задайте частоту дискретизации, пиковое усиление, наклонный коэффициент и нормированную частоту среза для трех отлогих эквалайзеров. Частота дискретизации находится в Гц. Пиковое усиление находится в дБ.
Fs = 44.1e3; gain1 = -6; gain2 = 6; gain3 = 12; slope = 0.8; Fc = 18000/(Fs/2);
Разработайте коэффициенты фильтра с помощью заданных параметров.
[B1,A1] = designShelvingEQ(gain1,slope,Fc,'hi'); [B2,A2] = designShelvingEQ(gain2,slope,Fc,'hi'); [B3,A3] = designShelvingEQ(gain3,slope,Fc,'hi');
Создайте матрицы фильтра, совместимые с fvtool.
SOS1 = [B1',[1,A1']]; SOS2 = [B2',[1,A2']]; SOS3 = [B3',[1,A3']];
Визуализируйте свой проект фильтра.
fvtool(dsp.BiquadFilter('SOSMatrix',SOS1),... dsp.BiquadFilter('SOSMatrix',SOS2),... dsp.BiquadFilter('SOSMatrix',SOS3),... 'Fs',Fs); legend('gain = -6 dB',... 'gain = 6 dB',... 'gain = 12 dB',... 'Location','NorthWest')
