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; %预测时
2、段内的时间，s；m=5; %T时段内通过的列车数，列；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:mt(i)=(l/v(i)*(1+0.8*(d/l);end%速度修正量的计算 Cv，假设列车运行参考速度 v=100km/hfor i=1:mCv(i)=30*log10(V/
3、v(i);%i 列车的速度修正量end%线路和轨道结构的修正量 Ct,在这里假设不考虑Ct=0;%垂直指向性修正量 Cthetaif theta=-10end%几何发散衰减的计算 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);%屏障衰减的计算(无限长)Abarf=500;%声波频率，Hzdeta=10;%声程差，mc=340;%声速，m/stt=40*f*deta/(3*c);if tt=1Abar=10*log10(3*
4、pi*sqrt(1-tt2)/(4*atan(sqrt(1-tt)/(1+tt);elseAbar=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)*(17+(300/d);if Agr10Agr=0;else Agr=Agr1;end%声波传播途径引起的衰减量 AA=Adiv+Abar+Aatm+Agr;%j 列车噪声修正量 Cjfor i=1:mC(i)=C1(i)-A;end%预测点 i 列车通过时段内的等效声级 LpiLp0=87,85,84,86,87;%参考点 j 列车通过时段内的最大垂直指向性方向上的噪声辐射源强%Lp0=89;for i=1:mLp(i)=Lp0(i)+C(i);%预测点 i 列车通过时段内的等效声级end%预测点列车运行噪声等效声级，dB。sum=0;for i=1:msum=sum+t(i)*10(0.1*Lp(i);endLeql=10*log10(sum/T)
