Pref=20; % 20microPascal
Y = fft(y)
db = 20*log10((abs(Y)/Pref))

3楼: Originally posted by 沉默18 at 2018-05-23 08:34:24
非常感谢您！Y=fft（y），Y不是是声压才可以用您下面那个分贝计算公式？您知道怎么才能得到声压吗？
...

AWeighting = weightingFilter('A-weighting', 48000)

[wav,fs_Hz]=audioread( "audiofile.wav"  %load the WAV file
my_cal_factor = 1.0;  %the value for your system to convert the WAV into Pascals
wav_Pa = wav * my_cal_factor;

%extract the envelope
smooth_sec = 0.125;  %"FAST" SPL is 1/8th of second.  "SLOW" is 1 second;
smooth_Hz = 1/smooth_sec;
[b,a]=butter(1,smooth_Hz/(fs_Hz/2),'low');  %design a Low-pass filter
wav_env_Pa = sqrt(filter(b,a,wav_Pa.^2));  %rectify, by squaring, and low-pass filter

%compute SPL
Pa_ref = 20e-6;  %reference pressure for SPL in Air
SPL_dB = 10.0*log10( (wav_env_Pa ./ Pa_ref).^2 ); % 10*log10 because signal is squared

%plot results
figure;
subplot(2,1,1);
t_sec = ([1:size(wav_Pa)]-1)/fs_Hz;
plot(t_sec,wav_Pa);
xlabel('Time (sec)');
ylabel('Pressure (Pa)');

subplot(2,1,2)
plot(t_sec,SPL_dB);
xlabel('Time (sec)');
ylabel('SPL (dB)');
yl=ylim;ylim(yl(2)+[-80 0]);

6楼: Originally posted by 沉默18 at 2018-05-24 17:44:43
请问，怎么能得到横轴为频率，纵轴为db呢？
...