Method and apparatus for tuning finite impulse response filter in in-band full duplex transceiver

Abstract

A method and an apparatus for tuning an FIR filter in an in-band full duplex transceiver. The method for tuning an FIR filter may include: setting attenuation of the FIR filter to be a first value and then estimating input information of the FIR filter; estimating a delta response using the estimated input information of the FIR filter; and updating the attenuation of the FIR filter to a second value using the estimated delta response.

Claims

1. A method for tuning a finite impulse response (FIR) filter in an in-band full duplex transceiver, the FIR filter receiving a transmission signal to cancel a self-transmitted interference signal, the method comprising: setting, by the in-band full duplex transceiver, attenuation of the FIR filter to be a predetermined first value and then estimating input information of the FIR filter; estimating, by the in-band full duplex transceiver, a delta response using the input information of the FIR filter; updating, by the in-band full duplex transceiver, the attenuation to a second value using the delta response; and applying, by the in-band full duplex transceiver, the second value to the FIR filter, wherein the estimating of the input information of the FIR filter includes: setting, by the in-band full duplex transceiver, the attenuation to be the first value and then generating and transmitting a first training signal, in a first period; measuring, by the in-band full duplex transceiver, a first signal that is a baseband equivalent signal for a signal obtained by subtracting an output signal of the FIR filter from the self-transmitted interference signal, in the first period; generating and transmitting, by the in-band full duplex transceiver, a second training signal in a second period; measuring, by the in-band full duplex transceiver, a second signal that is a baseband equivalent signal for the self-transmitted interference signal, in the second periods; and estimating, by the in-band full duplex transceiver, the input information of the FIR filter using the first signal and the second signal.

2. The method of claim 1, wherein the delta response corresponds to a value obtained by subtracting a baseband equivalent time impulse response for the FIR filter set as the first value from a baseband equivalent time impulse response for a radio channel.

3. The method of claim 1, wherein the updating includes: calculating, by the in-band full duplex transceiver, a third value that is the attenuation of the FIR filter to be updated by using the delta response; and calculating, by the in-band full duplex transceiver, the second value by adding the third value to the first value.

4. The method of claim 1, wherein the estimating of the delta response includes estimating, by the in-band full duplex transceiver, the delta response using the input information of the FIR filter and the first signal.

5. The method of claim 1, wherein the estimating of the input information of the FIR filter includes: subtracting, by the in-band full duplex transceiver, the first signal from the second signal to generate a third signal that is a baseband equivalent signal for the output signal of the FIR filter; and calculating, by the in-band full duplex transceiver, the input information of the FIR filter using the third signal.

6. The method of claim 1, wherein the measuring of the second signal includes suppressing, by the in-band full duplex transceiver, the output signal of the FIR filter from being generated and then measuring the second signal.

7. The method of claim 1, wherein the first training signal and the second training signal are the same signal.

8. The method of claim 1, wherein the input information of the FIR filter includes a linear component and a nonlinear component.

9. The method of claim 1, wherein the first period and the second period are included in a field for transmitting the first training signal and the second training signal in a time domain frame.

10. An in-band full duplex transceiver comprising: a transmitter generating and transmitting a transmitted signal; a receiver converting a received signal into a baseband frequency signal and demodulating the converted signal; a finite impulse response (FIR) filter receiving the transmitted signal to cancel a self-transmitted interference signal; and an estimator estimating input information of the FIR filter by performing fast Fourier transform (FFT) on output information of an analog digital converter included in the receiver and setting attenuation of the FIR filter using the input information of the FIR filter, wherein the estimator sets the attenuation to be a first value and then estimates the input information of the FIR filter, estimates a delta frequency response using the input information of the FIR filter, updates the attenuation to a second value using the delta frequency response, and applies the second value to the FIR filter, wherein the transmitter generates and transmits a first training signal in a first period and generates and transmits a second training signal in a second period, wherein in the first period, the estimator sets the attenuation to be the first value and then performs the FFT on the output information of the analog digital converter to acquire a first signal, and wherein in the second period, the estimator performs the FFT on baseband equivalent information for the self-transmitted interference signal to acquire a second signal.

11. The in-band full duplex transceiver of claim 10, wherein the delta frequency response corresponds to a value obtained by subtracting a frequency response for the FIR filter set as the first value from a frequency response for a radio channel.

12. The in-band full duplex transceiver of claim 10, wherein the estimator subtracts the first signal from the second signal to calculate a third signal and uses the third signal to calculate the input information of the FIR filter.

13. The in-band full duplex transceiver of claim 10, wherein the estimator sets the attenuation of the FIR filter in the second period to suppress an output signal of the FIR filter.

14. The in-band full duplex transceiver of claim 10, further comprising: an antenna transmitting the transmitted signal to receive the received signal; and a distributor transmitting the transmitted signal to the antenna and transmitting the received signal to the receiver.

15. The in-band full duplex transceiver of claim 10, wherein the first training signal and the second training signal are the same signal.

16. A method for tuning a finite impulse response (FIR) filter in an in-band full duplex transceiver, the in-band full duplex transceiver including a transmitter generating and transmitting a transmitted signal, a receiver converting a received signal into a baseband frequency signal, and the FIR filter receiving the transmitted signal to cancel a self-transmitted interference signal, the method comprising: setting, by the in-band full duplex transceiver, attenuation of the FIR filter to be a first value and then generating and transmitting a first training signal, in a first period; performing, by the in-band full duplex transceiver, FFT on output information of an analog digital converter included in the receiver to acquire a first signal, in the first period; generating and transmitting, by the in-band full duplex transceiver, a second training signal in a state in which an output signal of the FIR filter is suppressed, in a second period; performing, by the in-band full duplex transceiver, the FFT on the output information of the analog digital converter to acquire a second signal, in the second period; estimating, by the in-band full duplex transceiver, input information of the FIR filter using the first signal and the second signal; estimating, by the in-band full duplex transceiver, a delta frequency response using the input information of the FIR filter, updating, by the in-band full duplex transceiver, the attenuation to a second value using the delta frequency response; and applying, by the in-band full duplex transceiver, the second value to the FIR filter.

17. The method of claim 16, wherein the delta frequency response corresponds to a value obtained by subtracting a frequency response for the FIR filter set as the first value from a frequency response for a radio channel.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a diagram illustrating an in-band full duplex transceiver according to an exemplary embodiment of the present invention.

(2) FIG. 2 is a diagram illustrating an in-band full duplex transceiver according to another exemplary embodiment of the present invention.

(3) FIG. 3 is a diagram illustrating a time domain frame structure according an exemplary embodiment of the present invention.

(4) FIG. 4 is a flowchart illustrating a method for processing input information of an FIR filter in a time domain to estimate the input information according to an exemplary embodiment of the present invention.

(5) FIG. 5 is a flowchart illustrating a method for tuning, by an estimator, an FIR filter according to an exemplary embodiment of the present invention.

(6) FIG. 6 is a flowchart illustrating a method for processing input information of an FIR filter in a frequency domain to estimate the input information according to another exemplary embodiment of the present invention.

(7) FIG. 7 is a flowchart illustrating a method for tuning, by an estimator, an FIR filter according to another exemplary embodiment of the present invention.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(8) In the following detailed description, only certain exemplary embodiments of the present invention have been shown and described, simply by way of illustration. As those skilled in the art would realize, the described embodiments may be modified in various different ways, all without departing from the spirit or scope of the present invention. Accordingly, the drawings and description are to be regarded as illustrative in nature and not restrictive. Like reference numerals designate like elements throughout the specification.

(9) Throughout the specification, a terminal may be called a mobile terminal (MT), a mobile station (MS), an advanced mobile station (AMS), a high reliability mobile station (HR-MS), a subscriber station (SS), a portable subscriber station (PSS), an access terminal (AT), user equipment (UE), and the like and may include functions of all or some of the MT, the MS, the AMS, the HR-MS, the SS, the PSS, the AT, the UE, and the like.

(10) Further, a base station (BS) may be called an advanced base station (ABS), a high reliability base station (HR-BS), a nodeB, an evolved node B (eNodeB), an access point (AP), a radio access station (RAS), a base transceiver station (BTS), a mobile multihop relay (MMR)-BS, a relay station (RS) serving as a base station, a high reliability relay station (HR-RS) serving as a base station, and the like and may also include functions of all or some of the ABS, the HR-BS, the nodeB, the eNodeB, the AP, the RAS, the BTS, the MMR-BS, the RS, the HR-RS, and the like.

(11) Throughout the specification, a transceiver may refer to a terminal, a mobile terminal (MT), a mobile station (MS), an advanced mobile station (AMS), a high reliability mobile station (HR-MS), a subscriber station (SS), a portable subscriber station (PSS), an access terminal (AT), user equipment (UE), and the like and may include functions of all or some of the terminal, the MT, the MS, the AMS, the HR-MS, the SS, the PSS, the AT, the UE, and the like.

(12) Further, a transceiver may refer to a base station (BS), an advanced base station (ABS), a high reliability base station (HR-BS), a nodeB, an evolved node B (eNodeB), an access point (AP), a radio access station (RAS), a base transceiver station (BTS), a mobile multihop relay (MMR)-BS, a relay station (RS) serving as a base station, a high reliability relay station (HR-RS) serving as a base station, and the like and may also include functions of all or some of the Bs, the ABS, the HR-BS, the nodeB, the eNodeB, the AP, the RAS, the BTS, the MMR-BS, the RS, the HR-RS, and the like.

(13) FIG. 1 is a diagram illustrating an in-band full duplex transceiver 1000 according to an exemplary embodiment of the present invention.

(14) As illustrated in FIG. 1, the in-band full duplex transceiver 1000 according to the exemplary embodiment of the present invention includes a baseband digital region transmitting/receiving unit 100, an analog circuit region transmitting/receiving unit 200, and an antenna unit 300.

(15) The baseband digital region transmitting/receiving unit 100 includes a transmitting data generator 110 and an estimator 120. Although not illustrated in FIG. 1, the baseband digital region transmitting/receiving unit 100 may further include a receiving data demodulation generator. The transmitting data generator 110 performs encoding and modulation to generate a digital signal corresponding to data to be transmitted. In FIG. 1, an output signal of the transmitting data generator 110, that is, a baseband time domain signal is represented by x[m]. A receiving data demodulation generator (not illustrated in FIG. 1) performs demodulation and decoding on receiving data. The estimator 120 sets attenuations for each tab of a finite impulse response (FIR) filter 240. An estimator 120 according to an exemplary embodiment of the present invention serves to estimate input information of the FIR filter 240 to set the attenuations for each tab of the FIR filter 240. Further, the estimator 120 according to the exemplary embodiment of the present invention may use the estimated input information of the FIR filter 240 to set the attenuations for each tab of the FIR filter 240, thereby canceling SIC.

(16) The analog circuit region transmitting/receiving unit 200 includes a digital analog converter (DAC) 210, a mixer 220, a power amplifier (PA) 230, the FIR filter 240, a signal coupler 250, a low noise amplifier (LNA) 260, an integrator 270, and an analog digital converter (ADC) 280. The transmitting data generator 110, the DAC 210, the mixer 220, and the PA 230 form a transmitter. Further, the signal coupler 250, the LNA 260, the integrator 270, the ADC 280, and the receiving data demodulation generator form a receiver.

(17) The DAC 210 converts a digital signal into an analog signal, the mixer 220 uses a carrier frequency f.sub.C to convert a baseband frequency signal into a radio frequency band signal. The PA 230 amplifies and outputs the RF band signal. In FIG. 1, an output signal of the DAC 210 is represented by x.sub.b(t) and an output signal of the PA 230, that is, a final transmitted signal is represented by s(t). The s(t) that is the final transmitted signal includes a nonlinear component (harmonics component and phase noise component) due to the mixer 220 and the PA 230. As described below, the estimator 120 according to the exemplary embodiment of the present invention estimates the baseband equivalent information for the s(t) including the nonlinear component.

(18) The antenna unit 300 includes a transmitting antenna 310 and a receiving antenna 320. The transmitting antenna 310 transmits the transmitted signal s(t) to the outside and the receiving antenna 320 receives a received signal. The transmitting antenna 310 and the receiving antenna 320 physically approach each other, and therefore the transmitted signal s(t) is introduced through the receiving antenna 320 to act as an interference signal, in which the interference signal is a self-transmitted interference signal. The reason is that the transceiver 1000 according to the exemplary embodiment of the present invention is operated by an in-band full duplex type. The signal received through the receiving antenna 320 includes a self-transmitted interference signal as well as a normal received signal. In FIG. 1, the self-transmitted interference signal is represented by y(t). As described below, according to the exemplary embodiment of the present invention, the self-transmitted interference signal y(t) may be canceled using the FIR filter 240.

(19) The FIR filter 240 receives the transmitted signal s(t) and generates and outputs a signal minimizing the self-transmitted interference signal y(t). As illustrated in FIG. 1, the FIR filter 240 according to the exemplary embodiment of the present invention includes a plurality of delayers d.sub.0-d.sub.N-1, a plurality of attenuators a.sub.0-a.sub.N-1, and a signal coupler 241.

(20) The plurality of delayers d.sub.0-d.sub.N-1 each have a fixed delay. All delay intervals between the respective delayers d.sub.i (i=0, 1, . . . , N1) may be the same or different and the delayers may be divided into a plurality of groups having the same delay interval. Here, N represents the total number of tabs and d.sub.i (i=0, 1, 2, . . . , N1) represents a delay applied to the tab.

(21) The plurality of attenuators a.sub.0-a.sub.N-1 are connected to the plurality of delayers d.sub.0-d.sub.N-1, respectively and attenuate signals. Attenuations of the respective attenuators a.sub.i (i=0, 1, 2 . . . , N1) vary and are set by the estimator 120. Values of the attenuations of the respective attenuators a.sub.i (i=0, 1, 2 . . . , N1) may have one of a size, a size and phase, a real value, and a complex value. For convenience, the following description will be made under the assumption that the attenuation has the size or the real value.

(22) The signal coupler 241 couples output signals of the plurality of attenuators a.sub.0-a.sub.N-1. The signal coupler 241 adds and couples all the output signals of the plurality of attenuators a.sub.0-a.sub.N-1. In FIG. 1, the output signal of the signal coupler 241, that is, the output signal of the FIR filter 240 is represented by z(t).

(23) The signal coupler 250 couples the received signal with the output signal of the FIR filter 240 and then outputs the coupled signal to the LNA 260. The signal coupler 250 subtracts the output signal of the FIR filter 240 from the received signal and then couples the two signals. In this case, the FIR filter 240 outputs the signal minimizing the self-transmitted interference signal y(t), and therefore the signal coupler 250 outputs a signal obtained by canceling the self-transmitted interference signal y(t) from the received signal to the LNA 260. In FIG. 1, the output signal of the signal coupler 250 is represented by v(t).

(24) The LNA 260 receives the received signal from which the self-transmitted interference signal y(t) is canceled from the signal coupler 250 and cancels and amplifies noise from the received signal. The integrator 270 uses the carrier frequency f.sub.c to convert an RF band signal into a baseband frequency signal. The ADC 260 converts the baseband frequency signal that is an analog signal into the digital signal.

(25) FIG. 2 is a diagram illustrating an in-band full duplex transceiver 1000 according to another exemplary embodiment of the present invention.

(26) As illustrated in FIG. 2, the in-band full duplex transceiver 1000 according to another exemplary embodiment of the present invention includes the baseband digital region transmitting/receiving unit 100, an analog circuit region transmitting/receiving unit 200, and an antenna unit 300. Except that the antenna unit 300 is configured of one transmitting/receiving antenna and the analog circuit region transmitting/receiving unit 200 further includes a distributor 290, the in-band full duplex transceiver 1000 is the same as FIG. 1, and therefore the description of the overlapping components with FIG. 1 will be omitted.

(27) For the in-band full duplex type, the antenna unit 300 is configured of a transmitting/receiving antenna 330 that simultaneously performs a transmitter function and a receiver function. That is, through the antenna 330, the transmitted signal is transmitted and the received signal is received.

(28) The analog circuit region transmitting/receiving unit 200 further includes the distributor 290. The distributor 290 is connected to the antenna 330 and transmits the transmitted signal s(t) to the antenna 330. Further, the distributor 290 transmits the received signal received from the antenna 330 to a receiving module (signal coupler 250, LNA 260, or the like of FIG. 1). That is, the distributor 290 according to the exemplary embodiment of the present invention serves to transmit the transmitted signal to the antenna 330 and transmit the received signal to the receiving module. The distributor 290 may be implemented as a circulator, an electrical balance duplexer (EBD), etc. The circulator and the EBD may be appreciated by a person having ordinary skill in the art to which the present invention pertains and thus the detailed description thereof will be omitted.

(29) First, a method for processing, by an estimator 120, input information of an FIR filter 240 in a time domain to estimate the input information will be described. Input information z(t) of the FIR filter 240 includes a nonlinear component and the estimator 120 according to the exemplary embodiment of the present invention estimates the baseband equivalent information of s(t) including the nonlinear component. The method for estimating baseband equivalent information of z(t) described below will be described based on FIG. 2 but may also be equally applied to FIG. 1.

(30) To describe the method for estimating, by an estimator 120, input information of an FIR filter 240, various signals of the in-band full duplex transceiver 1000 are mathematically modeled.

(31) x[m] means a baseband time domain signal that is oversampled by d times. If d=1, x[m] means a time domain signal that is baseband-sampled. If the sampled signal passes through the DAC 210, a baseband analog signal x.sub.b(t) is generated by the following Equation 1.
x.sub.b(t)=x[m] sin c(Wtm)(Equation 1)

(32) In the above Equation 1, x[m] is x.sub.b(m/W) and sin c(t) is defined by the following Equation 2.

(33) $\begin{array}{cc}\sin c\left(t\right)=\frac{\sin \left(t\right)}{t}& \left(\mathrm{Equation}2\right)\end{array}$

(34) For the representation of the above Equation 1, a sampling theorem is applied. That is, the sampling theorem that all the baseband waveforms band-limited in W/2 may be represented by a linear combination of coefficients (i.e., x[m]) given by samples with an orthogonal basis {sin c(Wtm)}.sub.m is applied to the above Equation 1. In the above Equation 1, W may represent a bandwidth and may be a system bandwidth of the baseband and a bandwidth oversampled by d times. For convenience, the following description will be made under the assumption that the W is a bandwidth oversampled by one time (d=1), but a value of d may be set to be another value.

(35) If x.sub.b(t) passes through the mixer 220, it becomes the RF signal, and if the RF signal passes through the PA 230, it becomes the time domain input signal s(t) of the FIR filter 240. If the s(t) is mathematically represented, it becomes the following Equation 3.
s(t)={square root over (2P)}Re{{tilde over (x)}.sub.b(t)e.sup.j2f.sup.c.sup.t}(Equation 3)

(36) In the following Equation 3, {tilde over (x)}.sub.b(t) represents the baseband equivalent signal for the RF signal including the nonlinear component generated while the x.sub.b(t) passes through the PA 230. P represents transmit power amplified by the PA 230 and f.sub.C represents the carrier frequency of the mixer 220. It may be appreciated based on the above Equation 3 that the s(t) includes both of the linear component and the nonlinear component. Hereinafter, a frequency domain expression for the s(t) in the above Equation 3 is defined by SM. In this case, f represents the RF frequency.

(37) The RF signal s(t) in the above Equation 3 passes through the FIR filter 240 and the time domain expression of the passed RF signal becomes the following Equation 4.
z(t)=h.sub.fir(t)*s(t)(Equation 4)

(38) In the above Equation 4, h.sub.fir(t) represents a time impulse response for the RF region for the FIR filter 240. The baseband equivalent time impulse response for the h.sub.fir(t) becomes the following Equation 5.

(39) $\begin{array}{cc}{h}_{\mathrm{fir}}^b\left[l\right]=\underset{i}{.Math.}{a}_i{e}^{-j2f_c{d}_i}\sin c\left[l-d_{i}W\right]& \left(\mathrm{Equation}5\right)\end{array}$

(40) In the above Equation 5, d.sub.i (i=0, 1, . . . , N1) and a.sub.i (i=0, 1, . . . , N1) each represent time delays and path attenuations for each tab of the FIR filter 240. Further, l represents an index of a filter tab. If the above Equation 5 is applied to a digital (baseband) region equivalent signal z[m] for z(t) that is an output of the FIR filter 240, the z[m] becomes the following Equation 6.

(41) $\begin{array}{cc}\begin{array}{c}z\left[m\right]=\underset{l}{.Math.}\tilde{x}\left[m-l\right]h_{\mathrm{fir}}^b\left[l\right]\\ =\underset{l}{.Math.}\tilde{x}\left[m-l\right]\underset{i=0}{\overset{N-1}{.Math.}}{a}_i{e}^{-j2f_c{d}_i}\sin c\left[l-d_{i}W\right],\\ m=0,1,\mathrm{.Math.},M-1\end{array}& \left(\mathrm{Equation}6\right)\end{array}$

(42) In the Equation 6, by the above Equation 3, {tilde over (x)}[m]={tilde over (x)}.sub.b(m/W), in which M represents a maxim sample number of a time domain.

(43) Using the expression of the above Equations 4 to 6, the equivalent time domain expression y[m] in the digital (baseband) region for the self-transmitted interference signal y(t) becomes the following Equation 7.

(44) $\begin{array}{cc}\begin{array}{c}y\left[m\right]=\underset{l}{.Math.}\tilde{x}\left[m-l\right]h_{\mathrm{ch}}^b\left[l\right]\\ =\underset{l}{.Math.}\tilde{x}\left[m-l\right]\underset{i=0}{\overset{N^{}-1}{.Math.}}{b}_i{e}^{-j2f_c{}_i}\sin c\left[l-{}_{i}W\right],\\ m=0,1,\mathrm{.Math.},M-1\end{array}& \left(\mathrm{Equation}7\right)\end{array}$

(45) In the Equation 7, h.sub.ch.sup.b[l] represents the baseband (digital region) equivalent time impulse response for the channel generated by characteristics of the distributor 290 and the antenna 330. Further, b.sub.i and .sub.i each depend on the path attenuation and the time delay for each multipath i for the channel generated by the characteristics of the distributor 290 and the antenna 330. N represents the total number of multipaths for the channel.

(46) Referring to FIGS. 3 and 4, the method for processing, by an estimator 120 according to the exemplary embodiment of the present invention, input information of an FIR filter 240 in a time domain to estimate the input information will be described. In the following description, it is assumed that in each expression of the above Equations 5 to 7, the LNA 260, the integrator 270, and the ADC 280 do not have hardware impairment.

(47) FIG. 3 is a diagram illustrating a time domain frame structure according to an exemplary embodiment of the present invention.

(48) As illustrated in FIG. 3, the time domain frame according to the exemplary embodiment of the present invention includes an in-band full duplex (IFD) training field 410 and an IFD communication field 420. The IFD training field 410 is a separate field for transmitting a training signal to estimate the input information of the FIR filter 240. In the IFD training field 410, a first training signal 411 and a second training signal 412 are transmitted. The first training signal 411 and the second training signal 412 are separately generated signals to estimate the input information of the FIR filter 240 and two signals may be allocated in the frequency domain and may also be allocated in the time domain. Further, the first training signal 411 and the second training signal 412 may each be transmitted over at least one orthogonal frequency division multiplexing (OFDM) symbol section. Meanwhile, the first training signal 411 and the second training signal 412 are the same signal. Further, the IFD communication field 420 is a section in which data for in-band full duplex transmission/reception are transmitted.

(49) FIG. 4 is a flowchart illustrating a method for estimating input information of an FIR filter 240 according to an exemplary embodiment of the present invention.

(50) First, the estimator 120 uses the first training signal 411 to measure v[m] after self-interference cancellation (SIC) is performed (S410). The estimator 120 sets the attenuation a.sub.i (i=0, 1, 2 . . . , N1) of the FIR filter 240 to be any initial set value (corresponding to a(u) of Equation 22). Further, the transmitting data generator 110 generates the first training signal 411 and the estimator 120 measures the output information v[m] of the signal coupler 250 for the first training signal 411. Here, v[m] is a signal measured using the output of the ADC 280 for the first training signal 411 and the estimator 120 receives the output of the ADC 280 to measure v[m]. Meanwhile, using the above Equations 6 and 7, the baseband equivalent signal v[m] for the output signal v(t) of the signal coupler 250 is as the following Equation 8.
v[m]=y[m]z[m](Equation 8)

(51) The estimator 120 uses the second training signal 412 to measure the baseband equivalent signal y[m] of the self-interference signal (S420). That is, the estimator 120 measures the above Equation 7 that is the baseband equivalent signal y[m] for the self-transmitted interference signal y(t).

(52) For the estimator 120 to measure y[m], there is a need to introduce only the self-transmitted interference signal y(t) into the estimator 120. By a method for introducing a self-transmitted interference signal y(t) to an estimator 120, the estimator 120 may set the attenuation a.sub.i (i=0, 1, 2 . . . , N1) of the FIR filter 240 to be a very high value. By doing so, it is possible to physically suppress the output signal of the FIR filter 240 from being generated and only the self-transmitted interference signal y(t) is input to the estimator 120 through the LNA 260, the integrator 270, and the ADC 280. In this case, the estimator 120 measures the signal input from the ADC 280, in which the measured signal corresponds to the y[m].

(53) Next, the estimator 120 estimates the input information of the FIR filter 240 (S430). If the above Equation 8 is subtracted from the above Equation 7, it becomes a baseband equivalent signal z[m] of an output signal z(t) of the FIR filter 240 represented like the above Equation 6. That is, the estimator 120 may estimate z[m] by subtracting the v[m] measured in step S410 from the y[m] measured in step S420. If a matrix expression for the above Equation 6 is obtained to estimate the input information of the FIR filter 240 from the z[m] estimated as described above, it becomes the following Equation 9.
z=A.sub.{tilde over (x)}h.sub.fir(Equation 9)

(54) In the above Equation 9, A.sub.{tilde over (x)} is a Toeplitz matrix of {tilde over (x)}[m] and is represented by the following Equation 10.

(55) $\begin{array}{cc}{A}_{\tilde{x}}=[\begin{array}{cccccc}\tilde{x}\left[c\right]& \mathrm{.Math.}& \tilde{x}\left[1\right]& \tilde{x}\left[0\right]& \mathrm{.Math.}& \tilde{x}\left[c-1\right]\\ \tilde{x}\left[1-c\right]& \mathrm{.Math.}& \tilde{x}\left[0\right]& \tilde{x}\left[1\right]& \mathrm{.Math.}& \tilde{x}\left[c\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ \begin{array}{c}\tilde{x}[M-1-\\ c]\end{array}& \mathrm{.Math.}& \tilde{x}\left[M-2\right]& \tilde{x}\left[M-1\right]& \mathrm{.Math.}& \begin{array}{c}\tilde{x}[M+c-\\ 2]\end{array}\end{array}]\mathrm{or}{A}_{\tilde{x}}=[\begin{array}{cccccc}\tilde{x}\left[c-1\right]& \mathrm{.Math.}& \tilde{x}\left[1\right]& \tilde{x}\left[0\right]& \mathrm{.Math.}& \tilde{x}\left[c\right]\\ \tilde{x}\left[c\right]& \mathrm{.Math.}& \tilde{x}\left[2\right]& \tilde{x}\left[1\right]& \mathrm{.Math.}& \tilde{x}\left[1-c\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ \tilde{x}\left[c+M-2\right]& \mathrm{.Math.}& \tilde{x}\left[M\right]& \tilde{x}\left[M-1\right]& \mathrm{.Math.}& \tilde{x}\left[M-c-1\right]\end{array}]& \left(\mathrm{Equation}10\right)\end{array}$

(56) In the above Equation 10, c represents the number of non-causal elements (i.e., samples).

(57) In the above Equation 9, h.sub.fir is a vector expression for the above Equation 5 and is represented by the following Equation 11.
h.sub.firh.sub.fir.sup.b[c]h.sub.fir.sup.b[c+1] . . . h.sub.fir.sup.b[0]h.sub.fir.sup.b[1] . . . h.sub.fir.sup.b[c1]].sup.T or
h.sub.firh.sub.fir.sup.b[c1] . . . h.sub.fir.sup.b[1]h.sub.fir.sup.b[0]h.sub.fir.sup.b[1] . . . h.sub.fir.sup.b[c]].sup.T(Equation 11)

(58) In the above Equation 11, symbol T represents a transpose of the vector. Time delay values d.sub.0-d.sub.N-1 and the attenuations a.sub.i (i=0, 1, 2 . . . , N1) for each tab of the FIR filter 240 as described in the step S410 are set to be any initial value, and therefore the baseband time impulse response vector represented by the above Equation 11 may be obtained by calculation. That is, the estimator 120 already knows the baseband time impulse response vector as represented by the above Equation 11.

(59) By using a vector z of the Equation 9 estimated as described above and the previously acquired Equation 11, the baseband equivalent matrix for the input information of the FIR filter 240 is estimated by the following Equation 12.
.sub.{tilde over (x)}=zh.sub.fir.sup. or .sub.{tilde over (x)}=zh.sub.fir.sup.H(h.sub.firh.sub.fir.sup.H).sup.1(Equation 11)

(60) In the above Equation 12, h.sub.fir.sup. represents a pseudo-inverse matrix of b.sub.fir and h.sub.fir.sup.H represents the Hermitian of h.sub.fir.

(61) Finally, the estimator 120 extracts the {tilde over (x)}[m] that is the baseband equivalent component for the input signal of the FIR filter 240 in the matrix estimated by the above Equation 12. The {tilde over (x)}[m] extracted as described above is the baseband equivalent signal for the input signal of the FIR filter 240 to be obtained by the estimator 120. Further, as described above, the {tilde over (x)}[m] includes information on the nonlinear component as well as the information on the linear component.

(62) Meanwhile, to minimize the self-interference signal using the {tilde over (x)}[m], the estimator 120 resets the attenuation a.sub.i (i=0, 1, 2 . . . , N1) of the FIR filter 240. That is, the estimator 120 tunes the attenuation set as the initial value in the step S410 to be a value for minimizing the self-interference signal. Hereinafter, a method for tuning, by an estimator 120, an FIR filter 240 according to an exemplary embodiment of the present invention will be described with reference to FIG. 5.

(63) FIG. 5 is a flowchart illustrating a method for tuning, by an estimator 120, an FIR filter 240 according to an exemplary embodiment of the present invention.

(64) First, the estimator 120 according to the exemplary embodiment of the present invention estimates a delta response (S510). The v[m] that is the above Equation 8 may be represented by the following Equation 13.

(65) $\begin{array}{cc}\begin{array}{c}v\left[m\right]=y\left[m\right]-z\left[m\right]\\ =\underset{l}{.Math.}\tilde{x}\left[m-l\right]\left(h_{\mathrm{ch}}^b\left[l\right]-h_{\mathrm{fir}}^b\left[l\right]\right)\\ =\underset{l}{.Math.}\tilde{x}\left[m-l\right]h_{\mathrm{fir}}^b\left[l\right].\end{array}& \left(\mathrm{Equation}13\right)\end{array}$

(66) In the above Equation 13, h.sub.fir.sup.b[l] becomes the delta response to be estimated by the estimator 120. h.sub.fir.sup.b[l] becomes h.sub.ch.sup.b[l]h.sub.fir.sup.b[l]. Referring to the above Equation 5, h.sub.ch.sup.b[l]h.sub.fir.sup.b[l] is represented by the following Equation 14. h.sub.fir.sup.b[l]h.sub.ch.sup.b[l] represents the baseband equivalent time impulse response for a radio channel and h.sub.fir.sup.b[l] is the same as the above Equation 5. h.sub.ch.sup.b[l]

(67) $\begin{array}{cc}{h}_{\mathrm{fir}}^b\left[l\right]=\underset{i}{.Math.}{a}_i{e}^{-j2f_c{d}_i}\sin c\left[l-d_{i}W\right]& \left(\mathrm{Equation}14\right)\end{array}$

(68) In the above Equation 14, a.sub.i represents update attenuation of an i-th attenuator. If a matrix expression for the above Equation 13 is obtained, it becomes the following Equation 15.
v=A.sub.{tilde over (x)}h.sub.fir(Equation 15)

(69) In the above Equation 15, is the same as the above Equation 10 and is a value already obtained in FIG. 4. That is, A.sub.{tilde over (x)} corresponds to the matrix estimated by the above Equation 12. In the above Equation 15, h.sub.fir represents the delta response vector added to the current baseband equivalent time impulse response to update the time impulse response characteristics of the FIR filter 240. h.sub.fir h.sub.fir is represented by the following Equation 16. h.sub.fir
h.sub.firh.sub.fir.sup.b[c]h.sub.fir.sup.b[c+1] . . . h.sub.fir.sup.b[0]h.sub.fir.sup.r[1] . . . h.sub.fir.sup.b[c1]].sup.T
or
h.sub.firh.sub.fir.sup.b[c]h.sub.fir.sup.b[c+1] . . . h.sub.fir.sup.b[0]h.sub.fir.sup.b[1] . . . h.sub.fir.sup.b[c1]].sup.T
h.sub.firh.sub.fir[c1] . . . h.sub.fir.sup.b[1]h.sub.fir.sup.b[0]h.sub.fir.sup.b[1] . . . h.sub.fir.sup.b[c]].sup.T(Equation 16)

(70) In the above Equation 15, if a pseudo-inverse matrix for the A.sub.{tilde over (x)} is taken, the h.sub.fir is estimated by the following Equation 17.
.sub.fir=A.sub.{tilde over (x)}.sup.v(Equation 17)

(71) In the above Equation 17, A.sub.{tilde over (x)}.sup. represents the pseudo-inverse matrix of A.sub.{tilde over (x)}. The A.sub.{tilde over (x)} is a matrix already known by being estimated like the above Equation 12, and therefore A.sub.{tilde over (x)}.sup. the may be obtained by calculation. A.sub.{tilde over (x)} Further, v is also a value already obtained by being measured based on the first training signal in the step S410. Therefore, the delta response h.sub.fir to be estimated is estimated by the above Equation 17.

(72) Next, the estimator 120 according to the exemplary embodiment of the present invention calculates the attenuation a to be updated (S520). To calculate the update attenuation a of the FIR filter 240 from the delta response expression of the above Equation 14, the a of the above Equation 17 may be represented by the following Equation 18.
.sub.fir=sa(Equation 18)

(73) In the above Equation 18, s is represented by the following Equation 19.

(74) $\begin{array}{cc}s=[\begin{array}{ccc}\begin{array}{c}{e}^{-j2f_c{d}_0}\\ \sin c\left(-c-d_{0}W\right)\end{array}& \mathrm{.Math.}& \begin{array}{c}{e}^{-j2f_c{d}_{N-1}}\\ \sin c\left(-c-d_{N-1}W\right)\end{array}\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ \begin{array}{c}{e}^{-j2f_c{d}_0}\\ \sin c\left(c-1-d_{0}W\right)\end{array}& \mathrm{.Math.}& \begin{array}{c}{e}^{-j2f_c{d}_{N-1}}\\ \sin c\left(c-1-d_{N-1}W\right)\end{array}\end{array}]\mathrm{or}\text{.Math.}s=\left[\begin{array}{ccc}\begin{array}{c}{e}^{-j2f_c{d}_0}\\ \sin c\left(c-1-d_{0}W\right)\end{array}& \mathrm{.Math.}& \begin{array}{c}{e}^{-j2f_c{d}_{N-1}}\\ \sin c\left(c-1-d_{N-1}W\right)\end{array}\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ \begin{array}{c}{e}^{-j2f_c{d}_0}\\ \sin c\left(-c-d_{0}W\right)\end{array}& \mathrm{.Math.}& \begin{array}{c}{e}^{-j2f_c{d}_{N-1}}\\ \sin c\left(-c-d_{N-1}W\right)\end{array}\end{array}\right]& \left(\mathrm{Equation}19\right)\end{array}$

(75) In the s like the above Equation 19, f.sub.c is the carrier frequency and is a value already known, d.sub.0-d.sub.N-1 are a fixed delay value and are a preset value, and W is a bandwidth and is a value already known. Therefore, the s may be obtained by the above Equation 19.

(76) In the above Equation 18, the a is represented by the following Equation 20.
a=[.sub.0.sub.1 . . . .sub.N-1].sup.T(Equation 20)

(77) In the above Equation 18, the .sub.fir is the vector obtained by the estimation in step S510 and s is also the value calculated by the above Equation 19, which is already known. Therefore, a is obtained by the following Equation 21.
a=s.sup..sub.fir or =(s.sup.Hs).sup.1s.sup.H.sub.fir
a=s.sup..sub.fir(Equation 21)

(78) Finally, the estimator 120 according to the exemplary embodiment of the present invention uses the a to tune the attenuation of the FIR filter 240 (S530). By using the a obtained by the above Equation 21, the estimator 120 determines attenuation a(u+1) to be updated by the following Equation 22.
a(u+1)=a(u)+a(Equation 22)

(79) In the above Equation 22, a(u) represents the vector of the current attenuation and corresponds to any initial set value described in the step S410 of FIG. 4. Further, a(u+1) represents the vector of the attenuation to be finally updated.

(80) The FIR filter 240 applies the attenuation a(u+1) to be updated determined by the above Equation 22 to the attenuator of the FIR filter 240. By doing so, the estimator 120 according to the exemplary embodiment of the present invention may minimize the self-interference signal.

(81) According to the exemplary embodiment of the present invention as described above, the attenuation of the FIR filter 240 is set using the baseband equivalent signal for the input signal of the FIR filter 240 including the linear component and the nonlinear component, thereby obtaining the high SIC gain.

(82) Next, a method for processing, by an estimator 120, input information of a FIR filter 240 in a frequency domain to estimate the input information will be described.

(83) Referring to FIGS. 1 and 2, the RF signal like the above Equation 3 passes through the FIR filter 240 and a frequency domain expression for a transfer function of the FIR filter 240 is as the following Equation 23.

(84) $\begin{array}{cc}{H}_{\mathrm{fir}}\left(f\right)=\underset{q=0}{\overset{N-1}{.Math.}}{a}_q{e}^{-j2{\mathrm{fd}}_q}& \left(\mathrm{Equation}23\right)\end{array}$

(85) In the above Equation 23, d.sub.q represents a time delay of a q-th tab and a.sub.q represents path attenuation of the q-th tab. By using the expressions of the above Equations 3 and 23, the frequency domain and time domain expressions for the output signal of the FIR filter 240 are as the following Equation 24.

(86) 0$\begin{array}{cc}Z\left(f\right)=H_{\mathrm{fir}}\left(f\right)S\left(f\right)=\left(\underset{q=0}{\overset{N-1}{.Math.}}{\hat{a}}_q{e}^{-j2{\mathrm{fd}}_q}\right)S\left(f\right),z\left(t\right)=h_{\mathrm{fir}}\left(t\right)*s\left(t\right)& \left(\mathrm{Equatio}n24\right)\end{array}$

(87) By using the expressions of the above Equations 3 and 23, the self-transmitted interference signal y(t) and the frequency domain expression thereof are as the following Equation 25.

(88) $\begin{array}{cc}Y\left(f\right)=H_{\mathrm{ch}}\left(f\right)S\left(f\right)=\left(\underset{i=0}{\overset{N^{}-1}{.Math.}}{b}_i{e}^{-j2f{}_i}\right)S\left(f\right),y\left(t\right)=h_{\mathrm{ch}}\left(t\right)*s\left(t\right)& \left(\mathrm{Equation}25\right)\end{array}$

(89) In the above Equation 25, h.sub.ch(t) represents the time impulse response to a channel generated by characteristics of the distributor 290 and the antenna 330 and H.sub.Ch(f) represents the frequency domain expression of the h.sub.ch(t). Further, b.sub.i and .sub.i each depend on the path attenuation and the time delay for each multipath i for the channel generated by the characteristics of the distributor 290 and the antenna 330. N represents the total number of multipaths for the channel.

(90) By using the expressions of the Equations 24 and 25, the frequency domain expression for the output signal v(t) of the signal coupler 250 is as the following Equation 26.

(91) $\begin{array}{cc}\begin{array}{c}V\left(f\right)=Y\left(f\right)-Z\left(f\right)\\ =H_{\mathrm{ch}}\left(f\right)S\left(f\right)-H_{\mathrm{fir}}\left(f\right)S\left(f\right)\\ =\left(H_{\mathrm{ch}}\left(f\right)-H_{\mathrm{fir}}\left(f\right)\right)S\left(f\right).\end{array}& \left(\mathrm{Equation}26\right)\end{array}$

(92) Next, referring to FIG. 6, a method for processing, by the estimator 120 according to another exemplary embodiment of the present invention, the input information of the FIR filter 240 in the frequency domain to estimate the input information will be described. In the following description, it is assumed that in each expression of the above Equations 24 to 26, the LNA 260, the integrator 270, and the ADC 280 do not have hardware impairment.

(93) FIG. 6 is a flowchart illustrating a method for processing input information of an FIR filter 240 in a frequency domain to estimate the input information according to another exemplary embodiment of the present invention.

(94) First, the estimator 120 uses a first training signal 411 to obtain V[k] after the self-interference cancellation (SIC) is performed (S610). The estimator 120 sets the attenuation a.sub.i (i=0, 1, 2 . . . , N1) of the FIR filter 240 to be any initial value. Further, the transmitting data generator 110 generates the first training signal 411 and the estimator 120 obtains the output information V(k) of the signal coupler 250 for the first training signal 411. Here, the V[k] is a signal obtained by performing fast Fourier transform (FFT) on the output of the ADC 280 due to the first training signal 411 and the estimator 120 acquires the V(k) by receiving the output of the ADC 280 and performing the FFT on the received output. That is, the V[k] is a value obtained by performing the FFT on a baseband version (equivalent) information (i.e., output information of the ADC 280) for the v(t). If the V(k) is mathematically represented, it becomes the following Equation 27.

(95) $\begin{array}{cc}\begin{array}{c}V\left(k\right)=Y\left(k\right)-Z\left(k\right)\\ =H_{\mathrm{ch}}\left(k\right)S\left(k\right)-H_{\mathrm{fir}}\left(k\right)S\left(k\right)\\ =\left(H_{\mathrm{ch}}\left(k\right)-H_{\mathrm{fir}}\left(k\right)\right)S\left(k\right),\end{array}k=0,1,2,\mathrm{.Math.},F-1& \left(\mathrm{Equation}27\right)\end{array}$

(96) In the above Equation 27, k represents a subscriber index and F represents a size of the FFT. The V(k) of the above Equation 27 may be a version that is baseband sampled for the above Equation 26. Meanwhile, in step S610, the attenuation a.sub.i (i=0, 1, 2 . . . , N1) of the FIR filter 240 is set to be any value and the time delay d.sub.0-d.sub.N-1 is a fixed value, such that H.sub.fir(k) of the above Equation 27 may be calculated by the following Equation 28.

(97) $\begin{array}{cc}{H}_{\mathrm{fir}}\left(k\right)=\underset{q=0}{\overset{N-1}{.Math.}}{a}_q{e}^{-j\frac{2f_s{\mathrm{kd}}_q}{F}}\mathrm{or}{H}_{\mathrm{fir}}\left(k\right)=\underset{q=0}{\overset{N-1}{.Math.}}{a}_q{e}^{-j\frac{2{\mathrm{kd}}_q}{F}}& \left(\mathrm{Equation}28\right)\end{array}$

(98) In the above Equation 28, f.sub.s represents a baseband sampling rate and H.sub.fir(k) may be a version that is baseband sampled for H.sub.fir(f) defined in the above Equation 23.

(99) The estimator 120 uses the second training signal 412 to acquire Y(k) corresponding to the self-interference signal (S620). That is, the estimator 120 performs the FFT on the baseband version (equivalent) information (i.e., output information of the ADC 280) for the self-transmitted interference signal y(t) to acquire the Y(k) by the following Equation 29.
Y(k)H.sub.ch(k)S(k)(Equation 29)

(100) The Y(k) like the above Equation 29 may be a version that is baseband sampled for Y(f) defined in the above Equation 25.

(101) For the estimator 120 to acquire the Y(k), there is a need to introduce only the self-transmitted interference signal y(t) into the estimator 120. By a method for introducing a self-transmitted interference signal y(t) to an estimator 120, the estimator 120 may set the attenuation a.sub.i (i=0, 1, 2 . . . , N1) of the FIR filter 240 to be a very high value. By doing so, it is possible to physically suppress the output signal of the FIR filter 240 from being generated and only the self-transmitted interference signal y(t) is input to the estimator 120 through the LNA 260, the integrator 270, and the ADC 280. In this case, the estimator 120 performs the FFT on the signal input from the ADC 280 to acquire the Y(k).

(102) Next, the estimator 120 estimates the input information of the FIR filter 240 (S630). If the above Equation 27 is subtracted from the above Equation 29, it becomes Z(k) by the following Equation 30. That is, the estimator 120 may calculate the Z(k) by subtracting the V(k) measured in the step S610 from the Y(k) measured in the step S620. The Z(k) calculated as described above is represented by the following Equation 30.
Y(k)(k)=Z(k)=H.sub.fir(k)S(k)(Equation 30)

(103) As described in the above Equation 28, the H.sub.fir(k) is a value already calculated, and therefore the estimator 120 may use the above Equation 30 to estimate S(k) that is the input information of the FIR filter 240. The S(k) corresponds to the baseband equivalent frequency domain input information of the FIR filter 240. The expression for the S(k) estimated as described above is as the following Equation 31.

(104) $\begin{array}{cc}\hat{S}\left(k\right)=\frac{Y\left(k\right)-V\left(k\right)}{H_{\mathrm{fir}}\left(k\right)}& \left(\mathrm{Equation}31\right)\end{array}$

(105) When the estimator 120 performs inverse FFT (IFFT) on the above Equation 31, the baseband equivalent time domain input information of the FIR filter 240 may be acquired.

(106) Meanwhile, the estimator 120 according to another exemplary embodiment of the present invention uses the S(k) estimated in FIG. 6 to minimize the self-interference signal, thereby resetting the attenuation a.sub.i (i=0, 1, 2 . . . , N1) of the FIR filter 240. That is, the estimator 120 tunes the attenuation set as the initial value in the step S610 to be a value for minimizing the self-interference signal. Hereinafter, a method for tuning, by the estimator 120 according to another exemplary embodiment of the present invention, the FIR filter 240 will be described with reference to FIG. 7.

(107) FIG. 7 is a flowchart illustrating a method for tuning, by an estimator 120, an FIR filter 240 according to another exemplary embodiment of the present invention.

(108) First, the estimator 120 according to the exemplary embodiment of the present invention estimates a delta frequency response based on de-patterning (S710). In connection with the above Equation 27, the de-patterning is performed on X(k) by the following Equation 32.

(109) $\begin{array}{cc}\frac{V\left(k\right)}{X\left(k\right)}=\left(H_{\mathrm{ch}}\left(k\right)-H_{\mathrm{fir}}\left(k\right)\right)\left[\frac{S\left(k\right)}{X\left(k\right)}\right]=H_{\mathrm{fir}}\left(k\right)\left[\frac{S\left(k\right)}{X\left(k\right)}\right],k=0,1,2,\mathrm{.Math.},F-1& \left(\mathrm{Equation}32\right)\end{array}$

(110) In the above Equation 32, the X(k) is an expression of performing the FFT on the x[m] given in the above Equation 1 and H.sub.fir(k) is H.sub.ch(k)H.sub.fir(k) and represents the delta frequency response.

(111) If the above Equation 32 is divided by S(k)/X(k), the delta frequency response may be obtained by the following Equation 33.

(112) $\begin{array}{cc}\frac{V\left(k\right)}{X\left(k\right)}\left(\frac{X\left(k\right)}{S\left(k\right)}\right)=H_{\mathrm{fir}}\left(k\right),k=0,1,2,\mathrm{.Math.},F-1& \left(\mathrm{Equation}33\right)\end{array}$

(113) In the above Equation 33, the V(k) value is a value calculated in the step S610 of FIG. 6, the X(k) is also a value already known, and the S(k) is a value already estimated in FIG. 6. Therefore, the delta frequency response to be estimated may use the value already calculated to be calculated like the above Equation 33.

(114) Meanwhile, in the above Equation 33, the H.sub.fir(k) may be represented by the following Equation 34.

(115) $\begin{array}{cc}{H}_{\mathrm{fir}}\left(k\right)=\underset{q=0}{\overset{N-1}{.Math.}}{a}_q{e}^{-j\frac{2f_s{\mathrm{kd}}_q}{F}}\mathrm{or}\text{.Math.}{H}_{\mathrm{fir}}\left(k\right)=\underset{q=0}{\overset{N-1}{.Math.}}{a}_q{e}^{-j\frac{2{\mathrm{kd}}_q}{F}}{H}_{\mathrm{fir}}\left(k\right)=\underset{q=0}{\overset{N-1}{.Math.}}{a}_q{e}^{-j\frac{2f_s{\mathrm{kd}}_q}{F}}& \left(\mathrm{Equation}34\right)\end{array}$

(116) In the above Equation 34, a.sub.q represents an update attenuation of q-th attenuator and f.sub.s represents a baseband sampling rate.

(117) Next, the estimator 120 according to another exemplary embodiment of the present invention calculates the attenuation a to be updated (S720). To calculate the update attenuation a of the FIR filter 240 from the delta frequency response expression of the above Equation 34, the above Equation 34 may be represented by a matrix expression by the following Equation 35.

(118) $\begin{array}{cc}{H}_{\mathrm{fir}}=Fa[\begin{array}{c}{H}_{\mathrm{fir}}\left(0\right)\\ H_{\mathrm{fir}}\left(1\right)\\ \mathrm{.Math.}\\ H_{\mathrm{fir}}\left(F-1\right)\end{array}]=[\begin{array}{cccc}1& 1& \mathrm{.Math.}& 1\\ e^{-j\frac{2f_s{d}_0}{F}}& e^{-j\frac{2f_s{d}_1}{F}}& \mathrm{.Math.}& e^{-j\frac{2f_s{d}_{N-1}}{F}}\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ e^{-j\frac{2f_s\left(F-1\right)d_0}{F}}& e^{-j\frac{2f_s\left(F-1\right)d_1}{F}}& \mathrm{.Math.}& e^{-j\frac{2f_s\left(F-1\right)d_{N-1}}{F}}\end{array}][\begin{array}{c}{a}_0\\ a_1\\ \mathrm{.Math.}\\ a_{N-1}\end{array}]\mathrm{or}{H}_{\mathrm{fir}}=Fa[\begin{array}{c}{H}_{\mathrm{fir}}\left(0\right)\\ H_{\mathrm{fir}}\left(1\right)\\ \mathrm{.Math.}\\ H_{\mathrm{fir}}\left(F-1\right)\end{array}]=[\begin{array}{cccc}1& 1& \mathrm{.Math.}& 1\\ e^{-j\frac{2d_0}{F}}& e^{-j\frac{2d_1}{F}}& \mathrm{.Math.}& e^{-j\frac{2d_{N-1}}{F}}\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ e^{-j\frac{2\left(F-1\right)d_0}{F}}& e^{-j\frac{2\left(F-1\right)d_1}{F}}& \mathrm{.Math.}& e^{-j\frac{2\left(F-1\right)d_{N-1}}{F}}\end{array}]\left[\begin{array}{c}{a}_0\\ a_1\\ \mathrm{.Math.}\\ a_{N-1}\end{array}\right].& \left(\mathrm{Equation}35\right)\end{array}$

(119) In the above Equation 35, f.sub.s is the baseband sampling rate and is a value already known, F is a size of the FFT and is a value already known, and d.sub.0-d.sub.N-1 are fixedly delayed values and are preset values. Therefore, a matrix F may be a value that may be obtained by calculation.

(120) In above Equation 35, if the inverse matrix or the pseudo-inverse matrix is performed on the F, the a is represented by the following Equation 36.
a=F.sup.H.sub.fir(Equation 36)

(121) In the above Equation 36, F.sup. represents an inverse matrix or a pseudo-inverse matrix of the F when N=F and represents the inverse matrix of the F when N<F. H.sub.fir is a vector obtained by the previous estimation in step S710 and F is also a value calculated by the above Equation 35 and is already known, such that the a is calculated based on the above Equation 36. H.sub.fir

(122) Finally, the estimator 120 according to another exemplary embodiment of the present invention uses the a to tune the attenuation of the FIR filter 240 (S730). By using the a by the above Equation 36, the estimator 120 determines attenuation a(u+1) to be updated by the following Equation 37.
a(u+1)=a(u)a(Equation 37)

(123) In the above Equation 37, a(u) represents the vector of the current attenuation and corresponds to any initial set value described in the step S610 of FIG. 6. Further, a(u+1) represents the vector of the attenuation to be finally updated.

(124) The FIR filter 240 applies the attenuation a(u+1) to be updated determined by the above Equation 37 to the attenuator of the FIR filter 240. By doing so, the estimator 120 according to another exemplary embodiment of the present invention may minimize the self-interference signal.

(125) Meanwhile, unlike the above Equation 32, the estimator 120 according to another exemplary embodiment of the present invention may perform the de-patterning on the S(k) in the above Equation 27 by the following Equation 38.

(126) 0$\begin{array}{cc}\frac{V\left(k\right)}{S\left(k\right)}=\left(H_{\mathrm{ch}}\left(k\right)-H_{\mathrm{fir}}\left(k\right)\right)=H_{\mathrm{fir}}\left(k\right),k=0,1,2,\mathrm{.Math.},F-1& \left(\mathrm{Equation}38\right)\end{array}$

(127) In the above Equation 38, the V(k) value is a value calculated in the step S610 of FIG. 6 and the S(k) is a value already estimated in FIG. 6. Therefore, the delta frequency response to be estimated may use the value already calculated to be calculated like the above Equation 38.

(128) According to another exemplary embodiment of the present invention as described above, the input signal of the FIR filter 240 including the linear component and the nonlinear component is processed and estimated in the frequency domain and the attenuation of the FIR filter 240 is set using the estimated input signal, thereby obtaining the high SIC gain.

(129) While this invention has been described in connection with what is presently considered to be practical exemplary embodiments, it is to be understood that the invention is not limited to the disclosed embodiments, but, on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims.