请教通过三分之一倍频程计算声压级
我已经得到了声压信号通过下面的程序进行三分之一倍频程滤波
function = oct3bank(noi);
pi = 3.14159265358979;
Fs = 500000; %1/(t(2)-t(1)); % Sampling Frequency
N = 3; % Order of analysis filters.
%f=
ff=
%ff = ((2^(1/3)).^);% Exact center freq.
p = zeros(1,30);
% Design filters and compute RMS powers in 1/3-oct. bands
for i = 1:30
= oct3dsgn(ff(i),Fs,N);
y = filter(B,A,noi); %<-----滤波器
p(i) = sqrt(var(y)); %得到rms
end
function = 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;
= butter(N,);
然后再算A计权的分贝值
data=csvread('F0068CH2.CSV',0,3);
%t=data(:,1);
%Fs=1/(t(1)-t(2));
noi=data(:,2);
=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，是不是表示的是功率啊？我对于信号处理的东西很陌生，希望大家能够指点啊，谢谢！ 没加窗吧 我也不知道什么意思啊 请问你的问题解决了吗？
页:
[1]