``
clc;
c = 2.997924562e8; %Speed of Light 光速
dtor = pi/180; %5角度转换为弧度
% GrAng=app.GrAngEditField.Value;
GrAng=app.GrAngSlider.Value;
app.currentNumEditField.Value=GrAng;
phi=app.phiSpinner.Value;
ThWind =app.ThWindSpinner.Value;
Pol=app.PolDropDown.Value;

            % GrAng =[0.3];
            NGrAng = length(GrAng);%定义擦射角的数组的长度
            % ThWind = 0;
            % phi = 0;
            % freq = (0.5:0.1:40);
            % freq = app.fromEditField.Value;
            fx=app.fromSpinner.Value;
            step=app.stepSpinner.Value;
            fm=app.ToSpinner.Value;
            freq = (fx:step:fm);%fx为起始值，step为步长，fm为序列的结束值
            fMin=freq(1);
            fMax=max(freq);
            fMax=freq(end);

            % Variables for computing average abs deviation 计算平均误差
            SumDevH = 0.0;
            NValH = 0;
            SumDevV = 0.0;
            NValV = 0;

            phi_rad = phi .* dtor;
            % GRAPHICS DEFINITIONS 图形定义
            ha=[];hl=[];hp=[];ht=[];htx=[];hty=[];
            set(0,'Units','pixels')


            % Set SCreen Size 定义屏幕尺寸
            scnsize=get(0,'ScreenSize');
            pos1 = [5+scnsize(1),0.02*scnsize(4),.99*scnsize(3),.75*scnsize(4)];

            FigFont = 28;
            CircSize = 14;
            LineSize = 7;
            AxWidth = 4;

            % Initialize figure number 初始化图号
            
            % Defines units for values returned by ScreenSize 定义ScreenSize返回值的单位

            % Define colors for curves and points 定义曲线和点的颜色
            Cols = 'bgrcmyk';

            ColsMat = [ 0.161 0.698 0.725; % magenta
                0.608 0.18 0.8; % violet
                1.0 0.0 0.0; % red
                0.541 0.416 0.361; % brown
                0.0 0.0 0.0; % black
                1 0 1;% dark green
                0.067 0.204 0.898 ];%  dark blue
            cla(app.UIAxes14, 'reset');
            cla(app.UIAxes_4, 'reset');
            cla(app.UIAxes14_2, 'reset');
            cla(app.UIAxes14_3, 'reset');
            cla(app.UIAxes14_4, 'reset');
            %########################################################################
            % ALL GRAZING ANGLES 所有的擦射角
                % Do all sea states for HH  H极化下的海况。
                %% GIT
                for SeaSt=0:6
                    SeaState=SeaSt;
                    SS_Plus1 = SeaState + 1;

                    for ifreq = 1 : length(freq)
                        GTI_SigZHor(ifreq) = GTI_SigmaSea(freq(ifreq),SeaState,Pol,alpha,ThWind);
                    end

                    GTI_SigmaHH = GTI_SigZHor;
            
                    % Plot empirical curve first
                    semilogx(app.UIAxes14,freq,GTI_SigmaHH,'Color',ColsMat(SeaSt+1,:));
                    legend(app.UIAxes14,'GrAng =0.3 SeaState=1','GrAng =0.3 SeaState=2','GrAng =0.3 SeaState=3','GrAng =0.3 SeaState=4','GrAng =0.3 SeaState=5','GrAng =0.3 SeaState=6',Location='southeast')
                    hold(app.UIAxes14,"on")
                    grid(app.UIAxes14,'on');
                    xlabel(app.UIAxes14,'频率(GHz) ','FontSize',14);
                    ylabel(app.UIAxes14,'散射系数 - \sigma_{HH}^o (dB)','FontSize',14);
                    title(app.UIAxes14,'GIT模型HH极化','FontSize',14)
                end
                %% NRL
                for SeaSt=0:6
                    SeaState=SeaSt;
                    SS_Plus1 = SeaState + 1;

                    NRL_SigZHor = NRL_SigmaSea(freq,SeaState,Pol,alpha);
                    NRL_SigmaHH=NRL_SigZHor;
                    
                    % Plot empirical curve first

                    hl(end+1) = semilogx(app.UIAxes14_2,freq,NRL_SigmaHH);
                    legend(app.UIAxes14_2,'GrAng =0.3 SeaState=1','GrAng =0.3 SeaState=2','GrAng =0.3 SeaState=3','GrAng =0.3 SeaState=4','GrAng =0.3 SeaState=5',Location='southeast')
                    hold(app.UIAxes14_2,"on")
                    xlabel(app.UIAxes14_2,'频率(GHz) ','FontSize',14);
                    ylabel(app.UIAxes14_2,'散射系数 - \sigma_{HH}^o (dB)','FontSize',14);
                    title(app.UIAxes14_2,'NRL模型HH极化','FontSize',14)
                    grid(app.UIAxes14_2,'on');
                end
                %% HYB
                for SeaSt=0:6
                    SeaState=SeaSt;
                    SS_Plus1 = SeaState + 1;

                    for ifreq = 1 : length(freq)
                        HYB_SigZHor(ifreq) = HYB_SigmaSea(freq(ifreq),SeaState,Pol,alpha,ThWind);
                    end
                    HYB_SigmaHH=HYB_SigZHor;
                    % Plot empirical curve first
                    hl(end+1) = semilogx(app.UIAxes14_3,freq,HYB_SigmaHH);
                    legend(app.UIAxes14_3,'GrAng =0.3 SeaState=1','GrAng =0.3 SeaState=2','GrAng =0.3 SeaState=3','GrAng =0.3 SeaState=4','GrAng =0.3 SeaState=5',Location='southeast')
                    hold(app.UIAxes14_3,"on")
                    xlabel(app.UIAxes14_3,'频率(GHz) ','FontSize',14);
                    ylabel(app.UIAxes14_3,'散射系数 - \sigma_{HH}^o (dB)','FontSize',14);
                    title(app.UIAxes14_3,'HYB模型HH极化','FontSize',14)
                    grid(app.UIAxes14_3,'on');
                end
                %% TSC
                    for SeaSt=0:6
                    SeaState=SeaSt;
                    SS_Plus1 = SeaState + 1;

                    for ifreq = 1 : length(freq)
                        TSC_SigZHor(ifreq) = TSC_SigmaSea(freq(ifreq),SeaState,Pol,alpha,ThWind);
                    end
                    TSC_SigmaHH=TSC_SigZHor;
                    % Plot empirical curve first
                    hl(end+1) = semilogx(app.UIAxes14_4,freq,TSC_SigmaHH);
                    legend(app.UIAxes14_4,'GrAng =0.3 SeaState=1','GrAng =0.3 SeaState=2','GrAng =0.3 SeaState=3','GrAng =0.3 SeaState=4','GrAng =0.3 SeaState=5',Location='southeast')
                    hold(app.UIAxes14_4,"on")
                    % set(hl(end),'Color',ColsMat(SeaSt+1,:));
                    xlabel(app.UIAxes14_4,'频率(GHz) ','FontSize',14);
                    ylabel(app.UIAxes14_4,'散射系数 - \sigma_{HH}^o (dB)','FontSize',14);
                    title(app.UIAxes14_4,'TSC模型HH极化','FontSize',14)
                    grid(app.UIAxes14_4,'on');
                end

``