music2.txt

阵列信号处理波束形成算法仿真 有实际的应用价值 请参考

阵列信号处理波束形成卡朋算法仿真：该程序阵列数为N，采样点数为K，入射信号为30度角，干扰信号为60度角。 功能：提取有用信号抑制干扰信号 

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);