时间频域内的等级路面不平度含有上下限MATLAB

function [Gqf_highGqfGqf_down]=roadpsd_f_highdown（roadtypenua）
%计算8种路面激励的统计特性上下限
%输入参数：
%roadtype 路面类型
%n 空间频率
%输出参数：
%Gqf_high 位移功率谱密度向量上限
%Gqf_down 位移功率谱密度向量下限
if roadtype==‘A‘||roadtype==‘a‘
   Gqn0_high=32*10^-6;
   Gqn0_down=0;
   Gqn0=16*10^-6;
elseif roadtype==‘B‘||roadtype==‘b‘
    Gqn0_high=128*10^-6;
    Gqn0_down=32*10^-6;
    Gqn0=64*10^-6;
elseif roadtype==‘C‘||roadtype==‘c‘
    Gqn0_high=512*10^-6;
    Gqn0_down=128*10^-6;
    Gqn0=256*10^-6;
elseif roadtype==‘D‘||roadtype==‘d‘
    Gqn0_high=2048*10^-6;
    Gqn0_down=512*10^-6;
    Gqn0=1024*10^-6;
elseif roadtype==‘E‘||r