how to plot bar overall 1/3 octave by Poctave() as show below(bar color red)?

97 views (last 30 days)

Accepted Answer

Mathieu NOE
Mathieu NOE on 26 Feb 2021
hello
see example below :
Fs = 44.1E3;
xInp = rand(1,2^16);
% Compute the octave spectrum of the same signal, but this
% time, use a A-Weighting filter. Visualize the spectrum.
[P,f] = poctave(xInp,Fs,'FrequencyLimits',[20 20E3],'FilterOrder',6,'BandsPerOctave',3,'Weighting','A');
f = round(f);
figure(1), bar(1:length(f),20*log10(P)+200)
xlabel('Hz');
ylabel('Amplitude (dB(A))');
set(gca,'XTick',1:length(f));
set(gca,'XTickLabel',f);set(xticklabel_rotate([],60),'color',get(gca,'xcolor'));
  10 Comments
ngoc quy hoang ngoc quy
ngoc quy hoang ngoc quy on 27 Feb 2021
Edited: ngoc quy hoang ngoc quy on 27 Feb 2021
ok
but i think the poctave function output is sound power not Sound Pressure so Sound Pressure Level (SPL) does not calculate instead it will calculate sound power level ( LwA) by LwA = 10*log10(P/P0) with P0 = 10e-12 w.
I want to hear your opinion on that matter??

Sign in to comment.

More Answers (0)

Categories

Find more on Measurements and Spatial Audio in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!