高速列车平稳性指标计算

clear;clc;

aa=outaCheHC;
%均方根值
Arms=sqrt（（sum（aa.*aa））/length（aa））; 
fprintf（‘被动横向加速度均方根值为: %5.4f\n‘Arms）;
%计算平稳性指标
fc=50;
t_buchang=0.001;%步长
f_fangzhen=1/t_buchang;%仿真频率
aa1=outaCheHC（1:20*f_fangzhen）;%1-20000个数据
bb=reshape（aa1f_fangzhen*210）;%重塑一个200*10的矩阵	
freq=（0.5:0.5:f_fangzhen/2）;	%0.5：0.5:500
Ff（1:10）=0.8.*freq（1:10）.^2;    %freq（1:10）=0.5：0.5：5.0
Ff（11:52）=650./freq（11:52）.^2;  %freq（11:52）=5.5：0.5:26
Ff（53:f_fangzhen/2）=1;	        %freq（53:500）=26.5——250
cc=abs（fft（bb/9.8））;
dd=mean（cc‘）./（f_fangzhen）;	    %频率内加速度平均值
dd（1）=[];
Wz=7.08*（sum（dd（1:fc/0.5）.^3.*Ff（1:fc/0.5）./freq（1:fc/0.5）））^0.1;%频率计算从0.5~fc HZ
fprintf（‘被动平稳性指标为: %5.4f\n‘Wz）;
%加速度峰值
Amax=max（abs（outaCheHC））;   
fprintf（‘被动横向加速度峰值为:%5.4f\n‘Amax）;

%