程序如下
clear;clc;
%% param
t=0:0.01:5;
m=length(t); %快拍数
fc=3e3; %期望信号频率
fc_j=2e3; %干扰信号频率
M=5; %阵元数目
n=(0:M-1)';
theta_d=pi/9; %期望信号波达方向
theta_j=linspace(-pi/3,0,3); %干扰信号波达方向
theta_j2=theta_j/pi*180;
theta=[theta_d,theta_j];
A=exp(-1j*pi.*n*sin(theta));
