一种基于fpga的目标旋转微多普勒信号产生方法
【技术领域】
[0001] 本发明涉及模拟微波信号产生领域，特别涉及旋转微多普勒调制信号模拟产生的 方法。
【背景技术】
[0002] 目前，雷达目标模拟器技术的发展已经比较成熟。现在的大多数雷达目标模拟器 可以模拟出目标的径向速度、RCS、距离、个数以及模拟假目标和欺骗目标等等。雷达目标模 拟器的发展与应用极大地方便了雷达系统的开发与研制。随着人们对目标微多普勒特性逐 渐深入研究，发现分析目标微多普勒特性可以成为一种新的雷达目标识别的手段，基于不 同目标的微运动特性不同，可以利用雷达对目标回波的微多普勒特性进行分析，进而辨别 不同目标。但是现在的雷达目标模拟器中都没有产生相应微多普勒信号的模块，使得对能 够进行微多普勒特性检测并分析的雷达系统的研制与开发造成了一定程度上的困难。
【发明内容】
[0003] 为了克服雷达目标模拟器在模拟目标旋转微多普勒特性的一些不足，本发明提出 一种能够模拟目标旋转微多普勒特性的方法，达到实现模块结构简单，易于实现，能够灵 活的调整参数，实时的模拟出目标旋转的微多普勒特性的目的。
[0004] 本发明解决其技术问题所采用的技术方案是一种基于FPGA的目标旋转微多普勒 信号产生方法，该方法首先根据需要模拟的目标标量角速度对产生需要变化相位，再对该 相位转换为正余弦信号，将该正余弦信号对以后的模拟雷达信号进行变化，获得加入了旋 转微多普勒效应的信号相位信息，最后将该相位信息转化为正余弦信号，从而实现发明目 的。因而本发明一种基于FPGA的目标旋转微多普勒信号产生方法包括：
[0005] 步骤1 :根据需要模拟目标的标量角速度通过累加器产生相位信息；
[0006] 步骤2 :将相位信息通过坐标旋转数字计算方法（CORDIC)将相位信息转换为正弦 和余弦信息；
[0007] 步骤3 :用步骤2获得正弦和余弦信号对模拟的雷达信号进行变化，获得具有旋转 微多普勒效应的信号相位信息；
[0008] 步骤4 :最后再通过坐标旋转数字计算方法将该相位信息转换为正弦和余弦信号 即具有旋转微多普勒效应的雷达信号。
[0009] 进一步的所述步骤3的具体方法为：
[0010] 步骤3. 1 :将步骤2获得余弦信息乘以B，再除以标量的角速度，这个为第1路数 据；
[0011] 其中
[0012] 其中 为目标的旋转角速度，rP表示目标在目标本地坐标系中自 身所处的位置，目标的旋转初始欧拉角为<，，#，y):
[0013]
[0014]
[0015]
[0016]
[0017]
[0018]
[0019]
[0020]
[0021]
[0022]
[0023] 步骤3. 2 :将步骤2获得正弦信号乘以C，再除以标量角速度，为第2路信号，其中
[0024] 步骤3. 3 :将获得第1路和第2路fg号进行相加，获得相加 fg号；
[0025] 步骤3. 4 :米用步骤3. 1、3. 2、3. 3相同的方法获得额外两路相加{目号；
[0026] 步骤3.5 :将获得三路相加信号，及模拟目标方位角向量进行乘累加运算；
[0027] 步骤3. 6 :将步骤3. 5的运算结果与A相乘，其中A = 2f Ω /c，f为载波频率，Ω为 ω的模值，ω表示目标的旋转角速度，c为光速；
[0028] 步骤3. 7 :将步骤3. 6的结果进行数据统筹转换，获得相位信息。
[0029] 进一步的所述步骤3中模拟雷达信号的数据包括：旋转目标相对于雷达的方位角 和俯仰角、旋转目标所处方位角和俯仰角、旋转初始欧拉角、角速度、雷达载波单频信号频 率、散射中心点相对于旋转中心的位置坐标。
[0030] 本发明的有益效果是，可以根据旋转微多普勒特性的不同灵活实时地调整参数， 模拟出目标在一定范围内不同的旋转微多普勒特性。
【附图说明】
[0031] 图1为本发明一种基于FPGA的刚体目标旋转微多普勒信号产生方法的原理框 图；
[0032] 图2为数据转换模块内部组成结构图；
[0033] 图3为雷达与旋转目标间的几何关系；
[0034] 图4为微多普勒调制信号实部时域MATLAB仿真图；
[0035] 图5为微多普勒调制信号虚部时域MATLAB仿真图；
[0036] 图6为微多普勒调制信号的MATLAB时频分析图。
【具体实施方式】
[0037] 如图3,雷达空间固定坐标系是（X，Y，Z)，目标本地坐标系（X，y，z)，参考坐标系 是（X'，Y'，Z')平行于雷达空间固定坐标系并位于目标本地坐标系的原点。
[0038] 目标在雷达坐标系（X，Y，Z)中的方位角和俯仰角分别是α和β，某散射中心旋 转初始欧拉角为（供，沒，V )，旋转角速度g ?
[0039] 则因旋转引起的多普勒频移为：
[0040]
(1)
[0041] f 为载波频率，Ω 为 ω 的模值，n = [cos(6T)*cos(々）sin(a)*cos(0) sin(#)]T，
[0042]
[0043]
[0044]
[0045]
[0046]
[0047]
[0048]
[0049]
[0050]
[0051]
[0052] 以目标本地坐标系为观察坐标系，当目标围绕轴x、y和z以角速度ω = (ωχ，ωγ，ωζ)τ旋转时，在目标本地坐标系中所表示的位于r p= (Xp，yp，Zp)τ处的点散射 P将移到参考坐标系中用描述的新位置，并且旋转单位矢量变为
[0053]
[0054]
[0055] 在时间t旋转矩阵变为
[0056]
CN 105158746 A m ~P 4/4 页
[0057]
[0058] 由于公式（I)中所得的是所要产生的调制信号的频率值，由相位与频率的关系以 及利用公式（1)的积分得到相位值表达式：
[0059] (2)
[0060]
[0061] 该公式即为该新发明方法中所要用到的基本理论公式。
[0062] 如图1，累加器用于根据输入的频率字产生相应的相位，CORDIC模块1从累加器接 收相位产生相应的以对正余弦值进行一系列转换，（旋转目标相对于雷达的方位角和俯仰 角、旋转目标的空间旋转速度、雷达载波单频信号频率在A、B、C、D因子中有所包含，后文会 有所详述）。另外CORDIC模块2接收经转换所得的值直接作为输入所需的相位值，产生微 多普勒调制信号。
[0063] 假设雷达载波单频信号频率为1G，且位于（U = 1000m，V = 5000m，W = 5000m)的 目标围绕着x、y和z轴以初始欧拉角（炉=_3:(f，θ = 20°，Φ = 20° )和角速度ω = [JT，2 π，π ]作旋转。假设其中一个散射点的位置为（χ = 5m，y = 5m，z = 5m)。
[0064] 在公式（I)，A = 2f Ω /c，. &，c=6/ , , 所以 A、 B、C、D的值可提前通过matlab计算得到，A、B、C、D作为该发明模块的输入值，继而得到相 应的调制信号。
[0065] 在该例中输入的Ω要经过（Ω/2 π ) X 232的处理再作为频率字输入模块中，这是 因为在表达式（2)中cos Ω t和sin Ω t是作为角频率值，范围大小为2 π，而该模块的频率 字输入口为32位，范围大小为232。所以要先对Ω进行转换，实现数据的归一化。经归一化 后的数据作为频率字输入给相位累加器，累加器输出相应的相位值出输给CORDIC模块1， CORDIC模块输出对应的sin Ω t值和cos Ω t
[0066] 该例中在输入值方面，为了在内部数据处理的精确性，对B、C值输入模块同时提 前进行了 232的放大，对D值也同时进行2 16倍的扩大，为了保证数据输出的正确性，在处理 过程中，对数据会进行相应的截位。
[0067] 将得到的cos Ω t与B相乘再除以Ω (在该处，Ω为原值）得到行向量b (维数为 3)，将sin Ω t与C相乘再除以Ω (在该处，Ω为原值）得到行向量c (维数为3)，将b与c 相加后再与D相乘，最后再与A相乘。但该处得到的相位值不能直接作为相位值传给CORDIC 模块2,因为与A相乘得到的值为真实的相位值，范围为2 π，而CORDIC模块2的相位值输 入端口为15位，范围为215。所以要将与A相乘得到的结果除以2JI后再乘以2 15。最后得 到的结果作为相位输入值传入CORDIC模块2,从而得到相应的旋转微多普勒调制信号。相 关结果请见附图。
【主权项】
1. 一种基于FPGA的目标旋转微多普勒信号产生方法，该方法包括： 步骤1 :根据需要模拟目标的标量角速度通过累加器产生相位信息； 步骤2 :将相位信息通过坐标旋转数字计算方法将相位信息转换为正弦和余弦信息； 步骤3 :用步骤2获得正弦和余弦信号对模拟的雷达信号进行变化，获得具有旋转微多 普勒效应的信号相位信息； 步骤4 :最后再通过坐标旋转数字计算方法将该相位信息转换为正弦和余弦信号即具 有旋转微多普勒效应的雷达信号。2. 如权利要求1所述的一种基于FPGA的目标旋转微多普勒信号产生方法，其特征在于 所述步骤3的具体方法为： 步骤3. 1 :将步骤2获得余弦信息乘以B，再除以标量的角速度，这个为第1路数据； 其中其中ω为目标的旋转角速度，^表示目标在目标本地坐标系中自身所 处的位置，目标的旋转初始欧拉角为(供，夂，H步骤3. 2 :将步骤2获得正弦信号乘以C，再除以标量角速度，为第2路信号，其中 C-?/· 9!..... · T ； 步骤3. 3 :将获得第1路和第2路fg号进行相加，获得相加 fg号； 步骤3. 4 :采用步骤3. 1、3. 2、3. 3相同的方法获得额外两路相加信号； 步骤3. 5 :将获得三路相加信号，及模拟目标方位角向量进行乘累加运算； 步骤3. 6 :将步骤3. 5的运算结果与A相乘，其中A = 2f Ω /c，f为载波频率，Ω为ω 的模值，ω表示目标的旋转角速度，c为光速； 步骤3. 7 :将步骤3. 6的结果进行数据统筹转换，获得相位信息。3.如权利要求1所述的一种基于FPGA的目标旋转微多普勒信号产生方法，其特征在于 所述步骤3中模拟雷达信号的数据包括：旋转目标相对于雷达的方位角和俯仰角、旋转目 标所处方位角和俯仰角、旋转初始欧拉角、角速度、雷达载波单频信号频率、散射中心点相 对于旋转中心的位置坐标。
【专利摘要】该发明公开了一种基于FPGA的目标旋转微多普勒信号产生方法，涉及模拟微波信号产生领域。该方法包括：根据需要模拟目标的标量角速度通过累加器产生相位信息；将相位信息通过坐标旋转数字计算方法将相位信息转换为正弦和余弦信息；用步骤2获得正弦和余弦信号对模拟的雷达信号进行变化，获得具有旋转微多普勒效应的信号相位信息；最后再通过坐标旋转数字计算方法将该相位信息转换为正弦和余弦信号即具有旋转微多普勒效应的雷达信号。从而具有可以根据振动微多普勒特性的不同灵活实时地调整参数模拟出一定范围内不同的振动微多普勒特性。
【IPC分类】G01S7/40
【公开号】CN105158746
【申请号】CN201510494811
【发明人】周云, 林杰, 李俊慧, 姜庆, 汪学刚
【申请人】电子科技大学
【公开日】2015年12月16日
【申请日】2015年8月13日