a计权声功率级计算公式_超高层建筑气动噪声总声压级和A计权声压级的计算

用FLUENT中的FW-H噪声板块计算建筑的气动噪声，现在利用大涡模拟得到了监测点的瞬时风压时程，接下来要计算这个监测点的总声压级和A计权声压级，在网上找了个MATLAB程序，自己修改整合了之后的程序如下：

%输入需要分析的时程数据

x=importdata('xx.txt');

fs=1000;%采样频率

p0=2*10^(-5);%基准声压

y=x(:,2);%声压

n=length(y);%输入数据的长度

t=0:1/fs:(n-1)/fs;%时间

%%%%%%绘制声压时程图

figure(1);

plot(t,y);

xlabel('时间(s)');

ylabel('声压(pa)');

title('声压时程图');

grid on;

%定义三分之一倍频程的中心频率

f=[1.00 1.25 1.600 2.00 2.50 3.15 4.00 5.00 6.30 8.0];

f1=[20.00 25.0 31.5 40.0 50.0 63.0 80];

fc=[f1,100*f,1000*f,10000,12500,16000,20000];

%A计权声压级各频率对应的计权系数

cf=[-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,-9.3];

%加汉宁窗<