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）;

%