function?standard_gai

%?三角級數法求路面時域模型
clear?all;clc;
%??輸入仿真參數?
in=file_name;?
[roadtypeuaflfudftltudt]=file_parameters（in）;
%??確定仿真需要的其它主要參數?
[ftunfreqnt]=freq_parameters（flfudftltudtua）;?
%??8等級路面
Gqf=roadpsd_f（roadtypefu）;?
%??生成隨機數列
theta=suiji（nfreq）;?
%??路面時域模型
qt=shiyu（dfGqfthetafnfreqtnt）;?


%?將三角級數法所得路面功率譜密度與標準路面進行對比
figure
Fs=200;
NFFT=500;
ww=hamming（NFFT）;
%三角級數法所得路面的功率譜密度
[Gfff]=pwelch（qtww[]NFFTFs）;
loglog（ffGf‘k:‘fGqf‘k-‘）
xlabel（‘時間頻率?f/Hz‘）;?
ylabel（‘路面功率譜密度Gq（f）‘）;?
legend（‘模擬值‘‘理論值‘）
legend?boxoff
box?off

%?路面均方根值的對比
%?理論值的均方根值
i=1:nfreq;?
dGqf（i）=Gqf（i）.*df;
rms1=sum（dGqf）;?
evp1=sqrt（rms1）;
%?模擬值的均方根值
j=1:round（NFFT/2）;
dff=2*max（ff）/NFFT;
dGf（j）=Gf（j）.*dff;
rms2=sum（dGf）;?
evp2=sqrt（rms2）;


function?in=file_name?

%?功能：獲得原始數據文件的名稱并建立相應的通道號?
%?輸出參數：?
%?in??輸入文件通道號?

clear?
infile=input（?‘請輸入原始數據文件的名稱（包括后綴）:?‘‘s‘）;?
in=fopen（infile‘r‘）;?
if?in<0?
????disp（‘原始數據文件不存在，無法打開‘）;?
end?


function?[roadtypeuaflfudftltudt]=file_parameters（unit）?
%?功能：從文件輸入所需參數?
%?輸入參數：?
%?unit?????輸入文件通道號?
%?輸出參數：?
%?roadtype?路面類型?
%?ua???????車速km/h?
%?fl???????時間頻率的下限Hz?
%?fu???????時間頻率的上限Hz?
%?df???????時間頻率的增量Hz?
%?tl???????時間的下限s?
%?tu???????時間的上限s
%?dt???????時間的增量s?
%?在此注意：dt<=1/（2*fu）??nt>=2*fu*T
roadtype=fscanf（unit‘%s‘1）;?
ua=fscanf（unit‘%f‘1）;?
fl=fscanf（unit‘%f‘1）;?
fu=fscanf（unit‘%f‘1）;?
df=fscanf（unit‘%f‘1）;?
tl=fscanf（unit‘%f‘1）;?
tu=fscanf（unit‘%f‘1）;?
dt=fscanf（unit‘%f‘1）;?


function?[ftunfreqnt]=freq_parameters（flfudftltudtua）?

%??