三角级数法路面时域模型MATLAB

function Zg= Random_sine（ roadtypeuanlnunnt ）
%随机正弦叠加法
%输入参数：
%roadtype 路面等级
%ua 车速，km/h
%nl 空间频率下限
%nu 空间频率上限
%nn 频率划分区间的总数
%t 总采样时间
%输出参数：
%Zg 路面不平度随时间变化的量
if roadtype==‘A‘||roadtype==‘a‘
   G0=16*10^-6;
elseif roadtype==‘B‘||roadtype==‘b‘
    G0=64*10^-6;
elseif roadtype==‘C‘||roadtype==‘c‘
    G0=256*10^-6;
elseif roadtype==‘D‘||roadtype==‘d‘
    G0=1024*10^-6;
elseif roadtype==‘E‘||roadtype==‘e‘
    G0=4096*10^-6;
elseif