CN112748401A - 一种实时权值产生方法 - Google Patents

一种实时权值产生方法 Download PDF

Info

Publication number
CN112748401A
CN112748401A CN202011585142.0A CN202011585142A CN112748401A CN 112748401 A CN112748401 A CN 112748401A CN 202011585142 A CN202011585142 A CN 202011585142A CN 112748401 A CN112748401 A CN 112748401A
Authority
CN
China
Prior art keywords
phase
weight
weight phase
value
cosval
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.)
Granted
Application number
CN202011585142.0A
Other languages
English (en)
Other versions
CN112748401B (zh
Inventor
杨元浩
王宇
杜江
王伟
王沛
禹卫东
邓云凯
赵庆超
陈圳
张华春
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Chinese Academy of Sciences
Aerospace Information Research Institute of CAS
Original Assignee
University of Chinese Academy of Sciences
Aerospace Information Research Institute of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Chinese Academy of Sciences, Aerospace Information Research Institute of CAS filed Critical University of Chinese Academy of Sciences
Priority to CN202011585142.0A priority Critical patent/CN112748401B/zh
Publication of CN112748401A publication Critical patent/CN112748401A/zh
Application granted granted Critical
Publication of CN112748401B publication Critical patent/CN112748401B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00—Programme-control systems
    • G05B19/02—Programme-control systems electric
    • G05B19/04—Programme control other than numerical control, i.e. in sequence controllers or logic controllers
    • G05B19/042—Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本公开提供了一种实时权值产生方法，应用于基于FPGA的中频数字波束形成过程，包括：获取外部辅助数据，基于辅助数据，计算权重相位的相位生成余弦值cosval和相位生成正弦值sinval，基于相位生成余弦值cosval和相位生成正弦值sinval计算权重相位，将权重相位输入vivado中的cordic ip核，计算权值。本公开提供的方法能准确、稳定地实时计算出基于FPGA的中频数字波束形成过程中各通道需要的数字合成权值。

Description

一种实时权值产生方法
技术领域
本公开涉及数字信号处理技术领域，尤其涉及一种基于FPGA的实时权值产生方法。
背景技术
俯仰向数字波束形成(Digital Beam Forming，简称DBF)技术是随着数字信号处理方法的发展而建立起来的一门新技术，这种技术不仅能充分保留阵列天线上收集的信息，而且能利用复杂的数字信号处理方法对信号进行处理，在基带数字端运用数字端运用数字信号处理技术形成灵活的发射和接收波束输出。
目前，主要有传统的低频数字波束形成和改进的中频数字波束形成的方法。基于FPGA的中频数字波束形成的硬件实现还未成熟开发，实时权值生成器是基于FPGA的中频数字波束形成的核心组件，基于FPGA的中频数字波束形成时钟频率较高，且计算复杂度高，快速稳定的实时权值生成器是基于FPGA实现中频数字波束形成的必要保障。
综上所述，该方案能够提供一种基于FPGA的中频数字波束形成的实时权值产生器，保障基于FPGA的中频数字波束形成的实现。
发明内容
鉴于上述问题，本发明提供了一种基于FPGA的实时权值产生方法，应用于基于FPGA的中频数字波束形成过程，以解决上述技术问题。
本公开的一个方面提供了一种实时权值产生方法，应用于基于FPGA的中频数字波束形成过程，包括：获取外部辅助数据；基于所述辅助数据，计算权重相位的相位生成余弦值cosval和相位生成正弦值sinval；基于所述相位生成余弦值cosval和相位生成正弦值sinval计算所述权重相位；将所述权重相位输入vivado中的cordic ip核，计算权值。
可选地，所述辅助数据为天线雷达发射中频信号的参数信息，至少包括：所述FPGA的中心通道与地面目标的距离、总通道数、所述天线雷达对所述中频信号的采样频率、光速、所述天线雷达的轨道高度、地球半径、所述FPGA中每个通道相对于其中心通道的延迟、所述中频信号从发射到反射回所述天线雷达的最短时间、时宽、所述中频信号的波长、所述FPGA编号为ic的通道与其中心通道的距离、所述天线雷达的安装角。
可选地，所述基于所述辅助数据，计算权重相位的相位生成余弦值cosval和相位生成正弦值sinval包括：
SRval＝(0：(Nrg-1))*A1+A2(ic)；
cosval1＝SRval/A4；
cosval2＝A5./SRval；
cosval＝cosval1+cosval2；
sinval1＝1-cosval；
sinval2＝1+cosval；
sinval＝sqrt(1-cosvaj2)；
其中，A1、A2、A3、A4、A5分别表示第一辅助参数、第二辅助参数、第三辅助参数、第四辅助参数和第五辅助参数，根据所述辅助数据计算得到，SRval表示所述FPGA的中心通道与地面目标的距离，Nrg表示总的通道数，cosval1、cosval2表示所述相位生成余弦值cosval的分量，sinval1、sinval2表示所述相位生成正弦值sinval的分量。
可选地，所述第一辅助参数、第二辅助参数、第三辅助参数、第四辅助参数和第五辅助参数的计算公式如下：
Al＝1/Fs*C/2；
A2(ic)＝(tmin-delay(ic))*C/2：
A3＝(Height+Re)2-Re2；
A4＝2*(Re+Height)；
A5＝A3/A4；
Figure BDA0002864879570000031
其中，Fs为所述天线雷达对所述中频信号的采样频率，C为光速，Height为所述天线雷达的轨道高度，Re为地球半径，delay(ic)表示所述FPGA中每个通道相对于其中心通道的延迟，ic表示通道编号，Tr表示所述中频信号从发射到反射回所述天线雷达的最短时间，τ表示时宽。
可选地，得到所述相位生成余弦值cosval和相位生成正弦值sinval后，将所述相位生成余弦值cosval和相位生成正弦值sinval时钟同步。
可选地，所述基于所述相位生成余弦值cosval和相位生成正弦值sinval计算所述权重相位包括：
A6(ic)＝-2/lamda*chandistance1(ic)*cos(-lookbias*pi/180)；
A7(ic)＝-2/lamda*chandistance1(ic)*sin(-lookbias*pi/180)；
wt＝sinval*A6(ic)+cosval*A7(ic)；
其中，lamda表示所述中频信号的波长，chandistancel(ic)表示所述FPGA编号为ic的通道与其中心通道的距离，lookbias是所述天线雷达的安装角，A6、A7分别为第六辅助参数和第七辅助参数。
可选地，将所述权重相位输入vivado中的cordic ip核之前，将所述权重相位标准化至[-1，1]的范围内，使所述权重相位满足cordic ip核的输入要求。
可选地，所述将所述权重相位标准化至[-1，1]的范围内包括：将所述权重相位转换为至少有2位小数的定点数；判断所述权值相位的正负和奇偶；根据所述权值相位的正负、奇偶状态，将所述权值相位标准化至[-1，1]的范围内。
可选地，所述判断所述权值相位的正负和奇偶包括：判断所述权重相位的符号位是否为0；若所述符号位为0，则所述权重相位为正；若所述符号位为1，所述权重相位为负；判断所述权重相位的最后一位是否为0；若所述权重相位的最后一位为0，则所述权重相位为偶数；若所述权重相位的最后一位不为0，则所述权重相位为奇数。
可选地，所述根据所述权值相位的正负、奇偶状态，将所述权值相位标准化至[-1，1]的范围内包括：若所述权重相位为正，且所述权重相位为偶数，则所述权重相位的小数部分表示所述权重相位的标准化值；若所述权重相位为正，且所述权重相位为奇数，则所述权重相位为奇数，取所述权重相位的小数部分减1得到的值表示所述权重相位的标准化值；若所述权重相位为负，且所述权重相位为偶数，则所述权重相位为偶数，取所述权重相位的小数部分表示所述权重相位的标准化值；若所述权重相位为负，且所述权重相位为奇数，则所述权重相位为奇数，取所述权重相位的小数部分加1得到的值表示所述权重相位的标准化值。
在本公开实施例采用的上述至少一个技术方案能够达到以下有益效果：
本公开提供了一种实时权值产生方法，应用于基于FPGA的中频数字波束形成过程中所需的权重值，该方法能准确、稳定地实时计算出通道需要的数字合成权值。
附图说明
为了更完整地理解本公开及其优势，现在将参考结合附图的以下描述，其中：
图1示意性示出了本公开实施例提供的一种实时权值产生方法的流程图；
图2示意性示出了本公开实施例提供的一种实时权值生成器的结构示意图；
图3示意性示出了本公开实施例提供的一种基于FPGA的中频数字波束形成方法的结构示意图；
图4示意性示出了本公开实施例提供的MATLAB程序生成的实时权值结果示意图；
图5示意性示出了本公开实施例提供的一种实时权值产生方法生成的实时权值结果示意图。
具体实施方式
以下，将参照附图来描述本公开的实施例。但是应该理解，这些描述只是示例性的，而并非要限制本公开的范围。在下面的详细描述中，为便于解释，阐述了许多具体的细节以提供对本公开实施例的全面理解。然而，明显地，一个或多个实施例在没有这些具体细节的情况下也可以被实施。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本公开的概念。
在此使用的术语仅仅是为了描述具体实施例，而并非意在限制本公开。在此使用的术语“包括”、“包含”等表明了所述特征、步骤、操作和/或部件的存在，但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。
在此使用的所有术语(包括技术和科学术语)具有本领域技术人员通常所理解的含义，除非另外定义。应注意，这里使用的术语应解释为具有与本说明书的上下文相一致的含义，而不应以理想化或过于刻板的方式来解释。
根据中频数字波束形成理论，权值为：
Figure BDA0002864879570000051
其中，θ(t)表示雷达在t时刻接收到的回波方向下视角。
本公开实施例提供了一种实时权值产生方法，在FPGA中实现，权值在FPGA需要分拆成加权系数SRVal、cosval和sinval，并且由此来生成wt1，wt2，最后利用cordic核来计算实时权值，具体计算过程在实施例中进行详细说明。
如图1所示，本公开提供的一种实时权值产生方法，应用于基于FPGA的中频数字波束形成过程，包括步骤S110～S140。
S110，获取外部辅助数据。
在本公开实施例中，辅助数据为天线雷达发射中频信号的参数信息，至少包括：FPGA的中心通道与地面目标的距离、总通道数、天线雷达对中频信号的采样频率、光速、天线雷达的轨道高度、地球半径、FPGA中每个通道相对于其中心通道的延迟、中频信号从发射到反射回天线雷达的最短时间、时宽、中频信号的波长、FPGA编号为ic的通道与其中心通道的距离、天线雷达的安装角。
S120，基于辅助数据，计算权重相位的相位生成余弦值cosval和相位生成正弦值sinval。
在本公开实施例中，S120的具体计算公式包括：
SRval＝(0：(Nrg-1))*A1+A2(ic)；
cosval1＝SRval/A4；
cosval2＝A5./SRval；
cosval＝cosval1+cosval2；
sinval1＝1-cosval；
sinval2＝1+cosval；
sinval＝sqrt(1-cosvaj2)；
其中，A1、A2、A3、A4、A5分别表示第一辅助参数、第二辅助参数、第三辅助参数、第四辅助参数和第五辅助参数，根据辅助数据计算得到，SRval表示FPGA的中心通道与地面目标的距离，Nrg表示总的通道数，cosval1、cosval2表示相位生成余弦值cosval的分量，sinval1、sinval2表示相位生成正弦值sinval的分量。
第一辅助参数A1、第二辅助参数A2、第三辅助参数A3、第四辅助参数A4和第五辅助参数A5的计算公式如下：
A1＝1/Fs*C/2；
A2(ic)＝(tmin-delay(ic))*C/2；
A3＝(Height+Re)2-Re2；
A4＝2*(Re+Height)；
A5＝A3/A4；
Figure BDA0002864879570000071
其中，Fs为天线雷达对中频信号的采样频率，C为光速，Height为天线雷达的轨道高度，Re为地球半径，delay(ic)表示FPGA中每个通道相对于其中心通道的延迟，ic表示通道编号，Tr表示中频信号从发射到反射回天线雷达的最短时间，τ表示时宽。
S130，基于相位生成余弦值cosval和相位生成正弦值sinval计算权重相位。
步骤S130的计算公式如下：
A6(ic)＝-2/lamda*chandistance1(ic)*cos(-lookbias*pi/180)；
A7(ic)＝-2/lamda*chandistance1(ic)*sin(-lookbias*pi/180)；
wt＝sinval*A6(ic)+cosval*A7(ic)；
其中，lamda表示中频信号的波长，chandistance1(ic)表示FPGA编号为ic的通道与其中心通道的距离，lookbias是天线雷达的安装角，A6、A7分别为第六辅助参数和第七辅助参数。
需要说明的是，得到相位生成余弦值cosval和相位生成正弦值sinval后，将相位生成余弦值cosval和相位生成正弦值sinval时钟同步。由于相位生成正弦值sinval是基于相位生成余弦值cosval计算得到的，这样一般会导致cosval先计算出来，而sinval会延时几个时钟单元生成，为获得正确的权值，必须保持cosval和sinval同步，因此必须将cosval延时相应的时间单元，这样才能产生正确的权值，延时单元的确定往往是通过FIFO单元进行控制，利用同步时钟确定延时单元的多少。
S140，将权重相位输入vivado中的cordic ip核，计算权值。
在本公开实施例中，将权重相位输入vivado中的cordic ip核之前，将权重相位标准化至[-1，1]的范围内，使权重相位满足cordic ip核的输入要求。
具体地，权值的生成需要经过sin和cos计算，在FPGA中实现sin和cos计算有两种方式，一种是查表法，另一种则是cordic法，我们需要进行实时运算，且精度要求非常高，因此我们采取cordic法，vivado中的cordic ip核经过标准化，对输入要求为[-1，1]，因此我们需要对经由辅助数据计算出权重wt判定其正负和奇偶，来确定相位的范围，使得cordic核计算出正确的权值。
步骤S140将权重相位标准化至[-1，1]的范围内包括步骤S141～S143。
S141，将权重相位转换为至少有2位小数的定点数。
S142，判断权值相位的正负和奇偶。
判断权值相位的正负和奇偶包括：判断权重相位的符号位是否为0；若符号位为0，则权重相位为正；若符号位为1，权重相位为负。
判断权重相位的最后一位是否为0；若权重相位的最后一位为0，则权重相位为偶数；若权重相位的最后一位不为0，则权重相位为奇数。
S143，根据权值相位的正负、奇偶状态，将权值相位标准化至[-1，1]的范围内。
具体的，若权重相位为正，且权重相位为偶数，则权重相位的小数部分表示权重相位的标准化值。
若权重相位为正，且权重相位为奇数，则权重相位为奇数，取权重相位的小数部分减1得到的值表示权重相位的标准化值。
若权重相位为负，且权重相位为偶数，则权重相位为偶数，取权重相位的小数部分表示权重相位的标准化值。
若权重相位为负，且权重相位为奇数，则权重相位为奇数，取权重相位的小数部分加1得到的值表示权重相位的标准化值。
图2示意性示出了本公开实施例提供的一种实时权值生成器的流程图，该流程图示意性示出了加权生成器根据权值系数生成权值的过程。
需要说明的是，在上述过程中，参数tmin、delay(ic)还可以是根据辅助数据中的其它数据计算得到的。
tmin的计算步骤如下：
theta＝looklim_bowei*pi/180；
Figure BDA0002864879570000091
RR＝Re*sin(α)/sin(θ)；
Tr＝2*RR/C；
Figure BDA0002864879570000092
其中，looklim_bowei表示波位，theta表示发射天线雷达的可视角，H表示轨道高度，RR表示发射天线到地面目标的距离，Tr表示电磁波来回的时间，τ表示时宽。由于波位包括一个范围，因此Tr也存在范围，选择最小的Tr确定tmin。
delay(ic)的计算公式如下：
delay(ic)＝chandistance*sin(compensatetheta-lookbias*pi/180)/C
-chandistance*Bpt/Kr/lamda
chandistance＝(ic-(N_e+1)/2)*d
其中，lookbias是天线安装角，d表示通道间距，Kr是调频率，lamda是波长，N-e表示总的通道数，chandistance表示通道到中心的距离，Compensatetheta表示中心波位所对应的视角，Bpt计算中心通道处的导数生成的。
chandistance＝(ic-(N_e+1)/2)*d；
compensatetheta＝[looklim_bowei_min+looklim_bowei_max]/2*pi/180；
其中，looklim_bowei_min和looklim_bowei_max分别为波位的最小值和最大值。
Bpt的计算过程包括：
St＝t_center*C/2；
At＝St/2/(Re+Height)
+(2*Re*Height+Height2)/St/2/(Re+Height)；
Bt＝acos(At)-lookbias*pi/180；
Yt＝sin(Bt)；
Spt＝C/2；
Apt＝Spt/2/(Re+Height)
-Spt*(2*Re*Height+Height2)/St/St/2/(Re+Height)；
Bpt＝-1*Apt/sqrt(1-At2)；
其中，t_center表示天线的中心通道到地面目标的来回时间。
t_center＝R_center*2/C。
R_center表示天线的中心通道到地面目标的距离，根据RR的计算公式得到。
图3示意性示出了本公开实施例提供的一种基于FPGA的中频数字波束形成方法的结构示意图。
如图3所示，该基于FPGA的中频数字波束形成方法，将中频信号分为16个波段的信号后分别输入两片FPGA中进行处理，降低了输入FPGA的信号的波段范围，通过两片FPGA共同处理不同波段的信号的方式实现了对该中频信号的数字波束形成，其中，该两片FPGA从外部辅助数据接收到所需要计算权值的参数，分别输入到实时权值生成器中生成权值，与接收到的中频回波信号进行加权处理，然后，将两片FPGA加权处理后的实部和虚部信号进行累加，并进行正交解调和低通滤波处理，完成了基于FPGA的中频数字波束的形成。
图4示意性示出了本公开实施例提供的MATLAB程序生成的实时权值结果示意图，图5示意性示出了本公开实施例提供的一种实时权值产生方法生成的实时权值结果示意图。从图4、图5中的数据可知，本公开提供的方法能准确、稳定地实时计算出通道需要的数字合成权值。
本领域技术人员可以理解，本公开的各个实施例和/或权利要求中记载的特征可以进行多种组合或/或结合，即使这样的组合或结合没有明确记载于本公开中。特别地，在不脱离本公开精神和教导的情况下，本公开的各个实施例和/或权利要求中记载的特征可以进行多种组合和/或结合。所有这些组合和/或结合均落入本公开的范围。
尽管已经参照本公开的特定示例性实施例示出并描述了本公开，但是本领域技术人员应该理解，在不背离所附权利要求及其等同物限定的本公开的精神和范围的情况下，可以对本公开进行形式和细节上的多种改变。因此，本公开的范围不应该限于上述实施例，而是应该不仅由所附权利要求来进行确定，还由所附权利要求的等同物来进行限定。

Claims (10)

1.一种实时权值产生方法，应用于基于FPGA的中频数字波束形成过程，其特征在于，包括：
获取外部辅助数据；
基于所述辅助数据，计算权重相位的相位生成余弦值cosval和相位生成正弦值sinval；
基于所述相位生成余弦值cosval和相位生成正弦值sinval计算所述权重相位；
将所述权重相位输入vivado中的cordic ip核，计算权值。
2.根据权利要求1所述的方法，其特征在于，所述辅助数据为天线雷达发射中频信号的参数信息，至少包括：所述FPGA的中心通道与地面目标的距离、总通道数、所述天线雷达对所述中频信号的采样频率、光速、所述天线雷达的轨道高度、地球半径、所述FPGA中每个通道相对于其中心通道的延迟、所述中频信号从发射到反射回所述天线雷达的最短时间、时宽、所述中频信号的波长、所述FPGA编号为ic的通道与其中心通道的距离、所述天线雷达的安装角。
3.根据权利要求1所述的方法，其特征在于，所述基于所述辅助数据，计算权重相位的相位生成余弦值cosval和相位生成正弦值sinval包括：
SRval＝(0：(Nrg-1))*A1+A2(ic)；
cosval1＝SRval/A4；
cosval2＝A5./SRval；
cosval＝cosval1+cosval2；
sinval1＝1-cosval；
sinval2＝1+cosval；
sinval＝sqrt(1-cosvaj2)；
其中，A1、A2、A3、A4、A5分别表示第一辅助参数、第二辅助参数、第三辅助参数、第四辅助参数和第五辅助参数，根据所述辅助数据计算得到，SRval表示所述FPGA的中心通道与地面目标的距离，Nrg表示总的通道数，cosval1、cosval2表示所述相位生成余弦值cosval的分量，sinval1、sinval2表示所述相位生成正弦值sinval的分量。
4.根据权利要求3所述的方法，其特征在于，所述第一辅助参数、第二辅助参数、第三辅助参数、第四辅助参数和第五辅助参数的计算公式如下：
A1＝1/Fs*C/2；
A2(ic)＝(tmin-delay(ic))*C/2；
A3＝(Height+Re)2-Re2；
A4＝2*(Re+Height)；
A5＝A3/A4；
Figure FDA0002864879560000021
其中，Fs为所述天线雷达对所述中频信号的采样频率，C为光速，Height为所述天线雷达的轨道高度，Re为地球半径，delay(ic)表示所述FPGA中每个通道相对于其中心通道的延迟，ic表示通道编号，Tr表示所述中频信号从发射到反射回所述天线雷达的最短时间，τ表示时宽。
5.根据权利要求3所述的方法，其特征在于，得到所述相位生成余弦值cosval和相位生成正弦值sinval后，将所述相位生成余弦值cosval和相位生成正弦值sinval时钟同步。
6.根据权利要求1所述的方法，其特征在于，所述基于所述相位生成余弦值cosval和相位生成正弦值sinval计算所述权重相位包括：
A6(ic)＝-2/lamda*chandistancel(ic)*cos(-lookbias*pi/180)；
A7(ic)＝-2/lamda*chandistancel(ic)*sin(-lookbias*pi/180)；
wt＝sinval*A6(ic)+cosval*A7(ic)；
其中，lamda表示所述中频信号的波长，chandistancel(ic)表示所述FPGA编号为ic的通道与其中心通道的距离，lookbias是所述天线雷达的安装角，A6、A7分别为第六辅助参数和第七辅助参数。
7.根据权利要求1所述的方法，其特征在于，将所述权重相位输入vivado中的cordicip核之前，将所述权重相位标准化至[-1，1]的范围内，使所述权重相位满足cordicip核的输入要求。
8.根据权利要求7所述的方法，其特征在于，所述将所述权重相位标准化至[-1，1]的范围内包括：
将所述权重相位转换为至少有2位小数的定点数；
判断所述权值相位的正负和奇偶；
根据所述权值相位的正负、奇偶状态，将所述权值相位标准化至[-1，1]的范围内。
9.根据权利要求8所述的方法，其特征在于，所述判断所述权值相位的正负和奇偶包括：
判断所述权重相位的符号位是否为0；
若所述符号位为0，则所述权重相位为正；
若所述符号位为1，所述权重相位为负；
判断所述权重相位的最后一位是否为0；
若所述权重相位的最后一位为0，则所述权重相位为偶数；
若所述权重相位的最后一位不为0，则所述权重相位为奇数。
10.根据权利要求8所述的方法，其特征在于，所述根据所述权值相位的正负、奇偶状态，将所述权值相位标准化至[-1，1]的范围内包括：
若所述权重相位为正，且所述权重相位为偶数，则所述权重相位的小数部分表示所述权重相位的标准化值；
若所述权重相位为正，且所述权重相位为奇数，则所述权重相位为奇数，取所述权重相位的小数部分减1得到的值表示所述权重相位的标准化值；
若所述权重相位为负，且所述权重相位为偶数，则所述权重相位为偶数，取所述权重相位的小数部分表示所述权重相位的标准化值；
若所述权重相位为负，且所述权重相位为奇数，则所述权重相位为奇数，取所述权重相位的小数部分加1得到的值表示所述权重相位的标准化值。
CN202011585142.0A 2020-12-28 2020-12-28 一种实时权值产生方法 Active CN112748401B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011585142.0A CN112748401B (zh) 2020-12-28 2020-12-28 一种实时权值产生方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011585142.0A CN112748401B (zh) 2020-12-28 2020-12-28 一种实时权值产生方法

Publications (2)

Publication Number Publication Date
CN112748401A true CN112748401A (zh) 2021-05-04
CN112748401B CN112748401B (zh) 2024-02-06

Family

ID=75646430

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011585142.0A Active CN112748401B (zh) 2020-12-28 2020-12-28 一种实时权值产生方法

Country Status (1)

Country Link
CN (1) CN112748401B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2003264825A1 (en) * 2002-11-04 2004-06-07 Koninklijke Philips Electronics N.V. Superposition of waves using the cordic algorithm
US20100182053A1 (en) * 2009-01-19 2010-07-22 Himax Media Solutions, Inc. Sine/cosine value generation apparatus
US20150293209A1 (en) * 2012-10-08 2015-10-15 Mbda Uk Limited Radar receivers
CN106125050A (zh) * 2016-06-12 2016-11-16 上海无线电设备研究所 一种基于cordic核的正弦空间内波控码计算方法
CN106814990A (zh) * 2017-01-17 2017-06-09 中国科学院上海高等研究院 Cordic算法的优化方法
CN111736124A (zh) * 2020-07-31 2020-10-02 中国科学院空天信息创新研究院 雷达信号通道误差处理方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2003264825A1 (en) * 2002-11-04 2004-06-07 Koninklijke Philips Electronics N.V. Superposition of waves using the cordic algorithm
US20100182053A1 (en) * 2009-01-19 2010-07-22 Himax Media Solutions, Inc. Sine/cosine value generation apparatus
US20150293209A1 (en) * 2012-10-08 2015-10-15 Mbda Uk Limited Radar receivers
CN106125050A (zh) * 2016-06-12 2016-11-16 上海无线电设备研究所 一种基于cordic核的正弦空间内波控码计算方法
CN106814990A (zh) * 2017-01-17 2017-06-09 中国科学院上海高等研究院 Cordic算法的优化方法
CN111736124A (zh) * 2020-07-31 2020-10-02 中国科学院空天信息创新研究院 雷达信号通道误差处理方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHENG-HAN LI,: "A 4× 64 MIMO Detector for Generalized Spatial Modulation Systems", IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS–I: REGULAR PAPERS, pages 1 - 13 *
刘志强;冯晓伟;: "基于高速DA转换的雷达中频信号产生方法", 火力与指挥控制, no. 12, pages 2146 - 2149 *
李耀南;罗丁利;: "基于FPGA的DBF导向矢量实时算法", 火控雷达技术, no. 03, pages 27 - 31 *

Also Published As

Publication number Publication date
CN112748401B (zh) 2024-02-06

Similar Documents

Publication Publication Date Title
CN109459733B (zh) 基于调相方式的防撞雷达目标速度模拟装置、系统及方法
CN107479038B (zh) 一种高精度雷达目标回波实时模拟方法
CN110488285B (zh) 基于编码信号的非中断的双基sar相位同步信号处理方法
CN112379341B (zh) 一种雷达接收机的数字化建模方法
CN104502912A (zh) 高速运动目标逆合成孔径雷达成像方法
JP2012032348A (ja) Ｓａｒ装置
CN110596701B (zh) 基于二次椭圆模型的非平飞双站sar频域fenlcs成像方法
CN107607952A (zh) 基于电磁涡旋波的三维合成孔径雷达成像方法
Shaokun et al. Deception jamming method for ISAR based on sub-Nyquist sampling technology
CN114185046B (zh) 基于发射波形参数可调的视频sar成像系统及方法
CN118131236A (zh) 一种多通道联合阈值单比特sar成像方法
CN112748401A (zh) 一种实时权值产生方法
CN102540153A (zh) 一种基于站间直达波干扰的阵列幅相误差校准方法
CN103513239B (zh) 一种采用编码调相信号实现直达波干扰抑制的方法
CN105974374A (zh) 产生Ka波段调频线性连续波的频率合成器及工作方法
CN115308703B (zh) 一种高频mimo雷达辅助校准源及校准方法
CN105372644B (zh) 一种基于动态重修正的自适应波束形成方法及系统
CN114626006B (zh) 雷达成像中cs算法补偿因子实时生成的fpga实现方法
CN108508437A (zh) 基于时间变标的合成孔径雷达成像方法
CN115113199A (zh) 基于数字可编程超表面的雷达通信一体化方法
Munoz-Ferreras et al. Extended envelope correlation for range bin alignment in ISAR
CN109307868B (zh) 一种适用于太赫兹波段的单脉冲成像系统及方法
CN112702098B (zh) 基于fpga的中频数字波束形成方法
US12092720B2 (en) Radar signal processing device, radar device, radar signal processing method and non-transitory computer-readable medium
US11940556B2 (en) Testing device for testing a distance sensor that operates using electromagnetic waves

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