认知无线电网络传输时间设计问题

问题遇到的现象和发生背景

代码解读

问题相关代码，请勿粘贴截图

检测时间可变
clear all
clc;
close all;

fs = 12000; %抽样频率
Ts = 0.001 : 0.001 : 0.1; %检测时间
Sampling_Num1 = fs * Ts; %抽样数量
Sampling_Num = round(Sampling_Num1);

Simulation_Num = 800; %仿真次数
SNR_dB = -5; %信噪比（分贝） %SNR_dB_Vector = ( -20 : 4 : 10 );
Noise_Power = 1; % 噪声功率

for i1 = 1 : length( Ts )
Decision_Index = zeros( 1, Simulation_Num );

SNR = 10^( SNR_dB/10 );
PUser_Signal_Power = SNR * Noise_Power;
Noise = sqrt( Noise_Power ) * randn( 1, Sampling_Num(i1) );                                     %噪声
PUser_Signal = sqrt( PUser_Signal_Power ) * randn( 1, Sampling_Num(i1) );                          %授权用户信号

for i2 = 1:Simulation_Num
    Noise = sqrt( Noise_Power ) * randn( 1, Sampling_Num(i1) );                                     %噪声
    PUser_Signal = sqrt( PUser_Signal_Power ) * randn( 1, Sampling_Num(i1) );                          %授权用户信号
    if i2 <= Simulation_Num/2 
        Received_Signal = PUser_Signal + Noise;                %%%% 用于计算检测概率
    else
        Received_Signal = Noise;                               %%%% 用于计算虚警概率
    end;

    Pd = 0.8;
    Energy_Channel( i2 ) = sum( abs( Received_Signal ) .^ 2 );  %求能量
    Energy_Threshold ( i2 ) = norminv ( 1 - Pd, ( Sampling_Num(i1) + SNR ), sqrt( 2 * ( Sampling_Num(i1) + 2 * SNR) ) );  %高斯分布函数求阈值
    if ( Energy_Channel( i2 ) > Energy_Threshold( i2 ) );
       Decision_Index( i2 ) = 1;
    end;
end;
计算正确检测概率与虚警概率
Dedex( 1 : Simulation_Num/2 );
Decision_Index_2 = Decision_Index( ( Simulation_Num/2 + 1 ) : Simulation_Num );
ED_Num_CorDet = sum( Decision_Index_1 );                            %检测个数
ED_Num_FalseDet = sum( Decision_Index_2 );                          %虚警个数

ED_Prob_CorDet( i1 ) = ED_Num_CorDet / ( Simulation_Num/2 );              %检测概率
ED_Prob_MisDet( i1 ) = 1 -  ED_Prob_CorDet( i1 )    ;          %漏检概率
ED_Prob_FalseDet( i1 ) = ED_Num_FalseDet / ( Simulation_Num/2 );          %虚警概率

end;

hold on
grid on
plot ( Ts, ED_Prob_CorDet, 'b' );hold on
plot ( Ts, ED_Prob_MisDet, '*' );hold on
plot ( Ts, ED_Prob_FalseDet, 'r' );hold on
title('检测时间改变时检测概率\漏检\虚警概率曲线')
xlabel('Ts = 0.00 : 0.001 : 0.01')
h=legend('检测','漏检','虚警',3);

运行结果及报错内容

如何得出此图

我的解答思路和尝试过的方法

我想要达到的结果

得出上图