Trainsperling_列车_sperlingmatlab平稳性_平稳性指标_

data=dlmread('I:\Protokoll\Schienenverkehr\实验一 车辆平稳性实验\chuixiang352.txt','');%读取实验数据 z=data(:,1); % 导入数据 N=length(z); %采样点数 figure(1) t=1/1000:1/1000:N/1000; %建立采样时间数组t plot(t,z); %画出以时间t为横坐标，垂向加速度为纵坐标的图像 xlabel('时间 单位：s'); ylabel('垂向加速度 单位：m/s^2'); title('垂向加速度-时间图像'); grid on; axis([0,l/1000,-2.5,2.5]); axis fill fs=1000; %采样频率fs为1000Hz fz=fft(z,N); %傅里叶变换 a=abs(fz); %各频率点模值 a=a*2/N; %该点的模值除以N/2就是对应该频率下的信号的幅度 n=0:N-1; f=n*fs/N; %各点频率,第一个点为频率为0的分量，故N值减1 figure(2) %画出垂向加速度与频率的关系图像 plot(f,a); xlabel('频率 单位：Hz'); ylabel('垂向加速度 单位：m/s^2'); title('垂向加速度与频率的关系'); axis([0,30,0,0.1]); axis fill for i=1:N %计算垂向平稳性 f(i)=(i-1)*fs/N; %各点代表的振动频率，第一个点为频率为0的分量，故i值减一 j(i)=(a(i)*100)^3; %公式中振动加速度单位为cm/s^2 end syms i1 i2 i3 %将频率进行分组 i1=round(0.5*N/fs); %频率<=0.5Hz的点 i2=round(5.9*N/fs); %频率<=5.9Hz的点 i3=round(20*N/fs); %频率<=20Hz的点 sum=0; for i=i1:i2 %分段计算平稳性指标 sum=sum+(0.896*(j(i)*0.325*f(i))^0.1)^10; end for i=i2+1:i3 sum=sum+(0.896*(j(i)*400/(f(i)^3))^0.1)^10; end for i=i3+1:N sum=sum+(0.896*(j(i)/f(i))^0.1)^10; end Wz=(sum)^(0.1) %总平稳性指标