function?[Roadelevationtimeinterval]=inverse_FFT_time_domain_standardroad（roadtypeTimelengthSampletimeVehiclespeed）
%?功能：基于逆虛擬激勵方法生成路面不平度時域模型（路面激勵）
%?輸入參數：
%?Roadlevel?路面等級
%?Timelength?路面仿真總時間，單位s
%?Sampletime?路面采樣時間間隔，單位?s
%?Vehiclespeed?車速，單位km/h
%?輸出參數：
%?Roadelevation?路面高程，單位為?m
%?timeinterval?時間間隔，單位為?s



n0=.1;%參考空間頻率
w=2;%設置頻率指數

%?（1）確定路面類型
if?roadtype?==‘A‘?|?roadtype?==‘a‘?
????Gqn0=16*10^（-6）;?
elseif?roadtype?==‘B‘?|?roadtype?==‘b‘?
????Gqn0=64*10^（-6）;?
elseif?roadtype?==‘C‘?|?roadtype?==‘c‘?
????Gqn0=256*10^（-6）;?
elseif?roadtype?==‘D‘?|?roadtype?==‘d‘?
????Gqn0=1024*10^（-6）;?
elseif?roadtype?==‘E‘?|?roadtype?==‘e‘?
????Gqn0=4096*10^（-6）;?
elseif?roadtype?==‘F‘?|?roadtype?==‘f‘?
????Gqn0=16384*10^（-6）;?
elseif?roadtype?==‘G‘?|?roadtype?==‘g‘
?Gqn0=65536*10^（-6）;?
elseif?roadtype?==‘H‘?|?roadtype?==‘h‘?
????Gqn0=262144*10^（-6）;?????
end?
%?轉化車速為speedu單位m/s
speedu=Vehiclespeed/3.6;


%?給定時間頻率上下限
%?timefreq1?時間頻率上限，單位Hz
%?timefreq2?時間頻率下限，單位Hz
timefreq1=0.25;
timefreq2=50;
%?%?（2）確定采樣時間間隔Sampletime和總的采樣時間T的合理性
%??if?Timelength>=1/timefreq1
%??????Timelength1=Timelength;
%??else?
%?%??????仿真時間過短，請重新輸入更長的仿真時間
%????