sperling指标计算实验报告.docx

sperling指标计算实验报告一、 计算原理：数据给出的信息为时域中的加速度信号。时间的单位为sec，加速度单位为m/s21、将时域的信号转化为频域下。由于matlab中，fft快速傅里叶变化计算的值就是离散的，直接能够用于公式w=10i=1nai3fi*f(f)进行计算。fft所取得的点数多少对于计算sperling指标的影响不大。2、另一种方法是将整个fft的频率区段分频，将频率区段中的能量集中到一点上来，再计算sperling指标。二、 计算程序：主程序：% 读入数据str=t6_lh2k_v120_rear.dat;h,v,t=sfileread(str);h=h(1:length(h)-1);%只取10000个点v=v(1:length(v)-1);t=t(1:length(t)-1);%加速度单位为cm/s2%其采样频率为500hz,用1s除以其时间间隔即可% 垂向sperling指标sperlingv,yv,fv,ampv,fstrv=constantsperling(v,t);%直接将yv计算sperling，yv就是每个频率下的幅值，点数变多，自然会调整sum=0;yv=yv(fv50);fv=fv(fv=50 break end if fsim(i)=0.5&fsim(i)=5.9&fsim(i)1e-3),ampv(ampv1e-3);xlabel(frequency/hz);ylabel(accelaration/m*s-2);title(分段垂向sperling： ,num2str(sperlingv,4)% 横向spering指标sperlingh,yh,fh,amph,fstrh=constantsperling(h,t);%计算频率分段的sperlingsum=0;for i=1:length(amph) if fstrh(i)=0.5&fstrh(i)=5.4&fstrh(i)26 fi=650/fstrh(i)2; else fi=1; end sum=sum+(amph(i)*100)3/fstrh(i)*fi;endsperlingh=0.896*sum(10(-1);%不分频，计算sperlingsum=0;yh=yh(fh50);fh=fh(fh50 break end if fsim(i)=0.5&fsim(i)=5.4&fsim(i)26 fi=650/fsim(i)2; else fi=1; end sum=sum+(ysim(i)*100)3/fsim(i)*fi;%*100因为输入单位为m/s2，而计算单位为cm/s2endsum=0.896*sum(10(-1);figure(2);subplot(3,1,1);plot(t,h);xlabel(time/sec);ylabel(accelaration/m*s-2);title(str)subplot(3,1,2);plot(fh,yh);xlabel(frequency/hz);ylabel(accelaration/m*s-2);title(直接计算横向sperling： ,num2str(sum,4)subplot(3,1,3);stem(fstrh,amph);xlabel(frequency/hz);ylabel(accelaration/m*s-2);title(分段横向sperling： ,num2str(sperlingh,4)子程序：数据读入程序：function h,v,t=sfileread(filestr)f=fopen(filestr);for i=1:9 fscanf(f,%s ,1);endfor i=1:3 fscanf(f,%s ,1); endi=1;while feof(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_frontt1_lh2k_v140_reart6_lh2k_v120_reart2_lh2z_v160_reart3_lh2z_v160_frontt4_lh2k_v160_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随着速度加快，sperling指标开始恶化四、 附录：t6_lh2k_v120_rear.datsperling：2.732sperling：2.602t7_lh2k_v140_front.datsperling：2.917sperling：2.753t1_lh2k_v140_rear.datsperling：2.851sperling：3.252t2_lh2z_v160_rear.datsperling：1.587sperling：2.396t3_lh2z_v160_front.datsper