傅里叶逆变换方法生成路面时域模型

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 
% %      仿真时间过短，请重新输入更长的仿真时间
%