1. Field of Invention
This invention relates to mobile communications field, and more particularly, to a scheduling method using imperfect CQI (Channel Quality Information) feedback and a scheduler using the scheduling method, especially for WiMAX.
2. Description of Prior Art
IEEE 802.16, also known as WiMAX, has emerged as a strong candidate standard for the future wireless systems, primarily because it offers the potentials for high spectral efficiency, flexible spectrum options, and scalable carrier bandwidth options, as well as the most promising feature, mobility. To achieve above goals, is the WiMAX physical layer (PHY) is based on OFDM (Orthogonal Frequency Division Multiplexing), a scheme that offers good resistance to multi-path and allows WiMAX to operate in NLOS (Non Line Of Sight) conditions. The OFDM technology is also widely recognized as the PHY method for the next generation communications, like 3G LTE (Long Term Evolution). WiMAX's high spectral efficiency is also obtained by using high order modulation and error correction coding scheme under very good signal conditions.
In WiMAX, the whole spectrum is divided into many sub-carriers, and a frame includes many symbols. The combination of carriers and symbols is the radio resource which could be allocated to MS for data transmission. To track the wireless channel, a pilot signal is inserted into sub-carriers every symbol. The channel quality could be estimated by the received pilot signal.
For example, in a typical mobile WiMAX system shown in
In a Cell 1, there are lots of MSs (41, 42 and 43). Some are static users (static user 41), some are moving slowly (pedestrian user 42) while communicating with others through their WiMAX handsets, and some are using handsets on vehicles (vehicular users 43). However, due to multi-path fading and mobility, the wireless channel is not stable with the time. It is changing from time to time. Deep fading can cause one user's data transmission failed. This wastes the wireless bandwidth. To solve the problem, the proportional fairness scheduling algorithm is proposed. The BS scheduler always picks the user who has the best channel quality for data transmission. This is also known as the multi-user diversity. For example, in
Basically, there are two problems that the BS packet scheduler should solve: to increase the spectrum efficiency (also known as the cell throughput) and to guarantee the fairness among multiple MSs.
On the one hand, the goal of the BS packet scheduler is to allocate the radio resource (sub-carriers and symbols in a frame) to an appropriate MS among multiple MSs whose channel conditions are various. For example, when an MS's channel is in a good condition, allocation of resource to such station would gain good spectrum efficiency and a high cell throughput. But if the BS allocates resource to an MS whose channel condition is bad, the spectrum efficiency and cell throughput are low. This problem should be solved by the BS packet scheduler.
On the other hand, if all the resource is allocated to the MSs whose channel conditions are good, the throughput would be very low for the MSs whose channel conditions are bad. In this case, the fairness among multiple MSs is deteriorated. The BS packet scheduler should also handle such a problem.
To solve the problems, a proportional fairness (PF) scheduling algorithm is first proposed in Reference [1]. Then it is adopted as the default working scheduling algorithm in IEEE 802.16m (cf. Reference [2]).
Referring to
The key idea of the PF scheduling algorithm can be described as:
out of all N users will receive transmission opportunity in the next frame (S207).
According to the PF scheduling algorithm, when two users have the same history throughput, the one with higher instantaneous rate (high CQI) would get the transmission opportunity in the next frame, which increases the system throughput (spectrum efficiency). When the two users have the same instantaneous rate, the one with the lower history average rate (throughput) will transmit its data in the next frame, which guarantees the fairness between users. So, a PF scheduler could solve the problem listed above.
Reference [3] provides a PF scheduling method for high speed downlink packet access (HSDPA) system. At the scheduling time point, the BS queries each user's transmission block size to decide the user's current transmission rate Ricurrent in Equation (1).
Reference [4] applies a PF scheduling algorithm into a wireless network with relay stations. The PF scheduling part is the same as the existing solution (Reference [1]).
Reference [5] is the same as the traditional PF scheduling algorithm, except that it uses BLER (Block Error Rate) to correct the scheduling metric. It assumes the measured feedback CQI information is reliable, however, in a high mobility scenario, it is not true. So the unreliable CQI could lead to a bad system performance.
The algorithms in the prior arts all assume that the feedback CQI is reliable. This is true for static users (low mobility) in the most cases, but does not hold for mobile users. When the MS is moving, the channel Doppler spread effect results in the quick channel condition variation. When the channel degrades to a bad condition, the transmitted packet is lost. The Doppler effect is proportional to the speed of the mobile station. This means that at higher moving speeds, the channel changes more quickly, and it is more difficult for the receiver to track it. In such a case, the reported CQI of the fast moving MS is not as reliable as the low mobility MS. So the packet scheduling algorithm based on the reported CQI would cause the system throughput deteriorated.
Although the short CQI report interval could alleviate such a problem in some extent, however, for WiMAX system, the 5 ms report duration makes the Doppler effect obvious. On the other hand, a shorter frame length could not solve the problem when the MS is moving fast.
When the MS is moving, its feedback CQI information is not reliable for the PF scheduler to make an accurate resource allocation scheme. But the history CQI variation could be used to estimate the user's transmission rate. So this invention provides a method to correct the MS's transmission rate based on its history observed CQI. On the other hand, since fast channel variation could lead to packet loss, more radio resource should be allocated to the low channel variation users. So this invention also provides a method to computing a weight index of channel variation based on history CQI. Based on the weight index, a weighted PF scheduler is also provided.
According to a first aspect of the present invention, there is provided a proportional fairness scheduler used in a base station, which is applicable to a high mobility environment, the scheduler comprising:
a correction factor estimator which, for each user equipment in a serving cell of the base station, estimates a correction factor for an instantaneous transmission rate for a next frame according to a statistical result on changes in a channel quality in a current frame;
a metric calculator which, for each user equipment in the serving cell of the base station, calculates a metric based on the instantaneous transmission rate for the user equipment, the estimated correction factor and a history throughput recorded for the user equipment; and
a transmission opportunity granter which grants a transmission opportunity in the next frame to a user equipment having the optimal metric out of all the user equipments.
According to one embodiment of the present invention, the correction factor for the instantaneous transmission rate for the next frame is an estimated symbol good rate, and the statistical result on changes in the channel quality in the current frame is a symbol bad rate in the current frame, and
the correction factor estimator comprises:
Preferably, the modulation coding scheme to be used in the next frame is determined based on the received channel quality information by following a threshold-based mapping method.
Preferably, a modulation coding scheme to threshold mapping table is used to determine the modulation coding scheme to be used in the next frame based on the received channel quality information.
Preferably, a signal-to-noise ratio for the current frame is firstly calculated based on the received channel quality information, and the calculated signal-to-noise ratio is compared with signal-to-noise ratio thresholds in the modulation coding scheme to threshold mapping table to determine the modulation coding scheme to be used in the next frame.
Preferably, the rate calculating unit counts a symbol having a channel quality information value below a threshold corresponding to the determined modulation coding scheme as a bad symbol, and calculates the symbol bad rate in the current frame as a ratio of the number of the bad symbols in the frame to the total number is of the symbols in the current frame.
Preferably, the rate estimating unit estimates the symbol good rate in the next frame based on both the calculated symbol bad rate in the current frame and an estimated symbol bad rate in a previous frame.
Preferably, the rate estimating unit calculates the symbol good rate {circumflex over (α)}i(n+1) in the next frame as:
{circumflex over (α)}i(n+1)=1−{circumflex over (β)}i(n+1)
where {circumflex over (β)}i(n+1) denotes an estimated symbol bad rate in the next frame and is obtained by the rate estimating unit by following:
{circumflex over (β)}i(n+1)=γ*βi(n)+(1−γ)*{circumflex over (β)}i(n−1); nεZ+
where βi(n) denotes the calculated symbol bad rate in the current frame, {circumflex over (β)}i(n−1) denotes the estimated symbol bad rate in a previous frame, γ is a smoothing factor, the subscript i denotes the ith user equipment in the serving cell of the base station, n denotes the sequence number of frame, and the initial value {circumflex over (β)}i(0) when n=1 is set into 0.
Preferably, the metric calculator comprises:
Preferably, the effective rate estimating unit estimates the effective transmission rate Rieffective by following:
Rieffective=Ri*{circumflex over (α)}i(n+1); nεZ+
where the subscript i denotes the ith user equipment in the serving cell of the base station, n denotes the sequence number of frame, Ri denotes the instantaneous transmission rate, and {circumflex over (α)}i(n+1) denotes the estimated symbol good rate.
According to another embodiment of the present invention, the correction factor for the instantaneous transmission rate for the next frame is a weight index, and the statistical result on changes in the channel quality in the current frame is a channel quality information standard variance in the current frame, and
the correction factor estimator comprises:
Preferably, the channel quality information standard variance in the current frame is obtained by statistically analyzing channel quality information for respective symbols in the current frame received from the user equipment.
Preferably, the weight index Wi to be allocated to the user equipment is determined by the weight index allocator by:
where ΔCQIi(n) denotes the channel quality information standard variance in the current frame, the subscript i denotes the ith user equipment in the serving cell of the base station, and n denotes the sequence number of frame.
Preferably, the metric calculator comprises:
Preferably, the effective rate estimating unit estimates the effective transmission rate Rieffective by following:
Rieffective=Ri*Wi
where the subscript i denotes the ith user equipment in the serving cell of the base station, Ri denotes the instantaneous transmission rate, and Wi denotes the weight index to be allocated to the user equipment.
Preferably, the metric calculating unit calculates the metric Mi by following:
where the subscript i denotes the ith user equipment in the serving cell of the base station, Rieffective denotes the estimated effective transmission rate, and Rihistory denotes the history throughput.
Preferably, the transmission opportunity granter grants the transmission opportunity in the next frame to a user equipment having the maximum metric
out of all the user equipments.
Preferably, the instantaneous transmission rate Ri for the user equipment is determined based on a modulation coding scheme to be used in the next frame by following:
According to a second aspect of the present invention, there is provided a base station comprising a receiver which, for each user equipment in the serving cell of the base station, receives channel quality information for respective symbols in the current frame from the user equipment; an adaptive modulation coding unit which, for each user equipment in the serving cell of the base station, determines a modulation coding scheme to be used in the next frame based on the received channel quality information, and thereby determines the instantaneous transmission rate for the user equipment based on the determined modulation coding scheme; and the proportional fairness scheduler according to the first aspect of the present invention.
According to a third aspect of the present invention, there is provided a proportional fairness scheduling method used in a base station, which is applicable to a high mobility environment, the method comprising:
for each user equipment in a serving cell of the base station
According to one embodiment of the present invention, the correction factor for the is instantaneous transmission rate for the next frame is an estimated symbol good rate, and the statistical result on changes in the channel quality in the current frame is a symbol bad rate in the current frame, and
the step of estimating a correction factor for an instantaneous transmission rate for a next frame according to a statistical result on changes in a channel quality in a current frame comprises sub-steps of:
for each user equipment in a serving cell of the base station
Preferably, the modulation coding scheme to be used in the next frame is determined based on the received channel quality information by following a threshold-based mapping method.
Preferably, a modulation coding scheme to threshold mapping table is used to determine the modulation coding scheme to be used in the next frame based on the received channel quality information.
Preferably, a signal-to-noise ratio for the current frame is firstly calculated based on the received channel quality information, and the calculated signal-to-noise ratio is compared with signal-to-noise ratio thresholds in the modulation coding scheme to threshold mapping table to determine the modulation coding scheme to be used in the next frame.
Preferably, a symbol having a channel quality information value below a threshold corresponding to the determined modulation coding scheme is counted as a bad symbol, and the symbol bad rate in the current frame is calculated as a ratio of the number of the bad symbols in the frame to the total number of the symbols in the current frame.
Preferably, the symbol good rate in the next frame is estimated based on both the calculated symbol bad rate in the current frame and an estimated symbol bad rate in a previous frame. Preferably, the symbol good rate {circumflex over (α)}i(n+1) in the next frame is calculated as:
{circumflex over (α)}i(n+1)=1−{circumflex over (β)}i(n+1)
where {circumflex over (β)}i(n+1) denotes an estimated symbol bad rate in the next frame and is obtained by following:
{circumflex over (β)}i(n+1)=γ*βi(n)+(1−γ)*{circumflex over (β)}i(n−1); nεZ+
where βi(n) denotes the calculated symbol bad rate in the current frame, {circumflex over (β)}i(n−1) denotes the estimated symbol bad rate in a previous frame, γ is a smoothing factor, the subscript i denotes the ith user equipment in the serving cell of the base station, n denotes the sequence number of frame, and the initial value {circumflex over (β)}i(0) when n=1 is set into 0.
Preferably, the step calculating a metric based on the instantaneous transmission rate for the user equipment, the estimated correction factor and a history throughput recorded for the user equipment comprises sub-step of:
Preferably, the estimated effective transmission rate Rieffective is obtained by following:
Rieffective=Ri*{circumflex over (α)}i(n+1); nεZ+
where the subscript i denotes the ith user equipment in the serving cell of the base station, n denotes the sequence number of frame, Ri denotes the instantaneous transmission rate, and {circumflex over (α)}i(n+1) denotes the estimated symbol good rate.
According to another embodiment of the present invention, the correction factor for the instantaneous transmission rate for the next frame is a weight index, and the statistical result on changes in the channel quality in the current frame is a channel quality information standard variance in the current frame, and
the step of estimating a correction factor for an instantaneous transmission rate for a next frame according to a statistical result on changes in a channel quality in a current frame comprises sub-steps of:
for each user equipment in a serving cell of the base station
Preferably, the channel quality information standard variance in the current frame is obtained by statistically analyzing channel quality information for respective symbols in the current frame received from the user equipment.
Preferably, the weight index W, to be allocated to the user equipment is determined by:
where ΔCQIi(n) denotes the channel quality information standard variance in the current frame, the subscript i denotes the ith user equipment in the serving cell of the base station, and n denotes the sequence number of frame.
Preferably, the step calculating a metric based on the instantaneous transmission rate for the user equipment, the estimated correction factor and a history throughput recorded for the user equipment comprises sub-step of:
Preferably, the estimated effective transmission rate Rieffective is obtained by following:
Rieffective=Ri*Wi
where the subscript/denotes the ith user equipment in the serving cell of the base station, Ri denotes the instantaneous transmission rate, and Wi denotes the weight index to be allocated to the user equipment.
Preferably, the metric Mi is calculated by following:
where the subscript i denotes the ith user equipment in the serving cell of the base station, Rieffective denotes the estimated effective transmission rate, and Rihistory denotes the history throughput.
Preferably, the transmission opportunity in the next frame is granted to a user equipment having the maximum metric
out of all the user equipments.
Preferably, the instantaneous transmission rate Ri for the user equipment is determined based on a modulation coding scheme to be used in the next frame by following:
According to a fourth aspect of the present invention, there is provided a proportional fairness scheduling method used in a base station, which is applicable to a high mobility environment, the method comprising:
for each user equipment in a serving cell of the base station
Compared with the prior arts, the invention could combat the throughput loss due to the Doppler effect caused by the user's mobility. On the one hand, the invention provides more accurate rate estimation for PF scheduler. It not only allocates resource based on the user's current CQI, but also based on its history CQI variations. This could avoid the situation where the user's CQI is high now, but becomes low in the next frame. The accurate PF decision could increase the system throughput.
On the other hand, the present invention may also compare the history channel variation of each user's CQI. The BS prefers to allocate resources to those mobile stations whose CQI variation is low. This method gives more priority to the low CQI variation users, which could also increase the system throughput.
The above and other objects, features and advantages of the present invention will be clearer from the following detailed description about the non-limited embodiments of the present invention taken in conjunction with the accompanied drawings, in which:
a and 7b are block diagrams showing the PF scheduler 700 according to two embodiments of the present invention.
Hereunder, the present invention will be described in accordance with the drawings. In the following description, some particular embodiments are used for the purpose of description only, which shall not be understood as any limitation to the present invention but the examples thereof. While it may blur the understanding of the present invention, the conventional structure or construction will be omitted.
Referring to
In S201, each MS is required to report its CQI information to the BS in each frame. The CQI can be measured by the pilot signal or data message, and could be transmitted to the BS through a signaling channel or dedicated data channel.
In S203, after receiving each user's 001 information, the link adaptation module could map the CQI information to its current transmission rate. Although various mapping methods are available, here we use a threshold based one. For example a typical mapping table is shown in Table 1.
Based on the measured CQI and the threshold, the appropriate modulation coding scheme (MCS) for the next frame is gotten from Table 1. So the instantaneous transmission rate Ri for the ith MS is
In S404, after the transmission rate is decided in S203, the MCS which the MS will use in the next frame is also decided. BS counts the number of symbols in the last frame whose CQI is below the MCS threshold. On these symbols, the MS's transmission will failed. The ith MS's bad rate βi is defined as:
Since βi may change from one frame to another, here we use a one-order filter to smooth it. We define {circumflex over (β)}i(n+1) is the estimated bad rate in the (n+1)th frame. γ is a smoothing factor,
{circumflex over (β)}i(n+1)=γ*βi(n)+(1−γ)*{circumflex over (β)}i(n−1); nεZ+ (4)
where βi(n) is the ith MS's bad rate in the nth frame, {circumflex over (β)}i(n−1) is the ith MS's estimated bad rate in the (n−1)th frame. That is to say, the estimated bad rate {circumflex over (β)}i(n+1) in the (n+1)th frame (the next frame) is dependent on the ith MS's bad rate βi(n) in the nth frame (the last frame) and the ith MS's estimated bad rate {circumflex over (β)}i(n−1) in the (n−1)th frame (frame immediately before the last frame). The initial value {circumflex over (β)}i(0) (when n=1) is set into 0.
So the estimated effective transmission rate Rieffective for the ith MS in the (n+1)th frame is computed as:
Rieffective=Ri*(1−{circumflex over (β)}i(n+1)) (5)
in Equation (5), an estimated symbol good rate {circumflex over (α)}i(n+1) in the (n+1)th frame can be defined as: {circumflex over (α)}i(n+1)=1−{circumflex over (β)}i(n+1), so that Equation (5) may be modified into:
Rieffective=Ri*{circumflex over (α)}i(n+1) (5′)
As shown in Equation (5′), the estimated symbol good rate {circumflex over (α)}i(n+1) in the (n+1)th frame can be regarded as a correction factor for the instantaneous transmission rate Ri for the (n+1)th frame with respect to the ith MS (Equation (2)).
In S405, the BS uses each MS's corrected transmission rate for PF scheduling. The style of PF scheduling is the same as the prior arts. BS records each MS's history throughput Rihistory, and picks the MS which has the maximum metric
in Equation (6) for data transmission in the next frame (the (n+1)th number frame).
In S407, the user with the highest metric
out of all the users which is granted the transmission opportunity in the next frame will transmit its data in the (n+1)th frame.
Referring to
In S201, each MS is required to report its CQI information to the BS in each frame. The CQI can be measured by the pilot signal or data message, and could be transmitted to the BS through a signaling channel or dedicated data channel.
In S203, after receiving each user's CQI information, the link adaptation module could map the CQI information to its current transmission rate. Although various mapping methods are available, here we use a threshold based one. The mapping table is shown in the above Table 1. Based on the CQI, the appropriate MCS is gotten from Table 1. So the instantaneous transmission rate Ri for the ith MS is
In S504, BS computes each MS's weigh index whose value means how static the user's channel is. The value “1” means the channel is the most static. We define ΔCQIi(n) is the ith MS's CQI standard variance during the last frame (the nth frame). So the weight index Wi for the ith MS is,
So the estimated effective transmission rate Rieffective for the ith MS in the (n+1)th frame is computed as:
Rieffective=Ri*Wi (9)
As shown in Equation (9), the weight index Wi in the (n+1)th frame can be regarded as a correction factor for the instantaneous transmission rate Ri for the (n+1)th frame with respect to the ith MS (Equation (7)).
In S505, the BS uses each MS's corrected transmission rate for PF scheduling. The style of PF scheduling is the same as the prior arts. The BS records each MS's history throughput Rihistory, and picks the MS which has the maximum metric
in Equation (10) for data transmission in the next frame (the (n+1)th frame).
In S507, the user with the highest metric
out of all the users which is granted the transmission opportunity in the next frame will transmit its data in the (n+1)th frame.
In the above two embodiments, the weight index Wi and/or the estimated symbol good rate {circumflex over (α)}i(n+1) in the (n+1)th frame can be regarded as a correction factor for the instantaneous transmission rate Ri for the (n+1)th frame with respect to the ith MS (Equations (2) and (7)), or can be regarded as a correction factor for the metric
Hardware Implementations
As shown in
As discussed above, the AMC module 620 may determine the MCS to be used in the (n+1)th frame based on the received CQI by following a threshold-based mapping method. In this situation, the AMC module 620 uses a MCS mapping table (Table 1) to determine the MCS to be used in the (n+1)th frame based on the received CQI. the AMC module 620 firstly calculates an SNR for the nth frame based on the received CQI, and then compares the calculated SNR with the SNR thresholds in the MCS mapping table (Table 1) to determine the MCS to be used in the (n+1)th frame.
Based on the measured CQI and the threshold, the appropriate MCS for the next frame is gotten from Table 1. So the AMC module 620 may determine the instantaneous transmission rate Ri for the ith MS as
As shown in
out of all the MSs.
In one embodiment (
The rate estimating unit 7104 estimates a symbol good rate {circumflex over (α)}i(n+1) in the (n+1)th frame based on the calculated symbol bad rate βi(n) in the nth frame. Particularly, the rate estimating unit 7104 may calculates the symbol good rate {circumflex over (α)}i(n+1) in the (n+1)th frame as:
{circumflex over (α)}i(n+1)=1−{circumflex over (β)}i(n+1) (13)
where {circumflex over (β)}i(n+1) denotes an estimated symbol bad rate in the (n+1)th frame and is obtained by the rate estimating unit 7104 by following:
{circumflex over (β)}i(n+1)=γ*βi(n)+(1−γ)*{circumflex over (β)}i(n−1); nεZ+ (14)
where βi(n) is the ith MS's bad rate in the nth frame, {circumflex over (β)}i(n−1) is the ith MS's estimated bad rate in the (n−1)th frame. That is to say, the estimated bad rate {circumflex over (β)}i(n+1) in the (n+1)th frame (the next frame) is dependent on the ith MS's bad rate βi(n) in the nth frame (the last frame) and the ith MS's estimated bad rate {circumflex over (β)}i(n−1) in the (n−1)th frame (frame immediately before the last frame). The initial value {circumflex over (β)}i(0) (when n=1) is set into 0.
The metric calculator 720 includes an effective rate estimating unit 7202 for estimating an effective transmission rate Rieffective for the ith MS in the (n+1)th frame as Rieffective=Ri*{circumflex over (α)}i(n+1); nεZ+; and a metric calculating unit 7204 for calculating
the metric Mi as
In another embodiment (
The metric calculator 720′ includes an effective rate estimating unit 7206 for estimating an effective transmission rate Rieffective for the ith MS in the (n+1)th frame as Rieffective=Ri*Wi; and a metric calculating unit 7208 for calculating the metric Mi as
As shown in both
out of all the MSs.
Referring back to
out of all the users which is granted the transmission opportunity in the next frame will transmit its data in the (n+1)th frame.
Although the hardware implementations of the present invention are shown in
Simulation Results
Table 2 gives link level simulation parameters for the performance comparison of the proposed scheduling processes in respective sending SNR settings.
The throughput for “PF without AMC” is the lowest because when the channel is becoming into a bad condition, it is still using the original MCS, which would make the user's packet transmission failed. PF scheduling with the AMC module's help provides a better performance (“PF with AMC”). But when the channel is varying, this scheme is worse than the “PF with estimated rate” algorithm. The improvements lie in that we estimate the user's number of failed symbols in the next frame before scheduling according to the 1st embodiment of the present invention.
Since the channel prediction is hard to be totally accurate, we could provide more resource to the less channel variation users in “Weighted PF” algorithm according the 2nd embodiment of the present invention. So the throughput of the “Weighted PF” scheduler is much higher than the “PF with AMC” and “PF with estimated rate”.
When we provide much more radio resources to the users whose channel variation is low, the fairness between the multiple users will be deteriorated. Here we use the Jain fairness index (cf. Reference [7]) f to investigate the fairness loss due to the proposed PF scheduling processes. f is defined as
where xi is the throughput of the ith user, n is the number of users in the cell. The more the index is close to 1, the fairer the system is.
When SNR is set to 30 dB, the fairness performance is shown in Table 3. From Table 3, it can be seen that all algorithms except the weighted PF have a good Jain fairness performance.
From
The foregoing description gives only the preferred embodiments of the present invention and is not intended to limit the present invention in any way. Thus, any modification, substitution, improvement or like made within the spirit and principle of the present invention should be encompassed by the scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
2008 1 0178836 | Dec 2008 | CN | national |
Number | Name | Date | Kind |
---|---|---|---|
7912490 | Pietraski | Mar 2011 | B2 |
20040203476 | Liu | Oct 2004 | A1 |
Number | Date | Country | |
---|---|---|---|
20100135236 A1 | Jun 2010 | US |