MIMO雷达BP算法

%% 螺旋型MIMO + MC + BP
tic
clc;
clear all;
close all;
% %
N1 = 46; % 竖直方位维
N2 = 46; % 水平方位维
F =401; % 取频点35GHz
% dir=‘C:\Users\EW\Desktop\测试数据\1008\正向\左下\‘;
% S1 = [];
% for data_i=1:N1
%     for data_j = 1:N2  
%         i = （data_i-1）*N2 + data_j; 
%         MK1（data_j::）=importdata（[dirnum2str（i）‘.txt‘]）;
%     end
%     j = sqrt（-1）;
%     A1 = MK1（:1:）;
%     A1 = power（10A1/20）;    
%     P1 = MK1（:2:）;
%     S1（data_i::） = （A1）.*exp（1j*P1*pi/180）;
% end
% dir=‘C:\Users\EW\Desktop\测试数据\1008\正向\左上\‘;
% S2 = [];
% for data_i=1:N1
%     for data_j = 1:N2
%         i = （data_i-1）*N2 + data_j;
%         MK2（data_j::）=importdata（[dirnum2str（i）‘.txt‘]）; 
%     end
%     j = sqrt（-1）;
%     A2 = MK2（:1:）;
%     A2 = power（10A2/20）;    
%     P2 = MK2（:2:）;
%     S2（data_i::） = （A2）.*exp（1j*P2*pi/180）;
% end
% 
% dir=‘C:\Users\EW\Desktop\测试数据\1008\正向\右下\‘;
% S3 = [];
% for data_i=1:N1
%     for data_j = 1:N2
%         i = （data_i-1）*N2 + data_j;
%         MK3（data_j::）=importdata（[dirnum2str（i）‘.txt‘]）;
%     end
%     j = sqrt（-1）;
%     A3 = MK3（:1:）;
%     A3 = power（10A3/20）;    
%     P3 = MK3（:2:）;
%     S3（data_i::） = （A3）.*exp（1j*P3*pi/180）;
% end
% 
% dir=‘C:\Users\EW\Desktop\测试数据\1008\正向\右上\‘;
% S4 = [];
% for data_i=1:N1
%     for data_j = 1:N2
%         i = （data_i-1）*N2 + data_j;
%         MK4（data_j::）=importdata（[dirnum2str（i）‘.txt‘]）;
%     end
%     j = sqrt（-1）;
%     A4 = MK4（:1:）;
%     A4 = power（10A4/20）;    
%     P4 = MK4（:2:）;
%     S4（data_i::） = （A4）.*exp（1j*P4*pi/180）;
% end
% load（‘C:\Users\thz303\Desktop\WQL\数据\1021\左下数据\matlab.mat‘）;
% load（‘C:\Users\thz303\Desktop\WQL\数据\1021\左上数据\matlab.mat‘）;
% load（‘C:\Users\thz303\Desktop\WQL\数据\1021\右下数据\matlab.mat‘）;
% load（‘C:\Users\thz303\Desktop\WQL\数据\1021\右上数据\matlab.mat‘）;
% 
% % % % S1=flipud（S1）;%%左右
% % % % S1=fliplr（S1）;%%上下
AAA=1;
% % S1=fliplr（S1）;  
% % S2=flipud（S2）;
% % S3=flipud（S3）;
% % S4=fliplr（S4）;
% % S=S1+S2+S3+S4;
% S3=fliplr（S3）;
% S4=fliplr（S4）;

% Y_ch=S1;
%% 参数设置
c=3e8;
fc=500e9;%发射信号频率
fmin=490e9;
fmax=510e9;
B=20e9; %带宽
R=0.75;%探测平面与目标平面之间距离
lambda=c/fc;
w=2*pi*fc;
k=w/c;%波数
step=lambda/2;%收发机之间距离
Nx=20;%发射机个数
Ny=20;%接收机个数
P=801;
Nf=P;
f_step=B/（Nf-1）;            % 采样间隔
f=fmin+（0:Nf-1）*f_step;     % 采样频率
kmin=2*pi*fmin/c;           % 波数域频率
NB = 2*（fmax -fmin）/c;      % 带宽空域频点，单位距离由多少个频点来分辨（点数/米）
L_x=（Nx-1）*step;%总的发射机长度
L_y=（Ny-1）*step;%总的接收机长度

% theta_annt=20*pi/180;%方位角                       

%% X方位向分辨率

% tx = （-（Nx-1）/2:（Nx-1）/2）*step;  
% ty = zeros（1Nx）-（Nx-1）/2*step;  %左下，右下
% ty = zeros（1Nx）+（Nx-1）/2*step;  %左上，右上

tx = [0.4608  0.4656  0.4704  0.4752  0.4800  0.4848  0.4896  0.4944  0.4992  0.5040  -0.4608  -0.4656  -0.4704  -0.4752  -0.4800  -0.4848  -0.4896  -0.4944  -0.4992  -0.5040];
ty = [0.0425  -0.0425  0.0835  -0.0835  0.1245  -0.1245 0.1655  -0.1655  0