sperling指标计算实验报告
Sperling 指标计算实验报告 一、 计算原理: 数据给出的信息为时域中的加速度信号。时间的单位为 sec,加速度单位为 m/s^2 1、将时域的信号转化为频域下。由于 matlab 中,FFt 快速傅里叶变化计算的值就是离散的, 直接能够用于公式 进行计算。FFT 所取得的点数多少对于计算 𝑊 Σ = 10 ∑ 𝑛 𝑖 =1 [ ( 𝑎 𝑖 ) 3 𝑓 𝑖 ∗ 𝐹(𝑓)] 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/s^2 %其采样频率为 500Hz,用 1s 除以其时间间隔即可 %% 垂向 Sperling 指标 [SperlingV,yV,fV,ampV,fStrV]=constantSperling(V,T); %直接将 yV 计算 Sperling,yV 就是每个频率下的幅值,点数变多,自然会调整 Sum=0; yV=yV(fV=50breakendif fsim(i)=0.5 xlabel( Frequency/Hz ); ylabel( Accelaration/m*s^-2 ); title([ 分段垂向 Sperling : ,num2str(SperlingV,4)]) %% 横向 Spering 指标 [SperlingH,yH,fH,ampH,fStrH]=constantSperling(H,T); %计算频率分段的 Sperling Sum=0; for i=1:length(ampH)if fStrH(i)=0.5elseFi=1;endSum=Sum+(Ysim(i)*100)^3/fsim(i)*Fi;%*100 因为输入单位为 m/s^2 ,而计算单位为 cm/s^2 end Sum=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=(fileStr); for i=1:9fscanf(f, %s ,1); endfor i=1:3fscanf(f, %s ,1); end i=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; end fclose(f); end 三、 实验结果: 评价以客车 Sperling 指标限界为标准: 数据 文件 T7_LH2k_ v140_fron t T1_LH2k_ v140_rea r T6_LH2K_ v120_rea r T2_lh2z_ v160_rea r T3_lh2z_ v160_fron t T4_LH2K_ v160_rea r T5_LH2K_ v160_fron t 横向 2.753 3.252 2.602 2.396 2.46 3.516 2.973 垂向 2.917 2.851 2.732 1.587 1.5 2.928 3.228 合格 不合格 良好 舒适 舒适 不合格 不合格 LH2K 之间对比 120 140 160 横向(front) —— 2.753 2.973 垂向(front) —— 2.917 3.228 横向(rear) 2.602 3.252 3.516 垂向(rear) 2.732 2.851 2.928 随着速度加快,sperling 指标开始恶化