music2.m

来自「阵列信号处理波束形成算法仿真 有实际的应用价值 请参考」· M 代码 · 共 29 行

clc; 
clear all; 

N=12;
K=200;
M=2;
theta=[30,60]; 
Xt=zeros(N,K); 
SNR=0;
SIR=-20; 
Pd=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)));
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)); 

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);
plot(thetas,s);