首页 / 专利库 / 电信 / 波达方向 / 一种短波二维波达方向估计数字处理模块

一种短波二维波达方向估计数字处理模

阅读:938发布:2020-05-13

专利汇可以提供一种短波二维波达方向估计数字处理模专利检索,专利查询,专利分析的服务。并且本 发明 提供一种短波二维 波达方向 估计数字处理模 块 ,包括10片A/D转换芯片、3片数字下变频芯片、FPGA芯片以及OMAP芯片;A/D转换芯片将采集到的模拟 信号 转换为 数字信号 、并将转换后的数字信号分别送给3片数字下变频芯片;3片数字下变频芯片将接受到的数字信号进行数字下变频处理、并将处理后的10路变频信号送到FPGA芯片;FPGA芯片对接受到的10路变频信号进行滤波和二维波达方向估计运算,将估计得到的二维波达的 俯仰 角 信息送给OMAP芯片;OMAP芯片将FPGA高速串口送来的二维波达的俯仰角、方位角数据转换成 并行总线 数据送给CPCI 接口 芯片。采用处理能 力 较强的OMAP芯片和FPGA芯片,实现对短波信号进行二维波达方向估计,具有良好的 软件 可移植性;同时采用C PCI总线 ,通用性强。,下面是一种短波二维波达方向估计数字处理模专利的具体信息内容。

1.一种短波二维波达方向估计数字处理模,其特征在于:包括10片A/D转换芯片、3片数字下变频芯片、FPGA芯片以及OMAP芯片;
所述A/D转换芯片将采集到的模拟信号转换为数字信号、并将转换后的数字信号分别送给3片数字下变频芯片;
3片数字下变频芯片将接受到的数字信号进行数字下变频处理、并将处理后的10路变频信号送到FPGA芯片;
FPGA芯片对接受到的10路变频信号进行滤波和二维波达方向估计运算,将估计得到的二维波达的俯仰信息送给OMAP芯片;
OMAP芯片将FPGA高速串口送来的二维波达的俯仰角、方位角数据转换成并行总线数据送给CPCI接口芯片;
所述二维波达方向估计运算采用基于L型阵列的联合奇异值分解法,具体步骤如下:
设有K个短波信号从远场入射到平行阵列,其中K(K≤5),信号s(t)=[s1(t),…,sK(t)];阵元1-5收到的信号为X(t),阵元6-10收到的信号为Z(t),Nx(t)和Nz(t)分别是阵元1-
5和阵元6-10的噪声矢量,A(α)、A(θ),如下式:
X(t)=A(α)s(t)+Nx(t),
Z(t)=A(θ)s(t)+Nz(t);
其中:
A(α)=[α(φ1),α(φ2),...,α(φk),...,α(φK)],
A(θ)=[α(θ1),α(θ2),...,α(θk),...,α(θK)],
计算X(t)和Z(t)的互协方差矩阵:
Rxz=E[X(t)ZH(t)];
将Rxz分为两个5×4维的子矩阵:
Rxz,1=Rxz(:1,:4),
Rxz,2=Rxz(:,2:5);
将子矩阵Rxz,1和Rxz,2组成新的矩阵R:
对矩阵R进行奇异值分解,可以计算出左奇异矢量U1,将U1分为两个5×K维子阵U11和U12:
求取 #表示逆阵;然后可以计算出其第k(k=1,…,K)个特征值γk,从而可求得俯仰角和方位角:
俯仰角:
方位角:
2.根据权利要求1所述的一种短波二维波达方向估计数字处理模块,其特征在于:所述短波二维波达方向估计数字处理模块采用CPCI总线形式。
3.根据权利要求2所述的一种短波二维波达方向估计数字处理模块,其特征在于:所述短波二维波达方向估计数字处理模块由程序存储芯片EPCS16向FPGA芯片加载程序。
4.根据权利要求1所述的一种短波二维波达方向估计数字处理模块,其特征在于:所述二维波达方向估计运算采用基于双平行线性阵列的矩阵法,具体步骤如下:
设有K个短波信号从远场入射到平行阵列,其中K(K≤5),信号s(t)=[s1(t),…,sK(t)];阵元1-5收到的信号为X(t)=[x1(t),…,x5(t)],阵元6-10收到的信号计为Y(t)=[y6(t),…,y10(t)],Nx(t)和Ny(t)分别是阵元1-5和阵元6-10的噪声矢量,A为方向矢量,φ为KхK的对角矩阵;如下式所示:
X(t)=As(t)+Nx(t),
Y(t)=Aφs(t)+Ny(t);
先计算X(t)的自协方差矩阵:
Rxx=E[X(t)XH(t)];
再计算X(t)和Y(t)的互协方差矩阵:
Ryx=E[Y(t)XH(t)];
求取波达方向矩阵为:
其中#表示伪逆;接着对波达方向矩阵R进行特征分解,可以计算出其第k(k=1,…,K)个信号的非零特征值σk和相应的特征向量uk,进而可以求得方向角:
其中Ω0为频率,d为相邻阵元的间隔,D为两个平行阵列的间隔,最后可根据αk、βk,计算出俯仰角θk和方位角φk:

说明书全文

一种短波二维波达方向估计数字处理模

技术领域

[0001] 本发明属于电子电路设计,具体涉及一种短波二维波达方向估计数字处理模块。

背景技术

[0002] 智能天线作为移动通信系统中的核心技术,可以有效地对抗干扰信号,提高通信质量和通信的可靠性。短波通信由于存在时变衰落和干扰信号多的情况,常常接收信噪比低、通信效果差。将智能天线技术应用到短波通信中,可以改善通信的效果。波达方向估计是智能天线中的关键技术,通过判断信号的来波方向,自动调整发射和接收方向图,使它在干扰方向形成零陷,将干扰信号抵消,而在有用信号方向形成主波束,以达到提高信号增益、抑制干扰的效果。在实际的短波通信环境中,波达方向是二维的,需要采用二维波达方向估计算法估计来波的俯仰
[0003] 以往的数字处理模块常采用DSP+CPLD或DSP+FPGA的电路组合形式,DSP芯片种类很多,各芯片的语言也有所不同,给软件移植带来诸多不便。

发明内容

[0004] 针对现有技术存在的问题,本发明提供一种短波二维波达方向估计数字处理模块,采用处理能较强的OMAP芯片和FPGA芯片,实现对短波信号进行二维波达方向估计,具有良好的软件可移植性;通过对10路短波通道的中频信号采样和处理,采用新型的电路结构和高效的算法,可以快速准确的识别出短波信号二维波达方向;同时采用CPCI总线,通用性强。
[0005] 本发明的技术方案是:一种短波二维波达方向估计数字处理模块,包括10片A/D转换芯片、3片数字下变频芯片、FPGA芯片以及OMAP芯片;
[0006] 所述A/D转换芯片将采集到的模拟信号转换为数字信号、并将转换后的数字信号分别送给3片数字下变频芯片;
[0007] 3片数字下变频芯片将接受到的数字信号进行数字下变频处理、并将处理后的10路变频信号送到FPGA芯片;
[0008] FPGA芯片对接受到的10路变频信号进行滤波和二维波达方向估计运算,将估计得到的二维波达的俯仰角和方位角信息送给OMAP芯片;
[0009] OMAP芯片将FPGA高速串口送来的二维波达的俯仰角、方位角数据转换成并行总线数据送给CPCI接口芯片。
[0010] 进一步的,所述短波二维波达方向估计数字处理模块采用CPCI总线形式。
[0011] 进一步的,所述短波二维波达方向估计数字处理模块由程序存储芯片EPCS16向FPGA芯片加载程序。
[0012] 进一步的,所述二维波达方向估计运算基于双平行线性阵列的矩阵法,具体步骤如下:
[0013] 设有K个短波信号从远场入射到平行阵列,其中K(K≤5),信号s(t)=[s1(t),…,sK(t)];阵元1-5收到的信号为X(t)=[x1(t),…,x5(t)],阵元6-10收到的信号计为Y(t)=[y6(t),…,y10(t)],Nx(t)和Ny(t)分别是阵元1-5和阵元6-10的噪声矢量,A为方向矢量,φ为KхK的对角矩阵;如下式所示:
[0014] X(t)=As(t)+Nx(t),
[0015] Y(t)=Aφs(t)+Ny(t);
[0016] 先计算X(t)的自协方差矩阵:
[0017] Rxx=E[X(t)XH(t)];
[0018] 再计算X(t)和Y(t)的互协方差矩阵:
[0019] Ryx=E[Y(t)XH(t)];
[0020] 求取波达方向矩阵为:
[0021]
[0022] 其中#表示伪逆;接着对波达方向矩阵R进行特征分解,可以计算出其第k(k=1,…,K)个信号的非零特征值σk和相应的特征向量uk,进而可以求得方向角:
[0023]
[0024]
[0025] 其中Ω0为频率,d为相邻阵元的间隔,最后可根据αk、βk,计算出俯仰角θk和方位角φk:
[0026]
[0027]
[0028] 进一步的,所述二维波达方向估计运算采用基于L型阵列的联合奇异值分解法,具体步骤如下:
[0029] 设有K个短波信号从远场入射到平行阵列,其中K(K≤5),信号s(t)=[s1(t),…,sK(t)];阵元1-5收到的信号为X(t),阵元6-10收到的信号为Z(t),Nx(t)和Nz(t)分别是阵元1-5和阵元6-10的噪声矢量,A(α)、A(θ),如下式:
[0030] X(t)=A(α)s(t)+Nx(t),
[0031] Z(t)=A(θ)s(t)+Nz(t);
[0032] 其中:
[0033] A(α)=[α(φ1),α(φ2),...,α(φk),...,α(φK)];
[0034] A(θ)=[α(θ1),α(θ2),...,α(θk),...,α(θK)],
[0035]
[0036]
[0037] 计算X(t)和Z(t)的互协方差矩阵:
[0038] Rxz=E[X(t)ZH(t)];
[0039] 将Rxz分为两个5×4维的子矩阵:
[0040] Rxz,1=Rxz(:1,:4),
[0041] Rxz,2=Rxz(:,2:5);
[0042] 将子矩阵Rxz,1和Rxz,2组成新的矩阵R:
[0043]
[0044] 对矩阵R进行奇异值分解,可以计算出左奇异矢量U1,将U1分为两个5×K维子阵U11和U12:
[0045]
[0046] 求取 #表示逆阵;然后可以计算出其第k(k=1,…,K)个特征值γk,从而可求得俯仰角和方位角:
[0047] 俯仰角:
[0048]
[0049]
[0050] 方位角:
[0051] 本发明的有益效果是:
[0052] (1)本发明采用OMAP+FPGA的电路形式来实现短波二维波达方向估计,随着FPGA的发展,FPGA具有强大的并行处理能力和拥有丰富的IP核,可以实现各种数字信号处理算法。同时新型的OMAP芯片具有超强的处理器性能,可以实现ARM+DSP的处理能力,并且开发语言统一,便于软件移植、软件继承性好。
[0053] (2)本发明短波二维波达方向估计数字处理模块采用CPCI总线形式,插槽标准,传输速率高,通信可靠性好。
[0054] (3)在短波通信环境中,存在着多径传输,容易产生相干信号,对波达方向估计造成识别模糊,本发明使用基于双平行线性阵列的矩阵法、基于L型阵列的联合奇异值分解法等来估计相干信号的波达方向,可以有效消除识别模糊。附图说明
[0055] 图1为短波二维波达方向估计数字处理模块电路示意图;
[0056] 图2为采用矩阵法的双平行阵列示意图;
[0057] 图3为采用联合奇异值分解法的L型阵列示意图。

具体实施方式

[0058] 下面结构附图对本发明作进一步的说明。
[0059] 如图1所示,短波二维波达方向估计数字处理模块可以接收处理10个短波通道中频信号,采用OMAP3530作为模块中的控制芯片,该芯片内含ARM Cortex-A8内核和TMS320C64x系列的DSP处理器,完全软件兼容ARM9和C64x,有利于平台移植,并采用TI的SmartReflex技术实现低功耗。FPGA采用Altera公司的Cyclone II FPGA,运算速度快,A/D选用10片AD9244,选用3片下变频芯片HSP50216,FPGA芯片作逻辑控制处理和二维波达方向估计运算,CPCI接口芯片选用PCI9054。
[0060] OMAP芯片通过CPCI总线读入外部控制命令,同时OMAP芯片又将估计得到的二维数据和状态信息回传给CPCI总线。电路通过10片AD9244分别对10路短波中频信号进行采样,将采集到的模拟信号转换为数字信号、并将转换后的数字信号分别送给3片HSP50216芯片(每片可支持四个通道的数据处理)做数字下变频处理;然后将处理后的10路变频信号送到FPGA芯片,通过FPGA芯片对10路信号作滤波和二维波达方向估计运算,FPGA芯片同时负责模块的内部逻辑控制。最后FPGA将估计得到的俯仰角和方位角信息通过CPCI总线送给OMAP芯片,由该芯片通过CPCI总线传给外部主控模块。FPGA程序存储在EPCS16芯片中,程序由Quartus II编译生成,OMAP芯片运行Windows操作系统
[0061] 短波二维波达方向估计数字处理模块可以通过CPCI总线加载不同的程序,实现不同的二维波达方向估计算法。在短波通信环境中,存在着多径传输,容易产生相干信号,对波达方向估计造成识别模糊。本发明使用基于双平行线性阵列的矩阵法、基于L型阵列的联合奇异值分解法等来估计相干信号的波达方向,可以有效消除识别模糊。
[0062] 如图2所为采用矩阵法的双平行阵列示意图,基于双平行线性阵列的矩阵法的具体步骤如下:
[0063] 设有K个短波信号从远场入射到平行阵列,其中K(K≤5),信号s(t)=[s1(t),…,sK(t)];阵元1-5收到的信号为X(t)=[x1(t),…,x5(t)],阵元6-10收到的信号计为Y(t)=[y6(t),…,y10(t)],Nx(t)和Ny(t)分别是阵元1-5和阵元6-10的噪声矢量,A为方向矢量,φ为KхK的对角矩阵;如下式所示:
[0064] X(t)=As(t)+Nx(t),
[0065] Y(t)=Aφs(t)+Ny(t);
[0066] 先计算X(t)的自协方差矩阵:
[0067] Rxx=E[X(t)XH(t)];
[0068] 再计算X(t)和Y(t)的互协方差矩阵:
[0069] Ryx=E[Y(t)XH(t)];
[0070] 求取波达方向矩阵为:
[0071]
[0072] 其中#表示伪逆;接着对波达方向矩阵R进行特征分解,可以计算出其第k(k=1,…,K)个信号的非零特征值σk和相应的特征向量uk,进而可以求得方向角:
[0073]
[0074]
[0075] 其中Ω0为频率,d为相邻阵元的间隔,最后可根据αk、βk,计算出俯仰角θk和方位角φk:
[0076]
[0077]
[0078] 如图3所为采用联合奇异值分解法的L型阵列示意图,基于L型阵列的联合奇异值分解法的具体步骤如下:
[0079] 设有K个短波信号从远场入射到平行阵列,其中K(K≤5),信号s(t)=[s1(t),…,sK(t)];阵元1-5收到的信号为X(t),阵元6-10收到的信号为Z(t),Nx(t)和Nz(t)分别是阵元1-5和阵元6-10的噪声矢量,A(α)、A(θ),如下式:
[0080] X(t)=A(α)s(t)+Nx(t),
[0081] Z(t)=A(θ)s(t)+Nz(t);
[0082] 其中:
[0083] A(α)=[α(φ1),α(φ2),...,α(φk),...,α(φK)],
[0084] A(θ)=[α(θ1),α(θ2),...,α(θk),...,α(θK)],
[0085]
[0086]
[0087] 计算X(t)和Z(t)的互协方差矩阵:
[0088] Rxz=E[X(t)ZH(t)];
[0089] 将Rxz分为两个5×4维的子矩阵:
[0090] Rxz,1=Rxz(:1,:4),
[0091] Rxz,2=Rxz(:,2:5);
[0092] 将子矩阵Rxz,1和Rxz,2组成新的矩阵R:
[0093]
[0094] 对矩阵R进行奇异值分解,可以计算出左奇异矢量U1,将U1分为两个5×K维子阵U11和U12:
[0095]
[0096] 求取 #表示逆阵;然后可以计算出其第k(k=1,…,K)个特征值γk,从而可求得俯仰角和方位角:
[0097] 俯仰角:
[0098]
[0099]
[0100] 方位角:
[0101] 本实施例通过对10路短波通道的中频信号采样和处理,采用新型的电路结构和高效的算法,可以识别出短波信号二维波达方向;该模块同时具有良好的软件可移植性和通用性。
[0102] 虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈