CBF算法.docx
- 文档编号:2486556
- 上传时间:2022-10-30
- 格式:DOCX
- 页数:2
- 大小:12.44KB
CBF算法.docx
《CBF算法.docx》由会员分享,可在线阅读,更多相关《CBF算法.docx(2页珍藏版)》请在冰豆网上搜索。
%%%%%%非相干信号,矢量阵常规波束形成
clc;
closeall;
clearall;
M=10;%阵元数
m=[0:
M-1]';
N=1000;%采样点数
n=0:
N-1;
f=1000;%频率
fs=10*f;%采样频率
c=1500;%声速
lamd=c/f;%波长
%d=lamd/2;%阵元间距
d=0.7*lamd;
sita0=80*pi/180;%期望方向
sita1=80*pi/180;%目标1方位
sita2=100*pi/180;%目标2方位
sita3=225*pi/180;%目标3方位
snr1=20;%信噪比
snr2=20;
snr3=20;
s1=exp(-j*(2*pi*f*n/fs+2*pi*rand(1,N)));
s2=exp(-j*(2*pi*f*n/fs+2*pi*rand(1,N)));
s3=exp(-j*(2*pi*f*n/fs+2*pi*rand(1,N)));
s=[s1;s2;s3];
a0=exp(-j*2*pi*f*m*d*cos(sita0)/c);
u0=[1cos(sita0)sin(sita0)].';
Av0=kron(a0,u0);
a1=exp(-j*2*pi*f*m*d*cos(sita1)/c);
a2=exp(-j*2*pi*f*m*d*cos(sita2)/c);
a3=exp(-j*2*pi*f*m*d*cos(sita3)/c);
u1=[1cos(sita1)sin(sita1)].';
u2=[1cos(sita2)sin(sita2)].';
u3=[1cos(sita3)sin(sita3)].';
Av=[kron(a1,u1)kron(a2,u2)kron(a3,u3)];
%%%有噪声
S1=awgn(Av(:
1)*s1,snr1,'measured','db');
S2=awgn(Av(:
2)*s2,snr2,'measured','db');
S3=awgn(Av(:
3)*s3,snr3,'measured','db');
S=S1+S2;%+S3;
R=S*S'/N;
step=1;
Ang=0:
step:
360;
se=zeros(N,length(Ang));
A=zeros(3*M,0);
forkk=1:
length(Ang)
u=[1cos(Ang(kk)*pi/180)sin(Ang(kk)*pi/180)]';
a=kron(exp(-i*(2*m*pi*f*d*cos(Ang(kk)*pi/180)/c)),u);
Pc(kk)=a'*R*a;
end
figure
(1);
plot(Ang,10*log10(abs(Pc)/max(abs(Pc))));%,'LineWidth',2);
gridon;
axis([0360-350]);
xlabel('方位(^o)');
ylabel('归一化功率(dB)');
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- CBF 算法