路面不平度构造以及功率谱密度生成检测.doc

《路面不平度构造以及功率谱密度生成检测.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、路面不平度构造以及功率谱密度生成检测一构造C级路面不平度构造，采用谐波叠加法（又称为三角级数法），常见空间频率范围为：0.011

2、fl)/(2*N);%中心频率faik=2*pi*rand(1,N);%相位角ak=zeros(1,N);%系数fori=1:Nak(i)=sqrt(2*G0*n0^2*v*(fh-fl)/N/(fl+(i-1)*(fh-fl)/N)/(fl+i*(fh-fl)/N));endt=0:1/fh/2:L/v;%时间采样频率为最高时间频率的两倍Lt=length(t);qk=zeros(N,Lt);fori=1:Nqk(i,:)=ak(i)*sin(2*pi*fk(i)*t+faik(i));endqt

3、=zeros(Lt,1);qt(:,1)=sum(qk)';plot(t,qt(:,1)),gridontitle('c级路面时域不平度q（t）')xlabel('时间t/s')ylabel('位移q/m')GC=sqrt(sum(qt.^2)/length(qt))%路面不平度均方值L=0:0.05/2.83:1000;Ll=length(L);qll=zeros(N,Ll);nk=0.011+(2.83-0.011)/(2*N):(2.83-0.011)/N:2.83-(2.83-0.011)/

4、(2*N);fori=1:Nqll(i,:)=ak(i)*sin(2*pi*nk(i)*L+faik(i));endql=zeros(Ll,1);ql(:,1)=sum(qll)';figure(2);plot(L,ql(:,1)),gridonGC=0.0143路面功率谱密度：figure(3)ns=10*nh;%空间采样频率[pxx,m]=psd(qt,512,ns);%512为快速傅里叶变换点数loglog(m,pxx)gridon;holdon;n=0.011:0.001:2.83;y1=1

5、28e-6./n.^2;y2=512e-6./n.^2;loglog(n,y1);gridon;loglog(n,y2);gridon;xlabel('空间频率/(1/m)')ylabel('路面公路车谱密度G(n)/[m^2/(1/m)]')二构造F级路面不平度，常见空间频率范围为：0.011

6、距离fl=n1*v;fh=nh*v;%时间频域fk=fl+(fh-fl)/(2*N):(fh-fl)/N:fh-(fh-fl)/(2*N);%中心频率faik=2*pi*rand(1,N);%相位角ak=zeros(1,N);%系数fori=1:Nak(i)=sqrt(2*G0*n0^2*v*(fh-fl)/N/(fl+(i-1)*(fh-fl)/N)/(fl+i*(fh-fl)/N));endt=0:1/fh/2:L/v;%时间采样频率为最高时间频率的两倍Lt=length(t);qk=zeros

7、(N,Lt);fori=1:Nqk(i,:)=ak(i)*sin(2*pi*fk(i)*t+faik(i));endqt=zeros(Lt,1);qt(:,1)=sum(qk)';plot(t,qt(:,1)),gridontitle('f级路面时域不平度q（t）')xlabel('时间t/s')ylabel('位移q/m')GC=sqrt(sum(qt.^2)/length(qt))%路面不平度均方值L=0:0.05/2.83:1000;Ll=length(L);qll=zeros(N,Ll);n

8、k=0.011+(2.83-0.011)/(2*N):(2.83-0.011)/N:2.83-(2.83-0.011)/(2*N);fori=1:Nqll(i,:)=ak(i)*sin(2*pi*nk(i)*L+faik(i));endql=zeros(Ll,1);ql(:,1)=sum(qll)';figure(2);plot(L,ql(:,1)),gridon运算结果：GC=0.1210路面功率谱密度：figure(3)ns=10*nh;%空间采样频率[pxx,m