《城市轨道交通运输噪声模式MATLAB程序.doc》由会员分享，可在线阅读，更多相关《城市轨道交通运输噪声模式MATLAB程序.doc（3页珍藏版）》请在金锄头文库上搜索。
1、%城市轨道交通运输噪声预测模式%个参数表示含义如下%ljj列车长度，m；%vjj列车运行速度，m/s；%Lpj预测点j列车通过时段内的等效声级，按式（A.24）计算，dB（A）；%Lp0j参考点j列车通过时段内最大垂向指向性方向上的噪声辐射源强，dB(A)；%Cjj列车噪声修正量，dB（A）；%C1j j 列车车辆、线路条件及轨道结构等修正量，dB（A）；%Cvj j 列车速度修正量，dB（A）；%Ct 线路和轨道结构的修正量，dB（A）；%Ctheta 垂向指向性修正量，dB（A）；%A声波传播途径引起的衰减量，dB。clear;clc;T=10; %预测时段内的时间，s；m=5; %T时段
2、内通过的列车数，列；theta=20; %声源到预测点方向与水平角的夹角；d=30; %预测点到轨道中心线的水平距离，m；d0=7.5; %源强的参考距离，ml=20; %列车长度，这里假设三辆列车长度相等，mV=60; %列车参考运行速度v=55,63,65,60,68; %i列车的行驶速度%v=60;%t=zeros(m,1)%tii列车通过时段的等效时间，s；for i=1:m t(i)=(l/v(i)*(1+0.8*(d/l);end%速度修正量的计算Cv，假设列车运行参考速度v=100km/hfor i=1:m Cv(i)=30*log10(V/v(i);%i列车的速度修正量end%
3、线路和轨道结构的修正量Ct,在这里假设不考虑Ct=0;%垂直指向性修正量Cthetaif theta=-10&theta24 Ctheta=-0.012*(24-theta)1.5else Ctheta=-0.075*(theta-24)1.5end%i列车车辆，线路条件反轨道结构的修正量for i=1:m C1(i)=Cv(i)+Ct+Ctheta;end%几何发散衰减的计算AdivAdiv1=d*atan(l/(2*d0)+(2*l2)/(4*d02+l2);Adiv2=d0*atan(l/(2*d)+(2*l2)/(4*d2+l2);Adiv=10*log10(Adiv1/Adiv2);
4、%屏障衰减的计算(无限长)Abarf=500;%声波频率，Hzdeta=10;%声程差，mc=340;%声速，m/stt=40*f*deta/(3*c);if tt=1 Abar=10*log10(3*pi*sqrt(1-tt2)/(4*atan(sqrt(1-tt)/(1+tt);else Abar=10*log10(3*pi*sqrt(tt2-1)/(2*log(tt+sqrt(tt2-1);end%空气吸收引起的衰减Aatma=2.8Aatm=a*(d-d0)/1000;%地面效应衰减Agr% S=200;%面积，m2% hm=S/d;hm=1.5Agr1=4.8-(2*hm/d)*(1
5、7+(300/d);if Agr10 Agr=0;else Agr=Agr1;end%声波传播途径引起的衰减量AA=Adiv+Abar+Aatm+Agr;%j列车噪声修正量Cjfor i=1:m C(i)=C1(i)-A;end%预测点i列车通过时段内的等效声级LpiLp0=87,85,84,86,87;%参考点j列车通过时段内的最大垂直指向性方向上的噪声辐射源强%Lp0=89;for i=1:m Lp(i)=Lp0(i)+C(i);%预测点i列车通过时段内的等效声级end%预测点列车运行噪声等效声级，dB。sum=0;for i=1:m sum=sum+t(i)*10(0.1*Lp(i);endLeql=10*log10(sum/T)
