sperling指标 matlab,sperling指标计算实验报告.docx

Sperling指标计算实验报告

一、 计算原理：

数据给出的信息为时域中的加速度信号。时间的单位为sec，加速度单位为m/s^2

1、将时域的信号转化为频域下。由于matlab中，FFt快速傅里叶变化计算的值就是离散的，直接能够用于公式WΣ=10i=1n[ai3fi*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/s^2

%其采样频率为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;

for i=1:length(Ysim)

if fsim(i)>