路面谱密度matlab代码

包括标准路面谱频域和时域代码，以及基于Simulink搭建的滤波白噪声模型，并对采集的数据进行拟合得到的路面谱频域和时域，并与标准对比的代码，还对matlab与excel数据导入导出作部分介绍

%%
x1=xlsread（‘路面随机激励.xls‘‘A1:A10001‘）;
x=x1‘;
dt=0.001;
Nt=length（x）;     %  采样点（可能要修改）
V=30;             %  km/h
T=（Nt-1）*dt;      %  采样周期
t=0:dt:T;         %  采样时间点
v=V*1000/3600;    %  m/s
df=1/（Nt*dt）;     %  采样频率间隔
f=0:df:1/（2*dt）;  %  采用频率一半（可修改）
nf=length（f）;     %  采样点的个数
F=f/v             %  采样空间频率
for i=1:nf
	if 0.001<=F（i） & F（i）<=100  
		P（i）=16*10^-6*（0.1/F（i））^2;   %  获得标准空间路面谱
	else F（i）==0
		P（i）=0;
	end
end
p=P/v;            %  转化为频谱
figure（1）
loglog（FP）;      %  1输出标准路面空间频率谱G（n）
xlabel（‘空间频率n[1/m]‘）ylabel（‘路面不平度位移功率谱密度Gq（n）[m^3]‘）
figure（2）
loglog（fp）xlabel（‘时间频率[Hz]‘）ylabel（‘G（f）[m^2/Hz]‘） %  2输出时间频率谱G（f）
figure（3）
plot（tx）xlabel（‘Time [sec.]‘）ylabel（‘路面随机激励Z [m]‘）% 3输出路面随机激励 
xlim（[0 T]）;
nfft=2^nextpow2（Nt）;  % 取最接近Nt的2的n次幂
Fs=1000;
noverlap=20;          % noverlap是分段平均时重合的部分
[pxxfsd]=pwelch（x[]noverlapnfftFs/2）; % 作傅里叶逆变换
n=fsd/v;% 化为空间频率n
Pxx=pxx*v;%化成类似Gq（n）
figure（4）
loglog（fpfsdpxx）;% 4输出标准路面时间频率功率谱和逆变换的时间频率功率谱
xlabel（‘Frequency[Hz]‘）ylabel（‘PSD[m^2/Hz]‘）
legend（‘base PSD‘‘PSD of Road Surface‘）xlim（[1e-1 1e+2]）
figure（5）
loglog（FPnPxx）;  % 5输出标准路面空间频率功率谱和逆变换的空间频率功率谱
xlabel（‘空间频率n[1/m]‘）ylabel（‘PSD[m^3]‘）
legend（‘base PSD‘‘PSD of Road Surface‘）xlim（[1e-1 1e+2]）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      14282  2017-05-17 22:01  路面谱\csshuju.slx

     文件       1484  2017-05-18 10:30  路面谱\lumianpu.m

     文件      12468  2017-05-18 10:52  路面谱\matlab数据导出和导入.docx

     文件        876  2017-05-15 14:27  路面谱\Untitled2.m

     文件        120  2017-05-18 10:37  路面谱\参数介绍.txt

     文件     431616  2017-05-17 22:01  路面谱\路面随机激励.xls

     目录          0  2017-05-19 21:09  路面谱

----------- ---------  ---------- -----  ----

               460846                    7