基于matlab的归一化阵列响应幅值极坐标图仿真,对比四阵元,八阵元以及十六阵元

1.算法仿真效果

matlab2022a仿真结果如下：

 

2.算法涉及理论知识概要

       直线阵:由N个相同的单元天线等间距地排列在一条直线上构成。均匀直线阵:若各单元上的馈电电流大小相同，而相位沿线均匀递增或递减，这样的直线阵称为均匀直线阵。

 

 

3.MATLAB核心程序 `m=4 ;% array阵元

p=4; %  signal number信号数

N=3000;% recursive number 迭代次数 或快拍数

A=zeros(m,p); % array pattern阵列流型

theta=[30 0 -45 60]*pi/180;% the signal from the direction of 30 degree is expected. DOA 30为期望信号方向

j=sqrt(-1);

w=[0.01 0.2 0.3 0.4]*pi; % frequency for each signal.各个信号的数字频率

 

%%

s=to_get_s(w,N,p);

s_rec=get_s_rec(s,m,p,theta);

S=s_rec; %  output date matrix  .m*N 的阵元输出数据矩阵

%%%%——————————————%% 自适应调节权

y=S;% input data of GSC;

ad=exp(-jpi[0:m-1]'*sin(theta(1))) %steering vector in the direction of expected. 期望信号方向导向矢量

c=10;%  condition 波束形成条件

C=ad';

Wc=C'inv(CC')*c; %main path weight 主通道固定权

wa(1:m-1,1)=0;  % auxiliary path  辅助通道自适应权

B=get_B(m,theta); % get Block Matrix 得到阻塞矩阵

u=0.000001;

.................................................................

%%%%------------

%%%main path 主通道

wop=Wc;

drawpp(m,wop);

%%%%auxiliary path 辅助通道

wop=B'*wa(:,N)

drawpp(m,wop);

%%array response 总的阵列响应

wop=Wc-B'*wa(:,N);

drawpp(m,wop);`