CN104009944B - Apparatus and method for estimating channel response - Google Patents

Abstract

The invention provides a device and a method for estimating channel response. The method of the invention comprises the following steps: a receiving module receives first data and first reference data arriving in a first time interval, second data and second reference data arriving in a second time interval, and third data and third reference data arriving in a third time interval. An estimation module estimates channel responses corresponding to the first data, the third data, the first reference data, the second reference data and the third reference data. A coefficient calculation module performs a one-dimensional nanofiltration filter coefficient calculation procedure according to channel responses corresponding to the first reference data, the second reference data and the third reference data to generate a group of time domain interpolation coefficients. An interpolation module interpolates the channel response corresponding to the second data according to the time-domain interpolation coefficients and by using the channel response corresponding to the first data and the third data respectively.

Description

Apparatus and method for estimating channel response

technical field

The present invention is related to digital signal broadcasting technology, and in particular how to estimate the channel effect of integrated services digital broadcasting-terrestrial (ISDB-T) signals.

Background technique

With the advancement of communication technology, the development of digital TV broadcasting is becoming more and more mature. In addition to transmission via cable lines, digital TV signals can also be transmitted in the form of wireless signals through equipment such as base stations or artificial satellites. Integrated Services Digital Terrestrial Broadcasting (ISDB-T) is one of the standards widely adopted in this field at present. During the transmission of wireless signals, it is inevitable that they will be affected and interfered by the transmission environment. The ISDB-T receiver must evaluate the corresponding channel response, and then eliminate the influence of the channel response on the signal content, so that the received data can be interpreted correctly.

Each data frame (frame) in the ISDB-T signal contains 204 orthogonal frequency division multiplexing (OFDM) symbols (symbols), and each symbol contains multiple subcarriers (subcarrier) bearing content bits. Figure 1 is an example of ISDB-T signal content configuration, the horizontal axis is the frequency, and the vertical axis is the symbol number (the larger the number, the later the transmission time). As shown in FIG. 1 , subcarriers whose frequency index is a multiple of 3 (for example, 0, 3, 6, 9, . . . ) carry a scattered pilot signal (scatter pilot, SP) every four symbols.

Since the original content of the scattered pilot signal is known by the receiving end, it can be used as a basis for the receiving end to determine the channel response. For example, the ISDB-T receiver can first find out the frequency-domain channel response H _(t=0,f=0) of the scattered pilot signal whose time index is 0 and frequency index is 0, and the time index is 4 and the frequency index is 4. The frequency-domain channel response H _{(t=4, f=0)} of the scattered pilot signal is 0. Subsequently, the ISDB- _T receiver can use time-domain interpolation to determine the frequency _- domain channel response H _{( t=1, f=0)} , H _{(t=2, f=0)} , H _{(t=3, f=0)} . Similarly, the ISDB-T receiver can use H _{(t=1, f=3)} and H _{(t=5, f=3)} to interpolate H _{(t=2, f=3)} between them , H _{(t=3, f=3)} , H _{(t=4, f=3)} .

Generally speaking, the ISDB-T receiver determines the time-domain interpolation coefficient according to the ratio of the time interval, for example, let H _{(t=1, f=0)} =H _{(t=0, f=0)} *0.25+H _{(t= 4, f=0)} *0.75, that is, let the time-domain interpolation coefficients corresponding to H _{(t=1, f=0)} be 0.25 and 0.75. And so on, H _(t=2,f=0) =H _(t=0,f=0) *0.5+H _(t=4,f=0) *0.5, and H _{(t=3,f= 0)} =H _(t=0,f=0) *0.75+H _(t=4,f=0) *0.25. In the case where the channel response changes rapidly over time (for example, in a propagation environment with Doppler effect), this method of selecting the interpolation coefficients in the time domain is likely to produce erroneous interpolation results.

Contents of the invention

In order to solve the above problems, the present invention proposes a new channel response estimation device and a channel response estimation method. By using the Wiener filter (Wiener filter) coefficient calculation program to dynamically determine the interpolation coefficient in the time domain, compared with the prior art that determines the interpolation coefficient at a fixed time interval ratio, the estimation device and estimation method according to the present invention can more accurately reflect The variation in channel response over time.

A specific embodiment of the present invention is a channel response estimating device, which includes a receiving module, an estimating module, a coefficient calculating module and an interpolating module. The receiving module is used for receiving a first data and a first reference data arriving in a first time period, a second data and a second reference data arriving in a second time period, and receiving in a first time period A third data and a third reference data arrived in three time segments. The first data, the second data and the third data are transmitted through a first subcarrier. The first reference data, the second reference data and the third reference data are transmitted through a second subcarrier. The estimation module is used for estimating a channel response respectively corresponding to the first data, the third data, the first reference data, the second reference data and the third reference data. The coefficient calculation module is used for performing a Wiener filter coefficient calculation procedure according to the channel responses corresponding to the first reference data, the second reference data and the third reference data, so as to generate a set of time-domain interpolation coefficients. The interpolation module is used for interpolating the channel responses corresponding to the first data and the third data to generate a channel response corresponding to the second data according to the set of time-domain interpolation coefficients.

Another specific embodiment according to the present invention is a channel response estimation method. The method firstly performs a receiving step, receiving a first data and a first reference data arriving in a first time period, a second data and a second reference data arriving in a second time period, and A third data and a third reference data arriving in a third time period. The first data, the second data and the third data are transmitted through a first subcarrier. The first reference data, the second reference data and the third reference data are transmitted through a second subcarrier. Then, the method executes an estimating step, estimating a channel response respectively corresponding to the first data, the third data, the first reference data, the second reference data and the third reference data. Subsequently, the method executes a calculation step, performing a Wiener filter coefficient calculation procedure according to the channel responses corresponding to the first reference data, the second reference data and the third reference data, to generate a set of time-domain interpolation coefficient. Next, the method executes an interpolation step, using the channel responses corresponding to the first data and the third data to interpolate to generate a channel response corresponding to the second data according to the set of time-domain interpolation coefficients.

The advantages and spirit of the present invention can be further understood through the following detailed description of the invention and the accompanying drawings.

In order to have a better understanding of the above-mentioned and other aspects of the present invention, the preferred embodiments are specifically cited below, together with the accompanying drawings, and are described in detail as follows:

Description of drawings

Figure 1 is an example of ISDB-T signal content configuration.

FIG. 2 is a functional block diagram of an apparatus for estimating channel response according to an embodiment of the invention.

FIG. 3 is a flowchart of a channel response estimation method according to an embodiment of the invention.

Symbol Description

200: channel response estimation device 22: receiving module

24: Estimation module 26: Coefficient calculation module

28: Interpolation module S31~S34: Process steps

detailed description

An embodiment of the present invention is a channel response estimating device, the functional block diagram of which is shown in FIG. 2 . The channel response estimating device 200 includes a receiving module 22 , an estimating module 24 , a coefficient calculating module 26 and an interpolating module 28 . In practical applications, the channel response estimating device 200 can be integrated in various ISDB-T wireless receiving systems, or exist independently. Through the following description, those with ordinary knowledge in the technical field of the present invention can understand that the application field of the channel response estimation device 200 is not limited to the ISDB-T system, but can be extended to various time-domain interpolation coefficients that need to determine the channel response Wireless receiving system. The following embodiments mainly take the case where the channel response estimating apparatus 200 is used to cooperate with an ISDB-T receiving end as an example.

First, the receiving module 22 is responsible for receiving and temporarily storing each symbol that arrives at the ISDB-T receiving end in sequence. Its content is shown in FIG. P's transmission and multiplexing configuration control (transmission and multiplexing configuration control, TMCC) signal. These symbols are continuous in time, and the smaller the symbol number, the earlier the arrival time. The TMCC signal is used to transmit information such as modulation type, coding rate, scrambling length, multiplexing scheme, etc., for the reference of the receiving end. In practice, as the ISDB-T transmission mode changes, the number of TMCC subcarriers in each symbol may also be different. The following description assumes that P is a specific value corresponding to a certain TMCC subcarrier.

In this embodiment, the channel response estimating apparatus 200 uses the TMCC signal and the scattered pilot signal as the basis for estimating the unknown channel response. More specifically, the channel response estimating apparatus 200 uses the TMCC signal as a reference signal to determine time-domain interpolation coefficients, and then interpolates the channel response of the scattered pilot signal to generate the channel response of the data signal. Taking the case of determining the frequency-domain channel response H _{(t=1, f=0)} of a data signal with a time index of 1 and a frequency index of 0 as an example, the channel response estimating apparatus 200 first uses the time indexes of 0, 1, and 4 respectively , and the three TMCC signals whose frequency index is P are used as reference signals to determine a set of time-domain interpolation coefficients. Subsequently, the channel response estimating apparatus 200 generates H (t _{= 1, f = 0)} . The operation of the channel response estimation device 200 will be described in detail below.

Estimation module 24 is responsible for estimating the frequency domain channel response of each scattered pilot signal and TMCC signal. As mentioned earlier, the original content of the scattered pilot signal is known at the receiver, so the estimation module 24 can directly compare the known content with the received content to generate its frequency domain channel response. Relatively, the original content of the TMCC signals is unknown to the receiving end, and the estimation module 24 can first demodulate and decode each TMCC signal, and then evaluate its frequency-domain channel response. Since the ISDB-T transmitter uses differential binary phase shift keying (DBPSK) to generate TMCC signals, that is to say, there are only two types of TMCC signals, +1 and -1, and the decoding result of TMCC signals is correct. Therefore quite high. Correspondingly, the frequency-domain channel response of the TMCC signal generated by the estimation module 24 also has quite high reliability. It should be noted that the method of determining the frequency-domain channel response of the known signal is known to those skilled in the art of the present invention, and will not be repeated here.

Next, according to the frequency domain channel responses H _TMCC0 , H _TMCC1 , H _TMCC4 generated by the estimation module 24 for the three TMCC signals whose time index is 0, 1, and 4 respectively, and the frequency index is P, the coefficient calculation module 26 performs a one-dimensional A Wiener filter coefficient calculation program to generate a set of time-domain interpolation coefficients. One of the characteristics of the Wiener filter coefficient calculation program is that the amount of error in the estimation result can be minimized. In this embodiment, the coefficient calculation module 26 generates the first coefficient W ₁ and the second coefficient W ₂ included in the set of time-domain interpolation coefficients according to the following equation:

(Formula 1)

Equation 1 is one of the forms of the Wiener-Hopf equation.

Subsequently, the interpolation module 28 utilizes H _{(t=0, f=0)} and H _{(t=4, f=0)} generated by the estimation module 24 to generate H (t=1, f=0) according to the set of time-domain interpolation coefficients _{. f=0)} , its operation program can be expressed as:

(Formula 2)

In one embodiment, the coefficient calculation module 26 may further include a smoothing unit (not shown), for performing a smoothing process on the following matrix before generating the first coefficient W ₁ and the second coefficient W ₂ :

(Formula 3)

The matrix is made into a Toeplitz matrix, thereby simplifying the computational complexity of Equation 1. For example, the smoothing procedure may be, but not limited to, calculating the average value of the upper left element and the lower right element in the matrix, and replacing these two elements with the average value.

If the matrix in Equation 3 is smoothed into a Teeplitz matrix, the coefficient calculation module 26 can correspondingly use a Levinson recursion algorithm to generate the first coefficient W ₁ and the second coefficient W ₂ . It should be noted that the method of using the Levinson recursive algorithm to simplify the solution of the Wiener-Hopper equation is a known technology, and will not be repeated here.

As shown in Figure 1, a data signal with a time index of 0 and a frequency index of 0 and a TMCC signal with a time index of 0 and a frequency index of P are included in the same symbol, which can be regarded as being sent to ISDB-T in the same time segment Receiving end. Similarly, a data signal with a time index of 1 and a frequency index of 0 and a TMCC signal with a time index of 1 and a frequency index of P are sent to the ISDB-T receiving end in the same time period, while the time index is 4 and the frequency index is The data signal of 0 and the TMCC signal of time index 4 and frequency index P are sent to the ISDB-T receiving end in the same time segment. Since the frequency response characteristics in the same time segment are possibly similar, the channel response estimating device 200 uses the time correlations of H _TMCC1 with respect to H _TMCC0 and H _TMCC4 to estimate the relationship between H _{(t=1, f=0)} and H _{( t=0, f=0)} and the temporal correlation of H _{(t=4, f=0)} . By analogy, according to the frequency-domain channel responses produced by the estimation module 24 for the three TMCC signals whose time index is 0, 2, and 4 respectively, and whose frequency index is P, the coefficient calculation module 26 can perform one-dimensional filter coefficient calculation A program for generating another set of time-domain interpolation coefficients, suitable for generating a frequency-domain channel response of a data signal with a time index of 2 and a frequency index of 0. Different from the prior art in which the interpolation coefficients are determined at a fixed time interval ratio, the method of the channel response estimation apparatus 200 can more accurately reflect the variation of the channel response over time on average.

It should be noted that although in the above embodiments, the channel responses used for interpolation to generate the unknown channel responses correspond to known signals (scattered pilot signals), the concept of the present invention is not limited thereto. For example, if the channel response estimation 200 has determined the channel response of a data signal with a time index of 2 and a frequency index of 3, and a channel response of a data signal with a time index of 4 and a frequency index of 3, it may be further determined according to TMCC The signal determines a set of interpolation coefficients for generating the channel response of the data signal with a time index of 3 and a frequency index of 3.

In addition, the reference signal used to determine the interpolation coefficient is not limited to the TMCC signal. For example, in the digital video broadcasting-terrestrial (DVB-T) system, the continuous pilot signal (CP) contained in each symbol can also be used as a reference signal for the receiving end to determine the interpolation coefficient .

Another specific embodiment according to the present invention is a channel response estimation method, the flowchart of which is shown in FIG. 3 . The method first executes step S31, receiving a first data (for example, the scattered pilot signal whose time index is 0 and frequency index is 0 in Fig. 1 ) and a first reference data (such as Fig. 1 TMCC signal whose time index is 0 and frequency index is P), a second data arriving in a second time segment (such as a data signal whose time index is 1 and frequency index is 0 in FIG. 1 ) and a second Reference data (such as the TMCC signal whose time index is 1 and whose frequency index is P in FIG. 1 ), and a third data arriving in a third time zone (such as a time index of 4 and a frequency index of 0 in FIG. 1 scattered pilot signal) and a third reference data (for example, the TMCC signal whose time index is 4 and frequency index is P in FIG. 1 ). The first data, the second data and the third data are transmitted through a first subcarrier. The first reference data, the second reference data and the third reference data are transmitted through a second subcarrier. Next, the method executes step S32 , estimating a channel response corresponding to the first data, the third data, the first reference data, the second reference data and the third reference data respectively. Subsequently, the method executes step S33, performing a Wiener filter coefficient calculation procedure according to the channel responses corresponding to the first reference data, the second reference data, and the third reference data to generate a set of time-domain interpolation coefficients . Next, the method executes step S34 , using the channel response interpolation corresponding to the first data and the third data to generate a channel response corresponding to the second data according to the set of time-domain interpolation coefficients.

The various circuit operation changes previously described when introducing the channel response estimating apparatus 200 (for example, the way of simplifying the Wiener filter coefficient calculation procedure) can also be applied to the channel response estimation method shown in FIG. 3 , and the details will not be repeated here .

As mentioned above, the present invention proposes a new channel response estimation device and a channel response estimation method. By using the Wiener filter coefficient calculation program to dynamically determine the interpolation coefficients in the time domain, compared with the prior art in which the interpolation coefficients are determined at a fixed time interval ratio, the estimation device and estimation method according to the present invention can more accurately reflect the channel response as the changes over time.

Through the detailed description of the specific embodiments above, it is hoped that the characteristics and spirit of the present invention can be described more clearly, rather than limiting the scope of the present invention by the specific embodiments disclosed above. On the contrary, the intention is to cover various changes and equivalent arrangements within the scope of the claimed patent scope of the present invention.

Claims (10)

1. A channel response estimation device, comprising: A receiving module, used to receive a first data and a first reference data arriving in a first time interval, a second data and a second reference data arriving in a second time interval, and in a A third data and a third reference data arriving in the third time segment, wherein the first data, the second data and the third data are transmitted through a first subcarrier, the first reference data, the first reference data, The second reference data and the third reference data are transmitted through a second subcarrier; an estimating module, configured to estimate a channel response respectively corresponding to the first data, the third data, the first reference data, the second reference data and the third reference data; A coefficient calculation module, used to perform a Wiener filter coefficient calculation procedure according to the channel responses corresponding to the first reference data, the second reference data and the third reference data, so as to generate a set of time-domain interpolation coefficients; as well as An interpolation module, used for interpolating the channel responses corresponding to the first data and the third data to generate a channel response corresponding to the second data according to the set of time-domain interpolation coefficients; Wherein, the first data, the second data, the third data, the first reference data, the second reference data and the third reference data are sent by an integrated service digital terrestrial broadcasting transmitter, and the first data Each of the third data and the third data is a scattered pilot signal, and the first reference data, the second reference data and the third reference data are included in a transmission and multiplexing configuration control signal. 2. The channel response estimation device according to claim 1, wherein the set of time-domain interpolation coefficients includes a first coefficient W ₁ and a second coefficient W ₂ , and the coefficient calculation module generates the first coefficient according to the following equation Coefficient W ₁ and the second coefficient W ₂ : $\left[\begin{array}{c}{\mathrm{<span\; class="notranslate">\; W</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}\\ {\mathrm{<span\; class="notranslate">\; W</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}\end{array}\right]<span\; class="notranslate">\; =</span>{\left[\begin{array}{cc}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}^{<span\; class="notranslate">\; *</span>}& {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}^{<span\; class="notranslate">\; *</span>}\\ {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}^{<span\; class="notranslate">\; *</span>}& {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}^{<span\; class="notranslate">\; *</span>}\end{array}\right]}^{<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; 1</span>}}\left[\begin{array}{c}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}^{<span\; class="notranslate">\; *</span>}\\ {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}^{<span\; class="notranslate">\; *</span>}\end{array}\right]<span\; class="notranslate">\; ,</span>$ T ₁ represents the channel response corresponding to the first reference data, T ₂ represents the channel response corresponding to the second reference data, and T ₃ represents the channel response corresponding to the third reference data. 3. The channel response estimating device according to claim 2, wherein the interpolation module generates the channel response _H2 corresponding to the second data according to the following equation: ${\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; =</span><span\; class="notranslate">\; \&lsqb;</span>\begin{array}{cc}{\mathrm{<span\; class="notranslate">\; W</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}& {\mathrm{<span\; class="notranslate">\; W</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}\end{array}<span\; class="notranslate">\; \&rsqb;</span><span\; class="notranslate">\; \&CenterDot;</span>\left[\begin{array}{c}{\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}\\ {\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}\end{array}\right]<span\; class="notranslate">\; ,</span>$ Wherein H1 represents the channel response corresponding to the _first data, _H2 represents the channel response corresponding to the second data, and _H3 represents the channel response corresponding to the third data. 4. The channel response estimation device according to claim 2, wherein the coefficient calculation module comprises a smoothing unit for smoothing the following matrix before generating the first coefficient and the second coefficient program: ${\left[\begin{array}{cc}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}^{<span\; class="notranslate">\; *</span>}& {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}^{<span\; class="notranslate">\; *</span>}\\ {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}^{<span\; class="notranslate">\; *</span>}& {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}^{<span\; class="notranslate">\; *</span>}\end{array}\right]}^{<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; ,</span>$ to make the matrix a Teeplitz matrix. 5. The channel response estimating device as claimed in claim 4, wherein the coefficient calculation module utilizes a Levinson recursive algorithm to generate the first coefficient and the second coefficient. 6. A channel response estimation method, comprising: (a) receiving a first data and a first reference data arriving at a first time period, a second data and a second reference data arriving at a second time period, and at a third time A third data and a third reference data arriving in the segment, wherein the first data, the second data and the third data are transmitted through a first subcarrier, the first reference data, the second reference data and the third reference data are transmitted through a second subcarrier; (b) estimating a channel response respectively corresponding to the first data, the third data, the first reference data, the second reference data and the third reference data; (c) performing a Wiener filter coefficient calculation procedure according to the channel responses corresponding to the first reference data, the second reference data and the third reference data to generate a set of time-domain interpolation coefficients; and (d) using the channel response interpolation corresponding to the first data and the third data according to the set of time-domain interpolation coefficients to generate a channel response corresponding to the second data; Wherein, the first data, the second data, the third data, the first reference data, the second reference data and the third reference data are sent by an integrated service digital terrestrial broadcasting transmitter, and the first data Each of the third data and the third data is a scattered pilot signal, and the first reference data, the second reference data and the third reference data are included in a transmission and multiplexing configuration control signal. 7. The channel response estimation method according to claim 6, wherein the set of time-domain interpolation coefficients includes a first coefficient W ₁ and a second coefficient W ₂ , and step (c) includes generating the first coefficient W 2 according to the following equation A coefficient W ₁ and the second coefficient W ₂ : $\left[\begin{array}{c}{\mathrm{<span\; class="notranslate">\; W</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}\\ {\mathrm{<span\; class="notranslate">\; W</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}\end{array}\right]<span\; class="notranslate">\; =</span>{\left[\begin{array}{cc}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}^{<span\; class="notranslate">\; *</span>}& {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}^{<span\; class="notranslate">\; *</span>}\\ {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}^{<span\; class="notranslate">\; *</span>}& {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}^{<span\; class="notranslate">\; *</span>}\end{array}\right]}^{<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; 1</span>}}\left[\begin{array}{c}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}^{<span\; class="notranslate">\; *</span>}\\ {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}^{<span\; class="notranslate">\; *</span>}\end{array}\right]<span\; class="notranslate">\; ,</span>$ T ₁ represents the channel response corresponding to the first reference data, T ₂ represents the channel response corresponding to the second reference data, and T ₃ represents the channel response corresponding to the third reference data. 8. The channel response estimation method according to claim 7, wherein step (d) comprises generating the channel response H ₂ corresponding to the second data according to the following equation: ${\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; =</span><span\; class="notranslate">\; \&lsqb;</span>\begin{array}{cc}{\mathrm{<span\; class="notranslate">\; W</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}& {\mathrm{<span\; class="notranslate">\; W</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}\end{array}<span\; class="notranslate">\; \&rsqb;</span><span\; class="notranslate">\; \&Center\; Dot;</span>\left[\begin{array}{c}{\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}\\ {\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}\end{array}\right]<span\; class="notranslate">\; ,</span>$ Wherein H1 represents the channel response corresponding to the _first data, _H2 represents the channel response corresponding to the second data, and _H3 represents the channel response corresponding to the third data. 9. The channel response estimation method according to claim 7, wherein step (c) includes performing a smoothing procedure on the following matrix before generating the first coefficient and the second coefficient: ${\left[\begin{array}{cc}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}^{<span\; class="notranslate">\; *</span>}& {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}^{<span\; class="notranslate">\; *</span>}\\ {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}^{<span\; class="notranslate">\; *</span>}& {\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; 3</span>}}^{<span\; class="notranslate">\; *</span>}\end{array}\right]}^{<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; ,</span>$ to make the matrix a Teeplitz matrix. 10. The channel response estimation method as claimed in claim 9, wherein step (c) comprises using a Levinson recursive algorithm to generate the first coefficient and the second coefficient.