%%?螺旋型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