capon_algorithm.m

阵列信号处理波束形成卡朋算法仿真：该程序阵列数为N

clc;
clear all;
% 模拟信号
N=12;   % 阵源数（可变）
K=200;  % 采样点数
M=2;    % 第一个为期望信号，第二个为干扰信号
theta=[30,60];
Xt=zeros(N,K);
SNR=0;  % 信噪比（可变）
SIR=-20;    % 信干比（可变）
Pd=1;   % 发射信号功率设为1
Pn=Pd/(10^(SNR/10));    % 由信噪比计算噪声功率
Pi=Pd/(10^(SIR/10));    % 由信干比干扰信号功率
theta = theta*pi/180;
Ad=exp(-j*pi*[0:(N-1)]'*sin(theta(1)));   % 设波长等于2d
Ai=exp(-j*pi*[0:(N-1)]'*sin(theta(2)));
Sd=sqrt(Pd/2)*(randn(1,K)+j*randn(1,K));
Si=sqrt(Pi/2)*(randn(1,K)+j*randn(1,K));
Xt=Ad*Sd+Ai*Si;
% 考虑期望信号、干扰信号和噪声的最终接受到的信号：
Xt=Xt+sqrt(Pn/2)*(randn(N,K)+j*randn(N,K)); 
% Capon Beamform
Rx=Xt*Xt'/K;
Wop=pinv(Rx)*Ad;
thetas=-90:90;
tm=thetas*pi/180;
am=exp(-j*pi*[0:N-1]'*sin(tm));
s=abs(Wop'*am);  % 波束形成后的Y(t)
plot(thetas,s);