[go: up one dir, main page]

CN114844553B - Single symbol rate sampling method based on prior filtering for high-speed transmission - Google Patents

Single symbol rate sampling method based on prior filtering for high-speed transmission Download PDF

Info

Publication number
CN114844553B
CN114844553B CN202210315187.9A CN202210315187A CN114844553B CN 114844553 B CN114844553 B CN 114844553B CN 202210315187 A CN202210315187 A CN 202210315187A CN 114844553 B CN114844553 B CN 114844553B
Authority
CN
China
Prior art keywords
sampling
signal
interpolation
analog
digital conversion
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
CN202210315187.9A
Other languages
Chinese (zh)
Other versions
CN114844553A (en
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.)
Beihang University
CETC 54 Research Institute
Original Assignee
Beihang University
CETC 54 Research Institute
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 Beihang University, CETC 54 Research Institute filed Critical Beihang University
Priority to CN202210315187.9A priority Critical patent/CN114844553B/en
Publication of CN114844553A publication Critical patent/CN114844553A/en
Application granted granted Critical
Publication of CN114844553B publication Critical patent/CN114844553B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/14Relay systems
    • H04B7/15Active relay systems
    • H04B7/185Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851Systems using a satellite or space-based relay
    • H04B7/18515Transmission equipment in satellites or space-based relays
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/14Relay systems
    • H04B7/15Active relay systems
    • H04B7/185Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851Systems using a satellite or space-based relay
    • H04B7/18513Transmission in a satellite or space-based system
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/14Relay systems
    • H04B7/15Active relay systems
    • H04B7/185Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851Systems using a satellite or space-based relay
    • H04B7/18517Transmission equipment in earth stations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Analogue/Digital Conversion (AREA)

Abstract

The utility model provides a single code element rate sampling method based on prior filtering for high-speed transmission, which is applied to the technical field of satellite communication and signal processing, wherein, the receiver device is satellite device or ground device, and the receiver device comprises an analog-to-digital conversion chip and a digital processing chip. The analog-to-digital conversion chip processes the received time domain continuous signal to obtain a first sampling signal; the first sampled signal corresponds to a first sampling frequency, the first sampling frequency being a single symbol rate; the analog-to-digital conversion chip sends the first sampling signal to the digital processing chip; the digital processing chip calculates an interpolation value at each interpolation time according to the first sampling signal, and combines the first sampling signal and the interpolation value into a second sampling signal according to the interpolation time; and the second sampling frequency corresponding to the second sampling signal is 2 times of the code element rate. The method can reduce the requirement on the sampling rate of the analog-digital conversion chip and can reduce the power consumption of the analog-digital conversion chip.

Description

应用于高速传输的基于先验滤波的单倍码元速率采样方法Single symbol rate sampling method based on prior filtering for high-speed transmission

技术领域technical field

本公开涉及卫星通信与信号处理技术领域,尤其涉及应用于高速传输的基于先验滤波的单倍码元速率采样方法。The present disclosure relates to the technical field of satellite communication and signal processing, and in particular to a single symbol rate sampling method based on prior filtering applied to high-speed transmission.

背景技术Background technique

奈奎斯特采样定理是通信领域和信号处理领域需要遵守的基本规律之一。奈奎斯特采样定理包括:采样频率fs大于时域信号对应的最高频率fmax的2倍,即fs>=2fmax时采样之后的离散信号可以完整地保留时域信号中的信息,且不会造成频谱混叠。在采样频率不满足奈奎斯特采样定理时,会产生频谱混叠的情况,影响后续的信道补偿、频谱显示、时钟恢复等功能。The Nyquist sampling theorem is one of the basic laws that need to be observed in the field of communication and signal processing. The Nyquist sampling theorem includes: the sampling frequency f s is greater than twice the highest frequency f max corresponding to the time-domain signal, that is, when f s >=2f max , the discrete signal after sampling can completely retain the information in the time-domain signal, and will not cause spectral aliasing. When the sampling frequency does not satisfy the Nyquist sampling theorem, spectrum aliasing will occur, affecting subsequent functions such as channel compensation, spectrum display, and clock recovery.

在需要达到的码元速率较高时,无法获得合适的模数转换芯片是需要解决的技术问题。When the symbol rate to be achieved is high, the inability to obtain a suitable analog-to-digital conversion chip is a technical problem that needs to be solved.

发明内容Contents of the invention

为克服相关技术中存在的问题,本公开提供应用于高速传输的基于先验滤波的单倍码元速率采样方法,应用于接收方设备,应用于接收方设备,所述接收方设备是卫星设备或地面设备,所述接收方设备包括模数转换芯片和数字处理芯片,包括:In order to overcome the problems existing in related technologies, the present disclosure provides a priori filtering-based single symbol rate sampling method applied to high-speed transmission, which is applied to receiver equipment, and is applied to receiver equipment. The receiver equipment is satellite equipment Or ground equipment, the receiver equipment includes an analog-to-digital conversion chip and a digital processing chip, including:

模数转换芯片对接收到的时域连续信号进行处理,获得第一采样信号;其中,所述第一采样信号对应于第一采样频率,所述第一采样频率是单倍码元速率;The analog-to-digital conversion chip processes the received time-domain continuous signal to obtain a first sampling signal; wherein, the first sampling signal corresponds to a first sampling frequency, and the first sampling frequency is a single symbol rate;

所述模数转换芯片将所述第一采样信号发送至数字处理芯片;The analog-to-digital conversion chip sends the first sampling signal to a digital processing chip;

所述数字处理芯片根据所述第一采样信号计算每个内插时刻的内插值,按所述内插时刻将所述第一采样信号和所述内插值组合为第二采样信号;其中,所述第二采样信号对应于第二采样频率,所述第二采样频率为2倍码元速率。The digital processing chip calculates an interpolation value at each interpolation moment according to the first sampling signal, and combines the first sampling signal and the interpolation value into a second sampling signal according to the interpolation moment; wherein, the The second sampling signal corresponds to a second sampling frequency, and the second sampling frequency is twice the symbol rate.

在一实施方式中,所述方法还包括:In one embodiment, the method also includes:

发送方设备使用M倍的第一采样频率对所述时域连续信号进行采样,获得第三采样信号;其中,所述M是大于1的整数;其中,所述第一采样频率是单倍码元速率;The sending device uses M times the first sampling frequency to sample the time-domain continuous signal to obtain a third sampling signal; wherein, the M is an integer greater than 1; wherein, the first sampling frequency is a single code meta rate;

将所述第三采样信号中每个码元对应的连续的M个采样点中相同位置的M-1个采样点设置为零,获得第四采样信号;Set M-1 sampling points at the same position in the consecutive M sampling points corresponding to each symbol in the third sampling signal to zero to obtain a fourth sampling signal;

发送与所述第四采样信号对应的时域连续信号。and sending a time-domain continuous signal corresponding to the fourth sampling signal.

在一实施方式中,所述模数转换芯片对接收到的时域连续信号进行处理,获得第一采样信号,包括:所述模数转换芯片使用所述第一采样频率对所述时域连续信号进行采样,获得第一采样信号;其中,所述第一采样频率是单倍码元速率。In an embodiment, the analog-to-digital conversion chip processes the received time-domain continuous signal to obtain a first sampling signal, including: the analog-to-digital conversion chip uses the first sampling frequency to process the time-domain continuous signal The signal is sampled to obtain a first sampling signal; wherein, the first sampling frequency is a single symbol rate.

在一实施方式中,所述内插时刻是所述第一采样信号中相应的相邻两个码元峰值点中间的过零位置。In an implementation manner, the interpolation time is a zero-crossing position between corresponding peak points of two adjacent symbols in the first sampling signal.

在一实施方式中,所述内插时刻的内插值是所述第一采样信号中所有采样时刻的冲击信号经过带通滤波器后的时域信号在所述内插时刻的叠加值。In an embodiment, the interpolation value at the interpolation time is a superposition value at the interpolation time of the time-domain signals after the impact signals at all sampling time in the first sampling signal pass through a band-pass filter.

在一实施方式中,所述方法还包括:In one embodiment, the method also includes:

使用接收方设备和发送方设备均已知的测试信号,对目标信道的频谱响应进行测试,确定所述目标信道对应的带通滤波器的频谱特性;Using a test signal known to both the receiver device and the sender device, to test the spectral response of the target channel, and determine the spectral characteristics of the bandpass filter corresponding to the target channel;

其中,测试信号包含目标信道的带内所有频谱分量;所述目标信道包括发送方设备的成形滤波器、用于传输无线数据的无线信道、接收方设备的匹配滤波器,所述目标信道是所述线性时不变信道。Wherein, the test signal includes all spectral components in the band of the target channel; the target channel includes a shaping filter of the sending device, a wireless channel for transmitting wireless data, and a matched filter of the receiving device, and the target channel is the the linear time-invariant channel.

在一实施方式中,所述带通滤波器的系统函数为升余弦函数,成形系数为0.5,相应的码元个数大于30。In one embodiment, the system function of the bandpass filter is a raised cosine function, the shaping coefficient is 0.5, and the corresponding number of symbols is greater than 30.

在一实施方式中,所述码元个数为31。In one embodiment, the number of symbols is 31.

本公开的实施例提供的技术方案可以包括以下有益效果:The technical solutions provided by the embodiments of the present disclosure may include the following beneficial effects:

1、降低了信号采样频率,从而降低了对模数转换芯片的采样速率的需求,有利于实现对超高速信号的采样处理。1. The signal sampling frequency is reduced, thereby reducing the demand for the sampling rate of the analog-to-digital conversion chip, which is conducive to the realization of ultra-high-speed signal sampling processing.

2、降低了模数转换芯片的功耗。2. The power consumption of the analog-to-digital conversion chip is reduced.

3、降低了模数转换芯片与数字处理芯片件的信号传输速率,解决了高速数字接口难题。3. The signal transmission rate of the analog-to-digital conversion chip and the digital processing chip is reduced, and the problem of high-speed digital interface is solved.

应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本公开。It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present disclosure.

附图说明Description of drawings

此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the invention and together with the description serve to explain the principles of the invention.

图1是根据一示例性实施例示出的一种应用于高速传输的基于先验滤波的单倍码元速率采样方法的流程图;Fig. 1 is a flowchart of a single symbol rate sampling method based on a priori filtering applied to high-speed transmission according to an exemplary embodiment;

图2是根据一示例性实施例示出的另一种应用于高速传输的基于先验滤波的单倍码元速率采样方法的流程图;FIG. 2 is a flow chart of another priori filtering-based single symbol rate sampling method applied to high-speed transmission according to an exemplary embodiment;

图3是根据一示例性实施例示出的冲击信号经过带限后的时域延拓情况的示意图。Fig. 3 is a schematic diagram showing time-domain extension of an impulse signal after being band-limited according to an exemplary embodiment.

具体实施方式Detailed ways

这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。Reference will now be made in detail to the exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numerals in different drawings refer to the same or similar elements unless otherwise indicated. The implementations described in the following exemplary examples do not represent all implementations consistent with the present invention. Rather, they are merely examples of apparatuses and methods consistent with aspects of the invention as recited in the appended claims.

申请人在研究中发现:Applicants have found in their research that:

1、在码元速率较高时,接收方设备需要模数转换芯片具有较高采样速率,尤其在码元速率高到一定程度时,有可能没有合适的模数转换芯片达到使用需求。1. When the symbol rate is high, the receiving device requires an analog-to-digital conversion chip with a high sampling rate, especially when the symbol rate is high to a certain extent, there may be no suitable analog-to-digital conversion chip to meet the usage requirements.

2、在需要达到的码元速率较高时,接收方设备的模数转换芯片将采样后的数据发送至数字处理芯片后,模数转换芯片与数字处理芯片之间数据传输速率较高,造成巨大的数据传输压力。2. When the symbol rate to be achieved is high, after the analog-to-digital conversion chip of the receiving device sends the sampled data to the digital processing chip, the data transmission rate between the analog-to-digital conversion chip and the digital processing chip is high, resulting in Huge data transmission pressure.

本公开实施例提供了一种应用于高速传输的基于先验滤波的单倍码元速率采样方法。应用于接收方设备,此接收方设备是卫星设备或地面设备。所述接收方设备包括模数转换芯片和数字处理芯片。Embodiments of the present disclosure provide a single symbol rate sampling method based on a priori filtering that is applied to high-speed transmission. Applied to a receiver device, which is a satellite device or a terrestrial device. The receiver device includes an analog-to-digital conversion chip and a digital processing chip.

如图1所示,应用于高速传输的基于先验滤波的单倍码元速率采样方法,包括:As shown in Figure 1, the prior filtering-based single symbol rate sampling method applied to high-speed transmission includes:

步骤S101,模数转换芯片对接收到的时域连续信号进行处理,获得第一采样信号;其中,所述第一采样信号对应于第一采样频率,所述第一采样频率是单倍码元速率;Step S101, the analog-to-digital conversion chip processes the received time-domain continuous signal to obtain a first sampling signal; wherein, the first sampling signal corresponds to a first sampling frequency, and the first sampling frequency is a single symbol rate;

步骤S102,所述模数转换芯片将所述第一采样信号发送至数字处理芯片;Step S102, the analog-to-digital conversion chip sends the first sampling signal to a digital processing chip;

步骤S103,所述数字处理芯片根据所述第一采样信号计算每个内插时刻的内插值,按所述内插时刻将所述第一采样信号和所述内插值组合为第二采样信号;其中,所述第二采样信号对应于第二采样频率,所述第二采样频率为2倍码元速率。Step S103, the digital processing chip calculates an interpolation value at each interpolation time according to the first sampling signal, and combines the first sampling signal and the interpolation value into a second sampling signal according to the interpolation time; Wherein, the second sampling signal corresponds to a second sampling frequency, and the second sampling frequency is twice the symbol rate.

在一种实施方式中,步骤S101中单倍码元速率是所述时域连续信号对应的最高频率。2倍码元速率是所述时域连续信号对应的最高频率的2倍。In one implementation manner, the single symbol rate in step S101 is the highest frequency corresponding to the time-domain continuous signal. The double symbol rate is twice the highest frequency corresponding to the time-domain continuous signal.

本公开实施例可以达到以下技术效果:Embodiments of the present disclosure can achieve the following technical effects:

1、降低了信号采样频率,从而降低了对模数转换芯片的采样速率的需求,有利于实现对超高速信号的采样处理。1. The signal sampling frequency is reduced, thereby reducing the demand for the sampling rate of the analog-to-digital conversion chip, which is conducive to the realization of ultra-high-speed signal sampling processing.

2、降低了模数转换芯片的功耗。2. The power consumption of the analog-to-digital conversion chip is reduced.

3、降低了模数转换芯片与数字处理芯片件的信号传输速率,解决了高速数字接口难题。3. The signal transmission rate of the analog-to-digital conversion chip and the digital processing chip is reduced, and the problem of high-speed digital interface is solved.

本公开实施例提供了一种应用于高速传输的基于先验滤波的单倍码元速率采样方法。应用于接收方设备,所述接收方设备是卫星设备或地面设备。所述接收方设备包括模数转换芯片和数字处理芯片。Embodiments of the present disclosure provide a single symbol rate sampling method based on a priori filtering that is applied to high-speed transmission. Applied to a receiver device, the receiver device is a satellite device or a ground device. The receiver device includes an analog-to-digital conversion chip and a digital processing chip.

此方法包括步骤S100-1所述步骤S101至步骤S103。The method includes step S101 to step S103 described in step S100-1.

其中,步骤S100-1具体包括:Wherein, step S100-1 specifically includes:

发送方设备使用M倍的第一采样频率对所述时域连续信号进行采样,获得第三采样信号;其中,所述M是大于1的整数;其中,所述第一采样频率是单倍码元速率;The sending device uses M times the first sampling frequency to sample the time-domain continuous signal to obtain a third sampling signal; wherein, the M is an integer greater than 1; wherein, the first sampling frequency is a single code meta rate;

为了信号能够在带宽受限信道内传输,需要把传输信号进行白化,将每个码元只保留1个采样点,其它采样点设置为零值,即:将所述第三采样信号中每个码元对应的连续的M个采样点中相同位置的M-1个采样点设置为零,获得第四采样信号;In order that the signal can be transmitted in a bandwidth-limited channel, the transmission signal needs to be whitened, and only one sampling point is reserved for each symbol, and the other sampling points are set to zero values, that is: each of the third sampling signals is M-1 sampling points at the same position in the consecutive M sampling points corresponding to the symbol are set to zero to obtain a fourth sampling signal;

发送与所述第四采样信号对应的时域连续信号。and sending a time-domain continuous signal corresponding to the fourth sampling signal.

在一示例中,M的值为3。步骤S201中,所述第三采样信号中每个码元均对应于3个连续的采样点。其中,将所述第三采样信号中每个码元对应的连续的M个采样点中相同位置的M-1个采样点设置为零是指:In one example, the value of M is 3. In step S201, each symbol in the third sampling signal corresponds to 3 consecutive sampling points. Wherein, setting the M-1 sampling points at the same position in the consecutive M sampling points corresponding to each symbol in the third sampling signal to zero means:

设置所述第三采样信号中每个码元对应的连续的3个采样点中前2个采样点的值为0,Set the value of the first 2 sampling points in the 3 consecutive sampling points corresponding to each symbol in the third sampling signal to 0,

或者,设置所述第三采样信号中每个码元对应的连续的3个采样点中后2个采样点的值为0,Or, set the value of the last 2 sampling points in the 3 consecutive sampling points corresponding to each symbol in the third sampling signal to 0,

或者,设置所述第三采样信号中每个码元对应的连续的3个采样点中第1个采样点和第3个采样点的值为0。Or, set the values of the first sampling point and the third sampling point among the three consecutive sampling points corresponding to each symbol in the third sampling signal to 0.

在一示例中,M的值为4。步骤S201中,所述第三采样信号中每个码元均对应于4个连续的采样点。其中,将所述第三采样信号中每个码元对应的连续的M个采样点中相同位置的M-1个采样点设置为零是指:In one example, the value of M is 4. In step S201, each symbol in the third sampling signal corresponds to 4 consecutive sampling points. Wherein, setting the M-1 sampling points at the same position in the consecutive M sampling points corresponding to each symbol in the third sampling signal to zero means:

设置所述第三采样信号中每个码元对应的连续的4个采样点中前3个采样点的值为0,Set the value of the first 3 sampling points in the 4 consecutive sampling points corresponding to each symbol in the third sampling signal to 0,

或者,设置所述第三采样信号中每个码元对应的连续的4个采样点中后3个采样点的值为0,Or, set the values of the last 3 sampling points in the 4 consecutive sampling points corresponding to each symbol in the third sampling signal to 0,

或者,设置所述第三采样信号中每个码元对应的连续的4个采样点中除第2个采样点之外的其它3个采样点的值为0;Or, setting the values of the other 3 sampling points in the 4 consecutive sampling points corresponding to each symbol in the third sampling signal except the 2nd sampling point are 0;

或者,设置所述第三采样信号中每个码元对应的连续的4个采样点中除第3个采样点之外的其它3个采样点的值为0。Or, set values of the other three sampling points except the third sampling point among the four consecutive sampling points corresponding to each symbol in the third sampling signal to be 0.

本公开实施例提供了应用于高速传输的基于先验滤波的单倍码元速率采样方法。应用于接收方设备,所述接收方设备是卫星设备或地面设备。所述接收方设备包括模数转换芯片和数字处理芯片。Embodiments of the present disclosure provide a priori filtering-based single symbol rate sampling method applied to high-speed transmission. Applied to a receiver device, the receiver device is a satellite device or a ground device. The receiver device includes an analog-to-digital conversion chip and a digital processing chip.

此方法包括所述步骤S101至步骤S103。This method includes the steps S101 to S103.

其中,步骤S101中所述模数转换芯片对接收到的时域连续信号进行处理,获得第一采样信号,包括:Wherein, the analog-to-digital conversion chip described in step S101 processes the received time-domain continuous signal to obtain the first sampling signal, including:

所述模数转换芯片使用所述第一采样频率对所述时域连续信号进行采样,获得第一采样信号;其中,所述第一采样频率是单倍码元速率。The analog-to-digital conversion chip uses the first sampling frequency to sample the time-domain continuous signal to obtain a first sampling signal; wherein, the first sampling frequency is a single symbol rate.

本公开实施例提供了应用于高速传输的基于先验滤波的单倍码元速率采样方法。应用于接收方设备,所述接收方设备是卫星设备或地面设备。所述接收方设备包括模数转换芯片和数字处理芯片。Embodiments of the present disclosure provide a priori filtering-based single symbol rate sampling method applied to high-speed transmission. Applied to a receiver device, the receiver device is a satellite device or a ground device. The receiver device includes an analog-to-digital conversion chip and a digital processing chip.

此方法包括所述步骤S101至步骤S103。This method includes the steps S101 to S103.

并且,所述内插时刻是所述第一采样信号中相应的相邻两个码元峰值点中间的过零位置。In addition, the interpolation time is a zero-crossing position between corresponding peak points of two adjacent symbols in the first sampling signal.

所述内插时刻的内插值是所述第一采样信号中所有采样时刻的冲击信号经过带通滤波器后的时域信号在所述内插时刻的叠加值。The interpolation value at the interpolation moment is a superposition value at the interpolation moment of time-domain signals after the impact signals at all sampling moments in the first sampling signal pass through a band-pass filter.

本公开实施例中,通过将所述叠加值作为内插时刻的内插值,可以抑制插值点频谱的第一旁瓣,降低内插点的噪声值。In the embodiment of the present disclosure, by using the superimposed value as the interpolation value at the interpolation time, the first side lobe of the spectrum of the interpolation point can be suppressed, and the noise value of the interpolation point can be reduced.

本公开实施例提供了应用于高速传输的基于先验滤波的单倍码元速率采样方法。应用于接收方设备,所述接收方设备是卫星设备或地面设备。所述接收方设备包括模数转换芯片和数字处理芯片。Embodiments of the present disclosure provide a priori filtering-based single symbol rate sampling method applied to high-speed transmission. Applied to a receiver device, the receiver device is a satellite device or a ground device. The receiver device includes an analog-to-digital conversion chip and a digital processing chip.

如图2所示,应用于高速传输的基于先验滤波的单倍码元速率采样方法,包括:步骤S100以及所述步骤S101至步骤S103。As shown in FIG. 2 , the priori filtering-based single symbol rate sampling method applied to high-speed transmission includes: step S100 and the steps S101 to S103.

其中,步骤S100包括:使用接收方设备和发送方设备均已知的测试信号,对目标信道的频谱响应进行测试,确定所述目标信道对应的带通滤波器的频谱特性。Wherein, step S100 includes: using a test signal known to both the receiving device and the sending device, to test the spectral response of the target channel, and determine the spectral characteristics of the bandpass filter corresponding to the target channel.

其中,测试信号包含目标信道的带内所有频谱分量;所述目标信道包括发送方设备的成形滤波器、用于传输无线数据的无线信道、接收方设备的匹配滤波器,所述目标信道是所述线性时不变信道。Wherein, the test signal includes all spectral components in the band of the target channel; the target channel includes a shaping filter of the sending device, a wireless channel for transmitting wireless data, and a matched filter of the receiving device, and the target channel is the the linear time-invariant channel.

下面通过一个具体实施例进行说明。 A specific example will be described below.

具体实施例specific embodiment

以卫星遥感数据向地面传输数据的场景为例进行说明。Take the scenario of transmitting data from satellite remote sensing data to the ground as an example.

使用测试信号对目标信道行测试,确定所述目标信道对应的带通滤波器的频谱特 性。确定此带通滤波器的频域系统函数为

Figure 111553DEST_PATH_IMAGE001
,且符合奈奎斯特无码间串扰的条件,其时 域的系统函数为h(t) 。 The test signal is used to test the target channel to determine the spectrum characteristics of the band-pass filter corresponding to the target channel. Determine the frequency domain system function of this bandpass filter as
Figure 111553DEST_PATH_IMAGE001
, and meet the Nyquist condition of no intersymbol interference, the system function in the time domain is h(t) .

在信号正交传输时,对于I路和Q路信号进行并行处理,为了便于表述,这里只描述对于I路信号的处理过程。During the orthogonal transmission of signals, parallel processing is performed on the I-channel and Q-channel signals. For ease of expression, only the processing process for the I-channel signal is described here.

传输前,I路信号的表达式为s(t),如式(1)所示:Before transmission, the expression of channel I signal is s(t) , as shown in formula (1):

Figure 294273DEST_PATH_IMAGE002
(1)
Figure 294273DEST_PATH_IMAGE002
(1)

其中T为码元周期,a(t)为携带信息的随机变量,其幅度值随着基带信息和调制映射关系而变化,在正交相移键控(Quadrature Phase Shift Keying,QPSK)体制方式下,a (t)的取值如式(2)所示:Where T is the symbol period, a(t) is a random variable carrying information, and its amplitude value changes with the baseband information and the modulation mapping relationship. Under the quadrature phase shift keying (Quadrature Phase Shift Keying, QPSK) system , the value of a (t) is shown in formula (2):

Figure 667486DEST_PATH_IMAGE003
(2)
Figure 667486DEST_PATH_IMAGE003
(2)

卫星载荷的FPGA对s(t)进行M倍的采样,被采样后的信号为如式(3)所示的

Figure 125012DEST_PATH_IMAGE004
: The FPGA of the satellite load samples s(t) M times, and the sampled signal is shown in formula (3)
Figure 125012DEST_PATH_IMAGE004
:

Figure 408225DEST_PATH_IMAGE005
(3)
Figure 408225DEST_PATH_IMAGE005
(3)

其中,δ(t)为冲击信号;m表示每个码元的采样点序号,取值从0到M-1; T S 为采样周期。Among them, δ(t) is the impact signal; m represents the sampling point number of each symbol, and the value ranges from 0 to M-1; T S is the sampling period.

为了信号能够在带宽受限信道内传输,把传输信号进行白化处理,即每个码元只保留1个采样点,将其它采样点变为零值。In order for the signal to be transmitted in the bandwidth-limited channel, the transmission signal is whitened, that is, only one sampling point is reserved for each symbol, and other sampling points are changed to zero values.

s(t)在时间区间[nT,(n+1)T]的取值定义为a n ,即如式(4)所示:Define the value of s(t) in the time interval [ nT,(n+1)T ] as a n , as shown in formula (4):

Figure 558584DEST_PATH_IMAGE006
(4)
Figure 558584DEST_PATH_IMAGE006
(4)

其中,a n 为随机变量,在每个码元周期内取值恒定。Among them, a n is a random variable, which takes a constant value in each symbol period.

白化处理后的表达式如式(5)所示:The expression after whitening processing is shown in formula (5):

Figure 294459DEST_PATH_IMAGE007
(5)
Figure 294459DEST_PATH_IMAGE007
(5)

信号

Figure 758938DEST_PATH_IMAGE008
通过信道
Figure 257178DEST_PATH_IMAGE009
后,频域带宽受到切割,信号在时域内进行扩展,如图3 所示。 Signal
Figure 758938DEST_PATH_IMAGE008
through the channel
Figure 257178DEST_PATH_IMAGE009
After that, the frequency domain bandwidth is cut, and the signal is expanded in the time domain, as shown in Figure 3.

冲击信号经过带限系统后的延拓的时域表达式为系统的冲击响应,其中,信道的系统函数和冲击响应为一对傅立叶变换对。The time-domain expression of the continuation of the impulse signal after passing through the band-limited system is the impulse response of the system, where the system function of the channel and the impulse response are a pair of Fourier transform pairs.

Figure 516121DEST_PATH_IMAGE008
的第N个码元的冲击样点值为
Figure 516121DEST_PATH_IMAGE008
The impact sample point value of the Nth symbol is

Figure 739292DEST_PATH_IMAGE010
(6)
Figure 739292DEST_PATH_IMAGE010
(6)

该冲击信号经过信道后的时域信号为y N (t)The time-domain signal of the shock signal after passing through the channel is y N (t) :

Figure 804200DEST_PATH_IMAGE011
(7)
Figure 804200DEST_PATH_IMAGE011
(7)

由图3可知,在两个码元峰值点中间的过零点位置(内插时刻)插入内插值,此内插 值由

Figure 796427DEST_PATH_IMAGE008
中的各采样时刻的冲击信号经过带通滤波器后的时域信号在所述内插时刻的叠 加值; It can be seen from Figure 3 that an interpolation value is inserted at the zero-crossing position (interpolation time) between the peak points of two symbols, and the interpolation value is determined by
Figure 796427DEST_PATH_IMAGE008
The superposition value of the time-domain signal at the interpolation moment of the impact signal at each sampling moment in the bandpass filter;

即内插时刻应插入的数值为:That is, the value to be inserted at the time of interpolation is:

Figure 491850DEST_PATH_IMAGE012
(8)
Figure 491850DEST_PATH_IMAGE012
(8)

接收方设备将通过式(8)计算出的不同内插时刻的内插值与

Figure 61372DEST_PATH_IMAGE008
组合后,使得采 样信号恢复成2倍码元速率,使此采样序列满足了奈奎斯特采样定理。 The receiving device combines the interpolation values at different interpolation moments calculated by formula (8) with
Figure 61372DEST_PATH_IMAGE008
After combination, the sampling signal is recovered to double the symbol rate, and the sampling sequence satisfies the Nyquist sampling theorem.

下述描述用于验证内插算法可以对信号的第一旁瓣具有良好的抑制作用。The following description is used to verify that the interpolation algorithm can have a good suppression effect on the first side lobe of the signal.

对基于先验滤波的内插算法进行仿真,仿真条件为:带限信道的系统函数为升余弦函数,成形系数为0.5,分别取码元个数1、3、9、17、25、31。FFT运算1024点。The interpolation algorithm based on prior filtering is simulated. The simulation conditions are: the system function of the band-limited channel is a raised cosine function, the shaping coefficient is 0.5, and the number of symbols is 1, 3, 9, 17, 25, and 31 respectively. FFT operation 1024 points.

选取相关运算码元个数分别为1、3、9、17、25、31内插后的频谱图,通过对频谱图中不同内插码元个数信号频谱的测量,采样3个码元内插时,对第一旁瓣的抑制只能达到20dB;采用31个码元内插时,对第一旁瓣的抑制可达到60dB;随着内插码元个数的增加,对第一旁瓣的抑制不断增强。Select the interpolated spectrograms with the number of correlation operation symbols being 1, 3, 9, 17, 25, and 31 respectively, and measure the signal spectrum of different numbers of interpolated symbols in the spectrogram, and sample 3 symbols When interpolating, the suppression of the first side lobe can only reach 20dB; when 31 symbols are used for interpolation, the suppression of the first side lobe can reach 60dB; with the increase of the number of interpolation symbols, the suppression of the first side lobe The inhibition of the valve increased continuously.

因此,当先验滤波内插算法选取相关运算码元的个数大于30个,内插算法可以对信号的第一旁瓣具有良好的抑制作用。Therefore, when the number of correlation operation symbols selected by the prior filtering interpolation algorithm is greater than 30, the interpolation algorithm can have a good suppression effect on the first side lobe of the signal.

下述内容用于说明先验滤波的内插算法对信噪比的改善。The following content is used to illustrate the improvement of the signal-to-noise ratio by the interpolation algorithm of prior filtering.

信道内存在高斯白噪声,在内插前,任意取样时刻(t=nT)的高斯白噪声信号ζ n ,为其均值为a,方差为σ 2 。如果线性系统的输入是高斯信号,则系统的输出也是高斯信号。There is Gaussian white noise in the channel. Before interpolation, the Gaussian white noise signal ζ n at any sampling time ( t=nT ) has an average value of a and a variance of σ 2 . If the input to a linear system is a Gaussian signal, the output of the system is also a Gaussian signal.

经过内插后,在内插时刻的高斯白噪声信号为:After interpolation, the Gaussian white noise signal at the interpolation moment is:

Figure 805337DEST_PATH_IMAGE013
Figure 805337DEST_PATH_IMAGE013

经过内插后,在内插点的高斯白噪声信号由其特征值决定。After interpolation, the Gaussian white noise signal at the interpolation point is determined by its eigenvalues.

内插时刻的均值为:The mean of the interpolation moments is:

Figure 714387DEST_PATH_IMAGE014
Figure 714387DEST_PATH_IMAGE014

因为E[ζ n ]=0 ,所以内插时刻的随机信号的均值也为零,即E[ζ NI ]=0。Because E[ζ n ] =0 , the mean value of the random signal at the time of interpolation is also zero, that is, E[ζ NI ] =0.

内插时刻的方差为:The variance of the interpolation time is:

Figure 643029DEST_PATH_IMAGE015
Figure 643029DEST_PATH_IMAGE015

即:Right now:

Figure 575213DEST_PATH_IMAGE016
Figure 575213DEST_PATH_IMAGE016

将上式展开,在展开式中的各元素分为相同元素的平方和不同元素的乘积两种,不同元素的乘积表示在不同时刻噪声信号的乘积,因为高斯白噪声信号各时刻之间相互独立,所以不同元素的乘积均值为零。Expanding the above formula, each element in the expansion formula is divided into two types: the square of the same element and the product of different elements. The product of different elements represents the product of the noise signal at different times, because the Gaussian white noise signal is independent of each other at each time , so the mean of the product of different elements is zero.

从而通过以上分析上式可以表示为:Through the above analysis, the above formula can be expressed as:

Figure 746038DEST_PATH_IMAGE017
Figure 746038DEST_PATH_IMAGE017

随机信号方差代表该信号的交流功率,在均值为零的情况下,内插时刻噪声信号的交流功率即为该信号的总功率。上式表明,内插时刻的信号方差与传输信道的系统函数紧密相关。The variance of the random signal represents the AC power of the signal. When the mean value is zero, the AC power of the noise signal at the interpolation moment is the total power of the signal. The above formula shows that the signal variance at the time of interpolation is closely related to the system function of the transmission channel.

内插时刻的噪声功率可等效为系统函数在不同内插时刻取样的平方和,即The noise power at the interpolation time can be equivalent to the sum of the squares of the system function samples at different interpolation time, that is,

Figure 509594DEST_PATH_IMAGE018
Figure 509594DEST_PATH_IMAGE018

采用MATALAB软件对上式进行计算比较。MATALAB software was used to calculate and compare the above formula.

在不同的成形系数的情况下,对QPSK调制信号进行差值,其信噪比恶化如表1所示。In the case of different shaping coefficients, the difference value is performed on the QPSK modulated signal, and the signal-to-noise ratio deteriorates as shown in Table 1.

表1两倍内插不同成形系数的噪声功率恶化情况Table 1 The deterioration of noise power with different shaping coefficients for double interpolation

Figure 750083DEST_PATH_IMAGE019
Figure 750083DEST_PATH_IMAGE019

Figure 294197DEST_PATH_IMAGE020
Figure 294197DEST_PATH_IMAGE020

在31个码元相关运算时,由于内插值为过零点值,且都是信号的非峰值点内插,因此在内插点的噪声值不但没有恶化,反而比直接采样会优化。信号的成形系数越大,滤波后的拖尾越小,响应的噪声功率也会小,在内插中对噪声的抑制越明显;在成形系数为0.1时,内插点的信噪比比直接采样改善0.39dB;成形系数为0.9时,内插点的信噪比改善2.6dB。During the 31-symbol correlation calculation, since the interpolation value is the zero-crossing point value and is interpolated at the non-peak point of the signal, the noise value at the interpolation point is not only not deteriorated, but is better than direct sampling. The larger the shaping coefficient of the signal, the smaller the smear after filtering, the smaller the noise power of the response, and the more obvious the suppression of noise during interpolation; when the shaping coefficient is 0.1, the signal-to-noise ratio of the interpolation point is directly The sampling improvement is 0.39dB; when the shaping factor is 0.9, the signal-to-noise ratio of the interpolation point is improved by 2.6dB.

所以,先验滤波内插算法对信噪比是有改善的。Therefore, the prior filtering interpolation algorithm can improve the signal-to-noise ratio.

本领域技术人员在考虑说明书及实践这里公开的发明后,将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由下面的权利要求指出。Other embodiments of the invention will be readily apparent to those skilled in the art from consideration of the specification and practice of the invention disclosed herein. This application is intended to cover any modification, use or adaptation of the present invention, these modifications, uses or adaptations follow the general principles of the present invention and include common knowledge or conventional technical means in the technical field not disclosed in this disclosure . The specification and examples are to be considered exemplary only, with a true scope and spirit of the invention being indicated by the following claims.

应当理解的是,本发明并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。It should be understood that the present invention is not limited to the precise constructions which have been described above and shown in the accompanying drawings, and various modifications and changes may be made without departing from the scope thereof. The scope of the invention is limited only by the appended claims.

Claims (7)

1. A single code element rate sampling method based on prior filtering applied to high-speed transmission is applied to receiver equipment, the receiver equipment is satellite equipment or ground equipment, the receiver equipment comprises an analog-to-digital conversion chip and a digital processing chip, and the method is characterized by comprising the following steps:
the analog-to-digital conversion chip processes the received time domain continuous signal to obtain a first sampling signal; wherein the first sampled signal corresponds to a first sampling frequency, the first sampling frequency being a single symbol rate;
the analog-to-digital conversion chip sends the first sampling signal to a digital processing chip;
the digital processing chip calculates an interpolation value at each interpolation time according to the first sampling signal, and combines the first sampling signal and the interpolation value into a second sampling signal according to the interpolation time; wherein the second sampled signal corresponds to a second sampling frequency, the second sampling frequency being 2 times a symbol rate;
the sender equipment samples the time domain continuous signal by using M times of first sampling frequency to obtain a third sampling signal; wherein M is an integer greater than 1; wherein the first sampling frequency is a single symbol rate;
setting M-1 sampling points at the same position in continuous M sampling points corresponding to each code element in the third sampling signal to be zero to obtain a fourth sampling signal;
and transmitting a time domain continuous signal corresponding to the fourth sampling signal.
2. The method of claim 1,
the analog-to-digital conversion chip processes the received time domain continuous signal to obtain a first sampling signal, and the method comprises the following steps:
the analog-to-digital conversion chip samples the time domain continuous signal by using a first sampling frequency to obtain a first sampling signal; wherein the first sampling frequency is a single symbol rate.
3. The method of claim 1,
the interpolation instant is a zero-crossing position in the middle of the peak points of the respective two adjacent symbols in the first sampled signal.
4. The method of claim 1,
the interpolation value at the interpolation time is the superposition value of the time domain signals of the impact signals of all sampling times in the first sampling signals after passing through the band-pass filter at the interpolation time.
5. The method of claim 1,
the method further comprises the following steps:
testing the frequency spectrum response of a target channel by using a test signal known by both receiver equipment and sender equipment, and determining the frequency spectrum characteristic of a band-pass filter corresponding to the target channel;
wherein the test signal contains all spectral components in-band of the target channel; the target channel comprises a shaping filter of a sending device, a wireless channel used for transmitting wireless data and a matched filter of a receiving device, and the target channel is a linear time-invariant channel.
6. The method of claim 5,
the system function of the band-pass filter is a raised cosine function, the forming coefficient is 0.5, and the number of corresponding code elements is more than 30.
7. The method of claim 6,
the number of symbols is 31.
CN202210315187.9A 2022-03-29 2022-03-29 Single symbol rate sampling method based on prior filtering for high-speed transmission Active CN114844553B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210315187.9A CN114844553B (en) 2022-03-29 2022-03-29 Single symbol rate sampling method based on prior filtering for high-speed transmission

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210315187.9A CN114844553B (en) 2022-03-29 2022-03-29 Single symbol rate sampling method based on prior filtering for high-speed transmission

Publications (2)

Publication Number Publication Date
CN114844553A CN114844553A (en) 2022-08-02
CN114844553B true CN114844553B (en) 2023-03-10

Family

ID=82563754

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210315187.9A Active CN114844553B (en) 2022-03-29 2022-03-29 Single symbol rate sampling method based on prior filtering for high-speed transmission

Country Status (1)

Country Link
CN (1) CN114844553B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119449154A (en) * 2025-01-08 2025-02-14 北京邮电大学 Satellite communication signal rapid acquisition method, device, electronic equipment and storage medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101378251A (en) * 2007-08-28 2009-03-04 深圳迈瑞生物医疗电子股份有限公司 Method for improving signal display quality under a condition of low sampling rate
CN102255864A (en) * 2011-08-30 2011-11-23 豪威科技(上海)有限公司 Low-complexity general sampling recovery method and device
US8582676B1 (en) * 2011-08-30 2013-11-12 Omnivision Technologies (Shanghai) Co., Ltd. Low complexity general sampling recovery method and apparatus
CN104901708A (en) * 2015-01-30 2015-09-09 哈尔滨工程大学 Compressive sampling broadband digital receiver and signal processing method thereof
CN106849950A (en) * 2016-12-29 2017-06-13 中国电子科技集团公司第五十研究所 Radiofrequency signal A/D conversion system and method based on multi tate parallel sampling
CN110401510A (en) * 2019-07-23 2019-11-01 中国电子科技集团公司第五十四研究所 A kind of demodulator that can be changed frame by frame suitable for carriers rate and modulation coding scheme

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3023036A1 (en) * 2014-06-27 2016-01-01 Orange RE-SAMPLING BY INTERPOLATION OF AUDIO SIGNAL FOR LOW-LATER CODING / DECODING
CN110460337A (en) * 2019-08-13 2019-11-15 刘涛 A kind of high-speed sampling and reconstructing method
CN112511477A (en) * 2020-11-16 2021-03-16 南京融星智联信息技术有限公司 Hybrid satellite communication modulation identification method and system based on constellation diagram and deep learning

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101378251A (en) * 2007-08-28 2009-03-04 深圳迈瑞生物医疗电子股份有限公司 Method for improving signal display quality under a condition of low sampling rate
CN102255864A (en) * 2011-08-30 2011-11-23 豪威科技(上海)有限公司 Low-complexity general sampling recovery method and device
US8582676B1 (en) * 2011-08-30 2013-11-12 Omnivision Technologies (Shanghai) Co., Ltd. Low complexity general sampling recovery method and apparatus
CN104901708A (en) * 2015-01-30 2015-09-09 哈尔滨工程大学 Compressive sampling broadband digital receiver and signal processing method thereof
CN106849950A (en) * 2016-12-29 2017-06-13 中国电子科技集团公司第五十研究所 Radiofrequency signal A/D conversion system and method based on multi tate parallel sampling
CN110401510A (en) * 2019-07-23 2019-11-01 中国电子科技集团公司第五十四研究所 A kind of demodulator that can be changed frame by frame suitable for carriers rate and modulation coding scheme

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
From Theory to Practice: Sub-Nyquist Sampling of Sparse Wideband Analog Signals;Moshe Mishali,etc.;《IEEE Journal of Selected Topics in Signal Processing》;20100222;全篇 *
一种针对宽带卫星通信信号的欠采样处理方法;文韬等;《舰船电子对抗》;20180425(第02期);全篇 *

Also Published As

Publication number Publication date
CN114844553A (en) 2022-08-02

Similar Documents

Publication Publication Date Title
US7502425B2 (en) Efficient reduced complexity windowed optimal time domain equalizer for discrete multitone-based DSL modems
JP4409395B2 (en) Propagation path estimation method and estimation apparatus
US8817898B2 (en) OFDM receivers
US7961806B2 (en) Power adaptive channel estimation for a multi-path receiving
US7020212B1 (en) Method and system for a multiple dimensional adaptive frequency domain noise canceler for DMT transceivers
US20040252772A1 (en) Filter bank based signal processing
JP4637498B2 (en) Communication apparatus and communication method
US6718019B1 (en) Method and apparatus for wireline characterization
CN111935046B (en) Low-complexity frequency shift keying signal symbol rate estimation method
CN114844553B (en) Single symbol rate sampling method based on prior filtering for high-speed transmission
US6845125B2 (en) xDSL transceiver
KR20180052003A (en) Method and Apparatus for Distortion Compensation of Subcarrier in Orthogonal Frequency Division Multiplexing System
US6993099B2 (en) Communications receiver architectures and algorithms permitting hardware adjustments for optimizing performance
CN105515711B (en) Joint compressed sensing and the rangefinder Interference Suppression System for receiving diversity
US7463678B2 (en) Equalization scheme for DSL receivers in presence of an under-sampled or over-sampled transmit IDFT
WO2012132185A1 (en) Transmission-reception device
CN101001223A (en) A Blind Method for Improving Channel Estimation Accuracy in Wireless Mobile Communications
JP2004201267A (en) Digital filter, its factor calculation apparatus and factor calculation method
KR100747889B1 (en) Channel Estimation Apparatus and Method Using Mutual Conversion in Frequency Domain and Time Domain
Shentu et al. Blind frequency offset estimation for PCC-OFDM with symbols overlapped in the time domain
KR100741349B1 (en) Receiving device for wireless communication
KR20030042103A (en) The aparatus and method for Channel Estimation For Wirelss Telecommnunication System
Qiu et al. Frequency-domain compensation method for CIC filter based on OFDM modulation
Leyonhjelm et al. Designing for low ISI in an OFDM modem
KR20030029235A (en) Feedforward symbol synchronizing apparatus and method thereof

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