designParamEQ
Разработайте параметрический эквалайзер
Синтаксис
[B,A] =
designParamEQ(N,gain,centerFreq,bandwidth)[B,A] =
designParamEQ(N,gain,centerFreq,bandwidth,mode)Описание
[ разрабатывает параметрический эквалайзер th-порядка B,A] =
designParamEQ(N,gain,centerFreq,bandwidth)N с заданным усилением, центральной частотой и пропускной способностью. B и A являются матрицами числителя и коэффициентов знаменателя со столбцами, соответствующими каскадным фильтрам раздела второго порядка (SOS).
Примеры
Разработайте параметрический эквалайзер 2D полосы
Укажите, что порядок фильтра, пиковое усиление в дБ, нормировал центральные частоты и нормировал пропускную способность полос вашего параметрического эквалайзера.
N = [2,4]; gain = [6,-4]; centerFreq = [0.25,0.75]; bandwidth = [0.12,0.10];
Сгенерируйте коэффициенты фильтра с помощью заданных параметров.
[B,A] = designParamEQ(N,gain,centerFreq,bandwidth);
Создайте матрицу фильтра, совместимую с fvtool.
SOS = [B',[ones(sum(N)/2,1),A']];
Визуализируйте свой проект фильтра.
fvtool(SOS)
Отфильтруйте аудио Используя параметрический эквалайзер SOS
Разработайте параметрический эквалайзер разделов второго порядка (SOS) с помощью designParamEQ и отфильтруйте аудиопоток.
Создайте считывающее устройство звукового файла и Системные объекты средства записи аудио устройства. Используйте частоту дискретизации читателя как частота дискретизации средства записи. Вызовите 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 = fileReader(); deviceWriter(audio); count = count+1; end reset(fileReader);
Разработайте параметрический эквалайзер SOS.
N = [4,4]; gain = [-25,35]; centerFreq = [0.01,0.5]; bandwidth = [0.35,0.5]; [B,A] = designParamEQ(N,gain,centerFreq,bandwidth);
Визуализируйте свой проект фильтра.
SOS = [B',[ones(4,1),A']]; fvtool(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 Signals',... '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(scope) release(deviceWriter) release(fileReader)
Отфильтруйте аудио Используя параметрический эквалайзер FOS
Разработайте параметрический эквалайзер разделов четвертого порядка (FOS) с помощью designParamEQ и отфильтруйте аудиопоток.
Создайте считывающее устройство звукового файла и Системные объекты средства записи аудио устройства. Используйте частоту дискретизации читателя как частота дискретизации средства записи. Вызовите 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 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);Создайте БИХ-фильтры FOS.
section1 = dsp.IIRFilter('Numerator',B(:,1)','Denominator',[1,A(:,1)']); section2 = dsp.IIRFilter('Numerator',B(:,2)','Denominator',[1,A(:,2)']);
Визуализируйте частотную характеристику своего параметрического эквалайзера.
[H1,w] = freqz(section1,8192,sampleRate); H2 = freqz(section2,8192,sampleRate); H = 20.*log10(abs(H1.*H2)); semilogx(w,H); title('Magnitude Response (dB)') xlabel('Frequency (Hz)') ylabel('Magnitude (dB)') grid on
Создайте спектр, анализатор, чтобы визуализировать исходный звуковой сигнал и звуковой сигнал прошел через ваш параметрический эквалайзер.
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'});
Проигрывайте отфильтрованный звуковой сигнал и визуализируйте исходные и отфильтрованные спектры.
setup(scope,ones(frameSize,2)); count = 0; while count < 2500 x = fileReader(); y = section1(x); z = section2(y); scope([x(:,1),z(:,1)]); deviceWriter(z); count = count + 1; end release(fileReader) release(deviceWriter) release(scope)