CBF算法Word下载.docx

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=ex

2、p（-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=1 cos（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=1 cos（sita1） sin（sita1）.u2=1 cos（sita2） sin（sita2）.u3=1 cos（sita3） sin

3、（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）;for kk=1:length（Ang） u=1 cos（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;endfigure（1）;plot（Ang,10*log10（abs（Pc）/max（abs（Pc）;%,LineWidth,2）;grid on;axis（0 360 -35 0）;xlabel（方位（o）ylabel（归一化功率（dB）