How to use matlab audioread Plot the spectrum with frequency on the x-axis and dBA /dBC on the y-axis？

Hi,everyone~

How to use matlab audioread to plot the spectrum with frequency on the x-axis and dBA/dBC sound pressure level (SPL) on the y-axis?

As shown in Figure 1.

The data that is currently available is a audio(wav).

I would be grateful if you have example code!

thanks

Attach a sample audio.(See the attachment)

fig.1

[ZHU z]

Could anyone help me? ddl will be here soon.

Answer 1

[Chunru]

doc audioread     % for reading audio file
doc pwelch        % for computing power spectral density

[ZHU z]

thanks! Do you have the code for the example shown in the figure1?

My difficulty is mainly in the dBA and dBC (A/C weighted sound pressure level)~

[Chunru]

doc weightingFilter     % from audio toolbox

% Below is the example in the document

% Create a dsp.AudioFileReader System object.
samplesPerFrame = 1024;
reader = dsp.AudioFileReader('Filename', ...
    'RockGuitar-16-44p1-stereo-72secs.wav', ...
    'SamplesPerFrame',samplesPerFrame, ...
    'PlayCount',Inf);

% Create a weightingFilter System object. Use the sample rate of the reader as the sample rate of the weighting filter.
Fs = reader.SampleRate;
weightFilt = weightingFilter('A-weighting',Fs);

% Create a spectrum analyzer to visualize the original audio signal and the audio signal after frequency-weighted filtering.
scope = dsp.SpectrumAnalyzer( ...
    'SampleRate',Fs, ...
    'PlotAsTwoSidedSpectrum',false, ...
    'FrequencyScale','Log', ...
    'FrequencyResolutionMethod','WindowLength', ...
    'WindowLength',samplesPerFrame, ...
    'Title','A-Weighted Filtering', ...
    'ShowLegend',true, ...
    'ChannelNames',{'Original signal','Filtered signal'});


x = reader();           % read data
y = weightFilt(x);      % filter the data
scope([x(:,1),y(:,1)])  % compute the spectrum after filtering

[ZHU z]

OMG~Thank you for your full answer, this is what I want~

Thank you!

[ZHU z]

hi Chunru~

Sorry to bother you again, the result after I run the code is shown in the figure 1.

The vertical axis is in dBm, is this the correct result? How do I convert to dBA? Thanks a lot~

fig 1