matlab 实现了一个基于阵列信号处理的空间角和极化参数估计系统

clear all
% close all
clc

twpi             = 2*pi;
rad              = pi/180;
deg              = 180/pi;
%%%%%---------------------搜索范围及步长确定---------------------------%%%%
Angle_range      = 50;
Search_step      = 1;
%% 信号参数生成
c                = 3e8;%光速
snr              = 20;%信噪比，单载频不大于23db
Center_Frequency = 6e9;%中心频率
lamda            = c/Center_Frequency;%接收信号波长
Noise_Power      = 2;%噪声功率
Sample_frequency = 1e9;%A/D采样频率
Sample_Number    = 200;%快拍数，单载频不超过300
Sample_Index     = [0:Sample_Number-1];%采样下标
Signal_Amplitude = sqrt(2*Noise_Power*10^(snr/10));
Array_Frequency1 = Center_Frequency;%信号1频率
S0               = exp(1j*2*pi*Array_Frequency1/Sample_frequency*Sample_Index);
%% 阵列参数生成
beishu           = 3.18;
R_6G             = 0.05*beishu;%对应6G波长半径
kelm             = 9;
%%%%%------------------------阵元位置生成------------------------------%%%%
Location_Y = [];
Location_Z = [];
for i = 1:kelm
    m = i-1;
    Location_y = R_6G*cos(twpi*m/kelm);
    Location_z = R_6G*sin(twpi*m/kelm);
    Location_Y = [Location_Y;Location_y];
    Location_Z = [Location_Z;Location_z];
end
Location = [Location_Y,Location_Z];

figure(1);              
scatter(Location(:,1),Location(:,2))
xlabel('Z');
ylabel('Y');
grid;
title('阵元分布图')
xlim([-0.5,0.5])
ylim([-0.5,0.5])
%%%%%---------------------电场增益矩阵生成-----------------------------%%%%
r_theta        = 90*ones(1,kelm)*rad;     %偶极子指向角theta （YOZ平面均为90°）
r_phi          = (0:360/kelm:360/kelm*(kelm-1))*rad;    %偶极子指向角phi (360°均匀分布)
r_b            =