usa_rail.zip_美国5级谱_美国六级谱_轨道不平顺_轨道不平顺谱_轨道谱

%% 美国轨道功率谱PSD % w=2*pi/L;空间频率定义 最关键的是单位转换问题！ % 参数设置 fs=1000; %采样频率 l=0.5:1/fs:300; %波长范围 m 0.5m-300m W=2*pi./l; Av=0.2095;%粗糙度常数cm2.(rad/m) 5级 Av_1=0.2095*100*(2*pi);%粗糙度常数mm2.(1/m) 5级 Aa=0.0762;%粗糙度常数cm2.(rad/m) 5级 Aa_1=0.0762*100*(2*pi);%粗糙度常数mm2.(1/m) 5级 Wc=0.8245;%rad/m 5级？ Ws=0.8209;%rad/m 5级，，。了，。看，；，；。；看 Ws_1=0.8209*2*pi;%1/m 5级 k=0.25;%安全系数 Lc=2*pi/Wc;% 7.6206m 特征波长1 Ls=2*pi/Ws;% 7.6540m 特征波长2 %美国六级线路 Av2=0.0339; Ws2=0.4380; Ls2=2*pi/Ws2; %单位转换 cm2/(rad/m) -> 100 mm2/(rad/m) 100/(2*pi) mm2/(1/m) %美国五级线路 %以空间频率w=2*pi/L(rad/m)作为x轴自变量 S单位为cm2/(rad/m) S1_W=k*Av*(Wc)^2./(W.^2.*(W.^2+Wc^2)); %高低不平顺 S2_W=k*Aa*(Wc)^2./(W.^2.*(W.^2+Wc^2)); %方向不平顺 S3_W=4*k*Av*Wc^2./((W.^2+Wc^2).*(W.^2+Ws^2)); %水平及轨距不平顺 %转换成统一的标准，以波长L(m)作为x轴自变量 S单位为cm2.(1/m) =100*mm2.(1/m) S1=100*k*Av./(2*pi*(2*pi./(Wc*l)).^2+2*pi).*l.^2; %高低不平顺 S2=100*k*Aa./(2*pi*(2*pi./(Wc*l)).^2+2*pi).*l.^2; %方向不平顺 S3=100*2*k*Av./(pi*l.^2*Lc.^2.*(1./l.^2+1/Lc.^2).*(1./l.^2+1/Ls2.^2)).*l.^2; %水平及轨距不平顺 S1_L=100*k*Av*(Wc)^2./(W.^2.*(W.^2+Wc^2))*2*pi./(l.^2); %高低不平顺 S2_L=100*k*Aa*(Wc)^2./(W.^2.*(W.^2+Wc^2))*2*pi./(l.^2); %方向不平顺 S3_L=100*4*k*Av*Wc^2./((W.^2+Wc^2).*(W.^2+Ws^2))*2*pi./(l.^2); %水平及轨距不平顺 %参数设置 14.3452m S4=100*k*Av2./(2*pi*(2*pi./(Wc*l)).^2+2*pi).*l.^2; %高低不平顺 S5=100*k*Aa./(2*pi*(2*pi./(Wc*l)).^2+2*pi).*l.^2; %方向不平顺 S6=100*2*k*Av2./(pi*l.^2*Lc.^2.*(1./l.^2+1/Lc.^2).*(1./l.^2+1/Ls2.^2)).*l.^2; S4_L=100*k*Av2*(Wc)^2./(W.^2.*(W.^2+Wc^2))*2*pi./(l.^2); %高低不平顺 S5_L=100*k*Aa*(Wc)^2./(W.^2.*(W.^2+Wc^2))*2*pi./(l.^2); %方向不平顺 S6_L=100*4*k*Av2*Wc^2./((W.^2+Wc^2).*(W.^2+Ws^2))*2*pi./(l.^2); %% 画图 figure(1); loglog(l,S1,'r'); %绘制波长与功率谱密度的双对数坐标图h title('美国5级高低不平顺'); %题名，并设置字号为25 ylabel('功率谱密度/mm2.m'); %y轴名称 grid on %加网格 xlabel('波长/m'); %x轴名称 set(gca,'xdir','reverse'); figure(2) loglog(l,S2,'r'); title('美国5级方向不平顺'); ylabel('功率谱密度/mm2.m'); grid on xlabel('波长/m'); set(gca,'xdir','reverse'); figure(3) loglog(l,S6,'r'); title('美国6级水平及轨距不平顺'); ylabel('功率谱密度/mm2.m'); grid on xlabel('波长/m'); set(gca,'xdir','reverse');