一种MIMO系统码型分类判决门限的简便求解方法技术方案

本发明专利技术提供了一种MIMO系统码型分类判决门限的简便求解方法，步骤包括：识别统计量构建、建立判决门限求解方程、计算归一化常数以及求解判决门限。该简便求解方法利用极值定理，针对不完全Gamma函数的求逆问题，用Gumbel分布函数近似独立同分布Gamma随机变量的最大值分布，通过计算Gumbel分布函数归一化常数bm得到方程的闭合解；利用该简便求解方法可以得到具有2根以上接收天线情形下判决门限的闭合解；该简便求解方法的计算效率高，且易于硬件实现，在一定条件下可以拓展到其它相关处理中。

A simple method for determining the decision threshold of code classification in MIMO system

The invention provides a simple method for solving the decision threshold of MIMO system code type classification. The steps include identifying statistics construction, establishing decision threshold solving equation, calculating normalized constant and solving decision threshold. This simple method uses the extremum theorem to solve the inverse problem of incomplete Gamma function. The Gumbel distribution function is used to approximate the maximum distribution of independent and identically distributed Gamma random variables. The closed solution of the equation is obtained by calculating the normalization constant BM of the Gumbel distribution function. In the case of receiving antenna, the closed solution of decision threshold is obtained. The simple solution method has high computational efficiency and is easy to be implemented by hardware. Under certain conditions, it can be extended to other related processing.

【技术实现步骤摘要】
一种MIMO系统码型分类判决门限的简便求解方法
本专利技术涉及一种判决门限的简便求解方法，尤其是一种MIMO系统码型分类判决门限的简便求解方法。
技术介绍
近年来，服从中心卡方分布的随机向量最大值作为一种重要的识别统计量被广泛应用于调制方式估计、编码类型识别以及GNSS接收机的卫星捕获中。对于MIMO系统中的编码类型，可通过计算任意两根天线互相关函数的最大值进行识别。针对接收天线为2的情形，其识别判决门限可以通过闭合解析式计算得到。然而，当接收天线的数量超过2时，则互相关函数最大值的累积分布函数将包含不完全Gamma函数，若要求取识别门限，必然需要涉及不完全Gamma函数逆的求解问题，但在接收天线数量超过2的情形下，很难给出其解析解，因此实际应用中通常利用牛顿迭代，Bisection等数值解法进行求取。尽管数值解法具有较高的计算精度，但其计算复杂度也较高，不便在DSP平台中实现，从而限制了其工程应用前景。
技术实现思路
本专利技术要解决的技术问题是现有的求解方法计算复杂度高，不便在DSP平台中实现，限制了其工程应用前景。为了解决上述技术问题，本专利技术提供了一种MIMO系统码型分类判决门限的简便求解方法，包括以下步骤：步骤1，识别统计量构建：针对MIMO系统中编码类型的分类识别问题，构建识别统计量γ，将MIMO系统中编码类型的分类识别问题转化为二元假设检验问题；步骤2，建立判决门限求解方程：利用Gumbel函数代替门限求解中的不完全Gamma函数，并对判决门限的解析式用一阶泰勒级数展开，得到判决门限求解的近似方程；步骤3，计算归一化常数：根据极值理论，计算归一化常数bm；步骤4，求解判决门限：利用步骤2中的判决门限求解方程以及步骤3中得到的归一化常数结果求解判决门限的闭合解。进一步地，步骤1中，将分类识别问题转化为二元假设检验问题的具体步骤为：针对MIMO系统中编码类型的分类识别问题，以任意两根接收天线之间互相关函数的最大值为统计量，统计量可以定义为：Υ＝maxFc(τ)式中，Fij(τ)为任意两根接收天线之间的互相关函数，Nr为接收天线数量，τ为延时量，且Fc(τ)服从自由度为2Nc＝Nr(Nr-1)的中心卡方分布；根据恒虚警准则，设虚警概率为pfa，则通过下式求得判决门限λ：(1-pfa)1/(N+v)＝P(Nc,λ/2)式中，P(Nc，λ/2)为不完整Gamma函数，N为接收信号的样本长度，v为循环前缀的长度，v值设置为N/4；将MIMO系统中编码类型的分类识别转化为如下的二元假设检验问题：H0：编码方案为空间复用码H1：编码方案为空时分组码若Υ＜λ，判为H0，反之，则判为H1。进一步地，步骤2中，得到判决门限求解的近似方程的具体步骤为：由于服从中心卡方分布的随机向量的累积分布函数和概率密度函数由下式求得：F(x；Nc,2)＝P(Nc,x/2)式中，F(x；Nc,2)是Gamma分布，其形状参数为Nc，尺度参数为2；因此有：(1-pfa)1/(N+v)＝F(λ；Nc,2)于是判决门限通过下式求得：λ＝F-1[(1-pfa)1/(N+v)]对于底分布函数为F(x；Nc,2)的随机向量，用Mn＝max(X1,X2,....Xn)表示该随机序列的最大值，若存在常数则式成立，式中G(x)＝exp{-e-x},x∈R为Gumbel分布的分布函数，而归一化常数bn为：由于pfa<<1，利用一阶泰勒级数将(1-pfa)1/(N+v)展开得到：显然，(N+v)/Pfa＞＞1，因此取其整数部分，令m≈int[(N+v)/Pfa]，同时式中的高阶分量忽略，上式进一步转化为：因此，判决门限求解的近似方程为：进一步地，步骤3中，计算归一化常数bm的具体步骤为：根据接收信号长度N、循环前缀的长度v、虚警概率为pfa、接收天线数量Nr以及等效整数m计算归一化常数bm为：bm＝2{ln[m/Γ(Nc)]+(Nc-1)lnBm+[(Nc-1)2lnBm-(Nc-1)2ln(Nc-1)+Nc-1]/Bm}式中，Bm＝ln[m/Γ(Nc)]+(Nc-1)ln(Nc-1)，Г(x)为Gamma函数。进一步地，步骤4中，求解判决门限的闭合解的公式为：将步骤3中所计算的归一化常数bm代入上式，即可得到判决门限λ的闭合解。本专利技术的有益效果在于：该简便求解方法利用极值定理，针对不完全Gamma函数的求逆问题，用Gumbel分布函数近似独立同分布Gamma随机变量的最大值分布，通过计算Gumbel分布函数归一化常数bm得到方程的闭合解；利用该简便求解方法可以得到具有2根以上接收天线情形下判决门限的闭合解；该简便求解方法的计算效率高，且易于硬件实现，在一定条件下可以拓展到其它相关处理中。附图说明图1为本专利技术的方法流程图；图2为本专利技术的不同接收天线数量条件下MIMO系统中空时分组码识别判决门限计算结果图；图3为本专利技术的不同样本数量时算法的性能比较。具体实施方式如图1所示，本专利技术公开了一种MIMO系统码型分类判决门限的简便求解方法，包括以下步骤：步骤1，识别统计量构建：针对MIMO系统中编码类型的分类识别问题，构建识别统计量γ，将MIMO系统中编码类型的分类识别问题转化为二元假设检验问题；步骤2，建立判决门限求解方程：利用Gumbel函数代替门限求解中的不完全Gamma函数，并对判决门限的解析式用一阶泰勒级数展开，得到判决门限求解的近似方程；步骤3，计算归一化常数：根据极值理论，计算归一化常数bm；步骤4，求解判决门限：利用步骤2中的判决门限求解方程以及步骤3中得到的归一化常数结果求解判决门限的闭合解。进一步地，步骤1中，将分类识别问题转化为二元假设检验问题的具体步骤为：针对MIMO系统中编码类型的分类识别问题，以任意两根接收天线之间互相关函数的最大值为统计量，统计量可以定义为：Υ＝maxFc(τ)式中，Fij(τ)为任意两根接收天线之间的互相关函数，Nr为接收天线数量，τ为延时量，且Fc(τ)服从自由度为2Nc＝Nr(Nr-1)的中心卡方分布；根据恒虚警准则，设虚警概率为pfa，则通过下式求得判决门限λ：(1-pfa)1/(N+v)＝P(Nc,λ/2)式中，P(Nc，λ/2)为不完整Gamma函数，N为接收信号的样本长度，v为循环前缀的长度，v值设置为N/4；将MIMO系统中编码类型的分类识别转化为如下的二元假设检验问题：H0：编码方案为空间复用码H1：编码方案为空时分组码若Υ＜λ，判为H0，反之，则判为H1。进一步地，根据步骤1可知，MIMO系统中码型分类识别的判决门限求解问题可以转化为对不完全Gamma函数逆的求解，为了提出不完全Gamma函数逆的闭合解求解方法，需对方程中所涉及的参数进行分析，在实际应用中，虚警概率pfa为小数且数值较小，如0.01、0.001等；Nc为整数，当Nr＞2时其数值大于3；N+v大于接收信号的长度，为了得到不完全Gamma逆函数的解析表达式，可以引入极值理论，将该问题转化为对服从中心卡方分布随机向量的最大值极限分布的归一化常数bm的计算，因此，步骤2中，得到判决门限求解的近似方程的具体步骤为：由于服从中心卡方分布的随机向量的累积分布函数和概率密度函数由下式求得：F(x；Nc,2)＝P(N本文档来自技高网...

【技术保护点】
1.一种MIMO系统码型分类判决门限的简便求解方法，其特征在于，包括以下步骤：步骤1，识别统计量构建：针对MIMO系统中编码类型的分类识别问题，构建识别统计量γ，将MIMO系统中编码类型的分类识别问题转化为二元假设检验问题；步骤2，建立判决门限求解方程：利用Gumbel函数代替门限求解中的不完全Gamma函数，并对判决门限的解析式用一阶泰勒级数展开，得到判决门限求解的近似方程；步骤3，计算归一化常数：根据极值理论，计算归一化常数bm；步骤4，求解判决门限：利用步骤2中的判决门限求解方程以及步骤3中得到的归一化常数结果求解判决门限的闭合解。

【技术特征摘要】
1.一种MIMO系统码型分类判决门限的简便求解方法，其特征在于，包括以下步骤：步骤1，识别统计量构建：针对MIMO系统中编码类型的分类识别问题，构建识别统计量γ，将MIMO系统中编码类型的分类识别问题转化为二元假设检验问题；步骤2，建立判决门限求解方程：利用Gumbel函数代替门限求解中的不完全Gamma函数，并对判决门限的解析式用一阶泰勒级数展开，得到判决门限求解的近似方程；步骤3，计算归一化常数：根据极值理论，计算归一化常数bm；步骤4，求解判决门限：利用步骤2中的判决门限求解方程以及步骤3中得到的归一化常数结果求解判决门限的闭合解。2.根据权利要求1所述的MIMO系统码型分类判决门限的简便求解方法，其特征在于，步骤1中，将分类识别问题转化为二元假设检验问题的具体步骤为：针对MIMO系统中编码类型的分类识别问题，以任意两根接收天线之间互相关函数的最大值为统计量，统计量可以定义为：Υ＝maxFc(τ)式中，Fij(τ)为任意两根接收天线之间的互相关函数，Nr为接收天线数量，τ为延时量，且Fc(τ)服从自由度为2Nc＝Nr(Nr-1)的中心卡方分布；根据恒虚警准则，设虚警概率为pfa，则通过下式求得判决门限λ：(1-pfa)1/(N+v)＝P(Nc,λ/2)式中，P(Nc，λ/2)为不完整Gamma函数，N为接收信号的样本长度，v为循环前缀的长度，v值设置为N/4；将MIMO系统中编码类型的分类识别转化为如下的二元假设检验问题：H0：编码方案为空间复用码H1：编码方案为空时分组码若Υ＜λ，判为H0，反之，则判为H1。3.根据权利要求2所述的MIMO系统码型分类判决门限的简便求解方法，其特征在于，步骤2中，得到判决门限求...

【专利技术属性】
技术研发人员：吴珊珊，胡国兵，王书旺，周波，聂佰玲，汤滟，宋冰，丁宁，
申请(专利权)人：南京信息职业技术学院，
类型：发明
国别省市：江苏,32