您正在关注此问题

您将在关注的内容源中看到更新。
您可能会收到电子邮件，具体取决于您的通信预设项。

how to calculate sound level (spL) audio file (*.wav) ?

67 次查看（过去 30 天）

显示更早的评论

ngoc quy hoang ngoc quy 2021-2-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav

⋮

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav

评论： ranavind 2023-1-31

calculate sound level

回答（1 个）

jibrahim 2021-3-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#answer_639151

⋮

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#answer_639151

在 MATLAB Online 中打开

Hi ngoc,

You can use splMeter from Audio Toolbox.

% EXAMPLE: Calculate sound pressure level of signal
audFileName = 'FunkyDrums-44p1-stereo-25secs.mp3';
scope = timescope('SampleRate',Fs,...
    'TimeSpanSource','property',...
    'YLimits',[20 110], 'AxesScaling','auto',...
    'ChannelNames',...
    {'LAF', 'LAeq', 'LApeak', 'LAmax'});
sm = splMeter('TimeWeighting', 'Fast',...
    'FrequencyWeighting', 'A-weighting',...
    'SampleRate', Fs, ...
    'TimeInterval',2,...
    'CalibrationFactor',1);
x = audioread(audFileName);
[LAF,LAeq,LApeak,LAmax] = sm(x(:,1));
scope([LAF,LAeq,LApeak,LAmax])

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_1368721

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_1368721

thank you, But problem is.

I want to calculate the sound pressure level third ooctave using poctave function

I used poctave function to convert pressure data (p) (time,pressure) to 1/3 octave spectrum

[P3rd, cf] = poctave(p,fs,'FrequencyLimits',[df,fs/2],'BandsPerOctave',3);

In order to calculate the SPL-1/3 Octave (sploct)

sploct = 10*log10(P3rd/pref^2);

Is the equation for calculating SPL-1/3 octave is correct or not?

can you help me?

Jimmy Lapierre 2021-3-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_1369146

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_1369146

One thing that appears to be missing is a calibration for your input (in addition to pref). The advantage of using splMeter is that it includes a method to calibrate your microphone level, so you don't have to do that by yourself.

ZHU z 2022-4-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2084114

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2084114

hi jibrahim~

thanks for your codes~

When I run it matlab says: "undefined function or variable 'Fs'.

Error Calculate_sound_pressure_level_of_signal (line 4)

scope = timescope('SampleRate',Fs,..."

How should I handle this issue? Can you tell me, thanks a lot~

jibrahim 2022-4-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2084769

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2084769

在 MATLAB Online 中打开

Sure. Use this call to audioread and move it up (before the definition of scope)

[x,Fs] = audioread(audFileName);

ZHU z 2022-4-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2085885

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2085885

在 MATLAB Online 中打开

Thanks jibrahim~

Sorry to bother you again,I moved the "[x,Fs] = audioread(audFileName);" up as you said, but still the same error,Is it because I am missing something?

% EXAMPLE: Calculate sound pressure level of signal
[x , Fs] = audioread('SAODIJI.wav');
scope = timescope('SampleRate',Fs,...
    'TimeSpanSource','property',...
    'YLimits',[20 110], 'AxesScaling','auto',...
    'ChannelNames',...
    {'LAF', 'LAeq', 'LApeak', 'LAmax'});
sm = splMeter('TimeWeighting', 'Fast',...
    'FrequencyWeighting', 'A-weighting',...
    'SampleRate', Fs, ...
    'TimeInterval',2,...
    'CalibrationFactor',1);
[LAF,LAeq,LApeak,LAmax] = sm(x(:,1));
scope([LAF,LAeq,LApeak,LAmax])

jibrahim 2022-4-6

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2086450

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2086450

That code works for me. what is the error?

ZHU z 2022-4-6

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2086985

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/757874-how-to-calculate-sound-level-spl-audio-file-wav#comment_2086985

在 MATLAB Online 中打开

thanks, jibrahim~

the error is ："Undefined function or variable 'timescope'."(line 3) scope = timescope('SampleRate',Fs,...

% EXAMPLE: Calculate sound pressure level of signal
[x , Fs] = audioread('SAODIJI.wav');
scope = timescope('SampleRate',Fs,...
    'TimeSpanSource','property',...
    'YLimits',[20 110], 'AxesScaling','auto',...
    'ChannelNames',...
    {'LAF', 'LAeq', 'LApeak', 'LAmax'});
sm = splMeter('TimeWeighting', 'Fast',...
    'FrequencyWeighting', 'A-weighting',...
    'SampleRate', Fs, ...
    'TimeInterval',2,...
    'CalibrationFactor',1);
[LAF,LAeq,LApeak,LAmax] = sm(x(:,1));
scope([LAF,LAeq,LApeak,LAmax])