汽车平顺性加速度均方根值.txt
《汽车平顺性加速度均方根值.txt》由会员分享,可在线阅读,更多相关《汽车平顺性加速度均方根值.txt(1页珍藏版)》请在冰豆网上搜索。
![汽车平顺性加速度均方根值.txt](https://file1.bdocx.com/fileroot1/2022-10/29/eaae6a2e-b37d-494b-8ac0-f9687545bc64/eaae6a2e-b37d-494b-8ac0-f9687545bc641.gif)
%�����������۵�����ƽ˳�Ե����۷��������û��������۷���������Ȩ���ٶȾ�����ֵ��������ƽ˳�ԣ���������
clc;
clearall;
closeall;
t=0:
0.001:
5;
y=20*sin(2*pi*10.*t);
figure
(1)
plot(t,y,'r');
gridon
%���ٶ�ʱ������ͨ����Ӧ��Ƶ�ʼ�Ȩ�����˲��õ���Ȩ���ٶ�ʱ�����̣���������
%�Լ��ٶ�ʱ�����̽���Ƶ�����õ��������ܶȣ���������������Ȩ���ٶȾ�����ֵ
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');
gridon
set(gca,'Xscale','log','Yscale','log');
acc_x=f;
acc_y=Px;
ww=zeros;
acc_num=0;
fori=1:
length(acc_x)%Ƶ�ʼ�Ȩ����
ifacc_x(i)>=0.5&&acc_x(i)<2
ww(i)=0.5;
elseifacc_x(i)>=2&&acc_x(i)<4
ww(i)=acc_x(i)/4;
elseifacc_x(i)>=4&&acc_x(i)<12.5
ww(i)=1;
elseifacc_x(i)>=12.5&&acc_x(i)<80
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;%��Ȩ���ٶȾ�����ֵ