1/3倍频程程序
function = octspectrum(yy,fs,fl,fh)%yy信号;fs采样率;fl为分析频率下限,fh为分析频率上限
f=;
fc0=;
oc6=2^(1/6);
index = find(fc0 <= fh & fc0 >= fl);
fc = fc0(index);
nc=length(fc);
n=length(yy);
nfft=2^nextpow2(n);
a=fft(yy,nfft);
for j=1:nc,
fl=fc(j)/oc6;
fu=fc(j)*oc6;
nl=round(fl*nfft/fs+1);
nu=round(fu*nfft/fs+1);
b=zeros(1,nfft);
b(nl:nu)=a(nl:nu);
b(nfft-nu+1:nfft-nl+1)=a(nfft-nu+1:nfft-nl+1);
c=ifft(b,nfft);
%%%计算对应每个中心频率段的有效值,与幅度谱相差3分贝
yc(j)=sqrt(var(real(c(1:n))));
end
return songzy41 发表于 2007-6-1 14:03
为求绝对声压值,并进一步计算出声压级,对于麦克风(传声器)当然必须标定,测量它的灵敏度。B&K公司有 ...
你好,可以将这个A计权程序发到我邮箱吗,我无法下载这个附件,谢谢!!
ml_tang1989@163.com 你好楼主,可以将这个A计权的程序代码发给我吗,我无法下载这个附件,谢谢你!!
ml_tang1989@163.com
谢谢!
页:
1
[2]