%%
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