From da829edda095ed0f00779701bc8db933f07f9263 Mon Sep 17 00:00:00 2001 From: Macrad <2244695967@qq.com> Date: Tue, 27 Apr 2021 23:42:50 +0800 Subject: [PATCH 1/3] =?UTF-8?q?CFAR=E5=A2=9E=E5=8A=A0=E6=B3=A8=E9=87=8A?= =?UTF-8?q?=E5=92=8C=E8=BE=93=E5=85=A5=E5=8A=9F=E8=83=BD=EF=BC=8C=E4=BF=AE?= =?UTF-8?q?=E6=94=B9=E7=AC=AC=E4=BA=8C=E6=AE=B5=E7=A8=8B=E5=BA=8F=E8=AE=A9?= =?UTF-8?q?=E7=9C=9F=E5=AE=9E=E7=9B=AE=E6=A0=87=E6=95=B0=E7=9B=AE=E5=8F=AF?= =?UTF-8?q?=E8=A7=86=E5=8C=96=EF=BC=8CV1.1=E7=89=88=E6=9C=AC?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- EchoProcess/CFAR-V1.1.m | 158 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 158 insertions(+) create mode 100644 EchoProcess/CFAR-V1.1.m diff --git a/EchoProcess/CFAR-V1.1.m b/EchoProcess/CFAR-V1.1.m new file mode 100644 index 0000000..5d254e6 --- /dev/null +++ b/EchoProcess/CFAR-V1.1.m @@ -0,0 +1,158 @@ +%% Target Detection Pd Pfa,CA-CFAR +close all;clear all;clc; +N = 1e4; %Points +n = 1:N; +SNR = 10; %SNR in dB +Pfa = 1e-3; +sigma_dB = 0; %Noise power in dB +sigma = db2pow(sigma_dB); %Noise Power=1w + +noise = sigma/sqrt(2)*(randn(1,N)+1i*randn(1,N)); +A = sqrt(sigma^2 * 10^(SNR/10)); %Target signal +T_fix = sqrt(-sigma^2*log(Pfa)); %fix thresholding 只与噪声功率有关 假设噪声功率为常数 +c = 500; %one side,ref cell = c-guardc +guardc = 10; +CUTIdx = N/2; %目标的个数 +%Simulation of target detection% +sample = rand(1,N) < 0.001; %可随机生成几个目标 对应位置返回1 +A0 = A .* sample; +%A1 = abs(A0 + noise.*sample); +signal1 = A0 + noise; %Receive signal + +figure(1); +idx_target = find(A0 > 0); %真实目标位置 +idx_detect = find(abs(signal1) > T_fix); %通过阈值检测到的目标 +plot(n,abs(signal1),'-o','MarkerIndices',idx_detect,'MarkerSize',10); +hold on; plot(n(idx_target),abs(signal1(1,idx_target)),'p','MarkerEdgeColor','r','MarkerFaceColor','r'); +hold on; plot(n,T_fix*ones(1,N)); +xlim([c+2,N-c-1]); + +%Evaluation of Pd & Pfa % +signal2 = A + noise; +disp('Detection Rate：') +Pd_est = sum(abs(signal2) > T_fix) / N +disp('False alarm Rate：') +Pfa_est = sum(abs(noise(1,c+2:N-c-1)) > T_fix) / N + +%Simulation of CFAR detection +%Ncells = 23;%Ntrails =1e4; + + +%signal = noise; +sample = rand(1,N) < 0.001; %可随机生成几个目标 +signal = A.*sample; +signal = signal+noise; +%signal(1,CUTIdx) = signal(1,CUTIdx) + A; % 5000个目标 +k = Pfa^(-1/(2*(c-guardc)))-1; %门限系数公式（在均值为0下的独立高斯噪声背景下） + +T_cfar = T_fix .* ones(1,N); +for j=c+1:N-c + pn = sum(abs(signal(1,(j-c:j-guardc-1))).^2) + sum(abs(signal(1,(j+guardc+1:j+c))).^2); + T_cfar(j) =k.*abs(pn); +end +figure(2); +idx2_target = find(sample > 0); %真实目标位置 +idxC = find(abs(signal).^2 > T_cfar); %大于检测阈值的目标 +plot(n,abs(signal).^2,'-o','MarkerIndices',idxC,'MarkerSize',10); +hold on; plot(n(idx2_target),abs(signal(1,idx2_target)).^2,'p','MarkerEdgeColor','r','MarkerFaceColor','r'); +%hold on;stem(N/2,abs(signal(1,N/2)).^2,'-p','MarkerEdgeColor','r','MarkerFaceColor','r'); +hold on;plot(n,T_cfar); +xlim([c+2,N-c-1]); +Pfa_cfar = sum(abs(noise(1,c+2:N-c-1)).^2 > T_cfar(1,c+2:N-c-1))./(N-2*c) +% %Pfa_np = sum(x(CUTIdx,:)>T_ideal)/Ntrials + +%% Pd +N_trial = 10000; +sum_fix = 0; +sum_cfar = 0; +T_fix = sqrt(-sigma^2*log(Pfa)); +for n = 1:N_trial + noise = sigma/sqrt(2)*(randn(1,N)+1i*randn(1,N)); + signal = noise; + signal(1,CUTIdx) = signal(1,CUTIdx) + A; %5000个目标 + sum_fix = sum_fix + (abs(signal(1,CUTIdx)) > T_fix); + + pn = sum(abs(signal(1,(CUTIdx-c:CUTIdx-guardc-1))).^2) + sum(abs(signal(1,(CUTIdx+guardc+1:CUTIdx+c))).^2); + T_cfar0 =k.*abs(pn); + sum_cfar = sum_cfar + (abs(signal(1,CUTIdx)).^2 > T_cfar0); +end +pd_fix = sum_fix / N_trial %阈值固定下的恒虚警率 +pd_cfar = sum_cfar / N_trial + +%% CFAR COMPARE +close all;clear all;clc; +N = 6e3; +n = 1:N; +SNR = 20; +Pfa = 1e-6; +sigma_dB = 0; +sigma = db2pow(sigma_dB); %Noise Power=1w +ndeta = 100; +scale = 5; +noise = sigma*[randn(1,N/2-ndeta), scale*randn(1,N/2+ndeta)]; +A = sqrt((scale*sigma)^2 * 10^(SNR/10)); %Target signal amplitude +T_fix = sqrt(-sigma^2*log(Pfa)); %fix thresholding +c = N/40; %one side,ref cell = c-guardc +guardc = 10; %guard cell number +CUTIdx = N/2; %position of target + +%CFAR mode choose%%%%%,0 CA;1 GO;2 SO;3:OS +disp("请输入你接下来的操作： 0代表 CA-CFAR") +disp(" 1代表 GO-CFAR") +disp(" 2代表 SO-CFAR") +disp(" 3代表 OS-CFAR") +disp(" 4代表 退出 ") +flag=input("请输入："); +while(flag==0||flag==1||flag==2||flag==3) + cfar=flag; + %%Simulation of CFAR detection + %signal model + signal = noise; + signal(1,CUTIdx) = signal(1,CUTIdx) + A; %5000个目标 + %k = Pfa^(-1/(2*(c-guardc)))-1; + T_cfar = ones(1,N); + ncell = 0; + pn = 0; + for j=c+1:N-c + a = sum(abs(signal(1,(j-c:j-guardc-1))).^2); + b = sum(abs(signal(1,(j+guardc+1:j+c))).^2); + if cfar==0 + pn = a + b; + ncell = 2*(c-guardc); + end + if cfar==1 + pn = max(a,b); + ncell = c-guardc; + end + if cfar==2 + pn = min(a,b); + ncell = c-guardc; + end + k = Pfa^(-1/ncell)-1; + T_cfar(j) =k.*abs(pn); + end + n_target = 5; %target numbers + if cfar==3 + s = sort(signal); + T_cfar = ones(1,N) .* s(1,N-n_target).^2; + end + + %traget detected + idxC = find(abs(signal).^2 > T_cfar); + + figure; + plot(n,10*log(abs(signal).^2),'-o','MarkerIndices',idxC,'MarkerSize',10); + hold on;stem(N/2,10*log(abs(signal(1,N/2)).^2),'-p','MarkerEdgeColor','r','MarkerFaceColor','r'); + hold on;plot(n,10*log(T_cfar)); + xlim([c+2,N-c-1]);ylim([0,100]); + Pfa_cfar = sum(abs(noise(1,c+2:N-c-1)).^2 > T_cfar(1,c+2:N-c-1))./(N-2*c); + fprintf('虚警率:%s\n',num2str(Pfa_cfar)); + disp("当前程序执行完毕") + disp("请输入你接下来的操作： 0代表 CA-CFAR") + disp(" 1代表 GO-CFAR") + disp(" 2代表 SO-CFAR") + disp(" 3代表 OS-CFAR") + disp(" t代表 退出 ") + flag=input('请输入：'); +end +disp("已退出程序") \ No newline at end of file From 8be87cea4152f905cc3b21a6b7e741af23b29d3a Mon Sep 17 00:00:00 2001 From: Macrad <2244695967@qq.com> Date: Wed, 28 Apr 2021 14:25:33 +0800 Subject: [PATCH 2/3] =?UTF-8?q?=E4=BF=AE=E6=94=B9CFAR=E6=96=87=E4=BB=B6?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- EchoProcess/CFARV2.m | 158 +++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 158 insertions(+) create mode 100644 EchoProcess/CFARV2.m diff --git a/EchoProcess/CFARV2.m b/EchoProcess/CFARV2.m new file mode 100644 index 0000000..5e65da8 --- /dev/null +++ b/EchoProcess/CFARV2.m @@ -0,0 +1,158 @@ +%% Target Detection Pd Pfa,CA-CFAR +close all;clear all;clc; +N = 1e4; %Points +n = 1:N; +SNR = 10; %SNR in dB +Pfa = 1e-3; +sigma_dB = 0; %Noise power in dB +sigma = db2pow(sigma_dB); %Noise Power=1w + +noise = sigma/sqrt(2)*(randn(1,N)+1i*randn(1,N)); +A = sqrt(sigma^2 * 10^(SNR/10)); %Target signal +T_fix = sqrt(-sigma^2*log(Pfa)); %fix thresholding 只与噪声功率有关 假设噪声功率为常数 +c = 500; %one side,ref cell = c-guardc +guardc = 10; +CUTIdx = N/2; %目标的个数 +%Simulation of target detection% +sample = rand(1,N) < 0.001; %可随机生成几个目标 对应位置返回1 +A0 = A .* sample; +%A1 = abs(A0 + noise.*sample); +signal1 = A0 + noise; %Receive signal + +figure(1); +idx_target = find(A0 > 0); %真实目标位置 +idx_detect = find(abs(signal1) > T_fix); %通过阈值检测到的目标 +plot(n,abs(signal1),'-o','MarkerIndices',idx_detect,'MarkerSize',10); +hold on; plot(n(idx_target),abs(signal1(1,idx_target)),'p','MarkerEdgeColor','r','MarkerFaceColor','r'); +hold on; plot(n,T_fix*ones(1,N)); +xlim([c+2,N-c-1]); + +%Evaluation of Pd & Pfa % +signal2 = A + noise; +disp('Detection Rate：') +Pd_est = sum(abs(signal2) > T_fix) / N +disp('False alarm Rate：') +Pfa_est = sum(abs(noise(1,c+2:N-c-1)) > T_fix) / N + +%Simulation of CFAR detection +%Ncells = 23;%Ntrails =1e4; + + +%signal = noise; +sample = rand(1,N) < 0.001; %可随机生成几个目标 +signal = A.*sample; +signal = signal+noise; +%signal(1,CUTIdx) = signal(1,CUTIdx) + A; % 5000个目标 +k = Pfa^(-1/(2*(c-guardc)))-1; %门限系数公式（在均值为0下的独立高斯噪声背景下） + +T_cfar = T_fix .* ones(1,N); +for j=c+1:N-c + pn = sum(abs(signal(1,(j-c:j-guardc-1))).^2) + sum(abs(signal(1,(j+guardc+1:j+c))).^2); + T_cfar(j) =k.*abs(pn); +end +figure(2); +idx2_target = find(sample > 0); %真实目标位置 +idxC = find(abs(signal).^2 > T_cfar); %大于检测阈值的目标 +plot(n,abs(signal).^2,'-o','MarkerIndices',idxC,'MarkerSize',10); +hold on; plot(n(idx2_target),abs(signal(1,idx2_target)).^2,'p','MarkerEdgeColor','r','MarkerFaceColor','r'); +%hold on;stem(N/2,abs(signal(1,N/2)).^2,'-p','MarkerEdgeColor','r','MarkerFaceColor','r'); +hold on;plot(n,T_cfar); +xlim([c+2,N-c-1]); +Pfa_cfar = sum(abs(noise(1,c+2:N-c-1)).^2 > T_cfar(1,c+2:N-c-1))./(N-2*c) +% %Pfa_np = sum(x(CUTIdx,:)>T_ideal)/Ntrials + +%% Pd +N_trial = 10000; +sum_fix = 0; +sum_cfar = 0; +T_fix = sqrt(-sigma^2*log(Pfa)); +for n = 1:N_trial + noise = sigma/sqrt(2)*(randn(1,N)+1i*randn(1,N)); + signal = noise; + signal(1,CUTIdx) = signal(1,CUTIdx) + A; %5000个目标 + sum_fix = sum_fix + (abs(signal(1,CUTIdx)) > T_fix); + + pn = sum(abs(signal(1,(CUTIdx-c:CUTIdx-guardc-1))).^2) + sum(abs(signal(1,(CUTIdx+guardc+1:CUTIdx+c))).^2); + T_cfar0 =k.*abs(pn); + sum_cfar = sum_cfar + (abs(signal(1,CUTIdx)).^2 > T_cfar0); +end +pd_fix = sum_fix / N_trial %阈值固定下的恒虚警率 +pd_cfar = sum_cfar / N_trial + +%% CFAR COMPARE +close all;clear all;clc; +N = 6e3; +n = 1:N; +SNR = 20; +Pfa = 1e-6; +sigma_dB = 0; +sigma = db2pow(sigma_dB); %Noise Power=1w +ndeta = 100; +scale = 5; +noise = sigma*[randn(1,N/2-ndeta), scale*randn(1,N/2+ndeta)]; +A = sqrt((scale*sigma)^2 * 10^(SNR/10)); %Target signal amplitude +T_fix = sqrt(-sigma^2*log(Pfa)); %fix thresholding +c = N/40; %one side,ref cell = c-guardc +guardc = 10; %guard cell number +CUTIdx = N/2; %position of target + +%CFAR mode choose%%%%%,0 CA;1 GO;2 SO;3:OS +disp("请输入你接下来的操作： 0代表 CA-CFAR") +disp(" 1代表 GO-CFAR") +disp(" 2代表 SO-CFAR") +disp(" 3代表 OS-CFAR") +disp(" 4代表 退出 ") +flag=input("请输入："); +while(flag==0||flag==1||flag==2||flag==3) + cfar=flag; + %%Simulation of CFAR detection + %signal model + signal = noise; + signal(1,CUTIdx) = signal(1,CUTIdx) + A; %5000个目标 + %k = Pfa^(-1/(2*(c-guardc)))-1; + T_cfar = ones(1,N); + ncell = 0; + pn = 0; + for j=c+1:N-c + a = sum(abs(signal(1,(j-c:j-guardc-1))).^2); + b = sum(abs(signal(1,(j+guardc+1:j+c))).^2); + if cfar==0 + pn = a + b; + ncell = 2*(c-guardc); + end + if cfar==1 + pn = max(a,b); + ncell = c-guardc; + end + if cfar==2 + pn = min(a,b); + ncell = c-guardc; + end + k = Pfa^(-1/ncell)-1; + T_cfar(j) =k.*abs(pn); + end + n_target = 5; %target numbers + if cfar==3 + s = sort(signal); + T_cfar = ones(1,N) .* s(1,N-n_target).^2; + end + + %traget detected + idxC = find(abs(signal).^2 > T_cfar); + + figure; + plot(n,10*log(abs(signal).^2),'-o','MarkerIndices',idxC,'MarkerSize',10); + hold on;stem(N/2,10*log(abs(signal(1,N/2)).^2),'-p','MarkerEdgeColor','r','MarkerFaceColor','r'); + hold on;plot(n,10*log(T_cfar)); + xlim([c+2,N-c-1]);ylim([0,100]); + Pfa_cfar = sum(abs(noise(1,c+2:N-c-1)).^2 > T_cfar(1,c+2:N-c-1))./(N-2*c); + fprintf('虚警率:%s\n',num2str(Pfa_cfar)); + disp("当前程序执行完毕") + disp("请输入你接下来的操作： 0代表 CA-CFAR") + disp(" 1代表 GO-CFAR") + disp(" 2代表 SO-CFAR") + disp(" 3代表 OS-CFAR") + disp(" 4代表 退出 ") + flag=input('请输入：'); +end +disp("已退出程序") \ No newline at end of file From f76f0d5a6f468b7306ea933dcb51811ae8b2c920 Mon Sep 17 00:00:00 2001 From: Macrad <2244695967@qq.com> Date: Wed, 28 Apr 2021 14:36:21 +0800 Subject: [PATCH 3/3] =?UTF-8?q?=E5=88=A0=E9=99=A4=E9=94=99=E8=AF=AF?= =?UTF-8?q?=E6=8F=90=E4=BA=A4?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- EchoProcess/CFAR-V1.1.m | 158 ---------------------------------------- 1 file changed, 158 deletions(-) delete mode 100644 EchoProcess/CFAR-V1.1.m diff --git a/EchoProcess/CFAR-V1.1.m b/EchoProcess/CFAR-V1.1.m deleted file mode 100644 index 5d254e6..0000000 --- a/EchoProcess/CFAR-V1.1.m +++ /dev/null @@ -1,158 +0,0 @@ -%% Target Detection Pd Pfa,CA-CFAR -close all;clear all;clc; -N = 1e4; %Points -n = 1:N; -SNR = 10; %SNR in dB -Pfa = 1e-3; -sigma_dB = 0; %Noise power in dB -sigma = db2pow(sigma_dB); %Noise Power=1w - -noise = sigma/sqrt(2)*(randn(1,N)+1i*randn(1,N)); -A = sqrt(sigma^2 * 10^(SNR/10)); %Target signal -T_fix = sqrt(-sigma^2*log(Pfa)); %fix thresholding 只与噪声功率有关 假设噪声功率为常数 -c = 500; %one side,ref cell = c-guardc -guardc = 10; -CUTIdx = N/2; %目标的个数 -%Simulation of target detection% -sample = rand(1,N) < 0.001; %可随机生成几个目标 对应位置返回1 -A0 = A .* sample; -%A1 = abs(A0 + noise.*sample); -signal1 = A0 + noise; %Receive signal - -figure(1); -idx_target = find(A0 > 0); %真实目标位置 -idx_detect = find(abs(signal1) > T_fix); %通过阈值检测到的目标 -plot(n,abs(signal1),'-o','MarkerIndices',idx_detect,'MarkerSize',10); -hold on; plot(n(idx_target),abs(signal1(1,idx_target)),'p','MarkerEdgeColor','r','MarkerFaceColor','r'); -hold on; plot(n,T_fix*ones(1,N)); -xlim([c+2,N-c-1]); - -%Evaluation of Pd & Pfa % -signal2 = A + noise; -disp('Detection Rate：') -Pd_est = sum(abs(signal2) > T_fix) / N -disp('False alarm Rate：') -Pfa_est = sum(abs(noise(1,c+2:N-c-1)) > T_fix) / N - -%Simulation of CFAR detection -%Ncells = 23;%Ntrails =1e4; - - -%signal = noise; -sample = rand(1,N) < 0.001; %可随机生成几个目标 -signal = A.*sample; -signal = signal+noise; -%signal(1,CUTIdx) = signal(1,CUTIdx) + A; % 5000个目标 -k = Pfa^(-1/(2*(c-guardc)))-1; %门限系数公式（在均值为0下的独立高斯噪声背景下） - -T_cfar = T_fix .* ones(1,N); -for j=c+1:N-c - pn = sum(abs(signal(1,(j-c:j-guardc-1))).^2) + sum(abs(signal(1,(j+guardc+1:j+c))).^2); - T_cfar(j) =k.*abs(pn); -end -figure(2); -idx2_target = find(sample > 0); %真实目标位置 -idxC = find(abs(signal).^2 > T_cfar); %大于检测阈值的目标 -plot(n,abs(signal).^2,'-o','MarkerIndices',idxC,'MarkerSize',10); -hold on; plot(n(idx2_target),abs(signal(1,idx2_target)).^2,'p','MarkerEdgeColor','r','MarkerFaceColor','r'); -%hold on;stem(N/2,abs(signal(1,N/2)).^2,'-p','MarkerEdgeColor','r','MarkerFaceColor','r'); -hold on;plot(n,T_cfar); -xlim([c+2,N-c-1]); -Pfa_cfar = sum(abs(noise(1,c+2:N-c-1)).^2 > T_cfar(1,c+2:N-c-1))./(N-2*c) -% %Pfa_np = sum(x(CUTIdx,:)>T_ideal)/Ntrials - -%% Pd -N_trial = 10000; -sum_fix = 0; -sum_cfar = 0; -T_fix = sqrt(-sigma^2*log(Pfa)); -for n = 1:N_trial - noise = sigma/sqrt(2)*(randn(1,N)+1i*randn(1,N)); - signal = noise; - signal(1,CUTIdx) = signal(1,CUTIdx) + A; %5000个目标 - sum_fix = sum_fix + (abs(signal(1,CUTIdx)) > T_fix); - - pn = sum(abs(signal(1,(CUTIdx-c:CUTIdx-guardc-1))).^2) + sum(abs(signal(1,(CUTIdx+guardc+1:CUTIdx+c))).^2); - T_cfar0 =k.*abs(pn); - sum_cfar = sum_cfar + (abs(signal(1,CUTIdx)).^2 > T_cfar0); -end -pd_fix = sum_fix / N_trial %阈值固定下的恒虚警率 -pd_cfar = sum_cfar / N_trial - -%% CFAR COMPARE -close all;clear all;clc; -N = 6e3; -n = 1:N; -SNR = 20; -Pfa = 1e-6; -sigma_dB = 0; -sigma = db2pow(sigma_dB); %Noise Power=1w -ndeta = 100; -scale = 5; -noise = sigma*[randn(1,N/2-ndeta), scale*randn(1,N/2+ndeta)]; -A = sqrt((scale*sigma)^2 * 10^(SNR/10)); %Target signal amplitude -T_fix = sqrt(-sigma^2*log(Pfa)); %fix thresholding -c = N/40; %one side,ref cell = c-guardc -guardc = 10; %guard cell number -CUTIdx = N/2; %position of target - -%CFAR mode choose%%%%%,0 CA;1 GO;2 SO;3:OS -disp("请输入你接下来的操作： 0代表 CA-CFAR") -disp(" 1代表 GO-CFAR") -disp(" 2代表 SO-CFAR") -disp(" 3代表 OS-CFAR") -disp(" 4代表 退出 ") -flag=input("请输入："); -while(flag==0||flag==1||flag==2||flag==3) - cfar=flag; - %%Simulation of CFAR detection - %signal model - signal = noise; - signal(1,CUTIdx) = signal(1,CUTIdx) + A; %5000个目标 - %k = Pfa^(-1/(2*(c-guardc)))-1; - T_cfar = ones(1,N); - ncell = 0; - pn = 0; - for j=c+1:N-c - a = sum(abs(signal(1,(j-c:j-guardc-1))).^2); - b = sum(abs(signal(1,(j+guardc+1:j+c))).^2); - if cfar==0 - pn = a + b; - ncell = 2*(c-guardc); - end - if cfar==1 - pn = max(a,b); - ncell = c-guardc; - end - if cfar==2 - pn = min(a,b); - ncell = c-guardc; - end - k = Pfa^(-1/ncell)-1; - T_cfar(j) =k.*abs(pn); - end - n_target = 5; %target numbers - if cfar==3 - s = sort(signal); - T_cfar = ones(1,N) .* s(1,N-n_target).^2; - end - - %traget detected - idxC = find(abs(signal).^2 > T_cfar); - - figure; - plot(n,10*log(abs(signal).^2),'-o','MarkerIndices',idxC,'MarkerSize',10); - hold on;stem(N/2,10*log(abs(signal(1,N/2)).^2),'-p','MarkerEdgeColor','r','MarkerFaceColor','r'); - hold on;plot(n,10*log(T_cfar)); - xlim([c+2,N-c-1]);ylim([0,100]); - Pfa_cfar = sum(abs(noise(1,c+2:N-c-1)).^2 > T_cfar(1,c+2:N-c-1))./(N-2*c); - fprintf('虚警率:%s\n',num2str(Pfa_cfar)); - disp("当前程序执行完毕") - disp("请输入你接下来的操作： 0代表 CA-CFAR") - disp(" 1代表 GO-CFAR") - disp(" 2代表 SO-CFAR") - disp(" 3代表 OS-CFAR") - disp(" t代表 退出 ") - flag=input('请输入：'); -end -disp("已退出程序") \ No newline at end of file