sperling指标计算实验报告

1、Sperling指标计算实验报告一、计算原理：数据给出的信息为时域中的加速度信号。时间的单位为sec,加速度单位为m/sA21、将时域的信号转化为频域下。由于matlab中，FFt快速傅里叶变化计算的值就是离散的，直接能够用于公式=151罕?F(f)进行计算。FFT所取得的点数多少对于计算fiSperling指标的影响不大。2、另一种方法是将整个FFT的频率区段分频，将频率区段中的能量集中到一点上来，再计算Sperling指标。二、计算程序：主程序：%读入数据stUT6_LH2K_v120_rear.dat'H,V,T=SFileRead(str);H=H(1:length(H)-1)

2、;%只取10000个点V=V(1:length(V)-1);T=T(1:length(T)-1);%加速度单位为cm/sA2%其采样频率为500Hz,用1s除以其时间间隔即可%垂向Sperling指标SperlingV,yV,fV,ampV,fStrV=constantSperling(V,T);%直接将yV计算Sperling,yV就是每个频率下的幅值，点数变多，自然会调整Sum=0;yV=yV(fV<50);fV=fV(fV<50);Ysim=yV;fsim=fV;fori=1:length(Ysim)iffsim(i)>=50breakendiffsim(i)<0

3、.5continueelseiffsim(i)>=0.5&&fsim(i)<5.9Fi=0.325*fsim(i)A2;elseiffsim(i)>=5.9&&fsim(i)<20Fi=400/fsim(i)A2;elseFi=1;endSum=Sum+(Ysim(i)*100)A3/fsim(i)*Fi;%*100因为输入单位为m/s2，而计算单位为cm/s2endSum=0.896*SumA(10A(-i);%分频计算得到的Sperling指标就是SperlingVfigure(1);subplot(3,1,1);plot(T,V);

4、xlabel('time/sec');ylabel('Accelaration/m*sA-2');title(str)subplot(3,1,2);plot(fV,yV);xlabel('Frequency/Hz');ylabel('Accelaration/m*sA-2');title('直接计算垂向Sperling：',num2str(Sum,4)subplot(3,1,3);stem(fStrV(ampV>1E-3),ampV(ampV>1E-3);xlabel('Frequency/Hz

5、');ylabel('Accelaration/m*sA-2');title('分段垂向Sperling:',num2str(SperlingV,4)%横向Spering指标SperlingH,yH,fH,ampH,fStrH=constantSperling(H,T);%计算频率分段的SperlingSum=0;fori=1:length(ampH)iffStrH(i)<0.5continueelseiffStrH(i)>=0.5&&fStrH(i)<5.4Fi=0.8*fStrH(i)A2;elseiffStrH(i

6、)>=5.4&&fStrH(i)<26Fi=650/fStrH(i)A2;elseFi=1;endSum=Sum+(ampH(i)*100)A3/fStrH(i)*Fi;endSperlingH=0.896*SumA(10A(-1);%不分频，计算sperlingSum=0;yH=yH(fH<50);fH=fH(fH<50);Ysim=yH;fsim=fH;fori=1:length(Ysim)ifYsim(i)>50breakendiffsim(i)<0.5continueelseiffsim(i)>=0.5&&fsi

7、m(i)<5.4Fi=0.8*fsim(i)A2;elseiffsim(i)>=5.4&&fsim(i)<26Fi=650/fsim(i)A2;elseFi=1;endSum=Sum+(Ysim(i)*100)A3/fsim(i)*Fi;%*100因为输入单位为m/sA2，而计算单位为cm/sA2endSum=0.896*SumA(10A(-1);figure(2);subplot(3,1,1);plot(T,H);xlabel('time/sec');ylabel('Accelaration/m*sA-2');title(st

8、r)subplot(3,1,2);plot(fH,yH);xlabel('Frequency/Hz');ylabel('Accelaration/m*sA-2');title('直接计算横向Sperling：',num2str(Sum,4)subplot(3,1,3);stem(fStrH,ampH);xlabel('Frequency/Hz');ylabel('Accelaration/m*sA-2');title('分段横向Sperling:',num2str(SperlingH,4)子程序：数

9、据读入程序：functionH,V,T=SFileRead(fileStr)f=fopen(fileStr);fori=1:9fscanf(f,'%s',1);endfori=1:3fscanf(f,'%s',1);endi=1;whilefeof(f)H(i)=fscanf(f,'%f',1);V(i)=fscanf(f,'%f',1);T(i)=fscanf(f,'%f',1);i=i+1;endfclose(f);end三、实验结果:评价以客车Sperling指标限界为标准:数据文件T7_LH2k_v140_

10、frontT1LH2kv140_rearT6LH2Kv120_rearT2lh2zv160_rearT3lh2zv160_frontT4LH2Kv160_rearT5_LH2K_v160_front横向2.7533.2522.6022.3962.463.5162.973垂向2.9172.8512.7321.5871.52.9283.228合格不合格良好舒适舒适不合格不合格LH2K之间对比120140160横向（front）一2.7532.973垂向（front）一2.9173.228横向(rear)2.6023.2523.516垂向(rear)2.7322.8512.928随着速度加快，spe

11、rling指标开始恶化四、附录：T6_LH2K_v120_rear.dat七脚J冲旨Al|1Q01D15JQK»ill屿DI3SaoF4.WIWWP-Wi.leeet-fSperling：2.732SIQ<口14114Rtmpwc：不升，汁可为峭间am?_l_nvEED>4-u4r,!_土EgmnvffiQN|iHI£.|UHZ目堀，ftK205USU45annlanSperling：2.602T7_LH2k_v140_front.datB1Q13Im/iwT-jJJIiifBB)5p4!lryiJMT16ItXIQ15133531FlMhM1E0421

12、3;d«5Sl<ZJ25nflMriMliirflli<3ffiS32>，FRJEEJi'al-rL#Sperling：2.917E.£nM>t4f。10EQ4：m株旬F闯qunc衲七4iii)ENIDI'S!5IDDsc26an202&3DSperling：2.753T1_LH2k_v140_rear.datYf1Ilfift鼻鼻L.藤尊4停尊XH3ll>与"FM«tMirrHvFrefUfirFnE寸讨羊何陷l.i】工5*WWf叶看臼|勤2S&1Sperling：2.8511b七柚卜血1s

13、i不时为皓同1SID15M无刃3E如临MISNEcaljT'J型村:M神斯何割*（忡m'.'r.te2p*rQl_m3MSperling：3.252T2_lh2z_v160_rear.datThv'fe=lriwdalfi时uIJfc#讷hirrtw_a.审M-3T-3-4-廿讶.-Apr旨d-lnwSperling：1.587il:.q.-'2LIa->-lw!l:.r-E-LriLEEBs-lq-1-一*网炒朴仙*W林JL1Lj11ID14Vfi省H-BWRSpmngrW沸111111li1：15心IE20若JD等40%至FmiiwitjUHt

14、向三的Msmnei争-fL±_T_i礼m中彳行母Q£011UJLr25»ISm若Ml3£|UftSu侦|«|岫|哪1X.16C/5IESperling：2.396T3_lh2z_v160_front.datiJJfi2-.<a4qEE3=W*34BIf16ifl1mW5?ET?t1tiif菱职网fIjSMufcfXG一-quk.L1-3ja75X)uS.UIM的Sperling：1.51I11.fM心"佻5；IIIIII&1Q>口MII不,Tf*VMR$p>idi|lHIirSperling:2.46T4_LH2K_v160_rear.datIBiinHS'3DSperling:2.928号岑5,一IKO13T抻1UiJbl&I-Aah1I1111111n-、_I4.耳食MlSqJS;型03aqjaLmq土1,任jnqjaLa4uaa&25Frpr|lffiKyAliMUI-计asmaQaQI1?F夏！#jzMSVSperling:3.516T5_LH2K_