Sperling平稳性计算

data=dlmread（‘I:\Protokoll\Schienenverkehr\实验一 车辆平稳性实验\chuixiang352.txt‘‘‘）;%读取实验数据
z=data（:1）;  % 导入数据
N=length（z）;  %采样点数
figure（1）
t=1/1000:1/1000:N/1000; %建立采样时间数组t
plot（tz）;  %画出以时间t为横坐标，垂向加速度为纵坐标的图像
xlabel（‘时间 单位：s‘）;
ylabel（‘垂向加速度 单位：m/s^2‘）;
title（‘垂向加速度-时间图像‘）;
grid on;
axis（[0l/1000-2.52.5]）;
axis fill
fs=1000;   %采样频率fs为1000Hz
fz=fft（zN）;    %傅里叶变换
a=abs（fz）;    %各频率点模值
a=a*2/N;    %该点的模值除以N/2就是对应该频率下的信号的幅度
n=0:N-1;        
f=n*fs/N;       %各点频率第一个点为频率为0的分量，故N值减1
figure（2）       %画出垂向加速度与频率的关系图像
plot（fa）;
xlabel（‘频率 单位：Hz‘）;
yla