CN109284528B - 基于流水线设计的相控阵天线波束控制码计算方法及系统 - Google Patents
基于流水线设计的相控阵天线波束控制码计算方法及系统 Download PDFInfo
- Publication number
- CN109284528B CN109284528B CN201810850177.9A CN201810850177A CN109284528B CN 109284528 B CN109284528 B CN 109284528B CN 201810850177 A CN201810850177 A CN 201810850177A CN 109284528 B CN109284528 B CN 109284528B
- Authority
- CN
- China
- Prior art keywords
- phase
- address
- control code
- frequency
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 79
- 238000013461 design Methods 0.000 title claims abstract description 14
- 230000006870 function Effects 0.000 claims abstract description 29
- 238000013500 data storage Methods 0.000 claims description 36
- 238000013139 quantization Methods 0.000 claims description 23
- 238000000034 method Methods 0.000 claims description 18
- 238000012937 correction Methods 0.000 claims description 14
- 230000003111 delayed effect Effects 0.000 claims description 5
- 238000004891 communication Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明以实现瓦片式相控阵天线扫描为目的,提出一种基于流水线设计的相控阵天线波束控制码计算方法,该计算方法包括指向误差修正模块、函数计算模块、地址计算模块、查表模块,其中地址计算模块采用13级流水线设计,适合采用FPGA实现,计算速度快,且有效降低了计算冗余和资源消耗。本发明实现了大型相控阵天线波束控制码的高效计算,功能全面,可扩展性强。
Description
技术领域
本发明涉及一种在使用VM调制器芯片的相控阵天线系统中对波束控制码 的计算方法,具体是一种基于流水线设计的相控阵天线波束控制码计算方法及 系统。
背景技术
相控阵天线是由许多相同的辐射单元按一定规律排布构成的阵列天线。相 控阵天线通过特定的馈电方式,智能化地控制和调整阵列每个单元的幅度和相 位,使波束赋形和无惯性地灵活扫描。
相控阵天线的一个重要应用是星间通信,星间通信距离长,链路损耗大, 随着通信速率需求的不断增长,发展大型Ka频段相控阵是必然的趋势。波控系 统主要是在嵌入式处理芯片的控制下,产生需要的波束指向控制码,并送至每 个单元组件,完成波束的定向控制。在某些应用场景中,比如Ka频段运载中继 测控,天线阵面规模大,波束窄,扫描角速度较大,对波束的指向精度和计算 速度提出了更高的要求。波束控制码的计算是整个波控系统的核心,计算的准 确性与效率直接决定波束指向的成败。
在相位和控制码地址计算时,有很多个步骤,包括理论相位计算、相位补 偿、温度补偿、幅度补偿和频率补偿,通常认为这些操作具有相对独立性,设 计师容易想到分模块进行计算,而不是统一起来采用流水线计算。分开来计算 虽然看上去直观,但是会产生几个缺点:第一,模块与模块之间的通信增加了 资源开销;第二,后面的模块需要等待前面的模块计算完才开始计算,整体计 算速度受影响;第三,采用传统的量化方法计算,得到的相位值有正有负,阵 面规模较大时,最大相位有几千度,除了数据位宽增加,还需增加求余运算模 块,不仅增加了设计难度,而且增加了资源消耗。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种波束控制码的 高效计算方法,实现了大规模相控阵天线波控码的计算,同时兼顾准确性、灵 活性、计算速度和资源消耗。
本发明的技术方案如下:
波束控制码计算包括:指向误差修正模块、三角函数计算模块、地址计算 模块、查表模块。
该模块采用一种基于FPGA的流水线计算方法,流水线共分13级,最大延 迟13个时钟周期。包括理论相位计算、相位补偿、温度补偿、频率首地址计算、 通道号首地址计算、幅度等级计算、幅度等级首地址计算、相位等级地址计算、 量化方式转化、控制码地址写入。
量化数据存储地址按下式计算:
量化数据存储地址=频率首地址+通道号首地址+幅度等级首地址+相位等级 地址。
其中,
频率首地址=每个频率的字节数/2*频率等级=921600*频率等级;
通道号首地址=每个通道的字节数/2*通道编号=3600*通道编号;
幅度等级首地址=幅度等级*360;
相位等级地址=相位等级*2。
频率等级与输入频率编号一一对应。
幅度等级是输入频率、通道号的函数,存储在幅度补偿数据存储表中。
相位等级由通过计算得到的相位确定,相位计算公式如下
相位=理论相位+相位补偿+温度补偿
其中理论相位计算公式如下:
式中x和y分别是阵元的横坐标和纵坐标,λ是波长,λ=c/f,其中c是光速, f是频率。
相位补偿数据是输入频率及通道号的函数,通过整机装配后的通道校准得 到,存储在相位补偿数据存储表中。
温度补偿数据是输入温度及通道号的函数,存储在温度补偿数据存储表中。
相位计算完成之后,需要乘以一个常数k,再截位,目的是将按整数次幂等 分量化的相位转化为180个相位等级,进而求出相位地址。
相位地址、幅度等级首地址、频率等级首地址三者相加之后得到量化数据 存储地址,然后依次存入RAM中,供后续模块使用。
查表模块完成波束控制码的flash读取。输入为存储在RAM中的波控码flash 存储地址,输出为波控码数据。
本发明的有益效果如下:
(1)控制码是离轴角、旋转角、阵元坐标、频率、温度、幅度补偿数据存 储表、位补偿数据存储表、温度补偿数据存储表等输入变量的函数,自由度高, 功能全面,兼容性和可扩展性强。
(2)地址计算模块采用13级流水线计算方法,只消耗9个加法器和6个 乘法器,从第一个数据输入到最后一个数据输出仅隔N+13个时钟周期,其中N 为阵元个数。若以40M时钟计算速度来算,整个计算时间仅为6.725us(N=256), 计算效率非常高。
(3)相位计算、相位补偿、温度补偿等数据采用2的整数次幂等分量化法, 前面采用不低于16位的位宽进行计算,相位计算全部完成后再截去高位和低位, 相当于求余和取整计算,保证了计算精度且无需再进行求余操作,降低了计算 复杂度和资源消耗。
附图说明
图1为本发明模块组成框图;
图2为本发明地址计算模块流水线计算框图;
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明包括指向误差修正模块、三角函数计算模块、地址计算模块和查表 模块。
其中mod(X,2N)表示X对2N求余,round()表示就近取整,N表示量化位宽取 N=16。
三角函数量化值=round(三角函数值*(2(N-1)-1))
其中round()表示就近取整,N是量化位宽,取N=16。本发明所述地址计算 模块完成控制码的存储地址计算,该模块输入为三角函数和sinθ,阵 元坐标,频率,温度,幅度补偿数据存储表,相位补偿数据存储表,温度补偿 数据存储表,输出全部通道的量化数据在flash中的存储地址。
量化数据在flash中的存储地址是输入参数的函数,根据量化数据存储格式 表,量化数据共分12个频率等级,256个通道,10个幅度等级,180个相位等 级,I、Q两路,每路2字节数据,一共12*256*10*180*2*2=22118400字节数据, flash存储位宽为2字节。
量化数据存储地址按下式计算:
量化数据存储地址=频率首地址+通道号首地址+幅度等级首地址+相位等级 地址。
其中,
频率首地址=每个频率的字节数/2*频率等级=921600*频率等级;
通道号首地址=每个通道的字节数/2*通道编号=3600*通道编号;
幅度等级首地址=幅度等级*360;
相位等级地址=相位等级*2。
频率等级与输入频率编号一一对应。
幅度等级是输入频率、通道号的函数,存储在幅度补偿数据存储表中,存 储格式见幅度补偿数据存储格式表,共256*12=3072个数据,每个数据位宽为4。
相位等级由通过计算得到的相位确定,相位计算公式如下
相位=理论相位+相位补偿+温度补偿
其中理论相位计算公式如下:
式中x和y分别是阵元的横坐标和纵坐标,λ是波长,λ=c/f,其中c是光速,f是频率。坐标x和y与其中一个频点的波长λ相除后一起进行量化,量化公式如 下(以x为例):
其中mod(X,2N)表示X对2N求余,round()表示就近取整,N表示量化位宽,取 N=16。得到的16位二进制值的最高位表示符号,第高2位到高4位表示‘整数’, 低12位表示‘小数’。
相位补偿数据是输入频率及通道号的函数,通过整机装配后的通道校准得 到,存储在相位补偿数据存储表中,存储格式见相位补偿数据存储格式表,共 12*256=3072个数据,每个数据位宽为16。先计算补偿相位存储地址,再查表 即可得出,相位补偿数据的量化方法如下:
其中mod(X,2N)表示X对2N求余,round()表示就近取整,N是量化位宽, 取N=16。
温度补偿数据是输入温度及通道号的函数,存储在温度补偿数据存储表中, 存储格式见温度补偿数据存储格式表,共16个温度等级。温度补偿数据量化方 法与相位补偿数据量化方法一致。
相位计算完成之后,需要乘以一个常数k,再截位,目的是将按整数次幂等 分量化的相位转化为180个相位等级,进而求出相位地址。
相位地址、幅度等级首地址、频率等级首地址三者相加之后得到量化数据 存储地址,然后依次存入RAM中,供后续模块使用。
地址计算模块采用一种基于FPGA的全流水线详细计算方法,框图如图2。 流水线共分13级,最大延迟13个时钟周期。
第1级:
a)根据输入的阵元坐标地址(0~255),查表得到16位的x坐标值和y 坐标值。
第2级:
c)读取当前温度数据。
第3级:
a)根据温度数据判决得到相应温度补偿数据表中的温度等级首地址;
b)取通道编号1(0~255,相比阵元坐标地址延迟3个时钟周期);
c)根据当前频率计算出频率系数,频率系数是由于频率不同导致计算结 果与根据参考频率计算的结果之比,为保精度,该系数用17位表示, 包括1个‘整数’位和16个‘小数’位,预先由matlab算出;
第4级:
a)将温度等级首地址与通道号2相加得到温度补偿地址;
b)将频率系数与sinθ相乘,截去低16位,高位保留;
d)根据当前频率计算相位补偿频率首地址;
e)取通道号2(比通道号1延迟1个时钟周期);
第5级:
a)根据温度补偿地址查表得到温度补偿数据;
c)将相位补偿频率首地址与通道号2相加得到相位补偿地址;
d)根据当前频率计算出幅度补偿频率首地址;
e)取通道号3(比通道号1延迟2个时钟周期)。
第6级:
a)将16位温度补偿数据与16位理论相位相加,输出16位相位数据, 高位自动溢出。
b)根据相位补偿地址查表得到16位相位补偿数据。
c)将幅度补偿频率首地址与通道号相加得到幅度等级在幅度补偿数据 存储表中的地址。
第7级:
a)将第6级计算出的16位相位数据与16位相位补偿数据相加,输出16 位相位数据,高位自动溢出;
b)根据幅度等级在幅度补偿数据存储表中的地址查表得到幅度等级。 第8级:
a)将前一级计算的相位与系数K相乘,目的是将按本文所提量化方法到 按相位度数的整数倍量化方法之间的转换。
b)根据当前频率计算频率首地址;
c)根据幅度等级计算幅度首地址;
d)取通道号4(比通道号1延迟5个时钟周期);
第9级:
a)将与系数K相乘后的数据进行截位,使得其值在在0~180之间。
b)将频率补偿首地址与幅度补偿首地址相加;
c)将通道号与3600相乘,得到通道号首地址;
第10级:
a)将第9级a)步骤计算所得的数值乘以2(左移1位),且360用0代 替,得到范围在(0,2,4,…358)的相位地址。
b)将第9级b)和c)结果相加;
第11级:
a)将第10级a)和b)计算结果相加得到控制码在flash中的存储地址。
b)计算控制码flash存储地址的RAM地址。
第12级:
a)将控制码flash存储地址写入RAM中。
第13级:
a)得到已经存储完毕的控制码flash存储地址。
从t=0时刻开始从输入端依次输入N个阵元的坐标,输出端延迟13个时钟 周期后依次输出N个阵元对应的控制码地址,从第一个数据输入到最后一个数 据输出仅隔N+13个时钟周期,其中N为阵元个数。若以40M时钟计算速度来 算,整个计算时间仅为6.725us(N=256),计算效率非常高。该模块一共用了9 个加法器和6个乘法器,无需除法运算和求余运算,适合用FPGA进行定点计 算,在保证精度和速度的前提下,资源消耗大大降低。
查表模块完成波束控制码的flash读取。输入为存储在RAM中的波控码flash 存储地址,输出为波控码数据。
需要说明的是,上文只是对本发明进行示意性说明和阐述,本领域的技术 人员应当明白,对本发明的任意修改和替换都属于本发明的保护范围。
Claims (9)
1.一种基于流水线设计的相控阵天线波束控制码计算方法,其特征在于,包括以下步骤:
步骤3,地址计算模块完成控制码的存储地址计算,输入为三角函数 和sinθ,阵元坐标,频率,温度,幅度补偿数据存储表,相位补偿数据存储表,温度补偿数据存储表,输出全部通道的量化数据在flash中的存储地址;
步骤4,查表模块完成波束控制码的flash读取,输入为存储在RAM中的波控码flash存储地址,输出为波控码数据;
所述地址计算模块采用基于FPGA的流水线计算方法,流水线共分13级,最大延迟13个时钟周期。
3.如权利要求2所述的一种基于流水线设计的相控阵天线波束控制码计算方法,其特征在于,地址计算模块中量化数据存储地址按下式计算:
量化数据存储地址=频率首地址+通道号首地址+幅度等级首地址+相位等级地址;
其中,
频率首地址=每个频率的字节数/2*频率等级=921600*频率等级;
通道号首地址=每个通道的字节数/2*通道编号=3600*通道编号;
幅度等级首地址=幅度等级*360;
相位等级地址=相位等级*2;
频率等级与输入频率编号一一对应;
幅度等级是输入频率、通道号的函数,存储在幅度补偿数据存储表中。
7.如权利要求4所述的一种基于流水线设计的相控阵天线波束控制码计算方法,其特征在于,相位计算完成之后,乘以常数k,再截位,将按整数次幂等分量化的相位转化为180个相位等级,进而求出相位地址。
8.如权利要求3所述的一种基于流水线设计的相控阵天线波束控制码计算方法,其特征在于,将相位地址、幅度等级首地址、频率等级首地址三者相加之后得到量化数据存储地址,然后依次存入RAM中,供查表模块使用。
9.一种基于流水线设计的相控阵天线波束控制码计算系统,其特征在于,包括指向误差修正模块、三角函数计算模块、地址计算模块和查表模块:
查表模块完成波束控制码的flash读取,输入为存储在RAM中的波控码flash存储地址,输出为波控码数据;所述地址计算模块采用基于FPGA的流水线计算方法,流水线共分13级,最大延迟13个时钟周期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810850177.9A CN109284528B (zh) | 2018-07-28 | 2018-07-28 | 基于流水线设计的相控阵天线波束控制码计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810850177.9A CN109284528B (zh) | 2018-07-28 | 2018-07-28 | 基于流水线设计的相控阵天线波束控制码计算方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109284528A CN109284528A (zh) | 2019-01-29 |
CN109284528B true CN109284528B (zh) | 2023-06-09 |
Family
ID=65182857
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810850177.9A Active CN109284528B (zh) | 2018-07-28 | 2018-07-28 | 基于流水线设计的相控阵天线波束控制码计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109284528B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110068809A (zh) * | 2019-05-06 | 2019-07-30 | 成都泰格微电子研究所有限责任公司 | 一种用于相控阵雷达的波控处理模块 |
CN110068811A (zh) * | 2019-05-06 | 2019-07-30 | 成都泰格微电子研究所有限责任公司 | 一种相控阵雷达波束控制系统 |
CN117233716B (zh) * | 2023-11-16 | 2024-03-15 | 中科亿海微电子科技(苏州)有限公司 | 一种星载相控阵天线波束控制及测试一体化方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101858972A (zh) * | 2010-03-23 | 2010-10-13 | 深圳市蓝韵实业有限公司 | 基于延时参数实时计算和流水线的多波束合成方法和装置 |
CN103810319A (zh) * | 2013-10-31 | 2014-05-21 | 西安电子工程研究所 | 一种基于fpga的相控阵天线迭代配相算法的实现方法 |
CN107590085A (zh) * | 2017-08-18 | 2018-01-16 | 浙江大学 | 一种具有多级缓存的动态可重构阵列数据通路及其控制方法 |
-
2018
- 2018-07-28 CN CN201810850177.9A patent/CN109284528B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101858972A (zh) * | 2010-03-23 | 2010-10-13 | 深圳市蓝韵实业有限公司 | 基于延时参数实时计算和流水线的多波束合成方法和装置 |
CN103810319A (zh) * | 2013-10-31 | 2014-05-21 | 西安电子工程研究所 | 一种基于fpga的相控阵天线迭代配相算法的实现方法 |
CN107590085A (zh) * | 2017-08-18 | 2018-01-16 | 浙江大学 | 一种具有多级缓存的动态可重构阵列数据通路及其控制方法 |
Non-Patent Citations (2)
Title |
---|
基于并发流水线的直升机雷达多波束调度管理研究;陈杰等;《火控雷达技术》;20120925;第41卷(第03期);全文 * |
基于并行技术和流水线的LMS自适应滤波算法;杜秀丽等;《数据采集与处理》;20170315;第32卷(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109284528A (zh) | 2019-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109284528B (zh) | 基于流水线设计的相控阵天线波束控制码计算方法及系统 | |
CN109541550B (zh) | 一种相控阵天线快速波束切换的方法 | |
CN111859269B (zh) | 相控阵天线波控角度解算方法、系统及存储介质 | |
CN111144014A (zh) | 一种基于fpga的二维平面相控阵天线配相算法实现方法 | |
CN110488228B (zh) | 线性调频信号生成方法、装置及存储介质 | |
CN109521992B (zh) | 一种无乘法器的基于cordic算法的线性调频信号生成方法 | |
JP2019140594A (ja) | アレーアンテナ装置及びビーム形成方法 | |
CN111967608A (zh) | 数据处理方法、装置、设备及存储介质 | |
CN109558108B (zh) | 基于fpga的相控阵雷达移相器方法 | |
WO2022170811A1 (zh) | 一种适用于混合精度神经网络的定点乘加运算单元及方法 | |
WO2011100387A1 (en) | Rom-based direct digital synthesizer with pipeline delay circuit | |
CN114610267A (zh) | 一种基于指数函数和softmax函数的优化方法、硬件系统及芯片 | |
CN111325332A (zh) | 卷积神经网络的处理方法和装置 | |
CN112382855B (zh) | 基于稀疏阵列的有源相控阵天线波束控制方法 | |
CN109787670A (zh) | 一种可绘制方向图的阵列天线波束赋形系统 | |
US10938404B1 (en) | Digital-to-analog converter, transmitter, base station and mobile device | |
CN113138749A (zh) | 基于cordic算法的三角函数计算装置及方法 | |
US7676529B2 (en) | System and method for efficient rectangular to polar signal conversion using cordic algorithm | |
CN110458277B (zh) | 适用于深度学习硬件加速器的可配置精度的卷积硬件结构 | |
CN116719499A (zh) | 一种应用于5g最小二乘定位的自适应伪逆计算方法 | |
CN110676588B (zh) | 一种移相波控码硬件在线调节系统 | |
CN102566963A (zh) | 一种fpga进行数据处理的方法 | |
US11531896B2 (en) | Neural network circuit providing for operations with bit shifting and rounded values of weight information | |
CN115826674A (zh) | 毫赫兹级别频率步进的时钟信号发生器实现方法及系统 | |
CN113759777A (zh) | 相控阵天线波束控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |