CN112995084B - 信号的处理方法及处理装置 - Google Patents

信号的处理方法及处理装置 Download PDF

Info

Publication number
CN112995084B
CN112995084B CN202110169706.0A CN202110169706A CN112995084B CN 112995084 B CN112995084 B CN 112995084B CN 202110169706 A CN202110169706 A CN 202110169706A CN 112995084 B CN112995084 B CN 112995084B
Authority
CN
China
Prior art keywords
domain channel
time domain
window
frequency domain
signal
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
CN202110169706.0A
Other languages
English (en)
Other versions
CN112995084A (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.)
Bicocq Microelectronics Hangzhou Co ltd
Original Assignee
Bicocq Microelectronics Hangzhou Co ltd
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 Bicocq Microelectronics Hangzhou Co ltd filed Critical Bicocq Microelectronics Hangzhou Co ltd
Priority to CN202110169706.0A priority Critical patent/CN112995084B/zh
Publication of CN112995084A publication Critical patent/CN112995084A/zh
Application granted granted Critical
Publication of CN112995084B publication Critical patent/CN112995084B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2695—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters
    • H04B17/345—Interference values
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2668—Details of algorithms
    • H04L27/2681—Details of algorithms characterised by constraints
    • H04L27/2688—Resistance to perturbation, e.g. noise, interference or fading
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2691—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation involving interference determination or cancellation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Quality & Reliability (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Noise Elimination (AREA)

Abstract

本发明公开了一种信号的处理方法及处理装置。其中，该方法包括：获取待处理信号；将待处理信号转换为对应的第一频域信道响应；对第一频域信道响应进行加窗处理，得到第二频域信道响应，其中，加窗处理用于对由第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值的绝对值进行压缩；对第二频域信道响应进行变换获得目标时域信道估计，其中，上述变换为IDFT或者IFFT。本发明解决了由于用户时域信道响应扩展的干扰导致的噪声功率估计不准确的技术问题。

Description

信号的处理方法及处理装置
技术领域
本发明涉及信号处理领域，具体而言，涉及一种信号的处理方法及处理装置。
背景技术
SRS(Sounding reference signal，探测参考信号)用于估计上行信道频域信息，在使用SRS的一个空的循环移位窗进行噪声功率估计时，由于用户的时域信道响应扩展导致噪声功率估计偏大，图1为使用SRS进行噪声功率估计时，时域信道响应输出的绝对值的示意图，如图1所示，时域信道响应输出对应8个时域窗，其中包括7个用户的SRS信号，每一个用户的时域信道响应对应一个时域窗，时域信道响应输出中还有一个空的时域窗10用于计算噪声功率。由图1可看出，噪声所在的时域窗内存在来自其他时域窗内用户的时域信道响应的干扰，导致基于时域窗10内的噪声信号计算得到的噪声功率偏大，进而导致使用SRS所获得的噪声功率估计不准确。
针对上述由于用户时域信道响应扩展的干扰导致的噪声功率估计不准确的问题，目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种信号的处理方法及处理装置，以至少解决由于用户时域信道响应扩展的干扰导致的噪声功率估计不准确的技术问题。
根据本发明实施例的一个方面，提供了一种信号的处理方法，包括：获取待处理信号；将待处理信号转换为对应的第一频域信道响应；对第一频域信道响应进行加窗处理，得到第二频域信道响应，其中，加窗处理用于对由第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值进行压缩；对第二频域信道响应进行变换获得目标时域信道估计，其中，上述变换为IDFT或者IFFT。
进一步地，在对第二频域信道响应进行变换获得目标时域信道估计之后，方法还包括：确定噪声在目标时域信道估计中对应的第一时域窗；提取第一时域窗内的第二时域信道估计；根据第二时域信道估计计算噪声功率。
进一步地，第一时域窗为目标时域信道估计中空的时域窗，根据第二时域信道估计计算噪声功率，包括：获取第一时域窗中的多个噪声信号；根据第二时域信道估计算得到的多个噪声信号的幅值平方的平均值；确定多个噪声信号的幅值平方的平均值为噪声功率。
进一步地，对第一频域信道响应进行加窗处理，得到第二频域信道响应，包括：将第一频域信道响应乘以窗函数，得到第二频域信道响应；其中，窗函数为Hanning窗、Hamming窗、Kaiser窗、Gaussian窗中的任意一个。
进一步地，在对第二频域信道响应进行变换获得目标时域信道估计之后，方法还包括：确定用户在目标时域信道估计中对应的第二时域窗；提取第二时域窗内的第三时域信道估计；对第三时域信道估计补零且进行FFT变换，获得第三频域信道响应；将第三频域信道响应乘以窗函数的倒数，得到目标频域信道估计。
进一步地，待处理信号为多个用户信号的累加信号，提取第二时域窗内的第三时域信道估计，包括：确定第二时域窗的起始位置和窗的长度；根据起始位置和窗的长度，提取第三时域信道估计。
进一步地，待处理信号为SRS信号，将待处理信号转换为对应的第一频域信道响应，包括：将SRS信号乘以对应的SRS序列的共轭，得到第一频域信道响应。
根据本发明实施例的另一方面，还提供了一种信号的处理装置，包括：获取模块，用于获取待处理信号；转换模块，用于将待处理信号转换为对应的第一频域信道响应；加窗模块，用于对第一频域信道响应进行加窗处理，得到第二频域信道响应，其中，加窗处理用于对由第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值进行压缩；变换模块，用于对第二频域信道响应进行变换获得目标时域信道估计，其中，上述变换为IDFT或者IFFT。
根据本发明实施例的另一方面，还提供了一种存储介质，存储介质包括存储的程序，其中，在程序运行时控制存储介质所在设备执行上述任意的信号的处理方法。
根据本发明实施例的另一方面，还提供了一种处理器，处理器用于运行程序，其中，程序运行时执行上述任意的信号的处理方法。
在本发明实施例中，通过获取待处理信号，将待处理信号转换为对应的第一频域信道响应；对第一频域信道响应进行加窗处理，得到第二频域信道响应，加窗处理用于对由第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值进行压缩；对第二频域信道响应进行变换获得目标时域信道估计。由于加窗后获得的目标时域信道估计的信号输出幅值被压缩，使得目标时域信道估计中各用户时域信道响应的时域扩展被有效抑制，进而减小待处理信号对应的时域信道估计中用户所在的时域窗内的时域信道响应对噪声所在时域窗的影响，根据目标时域信道估计中空时域窗的信号参数可以获得更准确的噪声功率估计，解决了由于用户时域信道响应扩展的干扰导致的噪声功率估计不准确的技术问题。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1是根据现有技术中采用SRS进行噪声功率估计时，时域信道响应输出绝对值的示意图；
图2是根据本发明实施例的一种信号的处理方法的流程图；
图3是根据本发明实施例的一种可选的信号的处理方法的流程图；
图4是根据本发明实施例的一种可选的窗函数的示意图；
图5是根据本发明实施例所获得的时域信道响应输出绝对值的示意图；
图6是根据本发明实施例的一种信号的处理装置的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
实施例1
根据本发明实施例，提供了一种信号的处理方法的实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
图2是根据本发明实施例的信号的处理方法的流程图，如图2所示，该方法包括如下步骤：
步骤S201，获取待处理信号。
上述待处理信号为需要进行噪声分析的原始信号。在一种可选的实施例中，待处理信号为4G LTE或者5G NR标准协议中的SRS信号，在一种可选的实施例中，待处理信号的原始信号为SRS信号，作为待处理信号的单用户u发送的SRS信号的表达式为：
Figure BDA0002938613300000041
其中，pn为SRS基本序列，αu为用户u的循环移位窗索引，Nw为总的循环移位窗个数。
在另一种可选的实施例中，通过信号接收机可获取到多个用户的SRS信号的累加信号，多个用户的SRS信号的累加可作为待处理信号。具体的，多个用户经过信道后频域信号的累加结果rn为：
Figure BDA0002938613300000042
其中，
Figure BDA0002938613300000043
为单个用户u的频域信道响应，n＝1，2...N,N为SRS频域子载波的个数，vn为噪声信号，
Figure BDA0002938613300000044
为用户u的频域信道响应。
步骤S202，将待处理信号转换为对应的第一频域信道响应。
第一频域信道响应为待处理信号的原始信号转换获得的对应频域信号，可以用于后续的信号处理过程中的计算和处理。
在一种可选的实施例中，信号接收机获取到多个用户经过信道后频域信号的累加结果rn，则通过解调获得第一频域信道响应
Figure BDA0002938613300000051
为：
Figure BDA0002938613300000052
其中，
Figure BDA0002938613300000053
为SRS基本序列pn的共轭，v′n为噪声信号，由于|pn|＝1，噪声信号vn和v′n的功率相同。
步骤S203，对第一频域信道响应进行加窗处理，得到第二频域信道响应，其中，加窗处理用于对由第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值进行压缩。
上述变换可以为IDFT(Inverse Discrete Fourier Transform，离散傅立叶逆变换)或者IFFT(Inverse Fast Fourier Transform，快速傅立叶逆变换)，第一频域信道响应通过IDFT或者IFFT变换得到对应的第一时域信道估计。
第一时域信道估计输出的信号幅值可以为第一时域信道估计输出的信号功率，即待处理信号对应的时域信道响应的功率，第一时域信道估计输出的信号幅值压缩可理解为减小了第一时域信道响应的功率。例如，图5根据本发明实施例对图1所示的SRS信号第一频域信道响应加窗处理后，获得的第二频域信道响应经IFFT变换后得到时域信道响应的示意图，图5的横坐标为时间，纵坐标表示信号的幅值，该幅值可为信号的功率，图5中的时域信号波形与图1中的信号波形一致，均具有7个用户时域窗和1个噪声时域窗(图1中噪声时域窗为11，图5中噪声时域窗为51)，通过对比图1和图5各时域信道响应的幅值，可以得到图5中各时域信道响应的幅值绝对值均小于图1中相对应的时域信道响应，因此，通过对第一频域信道响应进行加窗处理，可以减小待处理信号对应的各用户时域窗的时域信道响应的功率，进而使用户时域信道响应的时域扩展被有效抑制。
在一种可选的实施例中，对第一频域信道响应进行加窗处理可以为将第一频域信道响应
Figure BDA0002938613300000054
乘以窗函数Wn，获得第二频域信道响应
Figure BDA0002938613300000055
Figure BDA0002938613300000056
窗函数Wn包括但不限于Hanning窗，Wn为Hanning窗时，窗函数Wn的表达式为：
Wn＝(n，α)＝(1-α)-αcos(2*π*n/(N-1))；
其中，0≦n≦N-1，N为SRS频域子载波的个数，α的取值可以为0.46。
步骤S204，对第二频域信道响应进行变换获得目标时域信道估计，其中，上述变换为IDFT或者IFFT。
目标时域信道估计为第一时域信道响应的幅值压缩后对应的时域信道估计，因此，目标时域信道估计输出信号的幅值绝对值小于第一时域信道估计输出信号的幅值绝对值，即目标时域信道估计中各用户时域信道响应的时域扩展被有效抑制，进而减小待处理信号对应的时域信道估计中用户所在的时域窗内的时域信道响应对噪声所在时域窗的影响。
在获得目标时域信道估计后，提取目标时域信道估计中空时域窗(即噪声所在的时域窗，即图5中的时域窗51)中的噪声信号，根据噪声信号的参数可以获得噪声功率估计。
本实施例中，通过获取待处理信号，将待处理信号转换为对应的第一频域信道响应；对第一频域信道响应进行加窗处理，得到第二频域信道响应，加窗处理用于对由第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值进行压缩；对第二频域信道响应进行变换获得目标时域信道估计。由于加窗后获得的目标时域信道估计的信号输出幅值被压缩，使得目标时域信道估计中各用户时域信道响应的时域扩展被有效抑制，进而减小待处理信号对应的时域信道估计中用户所在的时域窗内的时域信道响应对噪声所在时域窗的影响，根据目标时域信道估计中空时域窗的信号参数可以获得更准确的噪声功率估计，解决了由于用户时域信道响应扩展的干扰导致的噪声功率估计不准确的技术问题。
作为一种可选的实施例，在对第二频域信道响应进行变换获得目标时域信道估计之后，上述方法还包括：确定噪声在目标时域信道估计中对应的第一时域窗；提取第一时域窗内的第二时域信道估计；根据第二时域信道估计计算噪声功率。
目标时域信道估计中包含多个时域窗，多个时域窗包括用户所在的时域窗和噪声所在的时域窗，上述第一时域窗为噪声所在的时域窗。
具体的，根据上述步骤S201-S204可以获得多用户的第二频域信道响应
Figure BDA0002938613300000061
第二频域信道响应
Figure BDA0002938613300000062
对进行IDFT变换可以获得目标时域信道估计
Figure BDA0002938613300000063
在确定第一时域窗在目标时域信道估计
Figure BDA0002938613300000064
的位置后，可以确定第二时域信道估计
Figure BDA0002938613300000065
的表达式为：
Figure BDA0002938613300000071
其中，iN为第一时域窗的起始位置；n＝0，1...DN-1；DN为第一时域窗的窗长。
根据第二时域信道估计
Figure BDA0002938613300000072
可以获得噪声对应的时域信道响应进而计算得到噪声功率。
作为一种可选的实施例，第一时域窗为目标时域信道估计中空的时域窗，根据第二时域信道估计计算噪声功率，包括：获取第一时域窗中的多个噪声信号；根据第二时域信道估计算得到的多个噪声信号的幅值平方的平均值；确定多个噪声信号的幅值平方的平均值为噪声功率。
如图5所示，目标时域信道估计包含多个时域窗，其中空的时域窗51为噪声所在的时域窗，即第一时域窗。在确定了第一时域窗的第二时域信道估计
Figure BDA0002938613300000073
后，提取第一时域窗的多个噪声信号xn，n＝1，2..DN，则噪声功率Pn可通过如下公式计算获得:
Figure BDA0002938613300000074
其中，|xn|表示第n个噪声信号的幅值，DN为第一时域窗的窗长。
上述步骤中，通过提取目标时域信道估计中第一时域窗内时域信道估计，利用噪声所在的第一时域窗内的第二时域信道估计计算出噪声功率，由于目标时域信道估计输出信号的幅值小于加窗处理前原始待处理信号对应的第一时域信道估计，降低了用户所在时域窗的时域信道响应对第一时域窗的影响，进而使本实施例所计算获得的噪声功率估计更为准确。
作为一种可选的实施例，对第一频域信道响应进行加窗处理，得到第二频域信道响应，包括：将第一频域信道响应乘以窗函数，得到第二频域信道响应；其中，窗函数为Hanning窗、Hamming窗、Kaiser窗、Gaussian窗中的任意一个。
需要说明的是，窗函数用于压缩第一频域信道响应IDFT变换得到的第一时域信道响应的幅值，窗函数包括但不限于Hanning窗、Hamming窗、Kaiser窗、Gaussian窗，其它具有上述压缩效果的窗函数也在本发明所提出的窗函数的范围之内。
具体的，第二频域信道响应
Figure BDA0002938613300000075
为第一频域信道响应
Figure BDA0002938613300000076
乘以窗函数Wn，即：
Figure BDA0002938613300000081
图4为Hanning窗函数的示意图，Hanning窗函数Wn的表达式为：
Wn＝(n，α)＝(1-α)-αcos(2*π*n/(N-1))；
其中，0≦n≦N-1，α的取值可以为0.46。
通过将第一频域信道响应
Figure BDA0002938613300000082
乘以窗函数Wn，使得第二频域信道响应
Figure BDA0002938613300000083
经IDFT或者IFFT变换后的目标时域信道估计输出信号的幅值绝对值小于第一时域信道估计输出信号的幅值绝对值，即目标时域信道估计中各用户时域信道响应的时域扩展被有效抑制，进而减小待处理信号对应的时域信道估计中用户所在的时域窗内的时域信道响应对噪声所在时域窗的影响。
作为一种可选的实施例，在对第二频域信道响应进行变换获得目标时域信道估计之后，方法还包括：确定用户在目标时域信道估计中对应的第二时域窗；提取第二时域窗内的第三时域信道估计；对第三时域信道估计补零且进行FFT变换，获得第三频域信道响应；将第三频域信道响应乘以窗函数的倒数，得到目标频域信道估计。
由于对待处理信号的原始信号对应的第一频域信道进行加窗处理，使得第一时域信道估计输出信号的幅值被压缩，虽然在噪声功率计算上，降低了来自用户时域窗的时域信道估计的影响，但是加窗处理后的目标时域信道估计并不是真实的用户的时域信道估计，因此还需要将窗函数对目标时域信道估计造成的影响进行补偿，以获得真实的用户的时域信道估计。
目标时域信道估计中包括用户所在的第二时域窗和噪声所在的第一时域窗，将第二时域窗的第三时域信道估计
Figure BDA0002938613300000084
提取出来，进行补零和FFT(Fast Fourier Transform，快速傅立叶变换)变换，得到第三频域信道响应
Figure BDA0002938613300000085
在一种可选的实施例中，由于上述第二频域信道响应
Figure BDA0002938613300000086
为第一频域信道响应
Figure BDA0002938613300000087
乘以窗函数Wn获得，第三频域信道响应
Figure BDA0002938613300000088
为根据加窗处理后的第三时域信道估计
Figure BDA0002938613300000089
变换获得，则目标频域信道估计
Figure BDA00029386133000000810
的表达式为：
Figure BDA00029386133000000811
其中，n＝0,1…N-1，通过对第三频域信道响应
Figure BDA00029386133000000812
乘以窗函数的倒数，补偿了窗函数对用户的频域信道响应的影响。
需要说明的是，目标频域信道估计为单用户的频域信道估计，应与信号接收机获取的待处理信号中的单用户的频域信道估计一致。
作为一种可选的实施例，待处理信号为多个用户信号的累加信号，提取第二时域窗内的第三时域信道估计，包括：确定第二时域窗的起始位置和窗的长度；根据起始位置和窗的长度，提取第三时域信道估计。
具体的，多个用户信号的累加信号的时域信道估计为
Figure BDA0002938613300000091
则单用户u所在的第二时域窗的第三时域信道估计
Figure BDA0002938613300000092
的表达式为：
Figure BDA0002938613300000093
其中，n＝0，1...D-1，iu为用户u的第二时域窗的起始位置，D为用户u的第二时域窗的窗长。
通过上述步骤，可以从目标时域信道估计中提取出用户所在时域窗的第三时域信道估计。
作为一种可选的实施例，待处理信号为SRS信号，将待处理信号转换为对应的第一频域信道响应，包括：将SRS信号乘以其对应的SRS序列的共轭，得到第一频域信道响应。
信号接收机可获取到多个用户的SRS信号的累加信号，经过信道后频域信号的累加结果rn，则通过解调获得第一频域信道响应
Figure BDA0002938613300000094
为：
Figure BDA0002938613300000095
其中，
Figure BDA0002938613300000096
为SRS序列pn的共轭，v′n为噪声信号，由于|pn|＝1，噪声信号v和v′的功率相同。
图3是根据本发明实施例的一种可选的信号的处理方法的流程图，信号接收机接收到多用户发送的SRS累加信号rn：
Figure BDA0002938613300000101
其中，pn为标准的SRS序列，αu为用户u的循环移位窗索引，Nw为总的循环移位窗个数，
Figure BDA0002938613300000102
为单个用户u的频域信道响应，n＝1，2...N,vn为噪声信号。
如图3所示，对上述SRS累加信号rn进行处理的步骤包括：
步骤S301，SRS累加信号rn乘以本地SRS序列的共轭P，
Figure BDA0002938613300000103
获得SRS累加信号的第一频域信道估计
Figure BDA0002938613300000104
Figure BDA0002938613300000105
其中，v′n为噪声信号，由于|pn|＝1，噪声信号vn和v′n的功率相同。
步骤S302，乘以窗函数Wn，获得加窗后的第二频域信道响应
Figure BDA0002938613300000106
步骤S303，将第二频域信道响应
Figure BDA0002938613300000107
进行IDFT变换，获得多用户的目标时域信道
估计
Figure BDA0002938613300000108
在获得目标时域信道估计
Figure BDA0002938613300000109
之后，根据获取噪声功率估计和待处理信号的频域信道估计的不同目的，分别进入步骤S304和S306：
步骤S304，提取噪声所在时域窗内第二时域信道估计
Figure BDA00029386133000001010
如图5所示，目标时域信道估计
Figure BDA00029386133000001011
包含多个时域窗，其中空的时域窗51为噪声所在的时域窗。在确定了第二时域信道估计
Figure BDA00029386133000001012
后，提取第一时域窗的多个噪声信号xn，n＝1，2..DN，DN为噪声所在时域窗的窗长。
步骤S305，利用噪声时域信道估计计算噪声功率，噪声功率Pn可通过如下公式计算获得:
Figure BDA00029386133000001013
其中，|xn|表示第n个噪声信号的幅值。
步骤S306，提取用户的时域窗内第三时域信道估计
Figure BDA00029386133000001014
Figure BDA0002938613300000111
其中，n＝0，1...D-1，iu为用户u的时域窗的起始位置，D为用户u的时域窗的窗长。
步骤S307，对第三时域信道估计
Figure BDA0002938613300000112
进行补零和FFT变换，得到对应的第三频域信道响应
Figure BDA0002938613300000113
步骤S308，乘以窗函数Wn的倒数得到目标频域信道估计
Figure BDA0002938613300000114
通过步骤S304-S305可以获得噪声功率估计，通过步骤S306-S308可以补偿上述步骤S302加窗处理对用户的频域信道响应的影响，以获得真实的用户的频域信道响应。
实施例2
根据本发明实施例，提供了一种信号的处理装置的实施例，图6是根据本发明实施例的信号的处理装置的示意图，如图6所示，该装置包括：
获取模块61，用于获取待处理信号；转换模块62，用于将待处理信号转换为对应的第一频域信道响应；加窗模块63，用于对第一频域信道响应进行加窗处理，得到第二频域信道响应，其中，加窗处理用于对由第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值进行压缩；变换模块64，用于对第二频域信道响应进行变换获得目标时域信道估计，其中，上述变换为IDFT或者IFFT。
作为一种可选的实施例，上述装置还包括：第一确定模块，用于确定噪声在目标时域信道估计中对应的第一时域窗；第一提取模块，用于提取第一时域窗内的第二时域信道估计；噪声计算模块，用于根据第二时域信道估计计算噪声功率。
作为一种可选的实施例，第一时域窗为目标时域信道估计中空的时域窗，噪声计算模块，包括：获取子模块，用于获取第一时域窗中的多个噪声信号；计算子模块，用于根据第二时域信道估计算得到的多个噪声信号的幅值平方的平均值；功率确定子模块，用于确定多个噪声信号的幅值平方的平均值为噪声功率。
作为一种可选的实施例，加窗模块包括：窗函数乘法子模块，用于将第一频域信道响应乘以窗函数，得到第二频域信道响应；其中，窗函数为Hanning窗、Hamming窗、Kaiser窗、Gaussian窗中的任意一个。
作为一种可选的实施例，上述装置还包括：第二确定模块，用于确定用户在目标时域信道估计中对应的第二时域窗；第二提取模块，用于提取第二时域窗内的第三时域信道估计；FFT变换模块，用于对第三时域信道估计补零且进行FFT变换，获得第三频域信道响应；补偿模块，用于将第三频域信道响应乘以窗函数的倒数，得到目标频域信道估计。
作为一种可选的实施例，待处理信号为多个用户信号的累加信号，第二提取模块包括：窗长确定子模块，用于确定第二时域窗的起始位置和窗的长度；提取子模块，用于根据起始位置和窗的长度，提取第三时域信道估计。
作为一种可选的实施例，待处理信号为SRS信号，转换模块包括：序列乘法子模块，用于将SRS信号乘以对应的SRS序列的共轭，得到第一频域信道响应。
本实施例中的装置还包括实现实施例1中其他方法步骤的模块。
实施例3
根据本发明实施例的另一方面，还提供了一种存储介质，存储介质包括存储的程序，其中，在程序运行时控制存储介质所在设备执行实施例1中的信号的处理方法。存储介质被设置为存储用于执行以下步骤的程序代码：获取待处理信号；将待处理信号转换为对应的第一频域信道响应；对第一频域信道响应进行加窗处理，得到第二频域信道响应，其中，加窗处理用于对由第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值进行压缩；对第二频域信道响应进行变换获得目标时域信道估计，其中，上述变换为IDFT或者IFFT。本实施例中，由于加窗后获得的目标时域信道估计的信号输出幅值被压缩，使得目标时域信道估计中各用户时域信道响应的时域扩展被有效抑制，进而减小待处理信号对应的时域信道估计中用户所在的时域窗内的时域信道响应对噪声所在时域窗的影响，根据目标时域信道估计中空时域窗的信号参数可以获得更准确的噪声功率估计，解决了由于用户时域信道响应扩展的干扰导致的噪声功率估计不准确的技术问题。
实施例4
根据本发明实施例的另一方面，还提供了一种处理器，处理器用于运行程序，其中，程序运行时执行实施例1中的信号的处理方法。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
在本发明的上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中，应该理解到，所揭露的技术内容，可通过其它的方式实现。其中，以上所描述的装置实施例仅仅是示意性的，例如所述单元的划分，可以为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，单元或模块的间接耦合或通信连接，可以是电性或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (9)

1.一种信号的处理方法，其特征在于，包括：
获取待处理信号；
将所述待处理信号转换为对应的第一频域信道响应；
对所述第一频域信道响应进行加窗处理，得到第二频域信道响应，其中，所述加窗处理用于对由所述第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值进行压缩；
对所述第二频域信道响应进行变换获得目标时域信道估计，其中，所述变换为IDFT或者IFFT；
其中，在对所述第二频域信道响应进行变换获得目标时域信道估计之后，还包括：确定噪声在所述目标时域信道估计中对应的第一时域窗；提取所述第一时域窗内的第二时域信道估计；根据所述第二时域信道估计计算噪声功率。
2.根据权利要求1所述的方法，其特征在于，所述第一时域窗为所述目标时域信道估计中空的时域窗，根据所述第二时域信道估计计算噪声功率，包括：
获取所述第一时域窗中的多个噪声信号；
根据所述第二时域信道估计算得到的所述多个噪声信号的幅值平方的平均值；
确定所述多个噪声信号的幅值平方的平均值为所述噪声功率。
3.根据权利要求1所述的方法，其特征在于，对所述第一频域信道响应进行加窗处理，得到第二频域信道响应，包括：
将所述第一频域信道响应乘以窗函数，得到所述第二频域信道响应；
其中，所述窗函数为Hanning窗、Hamming窗、Kaiser窗、Gaussian窗中的任意一个。
4.根据权利要求3所述的方法，其特征在于，在对所述第二频域信道响应进行变换获得目标时域信道估计之后，所述方法还包括：
确定用户在所述目标时域信道估计中对应的第二时域窗；
提取所述第二时域窗内的第三时域信道估计；
对所述第三时域信道估计补零且进行FFT变换，获得第三频域信道响应；
将所述第三频域信道响应乘以所述窗函数的倒数，得到目标频域信道估计。
5.根据权利要求4所述的方法，其特征在于，所述待处理信号为多个用户信号的累加信号，提取所述第二时域窗内的第三时域信道估计，包括：
确定所述第二时域窗的起始位置和窗的长度；
根据所述起始位置和窗的长度，提取所述第三时域信道估计。
6.根据权利要求1所述的方法，其特征在于，所述待处理信号为SRS信号，将所述待处理信号转换为对应的第一频域信道响应，包括：
将所述SRS信号乘以对应的SRS序列的共轭，得到所述第一频域信道响应。
7.一种信号的处理装置，其特征在于，包括：
获取模块，用于获取待处理信号；
转换模块，用于将所述待处理信号转换为对应的第一频域信道响应；
加窗模块，用于对所述第一频域信道响应进行加窗处理，得到第二频域信道响应，其中，所述加窗处理用于对由所述第一频域信道响应经变换得到的第一时域信道估计输出的信号幅值进行压缩；
变换模块，用于对所述第二频域信道响应进行变换获得目标时域信道估计，其中，所述变换为IDFT或者IFFT；
其中，在对所述第二频域信道响应进行变换获得目标时域信道估计之后，还包括：确定噪声在所述目标时域信道估计中对应的第一时域窗；提取所述第一时域窗内的第二时域信道估计；根据所述第二时域信道估计计算噪声功率。
8.一种存储介质，其特征在于，所述存储介质包括存储的程序，其中，在所述程序运行时控制所述存储介质所在设备执行权利要求1至6中任意一项所述信号的处理方法。
9.一种处理器，其特征在于，所述处理器用于运行程序，其中，所述程序运行时执行权利要求1至6中任意一项所述信号的处理方法。
CN202110169706.0A 2021-02-07 2021-02-07 信号的处理方法及处理装置 Active CN112995084B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110169706.0A CN112995084B (zh) 2021-02-07 2021-02-07 信号的处理方法及处理装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110169706.0A CN112995084B (zh) 2021-02-07 2021-02-07 信号的处理方法及处理装置

Publications (2)

Publication Number Publication Date
CN112995084A CN112995084A (zh) 2021-06-18
CN112995084B true CN112995084B (zh) 2023-02-24

Family

ID=76349110

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110169706.0A Active CN112995084B (zh) 2021-02-07 2021-02-07 信号的处理方法及处理装置

Country Status (1)

Country Link
CN (1) CN112995084B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101557378A (zh) * 2009-05-18 2009-10-14 普天信息技术研究院有限公司 Ofdm系统中导频发送、信道估计和噪声功率估计方法
CN101764627A (zh) * 2008-12-26 2010-06-30 株式会社Ntt都科摩 确定上行链路的解调导频序列的方法、终端以及上行链路系统
WO2012114769A1 (ja) * 2011-02-25 2012-08-30 京セラ株式会社 チャネル推定装置およびチャネル推定方法
CN103441965A (zh) * 2013-08-30 2013-12-11 电子科技大学 一种Sounding信号的信噪比估计方法
CN110868369A (zh) * 2019-11-26 2020-03-06 武汉烽合智达信息技术有限责任公司 基于5g nr系统的上行信道估计方法及装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101616104B (zh) * 2009-07-27 2011-12-07 北京天碁科技有限公司 正交频分复用系统的信道估计方法和装置
CN102075465A (zh) * 2011-02-14 2011-05-25 电信科学技术研究院 一种噪声功率的估计方法及设备
CN109150783B (zh) * 2017-06-28 2020-12-04 大唐移动通信设备有限公司 一种信道估计方法及装置
CN109639303B (zh) * 2018-12-29 2021-06-18 南京天际易达通信技术有限公司 一种基于加窗处理的干扰检测和抑制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101764627A (zh) * 2008-12-26 2010-06-30 株式会社Ntt都科摩 确定上行链路的解调导频序列的方法、终端以及上行链路系统
CN101557378A (zh) * 2009-05-18 2009-10-14 普天信息技术研究院有限公司 Ofdm系统中导频发送、信道估计和噪声功率估计方法
WO2012114769A1 (ja) * 2011-02-25 2012-08-30 京セラ株式会社 チャネル推定装置およびチャネル推定方法
CN103441965A (zh) * 2013-08-30 2013-12-11 电子科技大学 一种Sounding信号的信噪比估计方法
CN110868369A (zh) * 2019-11-26 2020-03-06 武汉烽合智达信息技术有限责任公司 基于5g nr系统的上行信道估计方法及装置

Also Published As

Publication number Publication date
CN112995084A (zh) 2021-06-18

Similar Documents

Publication Publication Date Title
CN104125190B (zh) 适于低信噪比信道环境的ofdm系统符号定时同步实现方法
CN102347926B (zh) 载波频率捕获方法及装置
KR101294681B1 (ko) 기상 신호 처리장치 및 그 처리방법
CN106506430B (zh) 一种基于压缩感知技术的补偿峰均比非线性失真的新算法
CN106597408B (zh) 基于时频分析和瞬时频率曲线拟合的高阶pps信号参数估计方法
CN109581052A (zh) 一种迭代插值的实复转换频率估计方法
CN102082754B (zh) 一种ofdm信道估计方法及装置
CN109729035A (zh) 一种正交频分复用技术通信系统中的脉冲噪声抑制方法
CN108075807B (zh) 一种小区标识的检测方法及装置
CN112995084B (zh) 信号的处理方法及处理装置
CN102238123A (zh) 测距码检测方法和装置
CN110213820A (zh) 一种定时同步的方法及装置
CN108020741B (zh) 一种双频谐波衰减信号阻尼特征识别方法及装置
CN112583753B (zh) 一种相位补偿方法及电子设备
CN113552543B (zh) 基于set-stiaa的空间微动目标时频分析方法
CN110430156A (zh) 突发ofdm数据传输的帧同步方法及系统
CN115857013A (zh) 一种利用改进welch法计算地震计自噪声的方法
US20170201397A1 (en) Method, apparatus and system for detecting random access signal
CN113314147B (zh) 音频处理模型的训练方法及装置、音频处理方法及装置
CN114062992A (zh) 一种电流互感器一次侧检测过程超声波提取方法及装置
CN115695094A (zh) 信道估计方法、装置及通信设备
Lin Ridges reconstruction based on inverse wavelet transform
CN113178204A (zh) 一种单通道降噪的低功耗方法、装置及存储介质
JP2004004286A (ja) 雑音除去システムおよびプログラム
CN109818885B (zh) 一种信道估计方法、装置、设备及计算机可读存储介质

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