买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

龙图腾网获悉电子科技大学申请的专利采用低精度量化的毫米波通信系统定位方法获国家发明授权专利权，本发明授权专利权由国家知识产权局授予，授权公告号为：CN115604815B 。

龙图腾网通过国家知识产权局官网在2025-09-12发布的发明授权授权公告中获悉：该发明授权的专利申请号/专利号为：202211122290.8，技术领域涉及：H04W64/00；该发明授权采用低精度量化的毫米波通信系统定位方法是由关雅静;成先涛设计研发完成，并于2022-09-15向国家知识产权局提交的专利申请。

本采用低精度量化的毫米波通信系统定位方法在说明书摘要公布了：本发明属于无线通信技术领域，具体涉及一种采用低精度量化的毫米波通信系统定位方法。本发明在信道估计阶段提出了一种交替迭代算法框架。首先通过广义Turbo算法恢复出未量化信道，并利用多任务稀疏贝叶斯学习进行信道粗估计，随后基于期望最大化进行信道细估计，迭代上述过程完成信道参数估计。位置估计阶段将位置信息得到的信道参数与算法估计的信道参数的残差作为优化函数，利用牛顿法求解该优化问题。实验表明，本发明所提出的毫米波通信系统定位方法在低精度量化的情况下，依然能够实现精确的用户定位，并达到对应的理论下界。

本发明授权采用低精度量化的毫米波通信系统定位方法在权利要求书中公布了：1.采用低精度量化的毫米波通信系统定位方法，用于大规模多输入多输出正交频分复用系统，定义系统包括单天线的用户端和具有N根天线的基站端，其特征在于，定位方法包括以下步骤： S1、定义在系统上行链路中，子载波总数为M，第m个子载波上的频域信道表示为： 其中L为多径数，cl和τl为第l条多径的复增益和时延，φl,m是相应的空间方向，定义为 φl,m＝1+fmfcdsinθlλc2 是第m个子载波的频率，W是系统带宽，fc是载波频率，λc是载波波长，θl是第l条路径的到达角，d是天线间距，设置d＝λc2；aφl,m是阵列响应向量，考虑均匀线性阵列，有 在基站端，第m个子载波的接收信号如下： ym＝hmsm+nm4 其中sm是训练符号，表示均值为0和方差为σ2的加性复高斯噪声，将sm设置为1，基于公式4，第n条天线处接收的时域信号为： 是归一化的离散傅里叶变换矩阵，第i行第j列元素为 其中hm,n是hm的第n个元素，是噪声向量； 通过低精度ADC进行采样，将量化为数字信号qn，分别用和qn,p表示和qn的第p个元素，可得： 其中，和分别表示的实部和虚部，复数量化器由两个实数量化器组成，通过量化精度为Qb比特的ADC将和映射到离散值之一： 其中-∞＝u0u1…uB＝∞是量化门限，v1v2…vB是输出电平，对于中平均匀量化器： 其中Δ是量化间隔； 得到基于低精度ADC的系统模型为： S2、利用Gturbo算法恢复出未量化的频域信道Gturbo算法包括两个模块： 模块A：计算zn的后验均值和方差，由于各元素计算方式相同，下述计算过程省略下标n: 计算的外部均值和方差: 模块B：计算的后验均值和方差: 计算zn的外部均值和方差: 将作为的估计值，根据其中是的第m个元素，则可获得hm； 将模块A和模块B迭代执行，直到收敛； S3、信道参数的估计：将信道hm表示为字典中原子的线性组合，P表示到达角度在[-1,1]的均匀采样网格点数，将公式4重写为： ym＝Dmxm+nm20 其中，xm∈CN×1是只有L个非零元素的稀疏向量， 令xm服从相同的均值为0，方差为αx -1的复高斯分布，αx＝[αx,1,αx,2,...,αx,P]T，则有：其中Λx＝diag{αx}，噪声nm中每个元素服从均值为0，方差为β-1的相同的复高斯分布，则xm的后验分布也为复高斯分布，其均值和方差为： 更新αx和β的公式为： 其中Vxmp,p表示矩阵Vxm的第p个对角元素，按照公式24和公式25所示的更新表达式进行迭代更新，得到μxm后，即可取K个最大非零元的位置作为进一步估计的候选径，根据该位置保留字典Dm相应的列和μxm中相应的行，得到降维的字典和xm的估计值 S4、信道参数的细估计：对真正的未知字典线性近似引入量化误差，定义θk为候选径在字典中对应的角度值，将接收信号重构为： 其中且 表示对θk求导，γ和η表示候选径的复增益和时延，且δk∈[-1P,1P]，初始值的计算如下 表示的第k个元素； 利用EM算法估计{γ,δ,η,α,β}，在EM算法的E步骤中，更新γ的后验均值和方差： 其中，Λ＝diagα，将后验均值μ作为γ的估计值； 在EM算法的M步骤中，更新参数集{δ,η,α,β}，完全似然函数的对数期望为： 其中，Vk,k表示矩阵V的第k个对角元素，tr·表示矩阵的迹，const代表常数项，参数的更新公式如下 其中 其中 Bm＝Φmdiag{μ}36 E步和M步迭代更新，直至收敛；角度θ更新方式为θt+1＝θt+δ，θt表示第t次更新时的角度，将恢复出来的信号的均值和方差传回到Gturbo中，作为hn的先验均值和方差； S5、位置参数的估计，根据几何关系得： 上式为一个位置参数矢量到信道参数矢量的映射，具体映射关系为其中κl＝[τl,θl,cl]T，ζ1＝[px,py]T，ζl＝[sl,x,sl,y]T， 根据信道参数矢量κl＝[τl,θl,cl]T，位置参数矢量其中ζ1＝[px,py]T，ζl＝[sl,x,sl,y]T，构造优化函数： 其中，Jκ是一个加权矩阵，用单位阵替代，满足下式： 位置参数的更新公式如下 ζnew=ζ+h43 其中，h是步进距离，计算公式如下 h＝-H+μI-1g44 其中Jjacobi是雅可比矩阵，计算公式如下 算法中的μ是阻尼系数，取值由增益ρ决定， 其中 如果ρ0，增大μ的大小，反之，则减小μ的大小；迭代结束，至此获得了位置参数的估计值。

如需购买、转让、实施、许可或投资类似专利技术，可联系本专利的申请人或专利权人电子科技大学，其通讯地址为：611731 四川省成都市高新西区西源大道2006号；或者联系龙图腾网官方客服，联系龙图腾网可拨打电话0551-65771310或微信搜索“龙图腾网”。