文档介绍:
(1)matlab程序
clc;
clearall;
closeall;
imag=sqrt(-1);
element_num=32;%阵元数为8
d_lamda=1/2;%车元间距d与波长lamda的关系rndBnorm=20*log10(patternmagnorm);
figure(1)
subplot(1,2,1);
plot(theta,patterndBnorm);
gridon;
xlabel('theta/radian');
ylabel('amplitude/dB');
axis([-500]);
subplot(1,2,2);
pfft=fftshift(fft(w,256));
pfftmag=abs(pfft);
pfftmagnorm=pfftmag/max(max(pfftmag));
pfftdB=20*log10(pfftmagnorm);
pfftdBnorm=20*log10(pfftmagnorm);
plot(linspace(-pi/2,pi/2,256),pfftdBnorm);
gridon;
xlabel('theta/radian');
ylabel('FFT_amplitude/dB');
axis([-500]);
2)仿真结果
最大信噪比准则方向图和功率谱
matlab程序
clc;
clearall;
closeall;
imag=sqrt(-1);
element_num=8;%阵元数为8
d_lambda=;%间距为半波长
theta=-90::90;%扫描范围
theta0=0;%来波方位
theta1=20;%干扰方向
L=512;%采样点数
fori=1:L
amp0=10*randn(1);
amp1=200*randn(1);
ampn=1;s(:,i)=amp0*exp(imag*2*pi**sin(theta0*pi/180)*[0:element_num-1]');
j(:,i)=amp1*exp(imag*2*pi**sin(theta1*pi/180)*[0:element_num-1]');
n(:,i)=ampn*exp(randn(element_num,1)+imag*randn(element_num,1));
end
Rs=1/L*s*s';%信号自相关矩阵
Rnj=1/L*(j*j'+n*n');%干扰+噪声的自相关矩阵
[V,D]=eig(Rs,Rnj);%(Rs,Rnj)的广义特征值和特征向量
[D,I]=sort(diag(D));%特征向量排序
Wopt=V(:,I(8));%最优权矢量
forj=1:length(theta)
a=exp(imag*2*pi*d_lambda*sin(theta(j)*pi/180)*[0:element_num-1]');
f(j)=Wopt'*a;
p(j)=a'*Rs*a+a'*Rnj*a;
end
F=20*log10(abs(f)/max(max(abs(f))));
P=20*log10(abs(p)/max(max(abs(p))));
subplot(1,2,1)
plot(theta,F);
gridon;
holdon;
plot(theta0,-80:0,'.');
plot(theta1,-80:0,'.');
xlabel('theta/0');
ylabel('FindB');
title('max-SNR方向图');
axis([-9090-800]);
holdon;
subplot(1,2,2);
plot(theta,P,'r');
gridon;
xlabel('theta/0');
ylabel('功率indB');
title('max-SNR功率谱');
gridon;
axis([-9090-800]);
2)仿真结果
ASC旁瓣相消----MSE准则
(1)matlab程序
clc;
closeall;
clearall;
imag=sqrt(-1);
M=32;%ffi助天线数目
d_lambda=;%阵元间距
theta0=-30;%来波方向
theta1=60;%干扰方向
L=512;%采样单元数
s=zeros(1,512);%预划分一个区域
forii=1:L
amp0=1*randn(1);%信号的幅度随机产生,保证信号之间是不相关的
amp1=200*randn(1);
ampn=1;
jam(:,ii)=amp1*ex