白噪声法生成路面不平度时域信号

滤波白噪声法生成随机路面1.原理200()()()1()()nGqnGqnnfunGqfGqnufun22002()()()uGqGqnn……1设截止频率为0，则1式变换成2200220()()()uGqGqnn……2假设2式是白噪声激励下一阶线性系统的响应，设频率响应函数()aHjbj（a,b为未知常数）……32()()GqHjS（S为白噪声信号，设1S）……4连理式2、3、4解得000()anGqnub该线性系统的频率响应函数为000()()nGqnuHjj系统的微分方程为000()()()()nGqnuwtqtqt（()wt为白噪声的时域信号，()qt为激励，即为路面谱时域信号）状态方程000()()()()qtqtnGqnuwt由于000nu得0000()()()()qtnuqtnGqnuwt2.simulink仿真输入参数u=10;%速度（单位m/s）Gqn0=4096e-6;%E级路面的路面功率谱密度down=2048;%E级路面不平度下限up=8192;%E级路面不平度上限白噪声信号属性注意：采样周期为0.001s，采样频率为f=1/0.001=1000Hz,路面频率f1=un1000/2,满足采样定理。采样频率过低，可能引起低频功率下降。时域路面信号生成模型验证y=simout;Fs=1000;%采样频率[pxx,f]=pwelch(y,[],[],[],Fs);%功率谱密度函数n=f/u;%时间频率转换成空间频率pxx=pxx*u;loglog(n,pxx);gridon;holdon;y1=down.*1e-8./n.^2;%注意单位y2=up.*1e-8./n.^2;loglog(n,y1);loglog(n,y2);z=std(y)验证结果可见路面谱大部分落在了E级路面的区域内！！！完！！