佛山网站设计哪里好,网站开发英文文献,app开发免费,门头沟网站开发0.随机路面
仿真时间200秒 车速15m/s 采样间隔0.25秒#xff08;4Hz#xff09; A级路面平均振幅#xff1a;0.0044米 B级路面平均振幅#xff1a;0.0087米 C级路面平均振幅#xff1a;0.0173米 D级路面平均振幅#xff1a;0.0346米
1.先上代码其中之一
clc
clear
clo…0.随机路面
仿真时间200秒 车速15m/s 采样间隔0.25秒4Hz A级路面平均振幅0.0044米 B级路面平均振幅0.0087米 C级路面平均振幅0.0173米 D级路面平均振幅0.0346米
1.先上代码其中之一
clc
clear
close all
%%
SimTime200;
dt0.01;%仿真步长
time0:dt:SimTime;
sim_step length(time);
Ntlength(time); % 采样点可能要修改
u10; % m/s
df1/(Nt*dt); % 采样频率间隔
f0:df:1/(2*dt); % 采用频率一半(可修改
nflength(f); % 采样点的个数
nf/u; % 采样空间频率,F是浮动的随车速的改变而改变
%获得标准空间路面谱A级
PaLMP(16,nf,n);
paPa/u; % 转化为频谱
%获得标准空间路面谱B级
PbLMP(64,nf,n);
pbPb/u; % 转化为频谱
%获得标准空间路面谱C级
PcLMP(256,nf,n);
pcPc/u; % 转化为频谱
%获得标准空间路面谱D级
PdLMP(1024,nf,n);
pdPd/u; % 转化为频谱
%%
%1输出标准路面空间频率谱G(n)结果与车速、SimTime、dt无关
figure(1)
loglog(n,Pa); %1输出标准路面空间频率谱G(n)
hold on
loglog(n,Pb);
hold on
loglog(n,Pc);
hold on
loglog(n,Pd);
xlabel(空间频率n[1/m])
ylabel(路面不平度位移功率谱密度Gq(n)[m^3])
legend(A,B,C,D)%2输出时间频率谱G(f)结果与车速有关
% figure(2)
% loglog(f,pa)
% hold on
% loglog(f,pb)
% hold on
% loglog(f,pc)
% hold on
% loglog(f,pd)
% xlabel(时间频率[Hz])
% ylabel(G(f)[m^2/Hz]) % 2输出时间频率谱G(f)
% legend(A,B,C,D)
2.参考链接
链接1 链接2 链接3
3.笔记
Sine_RoadRoughnessGradeClassCal_230713.m中
3.1. 降SimTime可以降等级升omiga会使改线条等级范围变大降车速可降等级
但是在时间-输入关系不变的情况下低车速对应高空间频率路面高车速对应低空间频率路面
代码一SimTime20/4;dt0.01;omiga2pi1.5;A0.05;V18;
%%
%1.正弦路面
SimTime20/4;
dt0.01;%仿真步长
time0:dt:SimTime;
sim_step length(time);
omiga2*pi*1.5;
A0.05;%振幅
Z1A*sin(time*omiga);%路面高程值
Z_p1 A*cos(time*omiga)*omiga;%路面高程值导数
xZ1;Ntlength(x); % 采样点可能要修改
V18; % km/h结果C级路面 注意此图中的标准等级路面PSD-(1/n)函数线条与汽车理论208页图6-5一致
代码一点五SimTime20;dt0.01;omiga2*pi;A0.05;V18;
%%
%1.正弦路面
SimTime20;
dt0.01;%仿真步长
time0:dt:SimTime;
sim_step length(time);
omiga2*pi;
A0.05;%振幅
Z1A*sin(time*omiga);%路面高程值
Z_p1 A*cos(time*omiga)*omiga;%路面高程值导数
xZ1;Ntlength(x); % 采样点可能要修改
V18; % km/h结果B级路面
代码二SimTime20/2;dt0.01;omiga2pi2;A0.019;V18;
SimTime20/2;
dt0.01;%仿真步长
time0:dt:SimTime;
sim_step length(time);
omiga2*pi*2;
A0.019;
Z1A*sin(time*omiga);%路面高程值
Z_p1 A*cos(time*omiga)*omiga;%路面高程值导数
xZ1;
dt0.01;
Ntlength(x); % 采样点可能要修改
V18; % km/h结果A级路面
代码二点五SimTime20/4;dt0.01;omiga2*pi;A0.019;V9;
SimTime20/4;
dt0.01;%仿真步长
time0:dt:SimTime;
sim_step length(time);
omiga2*pi;
A0.019;
Z1A*sin(time*omiga);%路面高程值
Z_p1 A*cos(time*omiga)*omiga;%路面高程值导数
xZ1;
dt0.01;
Ntlength(x); % 采样点可能要修改
V9; % km/h结果B级路面
