CN101806885A - 一种多通道阵列信号发生方法和装置 - Google Patents

一种多通道阵列信号发生方法和装置 Download PDF

Info

Publication number
CN101806885A
CN101806885A CN 201010131237 CN201010131237A CN101806885A CN 101806885 A CN101806885 A CN 101806885A CN 201010131237 CN201010131237 CN 201010131237 CN 201010131237 A CN201010131237 A CN 201010131237A CN 101806885 A CN101806885 A CN 101806885A
Authority
CN
China
Prior art keywords
signal
frequency
square
shift register
wave
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
CN 201010131237
Other languages
English (en)
Other versions
CN101806885B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2010101312375A priority Critical patent/CN101806885B/zh
Publication of CN101806885A publication Critical patent/CN101806885A/zh
Application granted granted Critical
Publication of CN101806885B publication Critical patent/CN101806885B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种多通道阵列信号发生装置及其发生方法,该装置包括:信号参数设置模块,用于实现信号频率/脉宽/阵元位置/波束形成方向的选择;移位寄存器和反相器网络模块,用于生成频率为f=fclk/2N的多路延时/初始相位等差的方波信号,其中,fclk为时钟信号频率,N为移位寄存器的个数;开关矩阵,用于根据所需的信号参数,从方波信号集合中选择出具有相应延时/初始相位的方波信号;低通滤波器,用于对所述的多路方波信号整形,得到多路正弦波信号;电平转换电路,用于对所述的多路正弦波信号做电平转换。本发明适用于单频脉冲信号波束形成发射/接收信号的模拟,其设计巧妙,具有实现方便、体积小、成本低、功耗低、通道数多的优点。

Description

一种多通道阵列信号发生方法和装置
技术领域
本发明涉及声纳、雷达及无线通信技术领域,特别涉及多通道换能器/天线阵列系统中接收/发射波束形成相关技术。
背景技术
波束形成通常与声波发射/接收换能器阵列或者电磁波发射/接收天线阵列协同工作,广泛应用于声纳、雷达、通信等领域。波束形成包括对接收信号做波束形成和对发射信号做波束形成两部分:
(1)对接收信号做波束形成:为了使换能器/天线阵列加强从某个或多个特定方向到达的信号,基本的波束形成方法对每一路阵元接收的信号按该波束方向进行延时,再把经时延对齐的各路信号相加,得到指向特定方向的波束输出。如图1所示。
(2)对发射信号做波束形成:为了使换能器/天线阵列发射指向一个或多个特定的方向,设计每一路阵元发射的信号按指定波束方向具备一定的延时关系。
用于成像、测绘、定位、目标检测的声纳或者雷达系统,其在被动工作时,通常包含接收-处理部分;在主动工作时,通常包含接收-处理部分和处理-发射部分。对于接收-处理部分,由于先期一般需在室内实验室环境下,对信号接收/处理系统进行调试,而此时接收的换能器/天线阵列信号往往不具备实际使用现场的测试特性,有些调试工作还需消除接收电路与前端换能器/天线阵列的耦合。因而,它需要一种多通道阵列接收信号模拟器来模拟换能器/天线阵列接收特定位置信源的过程,产生多通道的阵列信号,以供所述信号接收-处理系统进行测试、校准、算法验证。对于所述处理-发射部分,则需要一种多通道阵列发射信号生成器以产生一组激励信号,通过功率放大后,加到发射阵列的各个换能器/天线,使得信号能量朝特定方向发射。
为实现指向性波束发射/接收,多通道信号之间必须满足一定的幅度和延时关系。事实上,“多通道阵列接收信号模拟器”生成信号的过程与“多通道阵列发射信号生成器”生成信号的过程源于相同的物理机制,在执行结构上,二者具有技术上的共通性,因此,通常可以把多通道阵列接收信号模拟器和多通道阵列发射信号生成器合成为多通道阵列信号发生器。当接收阵列接收某一个方向信源信号时,它能够仿真产生该接收阵列各阵元的输出接收信号;当发射阵列发射某一个定向波束时,它能够产生该发射阵列各阵元的输入激励信号。
多通道阵列信号发生器在声纳、雷达、通信系统的应用中占有很重要的地位。海底锚系声纳和自主水下航行器/水下滑翔机观测声纳是近年来水声工程技术发展的热点,其中采用多通道阵列发射/接收的主动声纳系统在水底地形成像、目标主动探测、水体性质(典型如流速)测量等方面有着重要的应用。这些场合中由于多采用电池供电,对设备功耗有很高的要求,自主水下航行器/水下滑翔机平台的承载空间也非常有限。
多通道阵列信号发生器是信号发生器的一种。现有的商用信号发生器,如安捷伦、泰克等信号发生器,虽然精度较高、通用性较好,但是通道数少,其通道数普遍为4通道,且各通道独立,通道间信号的关系不便调节,不宜用作多通道阵列信号的发生。
目前,一般实验室自制的多通道信号发生器和为数不多的商用专业多通道信号发生器多采用各通道独立数字波形控制加数模转换技术,如采用单片机、FPGA、CPLD控制多个D-A(数模转换芯片)的方法,来产生多通道的阵列信号。采用该技术,一个通道需要一片D-A。这样,这种技术就存在如下的缺点:单片机、FPGA、CPLD的管脚数资源有限,D-A价格较贵,体积较大,功耗较大;并行D-A所需的数据/控制管脚较多,串行D-A的数据/控制管脚较少而转换速度较慢,且各D-A之间数据更新的同步不好控制,实现过程复杂。这些缺点在需要的阵列信号路数较多时,尤为明显。
综上所述,发展小体积、低功耗、低成本、实现方便、通道数较多的多通道阵列信号发生器有很迫切的需要和广泛的应用前景。
发明内容
为克服上述现有技术的缺陷,本发明针对单频脉冲信号波束形成接收/发射信号的产生,提出一种多通道阵列信号的发生方法和装置,它能在需要信号通道数较多的情况下,以较小的功耗、体积、成本代价方便地模拟声波换能器/电磁波天线接收阵列输出接收信号,也可以生成声波换能器/电磁波天线发射阵列所需的激励信号。
为实现上述目的,本发明提供了一种多通道阵列信号发生装置,它包括:
(1)信号参数设置模块,包含时钟管理器、定时器、信号频率/脉宽/阵元位置/波束形成方向设置子模块、通道选择系数存储子模块,用于实现信号频率/脉宽/阵元位置/波束形成方向的选择;
(2)移位寄存器和反相器网络模块,用于生成频率为f=fclk/2N的多路延时/初始相位等差的方波信号,所述方波信号集合作为备选信号集合,其中,fclk为时钟信号频率,N为移位寄存器的个数;
(3)开关矩阵,用于根据所需的信号参数,从所述的备选信号集合中选择出具有相应延时/初始相位的方波信号,作为各通道输出的多路方波信号;
(4)低通滤波器,用于对所述的多路方波信号整形,得到多路正弦波信号;
(5)电平转换电路,用于对所述的多路正弦波信号做电平转换。
所述信号频率/脉宽/阵列形状阵元位置/接收波束形成方向的选择过程如下:
(1)首先,设置信号频率、阵元位置、波束形成方向,根据所设置的信号频率、阵元位置、波束形成方向计算得到各路方波信号的延时/初始相位,再根据各路方波信号的延时/初始相位计算得到控制开关矩阵所需的系数,存储于通道选择系数存储子模块中;不同的通道选择系数对应不同的阵元位置/波束形成方向;
(2)设置信号脉冲宽度,根据所设置的信号脉宽设计定时器的时间长度,得到高低电平以该时间长度交替的使能信号,去控制移位寄存器和反相器网络的使能端,得到有限脉宽信号;
(3)根据所设置的信号频率,由时钟管理器输出不同频率的时钟信号送入移位寄存器和反相器网络模块,来得到相应频率的信号;或是在时钟管理器输出时钟信号频率不变的情况下,通过改变移位寄存器和反相器网络中的移位寄存器的数目,来得到相应频率的信号。
所述各路信号的延时/初始相位计算过程:一个特定的阵列布置,其在三维空间XYZ坐标系下的第i个阵元的位置的表达式如下:
r i = r xi r yi r zi
接收/发射波束方向的表达式:
u = cos φ cos θ sin φ cos θ sin θ , φ是水平角,θ是俯仰角
假设接收到的信号是一个远场的信号exp(jωt),把它近似为平面波信号,在介质中传播速度为c,角频率ω,频率f。因此第i个阵元接收到的信号的表达式为:
x i ( t ) = exp ( jωt - 2 πf r i T · u / c ) i = 1,2 , · · · M
延时表达式为:
τ i = r i T · u / c = ( cos φ cos θ r xi + sin φ cos θ r yi + sin θ r zi ) / c
相位表达式为:
Figure GSA00000057860700045
由于信号的周期为2π,可以对相位以2π取模,得到:
Figure GSA00000057860700046
所述通道选择系数的计算过程如下:第i个阵元输出信号对应的通道选择系数表示对四舍五入取整,其中,[ ]为向下取整符号,
Figure GSA00000057860700053
为所述各路信号的初始相位,N为移位寄存器的个数,
Figure GSA00000057860700054
为所述移位寄存器和反相器网络生成多路延时/初始相位等差的方波信号的最小相位差。
所述移位寄存器和反相器网络模块是由多个移位寄存器和反相器首尾级联形成的环路,移位寄存器的初始状态均为0。该移位寄存器和反相器网络模块对输入频率为fclk的时钟信号分频得到所需频率的方波信号,所述方波信号频率f=fclk/2N,其中,N为移位寄存器的个数。在移位寄存器和反相器网络中,每个移位寄存器输出方波之间的相位有着等间隔的差异,这个相位的间隔为
Figure GSA00000057860700055
N为移位寄存器的个数,推导过程如下:
各个移位寄存器输出(即移位寄存器状态值)之间的关系可以表示如下式:
y 1 ( t + T ) = y N ( t ) ‾ y i + 1 ( t + T ) = y i ( t ) , 1 ≤ i ≤ N - 1
y 1 ′ ( t + T ) = y N ′ ( t ) ‾ y i + 1 ′ ( t + T ) = y i ′ ( t ) , 1 ≤ i ≤ N - 1
yi′(t)=yi(t)
其中yi′(t)表示把第i个移位寄存器的输出接到反相器,反相器在t时刻的输出值。由以上式子可以得出,yi+1(t)和yi(t)之间的时延是T,换算成相位就是
Figure GSA00000057860700059
若假设y1(t)的初始相位为0,则y2(t)为π/N,y3(t)为2π/N,以此类推,yi(t)为(i-1)π/N,yN(t)为(N-1)π/N。yi+1′(t)和yi′(t)之间的相位差也是
Figure GSA000000578607000510
由于yi′(t)=yi(t),相位相差π,假设y1(t)的初始相位为0,则y1′(t)的初始相位为π,y2′(t)为(N+1)π/N,以此类推,yi′(t)为(N+i-1)π/N,yN′(t)为(2N-1)π/N。由上可知,y1(t),y1′(t),y2(t),y2′(t)…yN(t),yN′(t)的初始相位均匀的分布在0~2π之间。于是,y1(t),y1′(t),y2(t),y2′(t)…yN(t),yN′(t)可以作为具有不同初始相位的方波信号集合,当我们需要某个初始相位的信号时,我们就可以在这个方波信号集合中,选择初始相位最接近所需信号初始相位的方波信号作为输出。当N越大时,相位差
Figure GSA00000057860700061
越小,则信号集合中信号的初始相位就可以越接近我们所需信号的初始相位,误差就会越小。
我们以某个移位寄存器的输出信号为参考,假定其初始相位为0,这样就可得到各个移位寄存器输出信号相对的初始相位,也就得到了2N个延时/初始相位等差的方波信号,作为所述的备选信号集合。
所述低通滤波模块,是用于对所述的多路方波信号整形,得到多路正弦波信号,其过程如下:
周期方波信号的表达式为:
x ( t ) = 1 , | t | < T 1 0 , T 1 < | t | &le; T 2 , 且x(t)=x(t+T)
傅立叶变换为:
X ( j&omega; ) = &Sigma; k = - &infin; + &infin; 2 sin k &omega; 0 T 1 k &delta; ( &omega; - k &omega; 0 )
因为该分频得到的方波的占空比为50%,即 T 1 = T 4 , 代入 &omega; 0 = 2 &pi; T , 所以傅立叶表达式化简为:
X ( j&omega; ) = &Sigma; k = - &infin; + &infin; 2 sin k&pi; / 2 k &delta; ( &omega; - k &omega; 0 )
可以看出,当k为偶数时傅立叶变换系数值为零,只有为奇数时才有非零值。也即是说,该方波只包含了奇数次的谐波,用截止频率为fc(fc满足f<fc<3f,f为多路方波信号的频率)的滤波器,把方波的三次及其以上的谐波滤掉,得到所述的多路正弦波信号。
所述电平转换电路根据下一级电路输入电平标准的要求,将多路正弦波信号电平转换为相应的电平形式。
本发明还提供一种多通道接收阵列输出信号的模拟生成方法,它包括以下步骤:
1)通过信号参数设置模块,对信号频率/脉宽/阵元位置/接收波束方向进行选择;
2)通过移位寄存器/反相器网络模块,生成频率为f=fclk/2N的多路延时/初始相位等差的方波信号,作为备选信号集合,其中,fclk为时钟信号频率,N为移位寄存器的个数;
3)根据所需的信号参数,选择相应的通道选择系数去控制开关矩阵,再由开关矩阵从备选信号集合中选择出具有相应延时/初始相位的多路方波信号,作为各通道输出的多路方波信号;
4)低通滤波器对所述的多路方波信号低通滤波,得到多路的正弦波信号;
5)电平转换电路对所述的多路正弦波信号做与输出接口匹配的电平转换。
所述信号频率/脉宽/阵元位置/接收波束方向的选择过程如下:
(1)首先,设置信号频率、阵元位置、接收波束方向,根据所设置的信号频率、阵元位置、接收波束方向计算得到各路方波信号的延时/初始相位,再根据各路方波信号的延时/初始相位计算得到控制开关矩阵所需的系数,存储于通道选择系数存储子模块中;不同的通道选择系数对应不同的阵元位置/接收波束方向;
(2)设置信号脉冲宽度,根据所设置的信号脉宽设计定时器的时间长度,得到高低电平以该时间长度交替的使能信号,去控制移位寄存器和反相器网络的使能端,得到有限脉宽信号;
(3)根据所设置的信号频率,由时钟管理器输出不同频率的时钟信号送入移位寄存器和反相器网络模块,来得到相应频率的信号;或是在时钟管理器输出时钟信号频率不变的情况下,通过改变移位寄存器和反相器网络中的移位寄存器的数目,来得到相应频率的信号。
本发明还提供一种多通道发射阵列激励信号的生成方法,它包括以下步骤:
1)通过信号参数设置模块,对信号频率/脉宽/阵元位置/发射波束方向进行选择;
2)通过移位寄存器和反相器网络模块,生成频率为f=fclk/2N的多路延时/初始相位等差的方波信号,作为备选信号集合,其中,fclk为时钟信号频率,N为移位寄存器的个数;
3)根据所需的信号参数,选择相应的通道选择系数去控制开关矩阵,再由开关矩阵从备选信号集合中选择出具有相应延时/初始相位的多路方波信号,作为各通道输出的多路方波信号;
4)低通滤波器对所述的多路方波信号低通滤波,得到多路的正弦波信号,作为发射阵列的激励信号。
所述信号频率/脉宽/阵元位置/发射波束方向的选择过程如下:
(1)首先,设置信号频率、阵元位置、发射波束方向,根据所设置的信号频率、阵元位置、发射波束方向计算得到各路方波信号的延时/初始相位,再根据各路方波信号的延时/初始相位计算得到控制开关矩阵所需的系数,存储于通道选择系数存储子模块中;不同的通道选择系数对应不同的阵元位置/发射波束方向;
(2)设置信号脉冲宽度,根据所设置的信号脉宽设计定时器的时间长度,得到高低电平以该时间长度交替的使能信号,去控制移位寄存器和反相器网络的使能端,得到有限脉宽信号;
(3)根据所设置的信号频率,由时钟管理器输出不同频率的时钟信号送入移位寄存器和反相器网络模块,来得到相应频率的信号;或是在时钟管理器输出时钟信号频率不变的情况下,通过改变移位寄存器和反相器网络中的移位寄存器的数目,来得到相应频率的信号。
本发明适用于单频脉冲信号波束形成发射/接收信号的模拟/生成,相比现有技术,本发明更加巧妙,具有实现方便、体积小、成本低、功耗低、通道数较多等显著特点,可在ASIC/FPGA中用硬件描述语言编程后,再加适当的外围电路得以实现。
附图说明
图1是现有技术阵列接收信号的“延时求和”波束形成原理示意图;
图2是本发明实施例的多通道阵列信号发生装置结构示意图;
图3是本发明实施例的多通道接收阵列输出信号的模拟生成方法流程图;
图4是本发明实施例的阵元位置和波束入射方向的三维空间坐标系表示图;
图5是本发明实施例的U型阵阵元位置示意图;
图6是本发明实施例的移位寄存器/反相器网络模块连接方法示意图;
图7是本发明实施例的128路所需延时/相位的方波信号产生过程示意图;
图8是本发明实施例的低通滤波器电路原理图;
图9是本发明实施例的电平转换电路原理图。
具体实施方式
下面,结合附图和具体实施例详细说明本发明。
参见附图。本实施例的多通道接收阵列输出信号的模拟发生方法,是针对声纳系统的应用环境描述的,但其原理在雷达、通信系统中亦是相通的。本发明的具体实施实例只是起示范的作用,并不用以限制本发明的保护范围。
本实施例可产生128通道差分输出、各通道间相位一致性较好、频率165/195khz、阵元位置由U形/均匀直线(ULA)阵决定、接收波束方向为正负0/30/45/60/90度、脉宽10ms的声纳接收阵列输出信号。多通道阵列信号发生装置的模块图如图2所示,主要由外围电路模块、信号参数设置模块、移位寄存器和反相器网络模块、开关矩阵、低通滤波器、电平转换电路组成。其中外围电路模块用于为整个电路提供系统时钟、电源、通信接口、复位控制信号、拨码开关等,该外围电路具体包括:20Mhz的晶振作为整个电路的系统时钟,正负5伏电源供电,RS232/JTAG作为通信接口,触发开关提供复位控制信号,多个拨码开关用于选择不同的信号频率和接收波束方向。在此装置的基础上提出了多通道接收阵列信号的模拟发生方法。如图3所示,实施该方法的具体流程主要包括:信号参数设置、备选信号集合的获取、满足所需延时/初始相位的方波信号的选取、低通滤波、单端-差分电平转换、外围电路设计。
所需要模拟的声纳接收阵列信号为单频正弦波脉冲信号,其具体参数为:频率165/195khz、阵形为U形/均匀直线阵、接收波束方向为正负0/30/45/60/90度、脉宽10ms。
多通道接收阵列信号的模拟发生方法,包括如下步骤:
1、通过信号参数设置模块,对信号频率/脉宽/阵元位置/接收波束方向进行选择,具体包括如下步骤:
(1)首先,设置信号频率、阵元位置、接收波束方向,根据所设置的信号频率、阵元位置、接收波束方向计算得到各路方波信号的延时/初始相位,所述各路信号的延时/初始相位计算过程如下:
图4中,ri表示阵元位置和u表示波束入射方向,表1是在声源频率165khz、U形阵(该阵有128个阵元,各个阵元的坐标位置如图5所示)、入射角为0度时,按照下列表达式:
延时表达式
&tau; i = r i T &CenterDot; u / c = ( cos &phi; cos &theta; r xi + sin &phi; cos &theta; r yi + sin &theta; r zi ) / c
相位表达式
Figure GSA00000057860700102
以2π取模后相位表达式
Figure GSA00000057860700103
算出的各个阵元输出信号对应的初始相位。计算中声速假设为1500m/s。
表1各阵元输出信号对应的初始相位
  阵元号   1   2   3   4   5   6   7   8   9   10   11   12   13   14   15
  相位   327   47   126   206   286   6   86   166   246   327   47   127   208   289   10
  阵元号   16   17   18   19   20   21   22   23   24   25   26   27   28   29   30
  相位   91   172   253   334   56   138   219   302   24   106   189   272   355   78   161
  阵元号   31   32   33   34   35   36   37   38   39   40   41   42   43   44   45
  相位   245   329   54   138   223   309   35   121   208   295   23   111   200   290   21
  阵元号   46   47   48   49   50   51   52   53   54   55   56   57   58   59   60
  相位   112   205   299   34   130   229   329   71   177   286   39   157   282   57   204
  阵元号   61   62   63   64   65   66   67   68   69   70   71   72   73   74   75
  相位   9   200   71   0   0   71   200   9   204   57   282   157   39   286   177
  阵元号   76   77   78   79   80   81   82   83   84   85   86   87   88   89   90
  相位   71   329   229   130   34   299   205   112   21   290   200   111   23   295   208
  阵元号   91   92   93   94   95   96   97   98   99   100   101   102   103   104   105
  相位   121   35   309   223   138   54   329   245   161   78   355   272   189   106   24
  阵元号   106   107   108   109   110   111   112   113   114   115   116   117   118   119   120
  相位   302   219   138   56   334   253   172   91   10   289   208   127   47   327   246
  阵元号   1   2   3   4   5   6   7   8   9   10   11   12   13   14   15
  阵元号   121   122   123   124   125   126   127   128
  相位   166   86   6   286   206   126   47   327
对其他接收波束方向或者其他类型接收阵,各阵元输出信号初始相位的计算方法同上。
再根据各路方波信号的延时/初始相位计算得到控制开关矩阵所需的系数,所述通道选择系数的计算过程如下:
第i个阵元输出信号对应的通道选择系数
Figure GSA00000057860700111
表示对
Figure GSA00000057860700112
四舍五入取整,[ ]为向下取整符号,
Figure GSA00000057860700113
为所述移位寄存器和反相器网络生成多路延时/初始相位等差的方波信号的最小相位差。与表1相对应的各阵元输出信号对应的通道选择系数如表2所示:
表2各阵元输出信号对应的通道选择系数
  阵元号   1   2   3   4   5   6   7   8   9   10   11   12   13   14   15
  系数值   220   32   85   139   193   4   58   112   166   220   32   86   140   194   7
  阵元号   16   17   18   19   20   21   22   23   24   25   26   27   28   29   30
  系数值   61   116   170   225   38   93   148   203   16   72   127   183   239   53   109
  阵元号   31   32   33   34   35   36   37   38   39   40   41   42   43   44   45
  系数值   165   222   36   93   151   208   24   82   140   199   16   75   135   195   14
  阵元号   46   47   48   49   50   51   52   53   54   55   56   57   58   59   60
  系数值   76   138   201   23   88   154   221   48   119   192   26   106   190   39   137
  阵元号   61   62   63   64   65   66   67   68   69   70   71   72   73   74   75
  系数值   6   135   48   0   0   48   135   6   137   39   190   106   26   192   119
  阵元号   76   77   78   79   80   81   82   83   84   85   86   87   88   89   90
  系数值   48   221   154   88   23   201   138   76   14   195   135   75   16   199   140
  阵元号   91   92   93   94   95   96   97   98   99   100   101   102   103   104   105
  阵元号   1   2   3   4   5   6   7   8   9   10   11   12   13   14   15
  系数值   82   24   208   151   93   36   222   165   109   53   239   183   127   72   16
  阵元号   106   107   108   109   110   111   112   113   114   115   116   117   118   119   120
  系数值   203   148   93   38   225   170   116   61   7   194   140   86   32   220   166
  阵元号   121   122   123   124   125   126   127   128
  系数值   112   58   4   193   139   85   32   220
对其他接收波束方向或者其他类型接收阵,各阵元输出信号的通道选择系数计算方法同上。
将计算得到各组通道选择系数存储于通道选择系数存储子模块中,不同组的通道选择系数对应不同的阵元位置/接收波束方向;
(2)设置信号脉宽10ms,根据所设置的信号脉宽设计定时器的时间长度为10ms,得到高低电平以10ms交替的使能信号,去控制移位寄存器和反相器网络的使能端,得到脉宽为10ms的信号;
(3)系统输入时钟为20Mhz,根据所设置的信号165/195khz频率,由时钟管理器输出特定频率的时钟信号送入移位寄存器和反相器网络模块,来得到165/195khz频率的信号;或是在时钟管理器输出时钟信号频率不变的情况下,通过改变移位寄存器和反相器网络中的移位寄存器的数目N,来得到165/195khz频率的信号。
2、通过移位寄存器和反相器网络模块,生成频率为f=fclk/2N的多路延时/初始相位等差的方波信号,作为备选信号集合,其中,fclk为时钟信号频率,N为移位寄存器的个数。对所述的f=fclk/2N表达式,可以验证如下(如表3所示):
表3:各时钟周期时刻移位寄存器的状态真值表
  y1(t)   y2(t)   y3(t)   ……   yN(t)
  0   0   0   0   ……   0
  T   1   0   0   ……   0
  2T   1   1   0   ……   0
  ……   ……   ……   ……   ……   ……
  (N-1)T   1   1   1   ……   0
  NT   1   1   1   ……   1
  y1(t)   y2(t)   y3(t)   ……   yN(t)
  (N+1)T   0   1   1   ……   1
  (N+2)T   0   0   1   ……   1
  ……   ……   ……   ……   ……   ……
  (2N-2)T   0   0   0   ……   1
  (2N-1)T   0   0   0   ……   1
  2NT   0   0   0   ……   0
  (2N+1)T   1   0   0   ……   0
  (2N+2)T   1   1   0   ……   0
  ……   ……   ……   ……   ……   ……
表3中,T为时钟周期,yi(t)为t时刻时第i个移位寄存器的状态值。由表3可知,在周期为T的时钟边沿触发作用下,每个移位寄存器交替地输出NT时间长的高电平,NT时间长的低电平,该高低电平的交替就得到了方波。可以算出移位寄存器输出方波的频率为f=1/(2NT)=fclk/(2N)。
所述的多路延时/初始相位等差的方波信号具体的生成过程如下:
此处选用D触发器作为移位寄存器,按照图6所示的方式首尾级联起来成为环路,得到D触发器和反相器网络模块,并且触发器的初始状态都设为0,则可对时钟信号分频得到所需频率的方波信号,信号频率f=fclk/2N,N为触发器的个数。此处20Mhz的系统时钟经过时钟管理器2倍频后fclk=40Mhz,f=165khz,可以计算出N=(fclk/f)/2≈121。实际产生的信号频率为f′=fclk/(2N)=40Mhz/(2*121)≈165.3khz,有0.3khz的误差。对于165khz信号,可实现的最小相位差
Figure GSA00000057860700141
换算成角度是180/121=1.488度。也即是说相位的精度为1.488度。可以看出,提高时钟的频率,相应的增加D触发器的个数,可提高最小延时/相位的精度。
我们以某个D触发器的输出信号为参考,假定其初始相位为0,这样就可得到各个移位寄存器输出信号相对的初始相位,也就得到了242个相位/延时等差的方波信号,作为开关矩阵需要的备选信号集合。此处D触发器/反相器网络用VHDL硬件描述语言在FPGA中实现。
3、根据所需的信号参数,选择相应的通道选择系数去控制开关矩阵,再由开关矩阵从备选信号集合中选择出具有相应延时/初始相位的128路方波信号,作为各通道输出的128路方波信号,所述的从备选信号集合中选择出满足延时需要的128路方波信号的具体过程如下:
如图7所示,121个移位寄存器对应165khz信号,102对应195khz。以得到128路165khz方波信号为例,当接收波束方向为0度时,与0度对应的一组通道选择系数用于控制开关矩阵,在备选信号集合242个信号中,选择最接近于所需延时/初始相位的方波输出作为输出信号,这样就得到了与0度对应的128路方波信号。可以通过外围的拨码开关去选择各组通道选择系数来控制开关矩阵,也可以通过通信接口(如JTAG/RS232口)直接用指令去选择各组通道选择系数进而去控制开关矩阵。此处开关矩阵用VHDL硬件描述语言在FPGA中实现。
以拨码开关为例,拨码开关的不同组合代表了各个接收波束方向(也对应于各组控制系数),由拨码开关去控制开关矩阵进行输出信号的选择。
4、低通滤波器对所述的多路方波信号低通滤波,得到128路的正弦波信号,低通滤波器截止频率和芯片的选择过程如下:
占空比为50%的周期方波信号,其傅立叶表达式:
X ( j&omega; ) = &Sigma; k = - &infin; + &infin; 2 sin k&pi; / 2 k &delta; ( &omega; - k &omega; 0 )
可以看出,当k为偶数时傅立叶变换系数值为零,只有为奇数时才有非零值。方波信号频率若为f0=165khz,要将该方波信号变为正弦波,需加低通滤波器进行滤波,滤掉三次及其以上的谐波。于是,我们的低通滤波器的截止频率fc应该满足f0<fc<3f0,也即是165khz<fc<495khz。同时考虑到信号为195khz时滤波的需要,195khz<fc<585khz,我们选择截止频率为fc=200khz。
Max275是低通滤波芯片,其中心频率可达300khz,截止频率是通过外加电阻网络设定的,该滤波芯片的特点之一为可以通过改变信号输入处的电阻对增益进行一定的调节。各个电阻和芯片Max275的连接关系如图8所示,各个电阻值可通过运行Max275 Filter Design Software计算得到。需要注意的是,外接电阻最小不宜小于5kΩ,最大不宜超过4MΩ,因为电阻过小可能达不到运放的驱动能力,过大会产生寄生电容。最后在输出加1uF的电容去掉偏置,就可得到所需的双极性正弦波。
5、根据下一级电路差分输入的实际需要,通过电平转换电路,将多路正弦波信号电平转换为差分信号。此处采用AD8032为主芯片,加外围电路构成如图9所示的单端转差分电路。

Claims (10)

1.一种多通道阵列信号发生装置,其特征在于包括:
(1)信号参数设置模块,包含时钟管理器、定时器、信号频率/脉宽/阵元位置/波束形成方向设置子模块、通道选择系数存储子模块,用于实现信号频率/脉宽/阵元位置/波束形成方向的选择;
(2)移位寄存器和反相器网络模块,用于生成频率为f=fclk/2N的多路延时/初始相位等差的方波信号,所述方波信号集合作为备选信号集合,其中,fclk为时钟信号频率,N为移位寄存器的个数;
(3)开关矩阵,用于根据所需的信号参数,从所述的备选信号集合中选择出具有相应延时/初始相位的方波信号,作为各通道输出的多路方波信号;
(4)低通滤波器,用于对所述的多路方波信号整形,得到多路正弦波信号;
(5)电平转换电路,用于对所述的多路正弦波信号做电平转换。
2.根据权利要求1所述的多通道阵列信号发生装置,其特征在于所述信号频率/脉宽/阵元位置/波束形成方向的选择过程如下:
(1)设置信号频率、阵元位置、波束形成方向,根据所设置的信号频率、阵元位置、波束形成方向计算得到各路方波信号的延时/初始相位,再根据各路方波信号的延时/初始相位计算得到控制开关矩阵所需的系数,存储于通道选择系数存储子模块中;不同的通道选择系数对应不同的阵元位置/波束形成方向;
(2)设置信号脉冲宽度,根据所设置的信号脉宽设计定时器的时间长度,得到高低电平以该时间长度交替的使能信号,去控制移位寄存器和反相器网络的使能端,得到有限脉宽信号;
(3)根据所设置的信号频率,由时钟管理器输出不同频率的时钟信号送入移位寄存器和反相器网络模块,来得到相应频率的信号,或是在时钟管理器输出时钟信号频率不变的情况下,通过改变移位寄存器和反相器网络中的移位寄存器的数目,来得到相应频率的信号。
3.根据权利要求1或2所述的多通道阵列信号发生装置,其特征在于所述的移位寄存器和反相器网络模块包含多个移位寄存器和反相器,该多个移位寄存器和反相器首尾级联成环路,每个移位寄存器的初始状态均设为0。
4.根据权利要求3所述的多通道阵列信号发生装置,其特征在于:所述多路正弦波信号是通过选择截止频率为fc的低通滤波器,将占空比为50%的周期方波的三次及三次以上谐波滤掉而得到,其中,fc满足f<fc<3f,f为所述多路方波信号的频率。
5.一种多通道接收阵列输出信号的模拟生成方法,其特征在于包括以下步骤:
(1)通过信号参数设置模块,对信号频率/脉宽/阵元位置/接收波束方向进行选择;
(2)通过移位寄存器和反相器网络模块,生成频率为f=fclk/2N的多路延时/初始相位等差的方波信号,作为备选信号集合,其中,fclk为时钟信号频率,N为移位寄存器的个数;
(3)根据所需的信号参数,选择相应的通道选择系数去控制开关矩阵,再由开关矩阵从备选信号集合中选择出具有相应延时/初始相位的多路方波信号,作为各通道输出的多路方波信号;
(4)低通滤波器对所述的多路方波信号低通滤波,得到多路的正弦波信号;
(5)电平转换电路对所述的多路正弦波信号做与输出接口匹配的电平转换。
6.如权利要求5所述的多通道接收阵列输出信号的模拟生成方法,其特征在于所述信号频率/脉宽/阵元位置/接收波束方向的选择过程如下:
(1)设置信号频率、阵元位置、接收波束方向,根据所设置的信号频率、阵元位置、接收波束方向计算得到各路方波信号的延时/初始相位,再根据各路方波信号的延时/初始相位计算得到控制开关矩阵所需的系数,存储于通道选择系数存储子模块中;不同的通道选择系数对应不同的阵元位置/接收波束方向;
(2)设置信号脉冲宽度,根据所设置的信号脉宽设计定时器的时间长度,得到高低电平以该时间长度交替的使能信号,去控制移位寄存器和反相器网络的使能端,得到有限脉宽信号;
(3)根据所设置的信号频率,由时钟管理器输出不同频率的时钟信号送入移位寄存器和反相器网络模块,来得到相应频率的信号;或是在时钟管理器输出时钟信号频率不变的情况下,通过改变移位寄存器和反相器网络中的移位寄存器的数目,来得到相应频率的信号。
7.如权利要求5或6所述的多通道接收阵列输出信号的模拟生成方法,其特征在于:所述多路正弦波信号是通过选择截止频率为fc的低通滤波器,将占空比为50%的周期方波的三次及三次以上谐波滤掉而得到,其中,fc满足f<fc<3f,f为所述多路方波信号的频率。
8.一种多通道发射阵列激励信号的生成方法,其特征在于包括以下步骤:
(1)通过信号参数设置模块,对信号频率/脉宽/阵元位置/发射波束方向进行选择;
(2)通过移位寄存器和反相器网络模块,生成频率为f=fclk/2N的多路延时/初始相位等差的方波信号,作为备选信号集合,其中,fclk为时钟信号频率,N为移位寄存器的个数;
(3)根据所需的信号参数,选择相应的通道选择系数去控制开关矩阵,再由开关矩阵从备选信号集合中选择出具有相应延时/初始相位的多路方波信号,作为各通道输出的多路方波信号;
(4)低通滤波器对所述的多路方波信号低通滤波,得到多路的正弦波信号,作为发射阵列的激励信号。
9.如权利要求8所述的多通道发射阵列激励信号的生成方法,其特征在于:所述信号频率/脉宽/阵元位置/发射波束方向的选择过程如下:
(1)设置信号频率、阵元位置、发射波束方向,根据所设置的信号频率、阵元位置、发射波束方向计算得到各路方波信号的延时/初始相位,再根据各路方波信号的延时/初始相位计算得到控制开关矩阵所需的系数,存储于通道选择系数存储子模块中;不同的通道选择系数对应不同的阵元位置/发射波束方向;
(2)设置信号脉冲宽度,根据所设置的信号脉宽设计定时器的时间长度,得到高低电平以该时间长度交替的使能信号,去控制移位寄存器和反相器网络的使能端,得到有限脉宽信号;
(3)根据所设置的信号频率,由时钟管理器输出不同频率的时钟信号送入移位寄存器和反相器网络模块,来得到相应频率的信号;或是在时钟管理器输出时钟信号频率不变的情况下,通过改变移位寄存器和反相器网络中的移位寄存器的数目,来得到相应频率的信号。
10.如权利要求8或9所述的多通道发射阵列激励信号的生成方法,其特征在于:所述多路正弦波信号是通过选择截止频率为fc的低通滤波器,将占空比为50%的周期方波的三次及三次以上谐波滤掉而得到,其中,fc满足f<fc<3f,f为所述多路方波信号的频率。
CN2010101312375A 2010-03-24 2010-03-24 一种多通道阵列信号发生方法和装置 Expired - Fee Related CN101806885B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101312375A CN101806885B (zh) 2010-03-24 2010-03-24 一种多通道阵列信号发生方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101312375A CN101806885B (zh) 2010-03-24 2010-03-24 一种多通道阵列信号发生方法和装置

Publications (2)

Publication Number Publication Date
CN101806885A true CN101806885A (zh) 2010-08-18
CN101806885B CN101806885B (zh) 2012-07-04

Family

ID=42608763

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101312375A Expired - Fee Related CN101806885B (zh) 2010-03-24 2010-03-24 一种多通道阵列信号发生方法和装置

Country Status (1)

Country Link
CN (1) CN101806885B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103117732A (zh) * 2013-02-22 2013-05-22 哈尔滨工程大学 多路视频脉冲信号发生装置及方法
CN103368529A (zh) * 2013-07-29 2013-10-23 中国石油大学(华东) 一种多通道声波相控任意波形信号发生器
CN103558790A (zh) * 2013-11-05 2014-02-05 济南诺顿科技有限公司 一种基准正弦信号产生电路
CN105699961A (zh) * 2016-04-01 2016-06-22 北方通用电子集团有限公司 超宽带时控阵探测装置
CN107219524A (zh) * 2017-05-11 2017-09-29 西北工业大学 一种基于全局最小相位近似的sar成像优化方法
CN108226909A (zh) * 2017-12-29 2018-06-29 苏州联视泰电子信息技术有限公司 一种基于水下多通道信号采集传输阵列系统的模拟装置
CN109298080A (zh) * 2018-10-31 2019-02-01 江苏大学 基于特征导波的焊缝缺陷检测专用阵列式传感器的分时激励系统及方法
CN111383655A (zh) * 2018-12-29 2020-07-07 北京嘉楠捷思信息技术有限公司 一种波束形成方法、装置及计算机可读存储介质
CN112187278A (zh) * 2020-10-15 2021-01-05 西北工业大学 一种声纳信号模拟器延时控制的声纳信号模拟方法及装置
CN114759987A (zh) * 2022-06-14 2022-07-15 天津大学 同频方波作为载波的调制解调方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5374903A (en) * 1988-04-22 1994-12-20 Hughes Aircraft Company Generation of wideband linear frequency modulation signals
JPH10300844A (ja) * 1997-04-25 1998-11-13 Tokimec Inc パルスレーダ距離測定装置及びパルスレーダ距離測定方 法
JP2006053161A (ja) * 2005-10-17 2006-02-23 Toshiba Corp パルスデータ生成方法、シェイプデータ生成方法、パルスデータ生成装置、シェイプデータ生成装置及び送信パルス信号生成装置
CN101082667A (zh) * 2006-06-01 2007-12-05 北京航空航天大学 一种毫米波捷变频雷达目标模拟器
CN201266240Y (zh) * 2008-10-14 2009-07-01 武汉大学 一种高频段的多通道dds信号发生器
CN101470199A (zh) * 2007-12-28 2009-07-01 中国科学院声学研究所 一种基于干涉型光纤水听器的铺设式共形阵声纳装置
CN101504457A (zh) * 2007-12-28 2009-08-12 中国科学院声学研究所 一种基于光纤光栅型光纤水听器的铺设式共形阵声纳装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5374903A (en) * 1988-04-22 1994-12-20 Hughes Aircraft Company Generation of wideband linear frequency modulation signals
JPH10300844A (ja) * 1997-04-25 1998-11-13 Tokimec Inc パルスレーダ距離測定装置及びパルスレーダ距離測定方 法
JP2006053161A (ja) * 2005-10-17 2006-02-23 Toshiba Corp パルスデータ生成方法、シェイプデータ生成方法、パルスデータ生成装置、シェイプデータ生成装置及び送信パルス信号生成装置
CN101082667A (zh) * 2006-06-01 2007-12-05 北京航空航天大学 一种毫米波捷变频雷达目标模拟器
CN101470199A (zh) * 2007-12-28 2009-07-01 中国科学院声学研究所 一种基于干涉型光纤水听器的铺设式共形阵声纳装置
CN101504457A (zh) * 2007-12-28 2009-08-12 中国科学院声学研究所 一种基于光纤光栅型光纤水听器的铺设式共形阵声纳装置
CN201266240Y (zh) * 2008-10-14 2009-07-01 武汉大学 一种高频段的多通道dds信号发生器

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《武汉大学学报(理学版)》 20091231 陈泽宗; 徐鹏; 阮航; 汪志水 多通道多频高频地波雷达DDS数字信号源的设计 第55卷, 第6期 2 *
《电声技术》 20070430 王留全; 黄建国; 张群飞; 雷开卓 水下多目标模拟信号源的设计及实现 第31卷, 第4期 2 *
《电声技术》 20100131 翟江鹏; 杨益新; 张宇 基于NIPCI-6722的水下多传感器阵列信号硬件仿真 第34卷, 第1期 2 *
《计算机测量与控制》 20080430 赵慎; 乔纯捷; 王跃科; 马熠 基于多DSP并行处理的鱼雷仿真信号源 第16卷, 第4期 2 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103117732A (zh) * 2013-02-22 2013-05-22 哈尔滨工程大学 多路视频脉冲信号发生装置及方法
CN103117732B (zh) * 2013-02-22 2015-12-09 哈尔滨工程大学 多路视频脉冲信号发生装置及方法
CN103368529A (zh) * 2013-07-29 2013-10-23 中国石油大学(华东) 一种多通道声波相控任意波形信号发生器
CN103368529B (zh) * 2013-07-29 2014-12-24 中国石油大学(华东) 一种多通道声波相控任意波形信号发生器
CN103558790A (zh) * 2013-11-05 2014-02-05 济南诺顿科技有限公司 一种基准正弦信号产生电路
CN105699961A (zh) * 2016-04-01 2016-06-22 北方通用电子集团有限公司 超宽带时控阵探测装置
CN107219524A (zh) * 2017-05-11 2017-09-29 西北工业大学 一种基于全局最小相位近似的sar成像优化方法
CN108226909A (zh) * 2017-12-29 2018-06-29 苏州联视泰电子信息技术有限公司 一种基于水下多通道信号采集传输阵列系统的模拟装置
CN108226909B (zh) * 2017-12-29 2021-09-14 苏州联视泰电子信息技术有限公司 一种基于水下多通道信号采集传输阵列系统的模拟装置
CN109298080A (zh) * 2018-10-31 2019-02-01 江苏大学 基于特征导波的焊缝缺陷检测专用阵列式传感器的分时激励系统及方法
CN109298080B (zh) * 2018-10-31 2023-12-15 江苏大学 基于特征导波的焊缝缺陷检测专用阵列式传感器的分时激励系统及方法
CN111383655A (zh) * 2018-12-29 2020-07-07 北京嘉楠捷思信息技术有限公司 一种波束形成方法、装置及计算机可读存储介质
CN111383655B (zh) * 2018-12-29 2023-08-04 嘉楠明芯(北京)科技有限公司 一种波束形成方法、装置及计算机可读存储介质
CN112187278A (zh) * 2020-10-15 2021-01-05 西北工业大学 一种声纳信号模拟器延时控制的声纳信号模拟方法及装置
CN114759987A (zh) * 2022-06-14 2022-07-15 天津大学 同频方波作为载波的调制解调方法

Also Published As

Publication number Publication date
CN101806885B (zh) 2012-07-04

Similar Documents

Publication Publication Date Title
CN101806885B (zh) 一种多通道阵列信号发生方法和装置
CN107332629B (zh) 一种基于声学轨道角动量多路复用的信号传输方法
CN100583696C (zh) 用于仿真无线信道的方法和设备
CN106772224A (zh) 一种采用时频分析的l型阵列二维波达方向估计算法
CN102866401A (zh) 一种基于mimo技术的三维成像方法
CN105049142A (zh) 一种双通道静态基带信道模拟装置及方法
CN109039506A (zh) 一种水下移动信道仿真方法
CN105204004A (zh) 基于数字延时和相位补偿的发射数字波束形成方法
CN103245936A (zh) 一种基于drfm的通用无线电高度模拟器
CN103592647A (zh) 阵列三维sar数据获取方法
Mai et al. Beampattern optimization for frequency diverse array with sparse frequency waveforms
CN105137401A (zh) 一种快速细步进捷变频雷达信号产生装置
CN105300437B (zh) 一种vlbi基带信号小数时延仿真方法
CN104535987A (zh) 适用于均匀圆阵列声纳系统的幅相误差自校正方法
CN105187232A (zh) 一种动态信道模拟装置及方法
CN104038959A (zh) 一种信道模拟器及其建模方法
US10396909B1 (en) Multi-antenna test system and method
CN101912277B (zh) 一种基于流水线设计的实时数字正交解调的方法与装置
CN103323831B (zh) 一种基于czt变换和剪除分裂基快速傅里叶变换的三维摄像声纳波束形成方法
CN102803986A (zh) 雷达装置
Kebkal et al. Development of upper-layer protocols with S2CR acoustic modems emulator
CN109375197B (zh) 一种小尺寸矢量阵低频散射校正方法
CN104077479A (zh) 一种基于守恒迎风格式获取参量阵声场空间分布的方法
CN104655111A (zh) 一种激光探测器模拟方法及装置
CN111177988A (zh) 基于fpga实现的智能多波束形成系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120704

Termination date: 20150324

EXPY Termination of patent right or utility model