The present invention relates to a radio communication device and particularly relates to a radio communication technique that enables to precisely estimate a channel response even in the circumstance of strong influences of interference waves coming from adjacent cells and to correctly demodulate an information signal.
In recent years, many more users demand fast data transmission in a radio communication system as the volume of communication increases. The multicarrier transmission represented by OFDM (Orthogonal Frequency Division Multiplexing) gets attention as a way of communication to realize the fast and high-volume data transmission. The OFDM, which is used in IEEE 802.11a being a radio system of a 5 GHz-band or digital terrestrial broadcast, provides for simultaneous communication by arranging tens to thousands of carriers in a minimum frequency interval that does not induce interference theoretically. Generally in OFDM, these carriers are referred to as subcarriers which are digitally modulated with PSK, QAM or the like for communication. It is known that OFDM and error correction are combined to obtain strong tolerance to frequency selective fading.
The patent literature 1 discloses a configuration example of a receiving device used in a communication system using the above OFDM.
A guard interval is first removed in the GI removing unit 1001 from a received signal that has been received in the antenna unit 1000 of the receiver shown in
The pilot signal transmitted to the IFFT unit 1007 is converted from a frequency domain signal into a time domain signal and a delay profile of a channel response is estimated in the delay profile estimating unit 1008. Then, the delay profile of the channel response estimated from the known pilot signal as described above is converted from a time domain signal into a frequency domain signal in the FFT unit 1009. The processing enables to seek a channel response in frequency domain.
The channel response in frequency domain obtained by the above processing is transmitted to the demodulating unit 1003-1 to -N to which an information signal has been transmitted, for use in compensation of the channel response, thereby enabling demodulation by compensating for influence of multipath fading and accordingly enabling to correctly demodulate an information signal.
The above prior art enables to estimate variance of a channel response for an information signal and compensate for influence of the variance. However, the art has a problem in that influence of interference waves coming from adjacent cells significantly degrades the accuracy of the estimation processing of a channel response in the circumstance of another OFDM system nearby using the same frequency band, or in the circumstance of a terminal being around a cell edge in an OFDM cellular system in which all cells use the same frequency band; as a result, an information signal cannot be demodulated correctly.
The art has another problem in that if a base station specific code is used in an OFDM symbol for estimating a channel response, to identify the base station in a cellular system, the base station specific code needs to be detected to estimate the channel response. This is problematic in that a usual correlation method in a frequency domain does not correctly work in the environment with strong frequency selective fading.
Similarly, in the case that a transmitting device implements a transmit diversity technique using a plurality of transmitting antennas, codes to estimate channel responses use the same code. That is, a technique has not been disclosed for both estimating channel responses and identifying antennas, using the same OFDM symbol.
It is an object of the present invention to precisely estimate a channel response and correctly demodulate an information signal even in the above circumstances.
In a radio communication technique according to the present invention, a transmitting device simultaneously transmits OFDM symbols for estimating channel responses, modulated with different codes for transmitting antennas, while a receiver converts the received OFDM symbols for estimating channel responses into frequency domain data, selects the codes used by the transmitting device for making a complex conjugate, and multiplies the obtained frequency domain data by the conjugate to calculate a frequency response of the channel response. Then, the signal is converted into a time domain signal to calculate a delay profile of the channel response. Additionally, the calculated delay profile is multiplied by an appropriate time window and again converted into frequency domain data, thereby allowing for estimating a highly accurate a channel response in frequency domain with reduced influence of interference, as well as enabling to identify antennas.
The estimating way of a channel response for each antenna and the identifying way of antennas are applied to a base station in a cellular system, thereby enabling to identify a base station and to highly accurately estimate a frequency response of a signal transmitted from each base station.
In the environment to receive a plurality of signals, a replica is created for an OFDM signal for estimating a channel response, which has resulted in an interferential signal, and is subtracted from a received signal. This improves the accuracy in estimating a channel response in frequency domain of a desired signal.
The radio communication technique according to the present invention applies the above antenna identification technique to select a transmitting antenna in transmit diversity if a transmitting device comprises a plurality of antennas. That is, a symbol to estimate a channel response using different codes for the antennas is transmitted and transmit diversity is performed to switch to a transmitting antenna. If a result of seeking a delay profile using a pre-determined code exceeds a threshold of power, it is determined that an antenna using such a code is used as a transmitting antenna. Similarly, a receiving side estimates the number of transmitting antennas in a MIMO system.
As described in the above, according to the present invention, an OFDM symbol to estimate channel responses is received that has been modulated with different codes for transmitting antennas, converted into a frequency domain signal, multiplied by complex conjugate signals of the codes used in a transmitting device, and converted into a time domain signal to calculate a delay profile. Then, the obtained delay profile is multiplied by an appropriate time window and again converted into a frequency domain signal, thereby enabling to calculate a highly accurate frequency response of a channel response. Additionally, the receiving device selects a code for the multiplication, thereby enabling simultaneous identification of a transmitting antenna.
Further, the highly accurately calculating way of frequency responses for antennas and the identifying way of a transmitting antenna are applied to a cellular system, thereby enabling to identify base stations and estimate channel responses.
Furthermore, a replica is created for an OFDM symbol for estimating a channel response, which has resulted in an interference wave, and is subtracted from a received signal, improving estimating institution of channel responses.
In a transmit diversity system, it is possible to select or switch a transmitting antenna.
In the MIMO system, it is possible to estimate the number of transmitting antennas used in a transmitting side, so that the number of transmitting antennas does not need to be notified on a communication protocol, thereby improving throughput in an entire system.
First of all, a parameter of an OFDM (Orthogonal Frequency Division Multiplexing) signal used herein is defined as follows: the number of subcarriers used for OFDM is represented as N, the length of an OFDM symbol is represented as Tsym, and the length of a guard interval is represented as Tgi.
The equation (1) means in view of the frequency axis that a real axis component of a subcarrier corresponding to fk is represented as ak and an imaginary axis component is represent as bk: Ck=ak+jbk (j is a unit of an imaginary number, i.e. j×j=−1). As antenna specific values, ak, x, bk and x are decided.
The present invention is characterized in that respective transmitting antennas transmit different preambles while a receiving device estimates a delay profile from each of antennas as well as estimates a channel response or identifies a transmitting antenna. If codes of code length N to generate preambles are Ckx and Cky, correlation between Ckx and Cky is represented as the following equation (2):
In the above equation, the symbol “*” means a complex conjugate. It is preferable to select codes such that Cor gives a smaller value, in selecting different codes for antennas.
In the following, radio communication techniques according to the embodiments of the present invention will be described.
First, a radio communication technique according to a first embodiment of the present invention will be described with reference to the drawings. The radio communication technique according to this embodiment is characterized by defining symbols used for estimating channel responses, which are transmitted simultaneously from a plurality of antennas, as different series (preambles) for the antennas, thereby separating and calculating delay profiles of signals transmitted from the antennas.
A radio communication device according to this embodiment is an example in which a transmitting device performs transmission using a plurality of antennas, and particularly, the technique according to the present invention is applied to transmit diversity.
The radio communication technique according to the first embodiment of the present invention is directed to downlink transmission, provides a plurality of antennas for a transmitting (base station) side, and relates to an antenna selection technique in performing transmitting antenna selection diversity. However, the invention should not be limited by this example and is applicable to uplink or other communications. This embodiment simultaneously transmits OFDM symbols for channel response estimation modulated with specific codes, from a plurality of transmitting antennas of a transmitting device. An example is illustrated in which a receiving side detects channel response estimating symbols transmitted from respective antennas, estimates which antenna among a plurality of antennas used in a transmitting side has transmitted a signal of the best quality (with the highest received power), and estimates a frequency response of a channel response.
The radio communication technique according to the first embodiment of the present invention will be described with reference to the drawings.
In the transmitting device of the above base station, the preamble (A, B0, B1) storage selecting unit 11 assigns a code A, a code B0 and a code B1 to carriers, which are stored in waveforms subjected to OFDM signal processing. If the antenna 15-a is selected according to antenna selection information, then the preamble (A, B0, B1) storage selecting unit 11-a outputs the preambles A, B0 and the preamble (A, B0, B1) storage selecting unit 11-b outputs the preamble A, B1. If the antenna 15-b is selected according to antenna selection information, then the preamble (A, B0, B1) storage selecting unit 11-a outputs the preambles A, B1 and the preamble (A, B0, B1) storage selecting unit 11-b outputs the preambles A, B0. That is, a preamble (A, B0, B1) storage selecting unit connected to a selected antenna outputs the preambles A, B0.
The following will describe the case that the antenna 15-a is selected according to antenna selection information. First, preambles are transmitted in transmitting data; that is, the switch units 12-a and 12-b first select and transmit preambles. At that time, the preamble (A, B0, B1) storage selecting unit 11-a outputs the preambles A and B0, while the preamble (A, B0, B1) storage selecting unit 11-b selects the preambles A and B1.
After the preambles have been transmitted, the switch unit 12-a performs switching to select data whose information data has been subjected to OFDM transmission processing from an error correction coding unit 1 to a GI (guard interval) inserting unit 6. On the other hand, the switch unit 12-b does not transmit data any more after the preamble transmitting.
The same transmitting operation is repeated until antenna switching is requested according to antenna information. When antenna switching is requested (the antenna 15-b is selected), preamble patterns for the preamble (A, B0, B1) storage selecting units 11-a and 11-b and control methods of the switch units 12-a and 12-b are replaced by each other.
Next, a receiving device of a terminal according to this embodiment will be described. As shown in
As described above, the preambles A, B0 and A, B1 are simultaneously transmitted from different antennas in the transmitting device of the base station. In a terminal receiving device, the preambles go through different channel responses and are simultaneously received by the single antenna 41. Signals that went through different channel responses go through the radio receiving unit 42 and the A/D converting unit 43 and are inputted to the synchronizing unit 44. The synchronizing unit 44 uses the preamble A to establish symbol synchronization; later processing is conducted at appropriate time.
Next, after a guard interval attached on a transmitting side is removed in the GI removing unit 45, signals are converted into parallel signals in the S/P converting unit 46 and inputted to the DFT unit 47. Then, the DFT unit 47 converts a received time domain signal into a frequency domain signal. The preambles B0 and B1 are simultaneously received, so that a frequency domain signal added with the preambles is inputted to the buffer unit 54. The inputted data is labeled R×B.
Next, the code selecting unit 55 selects a complex conjugate signal of the code B0 (a preamble transmitted from an antenna to which data is transmitted) and the code multiplying unit 48 multiplies the signal by R×B. The IDFT unit 49 performs an IDFT operation on the signal. The signal can be handled as a delay profile from an antenna that transmitted the preamble B0 (described below with reference to
Following the above sequence of processing or in parallel to the above sequence of processing, the code selecting unit 55 selects the code B1 and the code multiplying unit 48 multiplies the code by R×B. The result is subjected to IDFT in the IDFT unit 49, such that it is possible to seek a delay profile from an antenna that has transmitted preamble B1.
In the example shown in
The power measuring unit 56 estimates power from the delay profiles and seeks which antenna has a propagation path with higher power. Then, the result is notified as antenna change information from a transmitter (not shown) to a base station. If a change is requested to use an antenna having a propagation path with higher power, the base station changes an antenna as described above for communication.
As described in the above, the preamble B1 to select an antenna is simultaneously transmitted for the preamble B0 needed to transmit data and the above device configuration and processing are accomplished, so that time to estimate an antenna is not further needed to configure a transmit diversity system. Additionally, according to the method, a code of an antenna for transmitting data is defined to be B0, so that even if there is no antenna selection information beforehand, a transmitting side can select an antenna arbitrarily.
Moreover, assuming that an antenna control system for a transmit diversity described herein is an antenna of each base station in a cellular system, it is possible to identify a base station. However, since a transmitting antenna of each base station transmits an OFDM symbol for estimating a channel response, which is modulated with a code specific to the antenna, codes of antennas are not replaced by each other like the transmitting device according to the first embodiment.
If a terminal is at a position to receive electric waves from a plurality of base stations (such as at a cell edge or at a sector edge), it receives OFDM signals for estimating channel responses, which are modulated with different codes and simultaneously transmitted. The terminal can identify the next base station to be connected by performing the processing described in the first embodiment when a signal from a current connected base station weakens or when the quality of a signal declines. However, it is not easy for a distant base station to completely synchronize to receive an OFDM signal for estimating a channel response. However, complete synchronization is not needed since a guard section is set for an OFDM signal, but the station only needs to synchronize to some extent to prevent a problem.
Although the example of two transmitting antennas has been described in relation to the first embodiment, the number of base stations from which simultaneous reception is possible is not limited to two in a cellular system. Moreover, codes used by the base stations are not known. In that case, a method of notifying information such as codes of surrounding base stations from a current connected base station is adopted. A terminal selects a code based on the information and selects the next base station to be connected.
Next, a radio communication technique according to a second embodiment of the present invention will be described with reference to the drawings. The radio communication technique according to this embodiment is an example in which the radio communication technique is applied to MIMO (Multi-Input Multi-Output; hereinafter, referred to as MIMO) using a plurality of transmitting/receiving antennas in communication.
In a MIMO system, a transmitting/receiving station must know the number of antennas in use. The radio communication technique according to this embodiment includes the used antenna estimating unit 60 shown in
As above, a receiving device can estimate the number of transmitting antennas using an OFDM symbol to estimate a channel response in a MIMO system. As such, an efficient communication system can be realized compared to a communication system which requires previous notification of the number of antennas to be used.
When the reception starts (k=0), OFDM synchronization is performed at step S201, and then frequency conversion is performed for a region of a preamble B using DFT at step S202 (R×B). The processing in these steps is the same processing as S101 and S102, respectively. The R×B is multiplied by a complex conjugate of a code Bk (an integer being 0≦k<M) at step S203, and then a delay profile is estimated using IDFT at step S204. Power P is estimated from the delay profile at step S205, and then it is determined whether or not the power is above a threshold at step S206. If the power is not above the threshold (no), the processing ends. If the power is above the threshold, a delay profile sought using the code Bk is time-filtered at step S207, and then a channel response is estimated using DFT at step S208. It is determined whether or not k equals M−1 at step S209. If it equals (yes), the processing ends. If it does not equal (no), k is incremented by one at step S210 and the processing returns to step S203, from which the processing repeats till the end of the processing.
With the above processing, it is possible to estimate antennas in use being not more than the total number of transmitting antennas. Further, a single OFDM symbol is enough to estimate channel responses among all the antennas.
The description of this embodiment does not describe antenna selection at a transmitting side. However, if a highly reliable antenna is picked up by priority at a transmitting side when an antenna is selected at a transmitting side, the antenna number selection method at the receiving side illustrated in the above flowchart works more effectively.
Though the above embodiment has been described as an example of a method of giving priorities of use to preambles, a method of checking all preambles is also possible.
If preambles are given the priorities of use, a transmitting and a receiving devices must keep the same priority beforehand; or this example cannot be realized. Therefore, a receiving device seeks power from delay profiles for all possible preambles and processes according to determination that data of only preambles whose power exceeds a threshold has been transmitted, so that the device can determine the number of antennas in a similar way.
Next, a radio communication technique according to a third embodiment of the present invention will be described with reference to the drawings. The radio communication technique according to this embodiment is an example in which the present invention is applied to a base station constituting a cell.
As shown in
At the spot B, the mobile station knows a channel response from the base station BS-1 since the mobile station communicates with the base station BS-1. As such, in estimating a channel response from the base station BS-2, the accuracy of estimation of the channel response from the base station BS-2 can be improved by subtracting a preamble B0 component for the base station BS-1, from a received signal. A received signal of the preamble B0 transmitted from the base station BS-1 can be obtained by multiplying channel response information of the base station BS-1 stored in the channel response storage unit 62 by a code of the preamble B0 of the base station BS-1 selected by the code selecting unit (2) 64. The subtracting unit 61 performs an operation on a received signal with the preamble B0, so that it is possible to estimate a channel response from the base station BS-2 without the influence of the preamble B of the base station BS-1. Further, similar operation is repeated to estimate a channel response from the base station BS-3, and the channel response from the base station BS-3 can be estimated without influence of the base stations BS-1 and BS-2.
While channel responses are estimated from the base stations BS-1, BS-2 and BS-3 in that order in the above embodiment, estimation of channel responses cannot be in any order, but the estimation order resulting in higher performance is preferable. The order resulting in higher performance is the order of reliability in communication, by which, a channel response for a desired base station is sought, and then a channel response for a base station with a channel response considered to be highly reliable is sought in a normal state (when connected to a base station). At initial connection, channel responses for all connectable base stations are sought, and then channel responses are sought similarly from higher to lower reliability in communication.
A specific parameter indicating the reliability of communication used herein can simply be the received intensity, or can be SINR (Signal to Interference and Noise Power Ratio).
The above example has been described on the assumption that a mobile station knows a code of an interference wave beforehand. However, in many actual radio communication systems, there is a plurality of codes and a code of an interference wave is unknown. In the following, a flowchart is illustrated showing the flow of processing in such a case. In the processing, it is assumed that a mobile station is connected to a desired base station. A symbol Pinf in the flowchart is a power ratio when a delay profile is estimated using a certain code and is defined by an equation (3):
wherein Ptk means power of the sought delay profile at time tk.
A threshold of Pinf is defined as 10×log Pinf>3 dB in the flowchart, which is only one example, and the threshold is not limited to 3 dB. A symbol tguard indicates a sample point equal to the length of a guard interval.
In
First, steps S1 to S6 will be described. Step S1 is a step for generating a waveform Fpre(k) as a result of performing DFT on a received preamble, wherein k is an integer satisfying 0≦k<N and is a subcarrier number.
Step S2 is a step for multiplying a signal sought at step S1 by Corg*(k) being a complex conjugate of Corg(k) being a desired base station specific code. As a result, the obtained signal is labeled Fdp(k). The Fdp(k) is a channel response from the desired base station containing many interference components.
At step S3, IDFT (the symbol “<=” in
At step S4, Timp(t) is multiplied by Twd(t) being a time window to calculate Trimp(t).
At step S5, DFT (the symbol “<<=” in
Next, processing at steps S11 to S16 will be described. At step S11 following step S6, the channel response from the desired base station is subtracted from a waveform of the first sought preamble (see the equation in the drawing). In the calculation, a subtracted signal component is sought by multiplying Frdp(k) by a code specific to the desired base station. The signal sought at the step is a signal generated only from an interference signal and is defined as Finf(k).
Step S12 is a step for multiplying it by a conjugate of a code in order to seek channel responses for all codes of base station to be possibly interference. Herein, a code specific to a base station is represented by C(k, x) and a complex conjugate of the code is C*(k, x). An indicator x indicates a base station. For example, if there are ten codes, x is numbered from 0 to 9. A channel response is labeled Finfdp(k, x) from each interference base station sought at the step.
Step S13 is a step for performing IDFT on Finfdp(k, x) for every x. As a result, Tinfdp (t, x) is calculated. The Tinfdp (t, x) is a delay profile from each interference base station. At step S15, a power ratio Pinf (x) defined in the equation (3) is computed for every x. At step S16, x's are sought that satisfies Pinf(x)>3 dB and x's are sequenced in descending order of the magnitude of Pinf(x), the sequence is labeled “y”. The number of x's satisfying the above condition equation is labeled “w”. For example, if x is numbered from 0 to 9 as described before, there are three x's satisfying the condition equation, and the order of magnitude is x=3, x=1, x=2, then y equals 3, 1, 2 and w=3.
Finally, steps S21 to S27 will be described. Step S27 following step S16 is a step for determining whether or not a channel response of an interference wave has been estimated for every possible code. Step S21 is a step for decrementing w by one each time of a new loop. Step S22 is a step for seeking a channel response of an interference wave with likely maximum power among the interference waves remaining as interference wave components, in which the processing is same as in steps S2 and S12. As a result, a channel response of an interference wave can be estimated that has maximum power at that time.
Step S23 is a step for performing the same processing as at steps S3 and S13, or IDFT processing. Step S24 is a step for performing the same processing as at step S4, or for applying the time window. At step S25, the same processing is performed as at step S5, to identify a channel response for an interference base station having a code C(k, y). The path is stored at step S31. Step S26 is a step for removing the interference wave component and again execute the loop, the same processing as at step S11.
By executing the loop until w becomes 0, channel responses from all interference base stations can be estimated. Consequently, there is an advantage that data of an interference base station can be previously demodulated at handover, making the handover processing smooth.
Next, a radio communication device according to a fourth embodiment of the present invention will be described with reference to the drawings.
This embodiment features a method of time filtering in the receiving device. The method illustrated herein is not only utilized for the radio communication techniques according to the above first to third embodiments, but also applicable to a general system for calculating a delay profile from a preamble using an antenna specific code to estimate a channel response.
In this embodiment, each base station forms a cell shown in
It is important for a receiving device of a mobile station used in the above system to accurately estimate a channel response. For this purpose, the time filtering function of the receiving function shown in
The time filtering has a problem in which if the time for filtering is shorter, the components decreases that deteriorate with noises and the filtering performance improves, but when a signal component is removed, the accuracy of the estimation degrades. This embodiment involves methods of adaptively performing the time filtering.
With a first method, it is considered that there is a base station with various capabilities in a radio communication system with cell configuration; the base station notifies a mobile station of information of a cell area such as transmitting power; and the mobile station accordingly performs setting for the time filtering. According to the method, in communication with a base station constituting the radio communication system with a maximum cell area, the mobile station sets the time filtering into the guard interval length, shortening time for the time filtering for smaller cell areas.
With a second method, the time filtering length is changed depending on determination by the receiving device. It is important not to lose a signal from a desired base station as possible. For this purpose, a method is provided of broadening time window to a pre-determined level of a power ratio. Power of a k-th sample after IDFT is Ptk, and a ratio of sum of first to m-th powers to total power is Pr, then the equation is obtained:
It is proposed that the window width is broadened till Pr reaches a pre-determined power ratio. In that case, Pr of about 90% would be valid, but any value is not optimal. A method is proposed of setting the value in such a way to reduce influence from other cells and improve the accuracy of channel response estimation. For example, if a window is broadened to the width m satisfying Pr=90%, the estimation of the channel response is more accurate. However, such a window poses the case that m exceeds the guard interval length such as the case of a lot of noises, so a maximum value of m is set to the guard interval length to possibly prevent the accuracy from deteriorating.
It is also possible to approximate noise power at average power outside of the guard interval, subtract the result from a signal component, and then seek a power ratio. In that case, the equation (4) is modified as follows:
wherein tguard<a1<a2, the tguard is a sample point with the same length as the guard interval length. Handling the Pr′ in a similar way to the above example improves the characteristics, whereby a maximum value does not need to be set.
The present invention is applicable to a radio communication system.
Number | Date | Country | Kind |
---|---|---|---|
2005-009407 | Jan 2005 | JP | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/JP2006/000421 | 1/16/2006 | WO | 00 | 7/16/2007 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2006/075732 | 7/20/2006 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
5307376 | Castelain et al. | Apr 1994 | A |
5886987 | Yoshida et al. | Mar 1999 | A |
6359864 | Yoshida et al. | Mar 2002 | B1 |
7406130 | Huang et al. | Jul 2008 | B2 |
7436757 | Wilson et al. | Oct 2008 | B1 |
20010053143 | Li et al. | Dec 2001 | A1 |
20020012333 | Yoshida et al. | Jan 2002 | A1 |
20020122381 | Wu et al. | Sep 2002 | A1 |
20030076908 | Huang et al. | Apr 2003 | A1 |
20040091057 | Yoshida | May 2004 | A1 |
20040105512 | Priotti | Jun 2004 | A1 |
20040116077 | Lee et al. | Jun 2004 | A1 |
20040235421 | Matsuoka et al. | Nov 2004 | A1 |
20060023772 | Mudulodu et al. | Feb 2006 | A1 |
20060120467 | Miyoshi et al. | Jun 2006 | A1 |
20060140303 | Egashira et al. | Jun 2006 | A1 |
20060184862 | Kim et al. | Aug 2006 | A1 |
20060209765 | Ll et al. | Sep 2006 | A1 |
20070014376 | Huang et al. | Jan 2007 | A1 |
20070133393 | Bocquet | Jun 2007 | A1 |
20070171960 | Zhang et al. | Jul 2007 | A1 |
20070195865 | Joetten et al. | Aug 2007 | A1 |
20080317158 | Ketchum et al. | Dec 2008 | A1 |
20090296861 | Sampath | Dec 2009 | A1 |
20120127886 | Kim et al. | May 2012 | A1 |
Number | Date | Country |
---|---|---|
5-75568 | Mar 1993 | JP |
10-101007 | Apr 1998 | JP |
10-190537 | Jul 1998 | JP |
2002-101007 | Apr 2002 | JP |
2002-261670 | Sep 2002 | JP |
2003-18054 | Jan 2003 | JP |
2003-32217 | Jan 2003 | JP |
2003-124854 | Apr 2003 | JP |
2003-283441 | Oct 2003 | JP |
2004-88767 | Mar 2004 | JP |
2004-208254 | Jul 2004 | JP |
Number | Date | Country | |
---|---|---|---|
20080112309 A1 | May 2008 | US |