怎么用实际雷达数据仿真MUSIC算法

（1） 不管测向天线阵列形状如何，也不管入射来波入射角的维数如何，假定阵列由M个阵元组成，则阵列输出模型的矩阵形式都可以表示为：Y(t)=AX(t)+N(t)

其中，Y是观测到的阵列输出数据复向量；X是未知的空间信号复向量；N是阵列输出向量中的加性噪声；A是阵列的方向矩阵；此处，A矩阵表达式由图册表示。

MUSIC算法的处理任务就是设法估计出入射到阵列的空间信号的个数D以及空间信号源的强度及其来波方向。

（2） 在实际处理中，Y得到的数据是有限时间段内的有限次数的样本（也称快拍或快摄），在这段时间内，假定来波方向不发生变化，且噪声为与信号不相关的白噪声，则定义阵列输出信号的二阶矩：Ry。

（3） MUSIC算法的核心就是对Ry进行特征值分解，利用特征向量构建两个正交的子空间，即信号子空间和噪声子空间。对Ry进行特征分解，即是使得图册中的公式成立。

（4） U是非负定的厄米特矩阵，所以特征分解得到的特征值均为非负实数，有D个大的特征值和M-D个小的特征值，大特征值对应的特征向量组成的空间Us为信号子空间，小特征值对应的特征向量组成的空间Un为噪声子空间。

（5） 将噪声特征向量作为列向量，组成噪声特征矩阵 ，并张成M-D维的噪声子空间Un，噪声子空间与信号子空间正交。而Us的列空间向量恰与信号子空间重合，所以Us的列向量与噪声子空间也是正交的，由此，可以构造空间谱函数。

（6） 在空间谱域求取谱函数最大值，其谱峰对应的角度即是来波方向角的估计值。

你确定你的程序可以运行出来？s_pc=s_pc.*s_doppler这一条不报错？改了一下：

clc

clear

fz=10e3%脉冲重复频率

tz=1/fz%脉冲重复周期

fs=1e5%采样频率

ts=1/fs%采样周期

f_doppler=2.5e3%多普勒频率

N=tz/ts%一个周期宽度

N=fix(N)%取接近且小于n的整数

%%%产生雷达回波%%%%

echo_mobj_pulse=[zeros(1,100),1,1,zeros(1,N-101)]

s_pc=repmat(echo_mobj_pulse,1,16)

n=1:length(s_pc)

s_doppler=cos(n*2*pi*f_doppler/fs)

s_pc=s_pc.*s_doppler

s_noise=0.1*rand(1,n)

s_pc=s_pc+s_noise

figure,

plot(n.*ts,s_pc),xlabel('t(单位：s)'),title('回波信号')

最后一句啥意思？没用过这个函数，帮助也没看懂。。。

雷达(radar)原是“无线电探测与定位”的英文缩写。雷达的基本任务是探测感兴趣的目标，测定有关目标的距离、方问、速度等状态参数。雷达主要由天线、发射机、接收机（包括信号处理机）和显示器等部分组成。 雷达发射机产生足够的电磁能量，经过收发转换开关传送给天线。天线将这些电磁能量辐射至大气中，集中在某一个很窄的方向上形成波束，向前传播。电磁波遇到波束内的目标后，将沿着各个方向产生反射，其中的一部分电磁能量反射回雷达的方向，被雷达天线获取。天线获取的能量经过收发转换开关送到接收机，形成雷达的回波信号。由于在传播过程中电磁波会随着传播距离而衰减，雷达回波信号非常微弱，几乎被噪声所淹没。接收机放大微弱的回波信号，经过信号处理机处理，提取出包含在回波中的信息，送到显示器，显示出目标的距离、方向、速度等。 为了测定目标的距离，雷达准确测量从电磁波发射时刻到接收到回波时刻的延迟时间，这个延迟时间是电磁波从发射机到目标，再由目标返回雷达接收机的传播时间。根据电磁波的传播速度，可以确定目标的距离为：S=CT/2 其中S：目标距离 T：电磁波从雷达到目标的往返传播时间 C：光速

求采纳

---