汽车平顺性仿真分析程序.doc

clear; %清空 clc; clf; q=0; %定义初始变量 j=1; f_l=0.05;%给最小频率赋值 f_u=50; %给最大频率赋值 ns=1000;%将0.5-50HZ的频率分为1000段 delta_f=(f_u-f_l)/ns;%计算频率变化量 Gq_n0=256*1e-6; %选取B级路面不平度等级 n0=0.1; %给参考空间频率赋值 u=40/3.6; %速度换算 md=2*pi*rand(1001,1);%得出1001个随机变量矩阵 for t=0:0.01:100 %计算前轮路面输入 i=1; for n=1:1:1000 f_k=f_l+(n-0.5)*delta_f; Gq_f=Gq_n0*n0^2*u/f_k^2; a_k=sqrt(2*Gq_f*delta_f); q=q+a_k*sin(f_k*2*pi*(t+md(i))); i=i+1; end a(j)=t; %将时间变量赋值给变量a p(j)=q; %将输入变量赋值给p j=j+1; q=0; end figure(1); plot(a,p); %作图 grid on; %给图形加上栅格线 xlim([0,100]); ylim(auto); xlabel(时间(s)); %标注 ylabel(路面不平度位移(m)); title(路面时域随机数据（三角级数法）);%对模拟随机输入数据进行功率谱分析，采用矩形窗 %fs——分析频率,nf——频域的采样点数----------- [Py0,f]=periodogram(p(1:500),hann(500),500,100,onesided); %利用周期图法求前500个点（即前5s）的功率谱 Pya=Py0; for i=2:(2*20-1) [Py0,f]=periodogram(p((i-1)*250+1:((i-1)*250+500)),hann(500),500,100,onesided); %依次每隔250个点（即2.5s）求一次功率谱 Pya=Pya+Py0; end Pya=Pya/(2*20-1); %求平均 figure(2); %建立图形窗口2 % subplot(212); loglog(f(1:500/2),Pya(1:500/2),r); %用对数坐标作图 grid on; xlim([0.1,50]); xlabel(频率(Hz)); %标注 ylabel(路面位移功率谱密度(m^2/Hz)); hold on %将下面的图形一起绘制在figure（2）中 for i=1:500 %计算路面随机输入理论上的功率谱密度 fb(i)=0.1*i; Sq(i)=Gq_n0*0.01*u/(fb(i)^2); end loglog(fb,Sq,-); %采用对数坐标作图 %根据轴距lw和车速V（m/s）生成后轮输入时域数据 lw=6.53; %给轴距赋值 ========= td=lw/u; %计算前后轮时间差 l=round(td/0.01); %对round里面的数取整 pr=p(1:(10000-l)); pf=p((l+1):10000); m1=775; m2=2140; m3=19702; m5=1653; m7=730; m9=51; i4=238000; i6=270; i8=380; k1=1.96E+6; k2=7.84E+6; k3=3.82E+5; k4=1.78E+6; k5=3.14E+6; k6=2.45E+6; k7=1.96E+4; k8=1.96E+4; k9=18121; c3=3.18E+4; c4=1.17E+5; c5=5.00E+3; c6=3.64E+3; c7=7.5E+3; c8=7.3E+3; c9=362*0.75; l1=5.42; l2=(6.53-5.42); l3=(5.42+1.18); l4=(5.42+0.32); l5=(5.42-0.65); l6=(5.42-0.64); a=(1.18-0.54); b=(0.64+0.54); c=(0.21+0.32); d=(0.65-0.21); m=zeros(9,9); m(1,1:9)=m1; m(2,1:9)=m2; m(3,1:9)=m3; m