CN113406583B - 一种云mimo雷达目标检测概率的近似计算方法 - Google Patents

一种云mimo雷达目标检测概率的近似计算方法 Download PDF

Info

Publication number
CN113406583B
CN113406583B CN202110689341.4A CN202110689341A CN113406583B CN 113406583 B CN113406583 B CN 113406583B CN 202110689341 A CN202110689341 A CN 202110689341A CN 113406583 B CN113406583 B CN 113406583B
Authority
CN
China
Prior art keywords
detection
probability
quantization
approximate
cloud
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
Application number
CN202110689341.4A
Other languages
English (en)
Other versions
CN113406583A (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.)
Yangtze River Delta Research Institute of UESTC Huzhou
Original Assignee
Yangtze River Delta Research Institute of UESTC Huzhou
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 Yangtze River Delta Research Institute of UESTC Huzhou filed Critical Yangtze River Delta Research Institute of UESTC Huzhou
Priority to CN202110689341.4A priority Critical patent/CN113406583B/zh
Publication of CN113406583A publication Critical patent/CN113406583A/zh
Application granted granted Critical
Publication of CN113406583B publication Critical patent/CN113406583B/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
    • G01S7/40—Means for monitoring or calibrating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Probability & Statistics with Applications (AREA)
  • Software Systems (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Databases & Information Systems (AREA)
  • Operations Research (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及雷达技术领域，其公开了一种云MIMO雷达目标检测概率的近似计算方法，减少接收器与融合中心之间的通信负担，并可以获得闭式的检测概率。该方法包括：计算接收机与发射机之间的各个路径的检测统计量；对检测统计量进行量化处理，获得量化结果，各接收机将量化数据传输给融合中心；确定H1假设和H0假设情况下，融合中心处的总检测统计量；计算融合中心处的总检测统计量的概率分布函数；计算目标检测概率；针对量化间隔为Δ的均匀量化器，计算当Δ→0时均匀量化器输出结果的特征函数以及概率密度函数；确定云MIMO雷达中采用均匀量化器输出量化结果的近似概率密度函数；确定总检测统计量的近似分布，获得量化近似情况下的检测门限，计算近似的目标检测概率。

Description

一种云MIMO雷达目标检测概率的近似计算方法
技术领域
本发明涉及雷达技术领域，具体涉及一种云MIMO雷达目标检测概率的近似计算方法。
背景技术
在云MIMO(Multiple Input Multiple Out)雷达中，从目标返回的信号由每个本地传感器接收，它们可以通过回程网络与融合中心(FC)通信。本地传感器数量可能很大，并且接收天线可能分布在没有可用有线回程网络的较大地理区域中，因此本地传感器通过无线网络与FC通信，并且接收器采集的样本通常在每个接收器被量化后再发送给FC，以减轻通信负担。
通常，量化输出都是离散值，但是对于某些特殊的量化器，例如均匀量化器，有时可以将输出建模为输入加上量化噪声的连续值。在文献1(A.Sripad and D.Snyder,“Anecessary and sufficient condition for quantization errors to be uniformand white,”IEEE Transactions on Acoustics,Speech,and Signal Processing,vol.25,no.5,pp.442–448,1977)中，给出了标量均匀量化器输出可建模为输入加上均匀白的量化噪声的充要条件，并指出当量化间隔足够小时，零均值高斯输入信号的量化误差可以近似为均匀分布的噪声。在文献2(S.Khalili,O.Simeone,and A.M.Haimovich,“Cloudradio-multistatic radar:joint optimization of code vector and backhaulquantization,”IEEE Signal Processing Letters,vol.22,no.4,pp.494–498,2015)和文献3(W.Jiang and A.M.Haimovich,“Joint optimization of waveform andquantization in spectral congestion conditions,”in 201852nd AsilomarConference on Signals,Systems,and Computers,pp.1894–1898,2018)中，将最佳矢量量化器的量化输出表示为输入矢量加高斯量化噪声矢量。高斯量化误差近似是引入了加性高斯噪声来讨论量化影响，即将量化输出建模为输入加上高斯噪声，从而使得输出仍旧是连续随机变量，更容易进一步进行理论分析。
目标检测作为雷达的关键功能之一，在文献4(S.Khalili,O.Simeone,andA.M.Haimovich,“Cloud radio-multistatic radar:joint optimization of codevector and backhaul quantization,”IEEESignal Processing Letters,vol.22,no.4,pp.494–498,April 2015)和文献5(S.Jeong,O.Simeone,A.Haimovich,and J.Kang,“Optimization ofmultistatic cloud radar with multiple-access wirelessbackhaul,”in 2015IEEE Radar Conference(RadarCon),May 2015,pp.1650–1655)中，针对具有单个发射器和多个接收器的云雷达，利用量化的接收信号研究了最大化检测性能的发射波形优化。在文献6(W.Jiang and A.M.Haimovich,“Waveform optimization incloud radar with spectral congestion constraints,”in 2019IEEE RadarConference(RadarConf),April 2019,pp.1–6)中，基于量化的接收信号研究了云雷达系统中优化雷达波形和中继增益以使平均互信息最大。在文献7(Z.Wang,Q.He,and R.S.Blum,“Sampling rate and bits per sample tradeofffor cloud MIMO radar targetdetection,”in 201927th European Signal Processing Conference(EUSIPCO),2019,pp.1–5)中，利用量化的接收信号得到了云MIMO雷达的检测性能，并提出了采样率和每个样本的比特率之间的折中。
由于云MIMO雷达中天线数目较多且距离一般相距甚远，所以本地接收机通过无线信道将数据传输给融合中心。然而无线资源越来越紧张，需要在本地接收器对数据进行量化后再传输到融合中心。但是基于量化的离散输出，往往很难得到闭式的检测性能。因此，在云MIMO雷达中讨论量化的近似输出并得到闭式的检测概率是十分必要的。在现有的大部分基于量化数据的分布式参数估计文献中，还没有关于均匀量化器的非零均值高斯输入时量化输出近似分布的详细讨论。
而且，在云MIMO雷达中现有的讨论都是基于量化接收信号的，但当接收样本较多时，也可能带来不可接受的通信负担，需在本地接收器对接收信号进行一定的压缩。因此，在云MIMO雷达中进一步讨论一些减轻通信负担的算法是十分有必要的。
发明内容
本发明所要解决的技术问题是：提出一种云MIMO雷达目标检测概率的近似计算方法，减少接收器与融合中心之间的通信负担，并可以获得闭式的检测概率，以便后续进行理论分析和指导系统设计。
本发明解决上述技术问题采用的技术方案是：
一种云MIMO雷达目标检测概率的近似计算方法，包括以下步骤：
S1、云MIMO雷达通过接收机接收来自发射机的信号；
S2、计算接收机与发射机之间的各个路径的检测统计量；
S3、对检测统计量进行量化处理，获得量化结果，各接收机将量化数据传输给融合中心；
S4、确定H1假设和H0假设情况下，融合中心处的总检测统计量，其中，H1假设表示目标出现，H0假设表示目标未出现；
S5、计算H1假设和H0假设情况下，融合中心处的总检测统计量的概率分布函数；
S6、计算目标检测概率；
S7、针对量化间隔为Δ的均匀量化器，计算当Δ→0时均匀量化器输出结果的特征函数以及概率密度函数；
S8、确定当Δ→0时云MIMO雷达中采用均匀量化器输出量化结果的近似概率密度函数；
S9、确定总检测统计量的近似分布，获得量化近似情况下的检测门限，根据所述量化近似情况下的检测门限计算近似的目标检测概率。
作为进一步优化，步骤S1具体包括：
云MIMO雷达kTs时刻第n个接收机接收到来自第m个发射机的信号rnm[k]为
Figure GDA0003665591510000031
其中，第m个发射机在kTs时刻的采样值为
Figure GDA0003665591510000032
E是发射总能量，Ts为采样间隔，k是采样序号,k＝1,…,K，K为总的样本个数，fs为采样率，
Figure GDA0003665591510000033
为确定已知的目标反射系数；τnm，fnm表示对应于第n个接收机接收到第m个发射机信号路径的时延和多普勒频率，wnm[k]是第n个接收机接收到第m个发射机信号路径的噪声，且
Figure GDA0003665591510000034
()*表示复数的共轭，δ(k-k′)为狄拉克函数，k＝k'时函数取值为1，其他情况为0，k和k′代表的是采样时刻；M表示发射机的总个数；sm表示第m个发射机的发射信号，
Figure GDA0003665591510000035
表示噪声wnm[k]的方差。
作为进一步优化，步骤S2具体包括：
第nm个传输路径的检测统计量Tnm计算如下：
Figure GDA0003665591510000036
其中n＝1,…,N，m＝1,…,M，N和M表示接收机和发射机的总个数，
Figure GDA0003665591510000037
Re{·}表示对一个复数取实部。
作为进一步优化，步骤S3具体包括：
对检测统计量Tnm进行量化处理，得到量化后的结果qnm为
Figure GDA0003665591510000038
其中
Figure GDA0003665591510000039
表示对一个实数的量化操作；
Figure GDA0003665591510000041
其中D＝2b是量化总数，b是量化比特数，γ0,…,γD为量化器门限屏；
将第n个接收器量化后的结果排为列向量
Figure GDA0003665591510000042
其中
Figure GDA0003665591510000043
表示转置；
所有接收器将量化后的数据传送到融合中心，则融合中心收到的数据为
Figure GDA0003665591510000044
作为进一步优化，步骤S4具体包括：
在H1假设和H0假设的情况下，融合中心处的总检测统计量为：
Figure GDA0003665591510000045
其中，α是由虚警门限PFA确定的检测门限，T∈{0,1，…，NM(D-1)}，sum(y)表示对向量y的所有元素求和。
作为进一步优化，步骤S5具体包括：
计算H1和H0假设下qnm的条件概率分布函数为：
Figure GDA0003665591510000046
Figure GDA0003665591510000047
其中，d＝1,…,D-1为可能的量化结果，
Figure GDA0003665591510000048
表示在Hi，i＝0,1假设下事件
Figure GDA0003665591510000049
发生的概率，Q(·)表示标准高斯分布的累积分布函数，定义为
Figure GDA00036655915100000410
计算H1和H0假设下融合中心处总的检测统计量T的概率分布函数为：
Figure GDA0003665591510000051
Figure GDA0003665591510000052
其中，z∈{0,1，…，NM(D-1)}，znm∈{0,1，…，(D-1)}，n＝1,…,N,m＝1,…,M。
作为进一步优化，步骤S6具体包括：
根据以下公式计算检测概率：
Figure GDA0003665591510000053
其中，α是由虚警门限PFA确定的检测门限，
Figure GDA0003665591510000054
Figure GDA0003665591510000055
作为进一步优化，步骤S7具体包括：
针对量化间隔为Δ的均匀量化器，量化器输入为均值μ方差σ2的高斯信号x时，
当Δ/σ→0时量化误差ε＝q-x的近似概率密度函数为
Figure GDA0003665591510000056
其中q为量化器的输出结果；
根据下式得到当Δ→0时均匀量化器输出q的特征函数：
Figure GDA0003665591510000057
其中，
Figure GDA0003665591510000058
为量化误差ε的特征函数，
Figure GDA0003665591510000059
为量化器的输入x的特征函数；
根据下式得到当Δ→0时均匀量化器输出q的概率密度函数：
Figure GDA00036655915100000510
作为进一步优化，步骤S8具体包括：
计算获得当Δ→0时云MIMO雷达中采用均匀量化器时量化器输出qnm的近似概率密度函数：
Figure GDA0003665591510000061
Figure GDA0003665591510000062
其中
Figure GDA0003665591510000063
作为进一步优化，步骤S9具体包括：
计算得到检测统计量T的近似分布为：
Figure GDA0003665591510000064
Figure GDA0003665591510000065
其中
Figure GDA0003665591510000066
量化近似情况下的检测门限为：
Figure GDA0003665591510000067
其中PFA为虚警概率，Q-1(·)为Q(·)的反函数；
根据下式得到近似的检测概率为：
Figure GDA0003665591510000068
本发明的有益效果是：
由于在本地传感器首先计算得到检测统计量以减少需要传输的样本，可大大降低通信负担，更利于工程实现。且可以得到闭式的检测概率，更方便进行后续的理论分析和系统设计，可在实际应用中指导系统设计。
附图说明
图1为本发明的云MIMO雷达目标检测概率的近似计算方法流程图；
图2为实施例中接收机和发射机数量分别为N＝4和M＝8时的系统参数设置图；
图3为实施例中不同量化比特下检测概率PD随虚警概率PFA变化的示意图。
具体实施方式
为了方便描述，首先进行如下定义：
Figure GDA0003665591510000071
为转置，()H为共轭转置，()*表示复数的共轭，
Figure GDA00036655915100000712
表示数学期望，Re{·}表示对一个复数取实部。
在云MIMO雷达中，采用空间分置的M个发射机和N个接收机，在笛卡尔坐标系中，第m(m＝1,…,M)个发射天线和第n(n＝1,…,N)个接收天线分别位于
Figure GDA0003665591510000072
和
Figure GDA0003665591510000073
第m个发射机在kTs时刻的采样值为
Figure GDA0003665591510000074
E是发射总能量，Ts为采样间隔，k(k＝1,…,K)是采样序号，K为总的样本个数，fs为采样率，假设不同发射机的发射信号是正交的。假设存在目标时，目标位于(x,y)，运动速度为(vx,vy)，所以在kTs时刻第n个接收机接收到第m个发射机的信号为，
Figure GDA0003665591510000075
其中噪声wnm[k]假定为零均值的复高斯白噪声，且对不同的采样点k和k′有
Figure GDA0003665591510000076
为简化分析，假定发射系数
Figure GDA0003665591510000077
是确定已知参数。τnm表示与第nm条路径相关的延时，fnm表示与第nm条路径相关的接收信号多普勒频移。
将接收机n接收到的来自发射机m的所有时间样本信号排为向量rnm
Figure GDA0003665591510000078
其中，
Figure GDA0003665591510000079
将所有接收机收到的样本信号排为向量r
Figure GDA00036655915100000710
其中，
Figure GDA00036655915100000711
在H1假设(目标出现在测试单元中)和H0假设(目标不存在)下，集中式处理的检测问题为：
Figure GDA0003665591510000081
因此，可以得到两种假设下接收信号r的概率密度函数为
Figure GDA0003665591510000082
进一步可得到对数似然比为
Figure GDA0003665591510000083
根据NP准则，在集中式处理下的最优检测器为
Figure GDA0003665591510000084
其中αO是由虚警门限PFA,O确定的检测门限，
Figure GDA0003665591510000085
Figure GDA0003665591510000086
Figure GDA0003665591510000087
Figure GDA0003665591510000088
表示均值为μ方差为σ2的高斯分布。
在云MIMO雷达系统中，本地传感器通过通信回路将本地数据发送到融合中心。然后，融合中心对接收到的数据进行处理，完成目标检测。为了减少通信负担，本地数据在发送到融合中心之前需要进行量化处理。因此，式(7)中用于集中处理的最优检测器是不可能实现的，从(7)式中，我们看到最佳检测器计算了每个发射到接收路径的检测统计量Tnm，对于m＝1,…,M在本地接收器n计算和量化Tnm是合理的。量化的检测统计量被发送到融合中心，如式(7)所示，在融合中心将所有的检测统计量相加并与阈值进行比较，以做出最终决定，即在融合中心用量化的检测统计量替代(7)式中原始Tnm，这种方法可以大大地降低通信负担。
在具体实施上，如图1所示，本发明中采用如下步骤来计算基于量化检测统计量的云MIMO雷达的近似检测性能：
S1、云MIMO雷达通过接收机接收来自发射机的信号；
本步骤中，云MIMO雷达kTs时刻第n个接收机接收到来自第m个发射机的信号rnm[k]为
Figure GDA0003665591510000091
其中第m个发射机在kTs时刻的采样值为
Figure GDA0003665591510000092
E是发射总能量，Ts为采样间隔，k是采样序号,k＝1,…,K，K为总的样本个数，fs为采样率，
Figure GDA0003665591510000093
为确定已知的目标反射系数；τnm，fnm表示对应于第n个接收机接收到第m个发射机信号路径的时延和多普勒频率，wnm[k]是第n个接收机接收到第m个发射机信号路径的噪声，且
Figure GDA0003665591510000094
()*表示复数的共轭，δ(k-k′)为狄拉克函数，k＝k'时函数取值为1，其他情况为0，k和k′代表的是采样时刻；M表示发射机的总个数；sm表示第m个发射机的发射信号，
Figure GDA0003665591510000095
表示噪声wnm[k]的方差。
S2、计算接收机与发射机之间的各个路径的检测统计量；
本步骤中，计算第nm个传输路径的检测统计量Tnm，
Figure GDA0003665591510000096
其中n＝1,…,N，m＝1,…,M，N和M表示接收机和发射机的总个数，
Figure GDA0003665591510000097
Re{·}表示对一个复数取实部。
S3、对检测统计量进行量化处理，获得量化结果，各接收机将量化数据传输给融合中心；本步骤中，对检测统计量Tnm进行量化处理，得到量化后的结果qnm为
Figure GDA0003665591510000098
其中
Figure GDA0003665591510000099
表示对一个实数的量化操作，为
Figure GDA0003665591510000101
其中D＝2b是量化总数，b是量化比特数，γ0,…,γD为量化器门限；
将第n个接收器量化后的结果排为列向量：
Figure GDA0003665591510000102
其中
Figure GDA0003665591510000103
表示转置。
所有接收器将量化后的数据传送到融合中心，则融合中心收到的数据为
Figure GDA0003665591510000104
S4、确定H1假设和H0假设情况下，融合中心处的总检测统计量：
本步骤中，得到在假设目标出现在测试单元中为H1和假设目标不存在为H0的情况下，融合中心处的检测问题为
Figure GDA0003665591510000105
其中，α是由虚警门限PFA确定的检测门限，T∈{0,1，…，NM(D-1)}，sum(y)表示对向量y的所有元素求和。
S5、计算H1假设和H0假设情况下，融合中心处的总检测统计量的概率分布函数；
本步骤中，计算H1和H0假设下qnm的条件概率分布函数为
Figure GDA0003665591510000106
其中，d＝1,…,D-1为可能的量化结果，
Figure GDA0003665591510000107
表示在Hi，i＝0,1假设下事件
Figure GDA0003665591510000108
发生的概率，Q(·)表示标准高斯分布的累积分布函数，定义为
Figure GDA0003665591510000111
计算H1和H0假设下融合中心处总的检测统计量T的概率分布函数为
Figure GDA0003665591510000112
其中z∈{0,1，…，NM(D-1)}，znm∈{0,1，…，(D-1)}，n＝1,…,N,m＝1,…,M。
S6、计算目标检测概率；
本步骤中，根据以下公式计算检测概率
Figure GDA0003665591510000113
其中α是由虚警门限PFA确定的检测门限，
Figure GDA0003665591510000114
Figure GDA0003665591510000115
S7、针对量化间隔为Δ的均匀量化器，计算当Δ→0时均匀量化器输出结果的特征函数以及概率密度函数；
本步骤中，针对量化间隔为Δ的均匀量化器，量化器输入为均值μ方差σ2的高斯信号x时，
当Δ/σ→0时量化误差ε＝q-x的近似概率密度函数为：
Figure GDA0003665591510000116
其中q为量化器的输出结果。
根据下式得到当Δ→0时均匀量化器输出q的特征函数：
Figure GDA0003665591510000117
其中
Figure GDA0003665591510000118
为量化误差ε的特征函数，
Figure GDA0003665591510000119
为量化器的输入x的特征函数。
根据下式得到当Δ→0时均匀量化器输出q的概率密度函数：
Figure GDA0003665591510000121
S8、确定当Δ→0时云MIMO雷达中采用均匀量化器输出量化结果的近似概率密度函数；
本步骤中，利用步骤S7的结果可得到当Δ→0时云MIMO雷达中采用均匀量化器时量化器输出qnm的近似概率密度函数为：
Figure GDA0003665591510000122
其中
Figure GDA0003665591510000123
S9、确定总检测统计量的近似分布，获得量化近似情况下的检测门限，根据所述量化近似情况下的检测门限计算近似的目标检测概率。
本步骤中，得到检测统计量T的近似分布为：
Figure GDA0003665591510000124
其中
Figure GDA0003665591510000125
得到量化近似情况下的检测门限为：
Figure GDA0003665591510000126
其中PFA为虚警概率，Q-1(·)为Q(·)的反函数，
根据下式得到近似的检测概率为：
Figure GDA0003665591510000127
本发明的上述云MIMO雷达目标检测概率的近似计算方法的工作原理如下：
在云MIMO雷达中，在H1假设(目标出现在测试单元中)和H0假设(目标不存在)下，集中式处理的检测问题为
Figure GDA0003665591510000131
因此，可以得到两种假设下接收信号r的概率密度函数为
Figure GDA0003665591510000132
进一步可得到对数似然比为
Figure GDA0003665591510000133
根据NP准则，在集中式处理下的最优检测器为
Figure GDA0003665591510000134
其中αO是由虚警门限PFA,O确定的检测门限，
Figure GDA0003665591510000135
Figure GDA0003665591510000136
Figure GDA0003665591510000137
在云MIMO雷达系统中，本地传感器通过通信回路将本地数据发送到融合中心。然后，融合中心对接收到的数据进行处理，完成目标检测。为了减少通信负担，本地数据在发送到融合中心之前需要进行量化处理。因此，基于(7)式，我们可在本地传感器计算每个发射到接收路径的检测统计量Tnm，并将Tnm量化后发送到融合中心。为了简化分析，假设本地接收机到融合中心的通信回路是理想的，如(7)式所示，在融合中心可将所有量化后的检测统计量相加并与阈值进行比较，以做出最终决定。因此，在融合中心处的检测问题变为
Figure GDA0003665591510000138
其中α是由虚警门限PFA确定的检测门限。
根据检测统计量的分布，可得到H1和H0假设下qnm的条件概率分布函数为
Figure GDA0003665591510000141
进一步可得到两种假设下融合中心处总的检测统计量T的概率分布函数为
Figure GDA0003665591510000142
其中z∈{0,1，…，NM(D-1)}，znm∈{0,1，…，(D-1)}，n＝1,…,N,m＝1,…,M。
因此，检测概率为
Figure GDA0003665591510000143
其中α是由虚警门限PFA确定的检测门限，
Figure GDA0003665591510000144
Figure GDA0003665591510000145
由(12)式可看出，当我们将Tnm输入到量化器，量化输出结果是一个离散随机变量，这使得它难以获得检测概率的闭式表达式。因此，我们接下来讨论量化器输出的近似概率分布，并使用该近似分布来获得封闭形式的检测概率。
假设采用量化间隔为Δ的均匀量化器，量化器输入为均值μ方差σ2的高斯信号x，根据文献(A.Sripad and D.Snyder,“Anecessary and sufficient condition forquantization errors to be uniform and white,”IEEE Transactions onAcoustics,Speech,and Signal Processing,vol.25,no.5,pp.442–448,Oct 1977),量化器的输出q可以被建模为输入和量化误差ε的和，因此
ε＝q-x (36)
其中，
Figure GDA0003665591510000146
据此可以得到量化误差ε的累积分布函数为
Figure GDA0003665591510000151
根据累积分布函数，可得到量化误差ε的概率密度函数为
Figure GDA0003665591510000152
定义函数g(ε)为
Figure GDA0003665591510000153
并将g(ε)从区间ε∈[-Δ/2,Δ/2]周期扩展到ε∈[-∞,∞]，则扩展后的周期为Δ的函数可用傅里叶级数表示为
Figure GDA0003665591510000154
其中
Figure GDA0003665591510000155
为傅里叶级数，表示为
Figure GDA0003665591510000156
其中
Figure GDA0003665591510000157
表示输入x的特征函数。将(42)式带入(41)式中可得
Figure GDA0003665591510000158
将(43)式和(40)式带入(39)式中可得量化误差ε的概率密度函数为
Figure GDA0003665591510000161
由于输入x为值μ方差σ2的高斯信号，则x的特征函数为
Figure GDA0003665591510000162
将(45)式带入(44)式中可得量化误差ε的概率密度函数为
Figure GDA0003665591510000163
根据(46)式，量化误差ε的均值和方差分别为
Figure GDA0003665591510000164
和
Figure GDA0003665591510000165
从(46)式可以看出，其中的第二项
Figure GDA0003665591510000166
类似地，(48)式中的第二项满足
Figure GDA0003665591510000167
接下来我们分析当Δ/σ→0时，(49)和(50)式右端的项
Figure GDA0003665591510000168
的取值。由于指数和是非负的，因此
Figure GDA0003665591510000171
进一步可得
Figure GDA0003665591510000172
其中，当Δ/σ→0时
Figure GDA0003665591510000173
根据(51)，(52)和(53)式可得，当Δ/σ→0
Figure GDA0003665591510000174
当Δ/σ→0时，根据(49)，(59)和(54)式可将式(46)和(48)式结果简化为
Figure GDA0003665591510000175
和
E[ε2]＝Δ2/12 (56)
因此，只要Δ/σ足够小，量化误差ε可近似服从概率密度函数为(55)式的均匀分布。
由于量化器输出q＝x+ε可表示为高斯输入和均匀分布的量化误差的和，根据文献(M.Bertocco,C.Narduzzi,P.Paglierani,and D.Petri,“A noise model for digitizeddata,”IEEE Transactions on Instrumentation and Measurement,vol.49,no.1,pp.83–86,Feb 2000)可得输入和量化误差是相互独立的。因此，量化器输出q的特征函数可写为
Figure GDA0003665591510000176
其中
Figure GDA0003665591510000181
为量化误差ε的特征函数，
Figure GDA0003665591510000182
为量化器的输入x的特征函数。当Δ→0时，利用泰勒级数展开可得
Figure GDA0003665591510000183
Figure GDA0003665591510000184
因此，当Δ→0时，我们可以得到
Figure GDA0003665591510000185
将(60)式带入(57)式中可知，当Δ→0时，q的特征函数满足
Figure GDA0003665591510000186
随着Δ趋近于0，q的特征函数φq(t)逐渐趋近于φ(t)，根据文献(B.V.GnedenkoandA.N.Kolmogorov,“Limit distributions for sums ofindependent randomvariables,”IEEE Transactions on Instrumentation and Measurement,ser.Addison-Wesley Mathematics Series.Cambridge,MA:Addison-Wesley,1954)可知，特征函数φq(t)对应的概率密度函数收敛到特征函数φ(t)对应的概率密度函数。因此，可以得到，当Δ→0时均匀量化器输出q的概率密度函数可近似为
Figure GDA0003665591510000187
利用公式(62)的结果可得当Δ→0时云MIMO雷达中采用均匀量化器时量化器输出qnm的近似概率密度函数为
Figure GDA0003665591510000188
其中
Figure GDA0003665591510000191
进一步可得到检测统计量T的近似分布为
Figure GDA0003665591510000192
其中
Figure GDA0003665591510000193
利用(64)式的近似分布，可得到(32)式中检测问题的检测门限为
Figure GDA0003665591510000194
其中PFA为虚警概率，Q-1(·)为Q(·)的反函数，
因此，近似的检测概率为
Figure GDA0003665591510000195
实施例：
本实施例中，假设云MIMO雷达系统有M＝8个发射机和N＝4个接收机，其天线位置如图2所示，假设检测目标位于(150,130)m处的以(25,20)m/s的速度移动。仿真中发射信号为频率扩展高斯单脉冲信号
Figure GDA0003665591510000196
取T＝0.01s，Δf＝150Hz，fs＝600Hz。
定义
Figure GDA0003665591510000197
且设置噪声方差为
Figure GDA0003665591510000198
SNR＝-5dB。
图3绘制了基于不同量化比特数b下的检测概率PD与虚警概率PFA的关系图。图中将最优的集中式处理的结果用unquantized标注以衡量量化系统的检测性能，将步骤S6得到的基于离散量化输出分析的结果用direct标注，将步骤S9中量化输出用高斯近似表示后的近似结果用approx标注。从图中可以看出，随着b的增加量化间隔逐渐减小，量化输出采用高斯近似得到的相应曲线逐渐接近直接分析量化输出得到的曲线，当b不小于2时，近似的结果几乎与直接分析时的结果相同。因此，当量化比特足够大时，高斯近似得到的近似检测概率几乎与真实的检测概率相同。当量化位数b大于5时，量化性能几乎与未量化时的性能一致，说明此时量化的性能损失几乎可忽略。

Claims (3)

1.一种云MIMO雷达目标检测概率的近似计算方法，其特征在于，包括以下步骤：
S1、云MIMO雷达通过接收机接收来自发射机的信号；
S2、计算接收机与发射机之间的各个路径的检测统计量；
S3、对检测统计量进行量化处理，获得量化结果，各接收机将量化数据传输给融合中心；
S4、确定H1假设和H0假设情况下，融合中心处的总检测统计量，其中，H1假设表示目标出现，H0假设表示目标未出现；
S5、计算H1假设和H0假设情况下，融合中心处的总检测统计量的概率分布函数；
S6、计算目标检测概率；
S7、针对量化间隔为Δ的均匀量化器，计算当Δ→0时均匀量化器输出结果的特征函数以及概率密度函数；
S8、确定当Δ→0时云MIMO雷达中采用均匀量化器输出量化结果的近似概率密度函数；
S9、确定总检测统计量的近似分布，获得量化近似情况下的检测门限，根据所述量化近似情况下的检测门限计算近似的目标检测概率；
步骤S1具体包括：
云MIMO雷达kTs时刻第n个接收机接收到来自第m个发射机的信号rnm[k]为
Figure FDA0003665591500000011
其中，第m个发射机在kTs时刻的采样值为
Figure FDA0003665591500000012
E是发射总能量，Ts为采样间隔，k是采样序号,k＝1,…,K，K为总的样本个数，fs为采样率，
Figure FDA0003665591500000013
为确定已知的目标反射系数；τnm，fnm表示对应于第n个接收机接收到第m个发射机信号路径的时延和多普勒频率，wnm[k]是第n个接收机接收到第m个发射机信号路径的噪声，且
Figure FDA0003665591500000014
()*表示复数的共轭，δ(k-k′)为狄拉克函数，k＝k'时函数取值为1，其他情况为0，k和k′代表的是采样时刻；M表示发射机的总个数；sm表示第m个发射机的发射信号，
Figure FDA0003665591500000015
表示噪声wnm[k]的方差；
步骤S2具体包括：
第nm个传输路径的检测统计量Tnm计算如下：
Figure FDA0003665591500000021
其中n＝1,…,N，m＝1,…,M，N和M表示接收机和发射机的总个数，
Figure FDA0003665591500000022
Re{·}表示对一个复数取实部；
步骤S3具体包括：
对检测统计量Tnm进行量化处理，得到量化后的结果qnm为
Figure FDA0003665591500000023
其中
Figure FDA0003665591500000024
表示对一个实数的量化操作；
Figure FDA0003665591500000025
其中D＝2b是量化总数，b是量化比特数，γ0,…,γD为量化器门限屏；
将第n个接收器量化后的结果排为列向量
Figure FDA0003665591500000026
其中上标
Figure FDA0003665591500000027
表示转置；
所有接收器将量化后的数据传送到融合中心，则融合中心收到的数据为
Figure FDA0003665591500000028
步骤S4具体包括：
在H1假设和H0假设的情况下，融合中心处的总检测统计量为：
Figure FDA0003665591500000029
其中，α是由虚警门限PFA确定的检测门限，T∈{0,1，…，NM(D-1)}，sum(y)表示对向量y的所有元素求和；
步骤S5具体包括：
计算H1和H0假设下qnm的条件概率分布函数为：
Figure FDA0003665591500000031
Figure FDA0003665591500000032
其中，d＝1,…,D-1为可能的量化结果
Figure FDA0003665591500000033
Figure FDA0003665591500000034
表示在Hi，i＝0,1假设下事件
Figure FDA0003665591500000035
发生的概率，Q(·)表示标准高斯分布的累积分布函数，定义为
Figure FDA0003665591500000036
计算H1和H0假设下融合中心处总的检测统计量T的概率分布函数为：
Figure FDA0003665591500000037
Figure FDA0003665591500000038
其中，z∈{0,1，…，NM(D-1)}，znm∈{0,1，…，(D-1)}，n＝1,…,N,m＝1,…,M；
步骤S6具体包括：
根据以下公式计算检测概率：
Figure FDA0003665591500000039
其中，α是由虚警门限PFA确定的检测门限，
Figure FDA00036655915000000310
Figure FDA00036655915000000311
步骤S7具体包括：
针对量化间隔为Δ的均匀量化器，量化器输入为均值μ方差σ2的高斯信号x时，
当Δ/σ→0时量化误差ε＝q-x的近似概率密度函数为
Figure FDA00036655915000000312
其中q为量化器的输出结果；
根据下式得到当Δ→0时均匀量化器输出q的特征函数：
Figure FDA0003665591500000041
其中，
Figure FDA0003665591500000042
为量化误差ε的特征函数，
Figure FDA0003665591500000043
为量化器的输入x的特征函数；
根据下式得到当Δ→0时均匀量化器输出q的概率密度函数：
Figure FDA0003665591500000044
2.如权利要求1所述的一种云MIMO雷达目标检测概率的近似计算方法，其特征在于，
步骤S8具体包括：
计算获得当Δ→0时云MIMO雷达中采用均匀量化器时量化器输出qnm的近似概率密度函数：
Figure FDA0003665591500000045
Figure FDA0003665591500000046
其中
Figure FDA0003665591500000047
3.如权利要求2所述的一种云MIMO雷达目标检测概率的近似计算方法，其特征在于，
步骤S9具体包括：
计算得到检测统计量T的近似分布为：
Figure FDA0003665591500000048
Figure FDA0003665591500000049
其中
Figure FDA00036655915000000410
量化近似情况下的检测门限为：
Figure FDA0003665591500000051
其中PFA为虚警概率，Q-1(·)为Q(·)的反函数；
根据下式得到近似的检测概率为：
Figure FDA0003665591500000052
CN202110689341.4A 2021-06-22 2021-06-22 一种云mimo雷达目标检测概率的近似计算方法 Active CN113406583B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110689341.4A CN113406583B (zh) 2021-06-22 2021-06-22 一种云mimo雷达目标检测概率的近似计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110689341.4A CN113406583B (zh) 2021-06-22 2021-06-22 一种云mimo雷达目标检测概率的近似计算方法

Publications (2)

Publication Number Publication Date
CN113406583A CN113406583A (zh) 2021-09-17
CN113406583B true CN113406583B (zh) 2022-08-02

Family

ID=77682131

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110689341.4A Active CN113406583B (zh) 2021-06-22 2021-06-22 一种云mimo雷达目标检测概率的近似计算方法

Country Status (1)

Country Link
CN (1) CN113406583B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102043143A (zh) * 2010-11-09 2011-05-04 河海大学 统计mimo雷达目标检测仿真方法
CN106707273A (zh) * 2017-01-23 2017-05-24 西安电子科技大学 基于奈曼皮尔逊准则量化的多站雷达信号融合检测方法
CN108226893A (zh) * 2018-04-02 2018-06-29 电子科技大学 一种用于mimo雷达的低复杂度接收机设计方法
CN108957422A (zh) * 2018-06-01 2018-12-07 电子科技大学 一种基于量化数据的云mimo雷达的根克拉美罗下界计算方法
CN109521412A (zh) * 2018-12-26 2019-03-26 西安电子科技大学 基于局部统计量融合的雷达组网空域目标检测方法
CN110161471A (zh) * 2019-04-03 2019-08-23 电子科技大学 一种针对云mimo雷达的采样率和量化比特的计算方法
CN112363151A (zh) * 2020-11-05 2021-02-12 西安电子科技大学 频率分集阵多输入多输出雷达的自适应目标检测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102043143A (zh) * 2010-11-09 2011-05-04 河海大学 统计mimo雷达目标检测仿真方法
CN106707273A (zh) * 2017-01-23 2017-05-24 西安电子科技大学 基于奈曼皮尔逊准则量化的多站雷达信号融合检测方法
CN108226893A (zh) * 2018-04-02 2018-06-29 电子科技大学 一种用于mimo雷达的低复杂度接收机设计方法
CN108957422A (zh) * 2018-06-01 2018-12-07 电子科技大学 一种基于量化数据的云mimo雷达的根克拉美罗下界计算方法
CN109521412A (zh) * 2018-12-26 2019-03-26 西安电子科技大学 基于局部统计量融合的雷达组网空域目标检测方法
CN110161471A (zh) * 2019-04-03 2019-08-23 电子科技大学 一种针对云mimo雷达的采样率和量化比特的计算方法
CN112363151A (zh) * 2020-11-05 2021-02-12 西安电子科技大学 频率分集阵多输入多输出雷达的自适应目标检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Parameter Estimation Using Quantized Cloud MIMO Radar Measurements";Zhen Wang;《2018 IEEE 10th Sensor Array and Multichannel Signal Processing Workshop (SAM)》;20181030;全文 *
"Quantized Time Delay for Target Localization in Cloud MIMO Radar";wang zhen等;《2021 IEEE Radar Conference (RadarConf21)》;20210618;全文 *
"合作式MIMO雷达通信一体化系统的检测与互信息增益";王珍等;《信号处理》;20201020;全文 *

Also Published As

Publication number Publication date
CN113406583A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
CN108957422B (zh) 一种基于量化数据的云mimo雷达的根克拉美罗下界计算方法
CN105517019B (zh) 采用集成回归系统检测lte网络性能的方法
CN109490826B (zh) 一种基于无线电波场强rssi的测距与位置定位方法
CN109507641A (zh) 一种合作的mimo雷达和mimo通信系统性能计算方法
WO2024021440A1 (zh) 一种迭代聚焦式毫米波一体化通信与感知方法
JP2020190544A (ja) 融合された閾値／位相変調検出を使用してパルスを検出するためのシステムおよび方法
CN110161471B (zh) 一种针对云mimo雷达的采样率和量化比特的计算方法
CN109061556B (zh) 一种基于弹性网络的稀疏迭代波达角估计方法
CN113406583B (zh) 一种云mimo雷达目标检测概率的近似计算方法
CN112272064B (zh) 一种合作式mimo雷达的检测概率和互信息计算方法
CN111446998B (zh) 基于深度学习的波达方向估计方法
CN114690175B (zh) 一种基于被动外辐射源雷达的目标直接检测与跟踪方法
EP2824476A1 (en) Temporal analysis for user speed estimation in wireless networks
Wang et al. Parameter estimation using quantized cloud MIMO radar measurements
Chen et al. Data compression for multi-parameter estimation for emitter location
Mostofi et al. Compressive cooperative obstacle mapping in mobile networks
Gadhgadhi et al. Distance estimation using polynomial approximation and neural network based on rssi technique
Hong et al. Channel capacity analysis of indoor environments for location-aware communications
US20110142116A1 (en) Method and apparatus for estimating channel parameter
CN110391820B (zh) 一种基于dft的规避同频干扰的通信接收方法
Dashti et al. Optimum threshold for ranging based on toa estimation error analysis
Cheng et al. Gibbs Sampling Based Sparse Bayesian Learning for Direction-of-Arrival Estimation with Impulse Noise Towards 6G
Gormley et al. A spectrum sensor for CubeSat radios
Yan et al. Optimal local sensor decision rule for target detection with channel fading statistics in multi-sensor networks
US20150103958A1 (en) Methods, systems, and media for determining whether a signal of interest is present

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