%% 11月29日完全极化MUSIC方法仿真 存在完全非圆信号 完全非圆信号个数为M1,部分非圆信号个数M2(共计M个)
clear;close;
%% 选取M个互不相关的入射信号,阵元数为N,阵列形状为线阵,阵元等间距分布在y轴上(可以用其他形状的阵列代替,只需要改变空间相位因子即可)
%% 认为满足远场和窄带假设(当然不满足的情况也是改变空间相位矩阵即可)
%% 本次仿真因发现如果采用sin*sin的空间相位角,会发生相位模糊,于是采用sin*cos的空间相位角
%% 参数设置
N=10;
M=4;
M1=3; %3个非圆信号
M2=1; %一个圆信号
theta=[1.3 1.7 1.5 0.2];
phi=[1.5 1.2 1.3 1.4];
A1=zeros(N,N);
A2=zeros(N,N);
A3=zeros(N,N);
A4=zeros(N,N);
c=3e8;
fc=1e9; %载波频率
lamda=c/fc;
d=lamda/2;
L=100*N; %快拍数
%Ai是空间相位因子
for i=1:N
A1(i,i)=exp(1i*2*pi*d*sin(theta(1,1))*cos(phi(1,1))*(i-1)/lamda); %进行修改测试
A2(i,i)=exp(1i*2*pi*d*sin(theta(1,2))*cos(phi(1,2))*(i-1)/lamda);
A3(i,i)=exp(1i*2*pi*d*sin(theta(1,3))*cos(phi(1,3))*(i-1)/lamda);
A4(i,i)=exp(1i*2*pi*d*sin(theta(1,4))*cos(phi(1,4))*(i-1)/lamda);
end
A=[A1 A2 A3 A4];
%考虑g1=g2=g3=g10=1 阵元1到10的指向角分别(在yoz平面内,磁环线阵)为(90°,30°;90°,30°;90°,30°;270°,30°;270°,30°;270°,30°;270°,30°;90°,60°;90°,60°;90°,60°)
%极化敏感矩阵为
dabeta=[0 0 0 0 sin(pi/6) cos(pi/6);
0 0 0 0 sin(pi/6) cos(pi/6);
0 0 0 0 sin(pi/6) cos(pi/6);
0 0 0 0 sin(-pi/6) cos(-pi/6);
0 0 0 0 sin(-pi/6) cos(-pi/6);
0 0 0 0 sin(-pi/6) cos(-pi/6);
0 0 0 0 sin(-pi/6) cos(-pi/6);
0 0 0 0 sin(pi/3) cos(pi/3);
0 0 0 0 sin(pi/3) cos(pi/3);
0 0 0 0 sin(pi/3) cos(pi/3)];
dakesi=zeros(6,2*M);
for kkkkk=1:M
dakesi(:,(2*kkkkk-1):2*kkkkk)=[-sin(theta(1,kkkkk)) cos(phi(1,kkkkk))*cos(theta(1,kkkkk));
cos(theta(1,kkkkk)) cos(phi(1,kkkkk))*sin(theta(1,kkkkk));
0 -sin(phi(1,kkkkk));
cos(phi(1,kkkkk))*cos(theta(1,kkkkk)) sin(theta(1,kkkkk));
cos(phi(1,kkkkk))*sin(theta(1,kkkkk)) -cos(theta(1,kkkkk));
-sin(phi(1,kkkkk)) 0];
end
gama=[pi/4 pi/4 pi/3 pi/3]; %极化辅助角
eta=[-pi/4 0 pi/4 pi/2]; %极化相位角
jihuah=zeros(2,M);
for jjj=1:M
jihuah(1,jjj)=cos(gama(jjj));
jihuah(2,jjj)=sin(gama(jjj))*exp(1i*eta(jjj));
end
daoxiangshiliangA=zeros(N,M);
%导向矢量A
for kk=1:M
daoxiangshiliangA(:,kk)=A(:,(1+N*(kk-1)):(N*kk))*dabeta*dakesi(:,(2*kk-1):(2*kk))*jihuah(:,kk);
end
%对信源矢量S进行定义,不防将信源矢量定义为正弦函数sin(mt),整数点采样,m表示第m个信源;
%L为快拍数
%sin(0.3pikkk*x)为不同的信源表达式
%% 信源生成
S=zeros(M,L);
%% 圆信源与非圆信源产生待修改
for kkkk=1:L
S(1,kkkk)=20*(-1)^kkkk*sin(0.7*pi*kkkk);
S(2,kkkk)=20*sin(0.1*pi*kkkk)*exp(1i*pi/4); %前三个为非圆信号,最后一个为圆信号
S(3,kkkk)=20*sin(0.2*pi*kkkk)*exp(1i*pi/6); %信源进行修改,改为非圆周期的信源
S(4,kkkk)=10*exp(1i*0.7*pi*(kkkk-1));
%S(kkk,kkkk)=1;
end
%%
AM=daoxiangshiliangA(:,1:M1);
AN=daoxiangshiliangA(:,(M1+1:M));
feiyuanxiangwei=[0 pi/2 pi/3]; % 给出非圆相位的定义,M1个非圆信号的非圆相位
zengguang_AM=zeros(size(AM,1),size(AM,2));
for jj=1:M1
zengguang_AM(:,jj)=conj(AM(:,jj))*exp(-1i*feiyuanxiangwei(1,jj));
end
AR=[AM;zengguang_AM]; %增广的非圆导向矢量
ACM=zeros(2*size(AN,1),2*size(AN,2));
for kk=1:M2
ACM(:,2*kk-1)=[AN(:,kk);zeros(size(AN(:,kk),1),size(AN(:,kk),2))];
ACM(:,2*kk)=[zeros(size(AN(:,kk),1),size(AN(:,kk),2));conj(AN(:,kk))];
end
daB=[AR,ACM];
Srt=zeros(M1,L);
for kkk=1:M1
for jjj=1:L
Srt(kkk,jjj)=S(kkk,jjj); %非圆信源
end
end
Sct=zeros(2*M2,L);
for kk=1:M2
Sct(2*kk-1,:)=S(M1+kk,:);
Sct(2*kk,:)=conj(S(M1+kk,:));
end
Rt=[Srt;Sct];
nt=wgn(N,L,-1000);
zengguang_nt=[nt;conj(nt)];
yt=daB*Rt+zengguang_nt; %阵列增广输出矢量N*L L为快拍数
%% 开始估计
gujiRyy=zeros(2*N,2*N);
for u=1:L
gujiRyy=gujiRyy+(1/L)*(yt(:,u)*yt(:,u)'); %% 估计阵列样本协方差矩阵
end
[eig_E,eig_V]=eig(gujiRyy); %特征值分解 eig_E是特征向量矩阵,eig_V是特征值矩阵
[eig_V_sort,index] = sort(diag(eig_V),'descend');
eig_E_sort = eig_E(:,index);
Un=eig_E_sort(:,(2*M2+M1+1):2*N); %次特征矢量矩阵
Un1=Un(1:N,:);Un2=Un(N+1:2*N,:);
Un_gouzao=[Un1,zeros(size(Un1,1),size(Un1,2));zeros(size(Un2,1),size(Un2,2)),Un2];
%% theta和phi的搜索 分为非圆搜索和圆搜索
%% 首先进行非圆搜索
sousuotheta=0.01:0.01:1.8;
sousuophi=0.01:0.01:1.7;
deH_feiyuan=zeros(size(sousuotheta,2),size(sousuophi,2));
for t=1:size(sousuotheta,2)
for tt=1:size(sousuophi,2)
Dthetaphi=zeros(N,2);
ceshiA=zeros(N,N);
for z=1:N
ceshiA(z,z)=exp(1i*2*pi*d*sin(sousuotheta(1,t))*cos(sousuophi(1,tt))*(z-1)/lamda); %修改测试
end
%不同的theta,phi给出不同的空间相位矩阵
ceshikesi=[-sin(sousuotheta(1,t)) cos(sousuophi(1,tt))*cos(sousuotheta(1,t));
cos(sousuotheta(1,t)) cos(sousuophi(1,tt))*sin(sousuotheta(1,t));
0 -sin(sousuophi(1,tt));
cos(sousuophi(1,tt))*cos(sousuotheta(1,t)) sin(sousuotheta(1,t));
cos(sousuophi(1,tt))*sin(sousuotheta(1,t)) -cos(sousuotheta(1,t));
-sin(sousuophi(1,tt)) 0];
Dthetaphi=ceshiA*dabeta*ceshikesi;
Dthetaphi_zengguang=[Dthetaphi zeros(size(Dthetaphi,1),size(Dthetaphi,2));zeros(size(Dthetaphi,1),size(Dthetaphi,2)) conj(Dthetaphi)];
deH_feiyuan(t,tt)=1/(det(Dthetaphi_zengguang'*(Un)*(Un')*Dthetaphi_zengguang));
end
end
figure(1)
[meshphi,meshtheta]=meshgrid(sousuophi,sousuotheta);
mesh(meshtheta,meshphi,10*log10(abs(deH_feiyuan)));
colorbar;
colormap jet;
caxis([0 350]);
grid minor;
xlabel('\theta');
ylabel('\phi');
axis([0 1.8 0 1.7 0 350])
title('对完全非圆信号DOA估计')
% 以上完成非圆信号估计 估计个数为M1
%% 进行部分非圆搜索
deH_yuan=zeros(size(sousuotheta,2),size(sousuophi,2));
for t=1:size(sousuotheta,2)
for tt=1:size(sousuophi,2)
Dthetaphi=zeros(N,2);
ceshiA=zeros(N,N);
for z=1:N
ceshiA(z,z)=exp(1i*2*pi*d*sin(sousuotheta(1,t))*cos(sousuophi(1,tt))*(z-1)/lamda); %修改测试
end
%不同的theta,phi给出不同的空间相位矩阵
ceshikesi=[-sin(sousuotheta(1,t)) cos(sousuophi(1,tt))*cos(sousuotheta(1,t));
cos(sousuotheta(1,t)) cos(sousuophi(1,tt))*sin(sousuotheta(1,t));
0 -sin(sousuophi(1,tt));
cos(sousuophi(1,tt))*cos(sousuotheta(1,t)) sin(sousuotheta(1,t));
cos(sousuophi(1,tt))*sin(sousuotheta(1,t)) -cos(sousuotheta(1,t));
-sin(sousuophi(1,tt)) 0];
Dthetaphi=ceshiA*dabeta*ceshikesi;
Dthetaphi_zengguang=[Dthetaphi zeros(size(Dthetaphi,1),size(Dthetaphi,2));zeros(size(Dthetaphi,1),size(Dthetaphi,2)) conj(Dthetaphi)];
deH_yuan(t,tt)=1/(det(Dthetaphi_zengguang'*(Un_gouzao)*(Un_gouzao')*Dthetaphi_zengguang));
end
end
figure(2)
mesh(meshtheta,meshphi,10*log10(abs(deH_yuan)));
colorbar;
colormap jet;
caxis([0 350]);
grid minor;
xlabel('\theta');
ylabel('\phi');
axis([0 1.8 0 1.7 0 350])
title('对部分非圆信号或者圆信号DOA估计')
%% 对搜索结果进行分类
%% 非圆极化参数估计
tezhengzhisousuo_feiyuan=zeros(size(sousuotheta,2),size(sousuophi,2));
guangyitezhengshiliang_feiyuan=[];
for t=1:size(sousuotheta,2)
for tt=1:size(sousuophi,2)
Dthetaphi=zeros(N,2);
ceshiA=zeros(N,N);
for z=1:N
ceshiA(z,z)=exp(1i*2*pi*d*sin(sousuotheta(1,t))*cos(sousuoph
array.rar_天线极化_极化 阵列_极化天线_极化敏感阵列_阵列
版权申诉
5星 · 超过95%的资源 41 浏览量
2022-07-15
14:15:08
上传
评论 3
收藏 3KB RAR 举报
JaniceLu
- 粉丝: 78
- 资源: 1万+
评论6