AGC.rar_AGC matlab_AGC算法_HRRP_HRRP 识别_雷达 目标库

clc clear all M = 30; N = 1000; pt = 0.3; for kk = 1:10 %对目标1进行帧划分 X_11 = load('E:\master\1.txt'); %读取回波数据 X_12 = ReShape(X_11,N); X_13 = X_12'; X_1 = Perp(X_13,pt); clear X_11 clear X_12 clear X_13 [K_A,A] = AGC_12_12(M,X_1); %目标1的训练样本库建立 A_t = TEST_M(K_A,A,M); clear A %对目标2，采样点取：201:700，角度取32度之内 X_21 = load('E:\master\2.txt'); %读取回波数据 X_22 = ReShape(X_21,N); X_23 = X_22'; X_2 = Perp(X_23,pt); clear X_21 clear X_22 clear X_23 [K_B,B] = AGC_12_12(M,X_2); %目标2的训练样本库建立 B_t = TEST_M(K_B,B,M); clear B %对目标3，采样点取：201:700，角度取32度之内 X_31 = load('E:\master\3.txt'); %读取回波数据 X_32 = ReShape(X_31,N); X_33 = X_32'; X_3 = Perp(X_33,pt); clear X_31 clear X_32 clear X_33 [K_C,C] = AGC_12_12(M,X_3); %目标1的训练样本库建立 C_t = TEST_M(K_C,C,M); clear C %对目标4，采样点取：201:700，角度取32度之内 X_41 = load('E:\master\4.txt'); %读取回波数据 X_42 = ReShape(X_41,N); X_43 = X_42'; X_4 = Perp(X_43,pt); clear X_41 clear X_42 clear X_43 [K_D,D] = AGC_12_12(M,X_4); %目标1的训练样本库建立 D_t = TEST_M(K_D,D,M); clear D SNR = 10; %建立测试样本库 I = 2:4:300; J = length(I); j = 1; for i = 1:J X_A(:,j) = X_1(:,I(i)); X_A(:,j+1) = X_1(:,I(i)+1); X_A(:,j+2) = X_1(:,I(i)+2); j = j+3; end j = 1; for i = 1:J X_B(:,j) = X_2(:,I(i)); X_B(:,j+1) = X_2(:,I(i)+1); X_B(:,j+2) = X_2(:,I(i)+2); j = j+3; end j = 1; for i = 1:J X_C(:,j) = X_3(:,I(i)); X_C(:,j+1) = X_3(:,I(i)+1); X_C(:,j+2) = X_3(:,I(i)+2); j = j+3; end j = 1; for i = 1:J X_D(:,j) = X_4(:,I(i)); X_D(:,j+1) = X_4(:,I(i)+1); X_D(:,j+2) = X_4(:,I(i)+2); j = j+3; end tic d_AA = test_result(X_A,A_t,K_A); d_AB = test_result(X_A,B_t,K_B); d_AC = test_result(X_A,C_t,K_C); d_AD = test_result(X_A,D_t,K_D); d_A = [d_AA;d_AB;d_AC;d_AD]; [S_A II_A] = min(d_A,[],1); RATE_A(kk) = length(find(II_A ==1))/length(II_A); d_BA = test_result(X_B,A_t,K_A); d_BB = test_result(X_B,B_t,K_B); d_BC = test_result(X_B,C_t,K_C); d_BD = test_result(X_B,D_t,K_D); d_B = [d_BA;d_BB;d_BC;d_BD]; [S_B II_B] = min(d_B,[],1); RATE_B(kk) = length(find(II_B ==2))/length(II_B); d_CA = test_result(X_C,A_t,K_A); d_CB = test_result(X_C,B_t,K_B); d_CC = test_result(X_C,C_t,K_C); d_CD = test_result(X_C,D_t,K_D); d_C = [d_CA;d_CB;d_CC;d_CD]; [S_C II_C] = min(d_C,[],1); RATE_C(kk) = length(find(II_C ==3))/length(II_C); d_DA = test_result(X_D,A_t,K_A); d_DB = test_result(X_D,B_t,K_B); d_DC = test_result(X_D,C_t,K_C); d_DD = test_result(X_D,D_t,K_D); d_D = [d_DA;d_DB;d_DC;d_DD]; [S_D II_D] = min(d_D,[],1); RATE_D(kk) = length(find(II_D ==4))/length(II_D); toc end rate_a = sum(RATE_A)/10; rate_b = sum(RATE_B)/10; rate_c = sum(RATE_C)/10; rate_d = sum(RATE_D)/10;