data=xlsread('C:\Users\sang.linchen\Desktop\重庆\10号线数据Y--Z-环山公园--长河路.xlsx'); %读取实验数据
z1=data(:,38); % 导入原始数据
n=1:2:2000; %每2点取一个值
z2=z1(n,:); % 读取第n行数据
N=length(z2); %采样点数
A=mean(z2); %求均值
z=z2-A; %消除零点漂移
B=std(z); %求标准差
Aac=2.2*B; %求加速度
figure(1)
t=1/250:1/250:N/250; %建立采样时间数组t
plot(t,z); %画出以时间t为横坐标,垂向加速度为纵坐标的图像
xlabel('时间 单位:s');
ylabel('垂向加速度 单位:g');
title('垂向加速度-时间图像');
set(gca,'Fontname','Monospaced');
grid on;
axis([0,5.5,-0.5,0.5]);
axis fill;
fs=250; %采样频率fs为250Hz
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('垂向加速度 单位:g');
title('垂向加速度与频率的关系');
set(gca,'Fontname','Monospaced');
axis([0,30,0,0.07]);
axis fill
for i=1:N %计算垂向平稳性
f(i)=(i-1)*fs/N; %各点代表的振动频率,第一个点为频率为0的分量,故i值减一
j(i)=a(i)^3; %公式中振动加速度单位为g
end
syms i1 i2 i3 i4 %将频率进行分组
i1=round(0.5*N/fs); %频率<=0.5Hz的点
i2=round(5.9*N/fs); %频率<=5.9Hz的点
i3=round(20*N/fs); %频率<=20Hz的点
i4=round(30*N/fs); %频率<=30Hz的点
sum=0;
for i=i1+1:i2 %分段计算平稳性指标
sum=sum+(3.57*(j(i)*0.325*f(i))^0.1)^10;
end
for i=i2+1:i3
sum=sum+(3.57*(j(i)*400/(f(i)^3))^0.1)^10;
end
for i=i3+1:i4
sum=sum+(3.57*(j(i)/f(i))^0.1)^10;
end
Wz=(sum)^(0.1); %总平稳性指标
xlswrite('C:\Users\sang.linchen\Desktop\重庆\重庆.xlsx',Wz,1,'M22'); %调出平稳性指标
评论0