基于cyclic-MUSIC的DOA估计Matlab仿真.docVIP

信息与通信工程学院 阵列信号处理实验报告 （基于cyclic-MUSIC的DOA估计Matlab仿真） 学 号：XXXXXX 专 业：XXXXXX 学生姓名：XXX 任课教师：XXX 2015年X月  题目：基于cyclic-MUSIC的DOA估计Matlab仿真 算法简述： DOA估计算法一般使用与平稳信号，而生活中测量的信号往往是非平稳的。在非平稳的信号中，有一类信号属于周期平稳信号，利用周期平稳性得到的算法称循环（共轭循环）平稳算法。循环平稳的循环自相关函数表示为 其中为循环角频率。 共轭循环平稳的循环自相关函数表示为 基于天线阵列协方差矩阵的特征分类DOA估计算法中，MUSIC算法具有普遍的实用性。只要已知天线阵的分布形式，无论直线阵还是圆阵，不管阵元是否等间隔分布，都可以得到高分辨率的估计结果。 实验内容与结果： 实验先仿真分析BPSK信号的循环平稳和共轭循环平稳性与循环频率和延时的的关系，然后使用12阵元均匀线阵，阵元间距为信号波长的一半，输入信号为1个BPSK信号、设置载波频率10MHz、采样频率50MHz、快拍数1000、信号方位角-4、信噪比为25dB情况下，得到DOA估计结果。 图1 BPSK 的循环平稳性 图2 BPSK共轭循环平稳性 图3 DOA仿真结果 仿真分析 由仿真结果可知，BPSK信号具有循环平稳特性，当循环频率为，信号具有循环平稳性，当循环频率为，信号具有共轭循环平稳性。由图一图二可得到合适的循环频率和延时，例如图一，循环频率为、延时为40Ts时其循环自相关函数值最大。实验最后针对cyclic-music算法的DOA估计，使用共轭循环平稳性设计算法，并在给定实验条件下仿真，有效估计出了信号的来向。 程序 1. 程序1： clear all; close all; clc; %--------------------------------- %线阵 波束形成 %--------------------------------- %参数设置 S_No = 1; sensor_No = 12; azimuth = [40*pi/180 ]; Fs = 40e6; Fc = 10e6; %干扰功率 RB = 1e6; M = 2;%二进制 wavelength = 3e8/Fc; d = wavelength/2; K = 1000; %快拍数 %--------------- %dB Ps = -0; SNR = 20; M_No = 40; %码速率等于信息速率 Data_No = M_No*Fs/RB; t = 1/Fs:1/Fs:Data_No/Fs; %--------------------------------- %功率转换 Ps_l = 10^(Ps/10); Pn_l = 10^((Ps-SNR)/10); %--------------------------------- %信号生成 bit = randint(1,M_No); %产生信息序列 bitstream = []; for i = 1:M_No if bit(i) == 1 bitstream = [bitstream, ones(1, Fs/RB)]; else bitstream = [bitstream, -ones(1, Fs/RB)]; end end Carrier_R = cos(2*pi*Fc*t); S_R = Carrier_R .* bitstream; Carrier_I = sin(2*pi*Fc*t); S_I = Carrier_I .* bitstream; Signal_R = sqrt(Ps_l)*S_R/sqrt((S_R*S_R/length(S_R))); Signal_I = sqrt(Ps_l)*S_I/sqrt((S_I*S_I/length(S_I)));%生成BPSK复信号 Signal(1,:) = complex(Signal_R, Signal_I); %--------------------------------- %{ %--------------------------------- %模拟天线接收 As = zeros(sensor_No,S_No); for i = 1:sensor_No for ii