请教通过三分之一倍频程计算声压级

wangwenting

我已经得到了声压信号
通过下面的程序进行三分之一倍频程滤波
function [p,ff] = oct3bank(noi);
pi = 3.14159265358979;
Fs = 500000;    %1/(t(2)-t(1));     % Sampling Frequency
N = 3;      % Order of analysis filters.
%f=[1.00 1.25 1.60 2.00 2.50 3.15 4.00 5.00 6.300 8.00]
ff=[20 25 31.5 40 50 63 80 100 125 160 200 250 315 400 500 630 800 1000 1250 1600 2000 2500 3150 4000 5000 6300 8000 10000 12500 16000]
%ff = ((2^(1/3)).^[0:19]);  % Exact center freq.  
p = zeros(1,30);
% Design filters and compute RMS powers in 1/3-oct. bands
for i = 1:30
   [B,A] = oct3dsgn(ff(i),Fs,N);
   y = filter(B,A,noi);              %<-----滤波器
   p(i) = sqrt(var(y));        %得到rms
end

function [B,A] = oct3dsgn(Fc,Fs,N);

pi = 3.14159265358979;
f1 = Fc/(2^(1/6));
f2 = Fc*(2^(1/6));
Qr = Fc/(f2-f1);
Qd = (pi/2/N)/(sin(pi/2/N))*Qr;
alpha = (1 + sqrt(1+4*Qd^2))/2/Qd;
W1 = Fc/(Fs/2)/alpha;
W2 = Fc/(Fs/2)*alpha;
[B,A] = butter(N,[W1,W2]);

然后再算A计权的分贝值
data=csvread('F0068CH2.CSV',0,3);
%t=data(:,1);
%Fs=1/(t(1)-t(2));
noi=data(:,2);
[p,ff]=oct3bank(noi);
plot(ff,p);
Cj=[-50.5 -44.7 -39.4 -34.6 -30.2 -26.2 -22.5 -19.1 -16.1 -13.4 -10.9 -8.6 -6.6 -4.8 -3.2 -1.9 -0.8 0 0.6 1.0 1.2 1.3 1.2 1.0 0.5 -0.1 -1.1 -2.5 -4.3 -6.6]
La=0;
Lj=zeros(1,30);
for n=1:30
Lj(n)=10^(0.1*(10*log10(p(n)/(10*exp(-12)))+Cj(n)));
end
La=10*log10(sum(Lj));
这个算出来的是A计权的声功率级吗？
为什么跟实际测得差很多啊？

wangwenting

另外，通过三分之一倍频程滤波之后得到的p，是不是表示的是功率啊？
我对于信号处理的东西很陌生，希望大家能够指点啊，谢谢！

zhrelng

没加窗吧

350890819

我也不知道什么意思啊

檀转霞

请问你的问题解决了吗？