CN1003746B - 最大熵可编程雷达回波模拟器 - Google Patents
最大熵可编程雷达回波模拟器 Download PDFInfo
- Publication number
- CN1003746B CN1003746B CN85105633.4A CN85105633A CN1003746B CN 1003746 B CN1003746 B CN 1003746B CN 85105633 A CN85105633 A CN 85105633A CN 1003746 B CN1003746 B CN 1003746B
- Authority
- CN
- China
- Prior art keywords
- clutter
- circuit
- simulation
- output
- noise
- 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.)
- Expired
Links
- 238000004088 simulation Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims abstract description 14
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 230000008569 process Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 abstract description 2
- 230000003595 spectral effect Effects 0.000 description 5
- 241000700199 Cavia porcellus Species 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 229920002457 flexible plastic Polymers 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000003278 mimic effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000005309 stochastic process Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 235000015220 hamburgers Nutrition 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012353 t test Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本模拟器是一种雷达模拟技术,解决在实验室中产生合适的雷达视频回波信号的问题,它由目标、噪声、杂波产生电路及时序输出电路组成。杂波的产生采用了以最大熵算法算出的系数构成的自回归滤波器,以白噪声作为输入,使产生的杂波满足指定的功率密度谱要求,且具有随机性。本模拟器的杂波模拟方法适用于任何一种有色噪声的模拟,并且是可编程的。本模拟器可用于雷达信号处理设备的测试和调整。
Description
本发明是雷达模拟技术中视频回波信号的模拟技术。为了调整各种雷达信号处理设备,测试这些设备的性能指标是否达到设计要求,采用雷达回波模拟器是较好的方案。它能够在实验室中模拟雷达站实地的地形条件和气象条件,产生合适的参数可变的视频回波信号,代替实际雷达的高、中频部分(发射机、接收机等),对检波输出后的信号处理设备进行测试和调整。因此,本发明是一种视频雷达回波信号源。
为了使模拟的回波信号接近于真实的回波,现有的模拟器通常由目标模拟电路、接收机热噪声模拟电路、地物杂波和气象杂波(统称为雷达杂波)模拟电路等部分组成。雷达杂波是一种有色噪声,对于雷达杂波的模拟遇到的困难是:(一)真实的雷达杂波回波是一个随机过程,对于一个随机过程的模拟,必须保持其随机性及各种特征;(二)对真实的雷达杂波的研究尚不完善,尤其是对于杂波谱的研究虽然已提出了各种数学模型,但没有一种模型能普遍适用。这就要求在雷达杂波的模拟中具有一定的灵活性,能对各种数学模型甚至实测结果进行模拟。现有的雷达回波模拟器,以与我们最接近的日刊《三菱电机技报》1980.7.P.35~37中介绍的“雷达回波模拟器的研制”一文为例,解决雷达杂波模拟的方案是:将所需概率分布的杂波数据存贮在软磁盘内,通过微型计算机控制将数据由软磁盘送入随机存取存贮器(RAM),然后在同步信号控制下逐一读出,再用带通滤波器将杂波分为10路,把每一路的输出用起伏信号调制,以此将起伏附加上去,这种方案可以称为“存贮型”方案。用这种方案模拟产生的杂波无随机性,也不能满足杂波过程的重要特征,如功率密度谱等,无法测试频域处理设备,如动目标检测器(MTD)等。另外,还需要很大的存贮单元。
本发明的目的在于将最大熵谱估计的方法(MEM)应用于有色噪声的模拟。作为特例,应用于雷达杂波的模拟。本发明是一种“实时计算”型方案,用这种方案模拟产生的杂波具有很好的随机性,能满足任意指定的功率密度谱要求,并且只需要极少的存贮单元。
本发明的内容,是由目标模拟电路、噪声模拟电路、杂波模拟电路、时序和输出电路组成雷达回波模拟器。雷达回波模拟器的构成如图1所示。图中杂波模拟电路的核心是一个自回归滤波器,其输入是白噪声,滤波器的加权系数决定了滤波器的频率响应。因为输入的白噪声,所以同时决定了输出信号的功率谱密度形状。如果算好一组加权系数,使得由这组系数构成的滤波器系统函数的平方具有所需要的有色噪声的情况,那么,该滤波器的白噪声响应将是满足功率密度谱要求的有色噪声序列。因此,本发明的关键在于自回归滤波器加权系数的计算方法。
为了能得到有色噪声功率谱密度与滤波器加权系数之间的关系,采用了最大熵谱估计的方法。理论分析表明,最大熵谱与自回归过程的谱是一致的。因此,可以用现有的最大熵谱估计算法对一个所需要的有色噪声过程进行谱估计,得到的系数就是滤波器的加权系数。加权系数的计算过程如图2所示。这里,先给定功率密度谱模型,开方得到幅度谱,迭加随机相位后,作反傅立叶变换(IFFT),得到一个随机序列,再进行最大熵谱估计,得到滤波器加权系数。如果有实测杂波序列,可直接做最大熵谱估计,得到加权系数。
自回归滤波器的基本结构如图3所示。它有M个延迟单元,延迟后的序列值与对应的加权系数在乘法器内相乘,通过加法器全部相加再加上白噪声,得到输出,同时反馈至第一个延迟单元,由于实际电路中仅能完成实数运算,因而输出的实部和虚部必须分开计算,计算式如下:
式中:ai为自回归滤波器加权系数
M为滤波器阶数
Xn为输出序列
εn为输入白噪声序列
Xn-1为输出序列经过i个延迟单元后的延迟值。
滤波器实际电路结构如图4所示,它全部采用数字器件,电路中采用一片数字乘法-累加器,完成时分的乘法和加法运算。序列延迟值Xn-1迟值Xn-1虚部、白噪声实部、白噪声虚部通过选择器(3)选择进入乘法累加器与存放在加权系数存贮器中的加权系数相乘累加,加权系数ai存放在可编程只读存贮器(EPROM)中,由地址电路控制取出进入乘法累加器,运算的结果分别打入实部输出寄存器和虚部输出寄存器,得到了实部输出和虚部输出。输出值同时反馈到选择器(1)和(2),这两个选择器完成延迟-循环选择。在运算过程中,该选择器和延迟-循环寄存器构成循环移位电路,将序列延迟值Xn-1依次送到选择器(3),在运算结束后,该选择器和延迟-循环寄存器构成延迟移位电路,吐旧值纳新值,整个运算过程用一微程序控制器作时序控制,由微程序地址电路控制取出微程序。
在模拟雷达杂波时,只要采用雷达杂波的功率密度谱模型,算出加权系数后写入加权系数存贮器,就可得到所要模拟的雷达杂波输出。加权系数一次可以存入数组,存放在存贮器的不同区域,改变存贮器的高位地址,可得到不同模型的杂波输出。
本发明与现有技术相比具有以下的优点:
1.采用的最大熵法是一种普遍的方法,可以适于任何有色噪声的模拟;
2.用最大熵法模拟的雷达杂波,能满足指定的功率密度谱模型或实测结果的模拟;
3.可以以编程方式改变杂波的特性;
4.以随机白噪声作为输入,使模拟的有色噪声具有随机性;
5.加权系数仅需要2M+1存贮单元,M为滤波器阶数;
6.全部采用数字器件,工作稳定。
为了实现本发明,可以按图4电路制作杂波模拟电路。加权系数的计算可以在任何型号的计算机上以任何高级语言编写程序实现。最大熵算法可以采用伯格(Burg)算法,阶数M推荐为6,将算得的加权系数在单板计算机上或专用编程器上写入可编程只读存贮器(EPROM)中,然后插入电路,即可完成模拟,如果要改变杂波特性,可以将算得的系数写在另一片存贮器中,置换后即可得到改变后的杂波特性。
Claims (4)
1、一种雷达回波模拟器,由目标模拟电路、噪声模拟电路、杂波模拟电路、时序和输出电路组成,其特征在于杂波模拟电路由微程序地址电路、微程序控制器、加权系数存贮器及地址电路、乘法累加器、3个选择器、延迟-循环寄存器和输出寄存器构成,其中微程序地址电路接至微程序控制器,微程序控制器用作整个运算过程的时序控制;乘法累加器的输入一路来自加权系数存贮器,其地址由地址电路提供;另一路来自选择器(3),它选择白噪声实、虚部和延迟-循环寄存器的输出之一进入乘法累加器;乘法累加器的输出接至输出寄存器,其输出又反馈至选择器(1)和(2),它们和延迟-循环寄存器一起完成延迟-循环移位操作。整个杂波模拟电路的输入为随机白噪声,加权系数由最大熵算法编程得到。
2、按照权利要求1所述的杂波模拟电路,其特征在于采用了最大熵算法编程得到的系数来作为该电路中自回归滤波器加权系数的技术。在计算系数时,可以指定任意杂波功率密度谱模型,也可以对实测杂波序列进行计算。
3、按照权利要求1所述的杂波模拟电路,其特征在于可以以编程方式改变所模拟的杂波特性;并且可以一次存入数组加权系数,通过改变存贮器的高位地址,可得到不同模型的杂波输出。
4、按照权利要求1所述的杂波模拟电路,其特征在于将最大熵算法应用于有色噪声的模拟。作为特例,应用于雷达杂波的模拟。因此,该电路适于任何有色噪声的模拟。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN85105633.4A CN1003746B (zh) | 1985-07-25 | 1985-07-25 | 最大熵可编程雷达回波模拟器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN85105633.4A CN1003746B (zh) | 1985-07-25 | 1985-07-25 | 最大熵可编程雷达回波模拟器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN85105633A CN85105633A (zh) | 1987-01-28 |
CN1003746B true CN1003746B (zh) | 1989-03-29 |
Family
ID=4794605
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN85105633.4A Expired CN1003746B (zh) | 1985-07-25 | 1985-07-25 | 最大熵可编程雷达回波模拟器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1003746B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102435984B (zh) * | 2011-09-09 | 2013-10-02 | 陕西长岭电子科技有限责任公司 | 一种多普勒雷达回波模拟装置及其实现方法 |
CN106802593B (zh) * | 2016-12-20 | 2019-03-26 | 上海交通大学 | 雷达回波模拟器高精度延时控制方法及雷达回波模拟器 |
CN113945893A (zh) * | 2021-12-20 | 2022-01-18 | 清华大学 | 雷达杂波信号的模拟装置及方法 |
-
1985
- 1985-07-25 CN CN85105633.4A patent/CN1003746B/zh not_active Expired
Also Published As
Publication number | Publication date |
---|---|
CN85105633A (zh) | 1987-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Vorländer | Simulation of the transient and steady‐state sound propagation in rooms using a new combined ray‐tracing/image‐source algorithm | |
US3610901A (en) | Digital modified discrete fourier transform doppler radar processor | |
CN110208760B (zh) | 一种基于时域上采样的雷达回波仿真方法 | |
CN102831895A (zh) | 一种用现场可编程门阵列实现mfcc参数提取的方法 | |
CN107608935A (zh) | 基于时间重排压缩变换的冲击类信号时频分析与重构方法 | |
CN111830493A (zh) | 水中目标收发分置中高频声目标强度预报系统及方法 | |
CN112733364B (zh) | 一种基于阻抗矩阵分块的箔条云散射快速计算方法 | |
CN106646403A (zh) | K分布雷达杂波实时模拟方法及系统 | |
Chen et al. | Development of a statistical procedure for detecting the number of signals in a radar measurement | |
CN106802409A (zh) | 基于多cpu处理的外辐射源雷达实时信号处理方法 | |
CN103760535A (zh) | 一种高分辨雷达目标回波信号产生方法 | |
CN1003746B (zh) | 最大熵可编程雷达回波模拟器 | |
US6002914A (en) | Method and apparatus for simulating reverberation in a multipath sonar system | |
CN109446478A (zh) | 一种基于迭代和可重构方式的复协方差矩阵计算系统 | |
US5983067A (en) | Method and apparatus for simulating cross-correlation coefficients in a multipath sonar system | |
CN111654264A (zh) | 一种信号数据模拟器生成信号脉冲序列方法及系统 | |
CN101727514B (zh) | 一种被动雷达接收机仿真方法 | |
Simić et al. | An FPGA based implementation of a CFAR processor applied to a pulse-compression radar system | |
US5978646A (en) | Method and apparatus for simulating a lofargram in a multipath sonar system | |
CN111190152B (zh) | 一种球面多通道空间分布模拟角度的设计方法 | |
CN113702923A (zh) | 一种雷达杂波模拟方法及电子设备 | |
US5995803A (en) | Method and apparatus for simulating a multipath sonar system | |
RU2818373C1 (ru) | Способ имитации сигнально-помеховой обстановки | |
CN115562969B (zh) | 神经网络处理器仿真评估方法、系统、电子设备及介质 | |
CN116520270B (zh) | 一种基于评估模型的雷达电子战测试方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C13 | Decision | ||
GR02 | Examined patent application | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |