汽车平顺性加速度均方根值
% 根据汽车理论第六章平顺性的评价方法，采用基本的评价方法——加权加速度均方根值法来评价平顺性，程序如下
clc;
clear all;
close all;
t=0:0.001:5;
y=20*sin(2*pi*10.*t);
figure(1)
plot(t,y,'r');
grid on
% 加速度时间历程通过相应的频率加权函数滤波得到加权加速度时间历程，计算如下
% 对加速度时间历程进行频谱分析得到功率谱密度，按下述程序计算加权加速度均方根值
Fs=1000; % 采样频率
Nfft=length(y); % 采样点数
range='onesided';
window=hamming(length(y)); % 窗函数
noverlap=0; % 重叠度
[Px,f]=pwelch(y,window,noverlap,Nfft,Fs,range); % FFT
figure(2)
plot(f,Px,'r');
grid on
set(gca,'Xscale','log','Yscale','log');
acc_x=f;
acc_y=Px;
ww=zeros;
acc_num=0;
for i=1:length(acc_x) % 频率加权函数
if acc_x(i)>=0.5 && acc_x(i)
ww(i)=0.5;
elseif acc_x(i)>=2 && acc_x(i)
ww(i)=acc_x(i)/4;
elseif acc_x(i)>=4 && acc_x(i)
ww(i)=1;
elseif acc_x(i)>=12.5 && acc_x(i)
ww(i)=12.5/acc_x(i);
else
ww(i)=0;
end
acc_num=acc_num+ww(i)^2*acc_y(i)*(acc_x(2)-acc_x(1)); % 计算加权均方根值
end
acc_Q=(acc_num)^0.5; % 加权加速度均方根值