CBF算法Word下载.docx

CBF算法Word下载.docx

《CBF算法Word下载.docx》由会员分享，可在线阅读，更多相关《CBF算法Word下载.docx（2页珍藏版）》请在冰豆网上搜索。

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,'

S3=awgn（Av（:

3）*s3,snr3,'

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）'