三分之一频域加权倍频程计算分频率计权加速度振级.docxVIP

计算 三分之一频域加权倍频程计算算法matlab程序 --分频率计权加速度振级 1. 背景 结合上海软土地区典型工程实例，现场实测地铁减振前后隧道内及邻近建筑物内的振动加速度，据此对地铁轨道减振效果进行分析评估；在此基础上，提出了有限元与无限元相结合的数值分析简化预估方法，预测地铁交通激发的周边建筑物振动响应，并将计算分析结果与实测结果进行对比，分析地铁交通引起邻近建筑物振动随楼层的变化规律，探讨不同因素对建筑物振动响应的影响。研究结果表明：1）采用 Vanguard 扣件较科隆蛋扣件隧道内道床振级降低约 16dB，建筑物楼板中的振级降低 10.9~21.1dB，减振效果明显；2）随着楼层的增高，楼板中振动的振级呈现先减小后增大的现象；3）采用提出的数值分析得到的建筑物楼板振级和实测结果吻合较好，表明分析方法合理可行，为预测地铁周边建筑物振动提供一种快速简便的分析方法。 2. 参考文献 《地铁交通引起的建筑物振动实测与数值模拟分析》 岳建勇 3. 重要算法 某点的三分之一频域加权倍频程计算算法。 4.程序设计代码，可以直接使用 （1）主函数main %% 某点的三分之一频域加权倍频程 close all; clear all; clc; %% da=xlsread(.\da.xlsx);%导入要分析的激励加速度信号 x=da(:,1); y=da(:,2)/100; figure; plot(x,y); xlabel(时间t(s)); ylabel(加速度a（m/s^2）); title(加速度时程曲线); L=length(x); Fs=1000; [B tem]=BPC_ZJ(x,Fs); Y=abs(tem(end/2+1:end))/length(tem); f=(0:length(Y)-1)*80/length(Y); figure; plot(f,Y); xlabel(频率/Hz); ylabel(幅值/m^s/Hz); set(gca,xlim,[1,80],xtick,[0,20,40,60,80]);%%/////// H=B; %%绘制三分之一倍频程数据 f=[1.00,1.25,1.60,2.00,2.50,3.15,4.00,5.00,6.30,8.00]; fc=[0.50,0.63,0.80,f,10*f,100*f];%%定义三分之一倍频程的中心频率 figure; plot(fc(1:24),H(1:24)); %%%%%%横坐标根据论文需要进行更改//////////////////////////////// set(gca,xlim,[1,80],xtick,[0,20,40,60,80]);%%/////// xlabel(频率/Hz); ylabel(幅值/m^s/Hz); grid on; （2）被调用的脚本函数BPC_ZJ.m function [yq a] =BPC_ZJ(x,Fs) n=length(x); %样点个数 f=[1.00,1.25,1.60,2.00,2.50,3.15,4.00,5.00,6.30,8.00]; fc=[0.50,0.63,0.80,f,10*f,100*f];%%定义三分之一倍频程的中心频率 oc6=2^(1/6); nc=length(fc); %%取中心频率总的长度 nfft=2^nextpow2(n); %%大于并最接近n的2的幂次方长度 a=fft(x,nfft); %%进行FFT变换 for j=1:nc fl=fc(j)/oc6; %下限频率 fu=fc(j)*oc6; %上限频率 nl=round(fl*nfft/Fs+1); %下限频率对应的序号 nu=round(fu*nfft/Fs+1); %上限频率对应的序号 %如果上限频率大于折叠频率则循环中断 if fuFs/2 break; end %%%以每个中心频率段为通带进行带通频域滤波 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); %计算对应每个中心频率段的有效值 ycc(j)=sqrt(var(real(c(1:n)))); yc(j)=20.*log10(ycc(j)./10.^(-6)); %%分频加速度振级 QdB = [-7.56,-6.77,-6.44,-6.33,-6.29,-6.13,-5.50,-3.97,-1.86,-0.31,0.33,0.46,0.32,-0.10,-0.93,-2.22,-3.91, ... -5.84,-7.89,-10.01,-12.21,-14.62,-17.47,-