Research Article  Open Access
Pingyan Shi, Xiaohui Liu, Xiaofeng Tao, Jianchao Ji, "Efficient DSUWB MUD Algorithm Using Code Mapping and RVM", Mathematical Problems in Engineering, vol. 2016, Article ID 1930541, 9 pages, 2016. https://doi.org/10.1155/2016/1930541
Efficient DSUWB MUD Algorithm Using Code Mapping and RVM
Abstract
A hybrid multiuser detection (MUD) using code mapping and a wrong code recognition based on relevance vector machine (RVM) for direct sequence ultra wide band (DSUWB) system is developed to cope with the multiple access interference (MAI) and the computational efficiency. A new MAI suppression mechanism is studied in the following steps: firstly, code mapping, an optimal decision function, is constructed and the output candidate code of the matched filter is mapped to a feature space by the function. In the feature space, simulation results show that the error codes caused by MAI and the single user mapped codes can be classified by a threshold which is related to SNR of the receiver. Then, on the base of code mapping, use RVM to distinguish the wrong codes from the right ones and finally correct them. Compared with the traditional MUD approaches, the proposed method can considerably improve the bit error ratio (BER) performance due to its special MAI suppression mechanism. Simulation results also show that the proposed method can approximately achieve the BER performance of optimal multiuser detection (OMD) and the computational complexity approximately equals the matched filter. Moreover, the proposed method is less sensitive to the number of users.
1. Introduction
Ultra wide band (UWB) is an attractive wireless communication technology for its characteristics of high data transmission rate, low power density, high interference resistance, strong multipath resolution, and so on [1–3]. The application range of UWB has been broadened to Wireless Local Area Networks (WLAN), Wireless Sensor Networks (WSN), radar detection, and highspeed communications of indoor and outdoor applications [4–6].
There are mainly two standard schemes of UWB formulated by IEEE 802.15.3a, that is, the multibandbased orthogonal frequency division multiplexing (MBOFDM) and impulseradiobased direct sequence UWB (DSUWB) [7]. Compared with MBOFDM scheme, DSUWB has many attractive advantages such as low peaktoaverage power ratio, significant ability of information hidden, and less sensitivity to multipath fading [8, 9]. In a sense, the multiple access scheme of DSUWB is similar to code division multiple access (CDMA) systems; both of them use pseudorandom (PN) codes to distinguish different users. However, as in conventional CDMA systems, DSUWB systems also suffer from the multiple access interference (MAI). The optimal multiuser detector (OMD) proposed by Verdu [10] could achieve the single user’s BER performance, but it had a very high computational complexity and was too expensive to handle [11]. Therefore, suboptimal detectors which may approximate OMD’s BER performance with an acceptable computational complexity have become a focus of research. In recent years, artificial swarm algorithms are widely used in multiuser detection (MUD). Literature [12] investigates an efficient multiuser detector by selection of initial states based on code mapping for the artificial bee colony algorithm. A complexityperformancebalanced MUD method based on artificial fish swarm algorithm for DSUWB is introduced in [13]; the BER performance of these methods can approximate the performance of OMD while these artificial swarm algorithms need iterated operation. A multiuser frequencydomain turbo detector was employed which combines FD turbo equalization schemes with soft interference cancellation [14]. A codeaided interference suppression method was introduced for narrow band interference restriction in DSUWB systems [15]. Adaptive MUD methods using the recursive least square (RLS) principles were studied in [16, 17]. In [18], a lowcomplexity approximate SISO MUD using soft interference cancellation and linear minimum mean square error (MMSE) filtering for coded CDMA system was introduced. However, few studies were reported that can approximate OMD’s BER performance with a computational complexity approximates to linear MUD.
In this paper, a novel MUD algorithm which combines a code mapping method and a wrong code recognition based on RVM is proposed to achieve a BER performance approximate to OMD method with a very low complexity which almost equals the complexity of matched filter.
The remainder of this paper is organized as follows. In Section 2, the DSUWB system models are introduced. And in Section 3, the principles of the proposed algorithm are described, respectively. In Section 4, simulation results that compare the performance of the proposed algorithm and other approaches are illustrated and discussed, followed by conclusions given in Section 5.
2. System Models
Consider a user asynchronous DSUWB system in additive white Gaussian noise (AWGN) channel, and assume that each user employs the binary phaseshift key (BPSK) modulation. In this paper, Scholtz’s monocycle is used as the UWB pulse waveform, which is approximated to the second derivation of Gaussian pulse. And the expression of Scholtz’s monocycle is [19]where and are the pulse center and the pulse shape parameter, respectively. At the transmitter (), BPSK symbols are spread with the specific PN codes , which are the binary bit stream valued only by −1 or 1, and is the length of bits per packet. The symbol duration is denoted by . We consider that each BPSK symbol can be divided into chips each with duration , where equals . In each chip, a monocycle is transmitted with the duration of to represent the sign of the chip. Practically, the duration of a chip is much longer than the duration of a UWB pulse; that is, . The th user’s transmission signal can be written aswhere represents the random delay of the th transmitter’s monocycle, .
Assume that each transmitter uses a timeinvariant multipath channel in the same band. Let represent the impulse response of the transmission channel. Furthermore, suppose that UWB signals reach the receiving antenna by different paths. The channel impulse response can be written as [20]
The received signal that is transmitted from the th user is given by where denotes convolution operator.
The total received signal can be written aswhere is zeromean additive white Gaussian noise (AWGN) with the unilateral power spectral density of .
In the receiver, the traditional receiver of a DSUWB system consists of a pulse demodulator and a set of matched filters corresponding to each user. Regard signal as the input of the group of matched filters. Furthermore, the inter symbol interference (ISI) can be ignored when the baseband signal transmission rate is much lower than the UWB pulse rate.
Let vector represent the output of the matched filter group, and let vector represent the output of sign detectors, so the output of the matched filters can be represented as follows:where the vector denotes the correct bits of each user and the vector denotes the output of the AWGN from each user’s corresponding matched filter, while the random variable is normally distributed and and , denotes the crosscorrelation matrix, where (what is more, ) and , in which the diagonal element (, ) represents the signal amplitude of the th user.
3. The Proposed Hybrid Multiuser Detection
The main purpose of the proposed method is to pick out the error codes among the received signals and correct them. Firstly, the received codes are mapped into a feature space to make the error codes and the right ones different in some properties which can be easily distinguished. Secondly, identify and pick out the wrong codes by using some approaches in the feature space. Figure 1 shows the diagram of the algorithm. The first stage is corresponding to the code mapping, the second stage is corresponding to features extraction for RVM, and the third stage is corresponding to classification of the codes based on RVM and correcting the wrong codes.
3.1. Code Mapping
Here, a mapping function is constructed to map the output codes set of the matched filter to a onedimensional feature space [21].
According to the theory of optimal multiuser detection described in [22], optimal multiuser detection computes the likehood functionand selects the sequence that minimizes Thus the output of OMD can be expressed as
It can be easily seen that (9) needs trials to get the expected . The computational complexity of OMD is too high. Thus other methods need to be explored to suppress the multiple access interference. Let
According to (9), it can be seen that when the value of function achieves the minimum, the optimal can be obtained. By expanding (10), we get
From (10) and (11), it is clear that function is a quadratic form of vector , which means this is a nonlinear equation with the order of two. And it is not easy to solve the minimum of and the corresponding value of . In addition, the function embraces crosscomponent (, ), and we cannot judge a code independently by . Hence, it is not appropriate for function to be the mapping function.
In this case, it is better to calculate the derivative of function to decrease the order of mapping function. Making the partial derivation of (10), we get
By expanding (12), we get thorder linear equations given as follows:
Because is a discrete equation, we cannot just set the partial derivative described above equal to 0 and get the stationary point of . Let , . Substituting the candidate sequence into (13) (the output of matched filter as shown in (7) or results of other conventional suboptimal multiuser detectors can be employed as the candidate sequence ), there will be two situations as discussed below.
(1) No Wrong Code in . Based on the theory of extreme value, if MAI is the only interference resource without AWGN and the elements in are all correct, the result of (13) strictly equals . And in the condition of high SNR, substituting (6) into , we can see that () when there is no error bit in the candidate codes set .
(2) Wrong Codes Exist in . Suppose that (, ) is the wrong code (in other words, is the correct code); the other codes are correct. Substituting it to the th equation of (13), we getAs for which is not equal to , we getAccording to (14) and (15), we can see when the th user’s code is wrong; then we get , , . Therefore, the function can obviously differentiate the wrong codes and the right codes through the absolute value of it. In addition, is thorder linear equations which can get the result of without complex computations. In conclusion, it is appropriate to set as the optimal decision function in the hybrid MUD algorithm.
It is accessible to map into a onedimension feature space to identify the wrong codes in the candidate set. Figure 2 shows an example of the feature space mapping in the scenario of 10 users, with 6 dB signaltonoise ratio (SNR). of 3200 codes are shown in this figure (candidate is selected from the result of matched filter as in (7)).
According to Figure 2, it is clear that the wrong codes and the right ones have a significant difference in the feature space mapped by .
3.2. Code Classification and Correction
3.2.1. Theory of RVM
Tipping proposed the relevance vector machine (RVM) in [23], a sparse Bayesian modeling approach, which enables sparse classification and regression functions to be obtained by linearly weighting a small number of fixed basis functions from a large dictionary of potential candidates. And a significant advantage to support vector machine is that the kernel function of RVM avoids satisfying Mercer’s condition [24–26].
In classification, the output function is defined bywhere and denotes the weight matrix.
Suppose is a Gaussian conditional probability [27], with the 0 expectation and variance . For two classes’ classification, the likelihood function is defined bywhere denote the target value.
Seeking the maximum posterior probability estimation is equivalent to seeking the mode point of the Gaussian function, namely, .
Due to the fact thatthe maximum posterior probability estimation according to is equivalent to maximizewhere and denotes a constant. Similarly, the marginal likelihood function can be given by
Suppose , the approximation of the Gaussian posterior distribution; that is, , with the variance . The logarithm of the approximate marginal likelihood function is given bywhere .
Reference [28] proposed fast marginal likelihood maximization for sparse Bayesian models, which can reduce the learning time of RVM effectively. To simplify forthcoming expressions, it is defined that
It is showed that (20) has a unique maximum with respect to :
3.2.2. Wrong Codes Recognition Using RVM
We use RVM clustering approach to classify the candidate codes into right bits and wrong ones. Firstly, get the output of the matched filters by (6) and (7), the output of the matched filter group is , and the output of sign detectors is . Secondly compared with the vector and the correct bits vector , we can classify the candidate codes into right set and wrong set ; let and be the training samples to train the RVM. Thirdly two useful feature parameters of the training samples can be extracted below: and . Fourthly use the training samples and to train the RVM. Fifthly use the trained RVM clustering approach to classify the candidate codes into right bits and wrong ones.
4. Results and Discussions
In this section, simulation results are presented in order to show the superiority of the hybrid multiuser detector using code mapping based on relevance vector machine with respect to both system performance improvement and computational complexity reduction.
We present simulation results for the performance of a DSUWB system with four types of multiuser detector: () the matched filter detector; () the MMSE multiuser detector; () the OMD; and () the proposed hybrid multiuser detector. Set the chip transmission rate as 2.55 GHz with the UWB transmission bandwidth 4~8 GHz and use length 255 Kasami spreading codes for the transmitters. Assume that the signal was transmitted in two typical channels: one is AWGN channel without ISI and the other is IEEE 802.15.3a multipath channel and the multipath signal is processed by a selectiveRAKE (SRAKE) receiver with 8 branches and using MRC (Maximal Ratio Combining) rule. Our main objective was to examine a receiver’s BER which directly bases the performance of the receiver on the condition of different SNR environments. In addition, we can check the recognition rate to identify the SNR in this scenario.
4.1. Results of Experiment 1: BER Performance versus SNR
First, consider the AWGN channel environment and IEEE 802.15.3a channel model (without loss of generality, CM2 channel is studied in the simulation experiment). There are four kinds of detectors to be simulated: matched filter detector, MMSE detector, OMD, and the proposed hybrid multiuser detector. The number of transmitters is set to 10, the spread spectrum gain equals 255, and the range of SNR output from the receiver is 0 dB~12 dB. The curves of BER performance versus SNR are depicted in Figure 3.
(a)
(b)
As shown in Figure 3, on the whole, the BER performance of the hybrid multiuser detector is superior to that of the matched filter detector and the MMSE detector in both AWGN channel and the IEEE 802.15.3a multipath CM2 channel. Especially in conditions of high SNR (greater than 8 dB), the detection performance of the hybrid multiuser detector is much better than that of the other two detectors (approximately 1~3 orders of magnitude higher). Hence, what is more important is that the BER performance of the hybrid multiuser detector is very close to that of OMD. In theory, the OMD can completely eliminate the MAI and can arrive at the lower limit of BER in the condition of multiple users. This means that the hybrid multiuser detector can correct almost all the wrong codes caused by MAI. Another experiment was carried out to study the BER performance of the proposed hybrid MUD method in different IEEE 802.15.3a channel models, which is presented in Figure 4. Seen from Figure 4, the proposed method is sensitive to channel environments; the performance trend worsened with the deterioration of the channel environment.
4.2. Results of Experiment 2: User Capacity of System
Furthermore, consider another situation with different number of transmitters under the same SNR condition. To make the comparison apparent, the algorithms of matched filter detector MMSE and hybrid multiuser detector are simulated. The simulation results are shown in Figure 5 with the SNR of 4 dB, 6 dB, 8 dB, and 10 dB; the other simulation environments yield the same results as BER performance.
(a)
(b)
(c)
(d)
Based on the results in Figure 5, it can be seen that the BER performance of the hybrid multiuser detector remains nearly unchanged with an increase in the number of transmitters compared with that of the matched filter receiver. This shows that the hybrid multiuser detection can suppress multiaccess interference. Therefore, under the stated specific SNR conditions, the hybrid multiuser detector is still able to give a superior detection performance when there are many transmitters occupying the same channel simultaneously, which indicates that the BER performance of the hybrid multiuser detector is not sensitive to the channel capacity. This is one of the great advantages of this type of detector.
4.3. Results of Experiment 3: The NearFar Effect (NFE) Resistant Capability
Due to the spatial distribution of terminals, NFE is one of the outcomes in Multiple Address (MA) systems, which is discussed in this simulation in both AWGN channel and IEEE 802.15.3a CM2 channel. In this experiment, the number of users was set as 10; the first user received signal with the stationary SNR of 6 dB, while in that of other users, SNR varies from 0 to 12 dB synchronously. And other parameters are the same as those in the former experiments.
The nearfar effect resistant abilities of different multiuser detector in both AWGN channel and IEEE 802.15.3a CM2 channel are shown in Figure 6. It is evident that the matched filter detector has the worst NFE resistance ability, and both the proposed hybrid multiuser detector and the MMSE have the stronger NFE resistant ability compared to matched filter detector. However, considering the BER performance of them, hybrid MUD is the best.
(a)
(b)
4.4. Results of Experiment 4: Computational Complexity
OMD can obtain the best BER performance. However, as is shown in (8), in OMD the value of () needs computations to find the maximum one ( is the number of users). The multiuser proposed in this paper computes the mapping function and employs RVM to distinguish the wrong codes and then correct them. Let the execution time of matched filter normalize to 1. In an environment of 10user AWGN channel, Table 1 lists the relative execution time of matched filter, MMSE, OMD, and the proposed hybrid MUD.

From Table 1, it can be seen that the computational complexity of the proposed hybrid multiuser detector is much lower than that of OMD. And it is in the same order of magnitude to matched filter and MMSE. Thus the proposed multiuser detector based on code mapping and RVM can get good BER performance with computational complexity.
5. Conclusions
In this paper, a new MAI suppression mechanism is proposed for multiuser detection of DSUWB system. The algorithm embraces three steps. First, as the output of matched filters, the candidate codes set is mapped into a onedimensional feature space by the optimal decision function. The purpose is to make the right codes and wrong codes easy to be distinguished. Second, calculate the statistical characteristics of the codes in the feature space. Third, use RVM classifier to classify the receiver multiuser function mapping characteristic value, and then correct the wrong codes by sign reversal. The simulation results show that the BER performance and the nearfar effect resistant abilities of the new MAI suppression mechanism are superior to those of the matched filter detector and the MMSE detector in both AWGN channel and IEEE 802.15.3a channel, especially when the output SNR is greater than 8 dB. Its performance is very close to that of OMD. Compared to the OMD, the new MAI suppression mechanism has lower computational complexity that increases linearly with the number of users. Hence, we conclude that the new MAI suppression mechanism can give an admirable detection performance under high SNR and large number of users, which is very useful for practical applications. Furthermore, the performance of the new MAI suppression mechanism can hardly be influenced by changes in numbers of users, which means that this hybrid approach for multiuser detection will be an efficient method to solve interference problems when there are large quantities of transmitters occupying the same channel simultaneously.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The research in this paper is supported by the National Natural Science Foundation of China (Grant no. 61102084), the Fundamental Research Funds for the Central Universities in China (Grant no. HIT.NSRIF.2010092), and the China Postdoctoral Science Foundation (Grant no. 2011M500665).
References
 S. J. Ahmadi and R. C. Qiu, “Ultrawideband (UWB) communications channel measurements—a tutorial review,” International Journal of Ultra Wideband Communications and Systems, vol. 1, no. 1, pp. 11–30, 2009. View at: Google Scholar
 X. Shen, M. Guizani, H. H. Chen, R. C. Qiu, A. F. Molisch, and L. B. Milstein, “Guest editorial ultrawideband wireless communications—theory and applications,” IEEE Journal on Selected Areas in Communications, vol. 24, no. 4, pp. 713–716, 2006. View at: Publisher Site  Google Scholar
 H. Zhang, S. Hong, and K. Chang, “Class 4 active RFID multihop relay system based on IEEE 802.15.4a lowrate UWB in sensor network,” KSII Transactions on Internet and Information Systems, vol. 4, no. 3, pp. 258–272, 2010. View at: Publisher Site  Google Scholar
 Z. L. Xiong and Q. S. Xiang, “A novel multiuser detector in UWB communication systems,” in Proceedings of the Sensor Array and Multichannel Signal Processing Workshop Proceedings, pp. 129–132, IEEE, Barcelona, Spain, July 2004. View at: Publisher Site  Google Scholar
 J. Chóliz, Á. Hernández, and A. Valdovinos, “A framework for UWBbased communication and location tracking systems for wireless sensor networks,” Sensors, vol. 11, no. 9, pp. 9045–9068, 2011. View at: Publisher Site  Google Scholar
 M. H. Kabir and R. Kohno, “A hybrid TOAfingerprinting based localization of mobile nodes using UWB signaling for non lineofsight conditions,” Sensors, vol. 12, no. 8, pp. 11187–11204, 2012. View at: Publisher Site  Google Scholar
 J. Maunu, T. Koivisto, M. Laiho, and A. Paasio, “An analog viterbi decoder array for DSUWB receiver,” in Proceedings of the IEEE International Conference on UltraWideband (ICUWB '06), pp. 31–36, Waltham, Mass, USA, September 2006. View at: Publisher Site  Google Scholar
 S.S. Tan, A. Nallanathan, and B. Kannan, “Performance of DSUWB multipleaccess systems with diversity reception in dense multipath environments,” IEEE Transactions on Vehicular Technology, vol. 55, no. 4, pp. 1269–1280, 2006. View at: Publisher Site  Google Scholar
 H. Sato and T. Ohtsuki, “Frequency domain channel estimation and equalization for direct sequence ultra wideband (DSUWB) system,” IEE Proceedings: Communications, vol. 153, no. 1, pp. 93–98, 2006. View at: Publisher Site  Google Scholar
 S. Verdu, “Optimum multiuser asymptotic efficiency,” IEEE Transactions on Communications, vol. 34, no. 9, pp. 890–897, 1986. View at: Publisher Site  Google Scholar  MathSciNet
 Y. C. Yoon and R. Kohno, “Optimum multiuser detection in ultrawideband (UWB) multipleaccess communication systems,” in Proceedings of the International Conference on Communications (ICC '02), pp. 812–816, New York, NY, USA, May 2002. View at: Google Scholar
 Z. Yin, X. Liu, and Z. Wu, “A multiuser detector based on artificial bee colony algorithm for dsuwb systems,” The Scientific World Journal, vol. 2013, Article ID 547656, 8 pages, 2013. View at: Publisher Site  Google Scholar
 Z. Yin, Z. Zong, H. Sun, Z. Wu, and Z. Yang, “A complexityperformancebalanced multiuser detector based on artificial fish swarm algorithm for DSUWB systems in the AWGN and multipath environments,” EURASIP Journal on Advances in Signal Processing, vol. 2012, no. 1, article 229, 13 pages, 2012. View at: Publisher Site  Google Scholar
 P. Kaligineedi and V. K. Bhargava, “Frequencydomain turbo equalization and multiuser detection for DSUWB systems,” IEEE Transactions on Wireless Communications, vol. 7, no. 9, pp. 3280–3284, 2008. View at: Publisher Site  Google Scholar
 C. Wang, M. Ma, R. Ying, and Y. Yang, “Narrowband interference mitigation in DSUWB systems,” IEEE Signal Processing Letters, vol. 17, no. 5, pp. 429–432, 2010. View at: Publisher Site  Google Scholar
 G. S. Biradar, S. N. Merchant, and U. B. Desai, “Performance of constrained blind adaptive DSCDMA UWB multiuser detector in multipath channel with narrowband interference,” in Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM '08), pp. 1–5, IEEE, New Orleans, La, USA, December 2008. View at: Publisher Site  Google Scholar
 Q. Z. Ahmed and L.L. Yang, “Reducedrank adaptive multiuser detection in hybrid directsequence timehopping ultrawide bandwidth systems,” IEEE Transactions on Wireless Communications, vol. 9, no. 1, pp. 156–167, 2010. View at: Publisher Site  Google Scholar
 X. Wang and H. Vincent Poor, “Iterative (Turbo) soft interference cancellation and decoding for coded CDMA,” IEEE Transactions on Communications, vol. 47, no. 7, pp. 1046–1061, 1999. View at: Publisher Site  Google Scholar
 R. M. Cramer, M. Z. Win, and R. A. Scholtz, “Impulse radio multipath characteristics and diversity reception,” in Proceedings of IEEE International Conference on Communications (ICC '98), pp. 1650–1654, Atlanta, Ga, USA, June 1998. View at: Publisher Site  Google Scholar
 C. Thejaswi, P. S. Manohar, S. V. Ganapathy, R. Patro, and M. Raina, “Simple multiuser detectors for DSUWB systems,” in Proceedings of the 17th Annual IEEE International Symposium on Personal, Indoor and Mobile Radio Communications, pp. 1–5, IEEE, Helsinki, Finland, September 2006. View at: Publisher Site  Google Scholar
 Z. Yin, Y. Kuang, Z. Wu, and W. Tang, “A hybrid multiuser detector for DSUWB systems with matched filter and errorbit recognizer,” Journal of the Chinese Institute of Engineers, vol. 35, no. 8, pp. 1003–1012, 2012. View at: Publisher Site  Google Scholar
 S. Verdu, “Minimum probability of error for asynchronous Gaussian multipleaccess channels,” IEEE Transactions on Information Theory, vol. 32, no. 1, pp. 85–96, 1986. View at: Publisher Site  Google Scholar  MathSciNet
 M. E. Tipping, “Sparse Bayesian learning and the relevance vector machine,” Journal of Machine Learning Research, vol. 1, no. 3, pp. 211–244, 2001. View at: Publisher Site  Google Scholar  MathSciNet
 M. E. Tipping, “Fast marginal likelihood maximisation for sparse bayesian models,” in Proceedings of the 9th International Workshop on Artificial Intelligence and Statistics, Key West, Fla, USA, January 2003. View at: Google Scholar
 P.K. Wong, Q. Xu, C.M. Vong, and H.C. Wong, “Ratedependent hysteresis modeling and control of a piezostage using online support vector machine and relevance vector machine,” IEEE Transactions on Industrial Electronics, vol. 59, no. 4, pp. 1988–2001, 2012. View at: Publisher Site  Google Scholar
 Q. Xu and P.K. Wong, “Hysteresis modeling and compensation of a piezostage using least squares support vector machines,” Mechatronics, vol. 21, no. 7, pp. 1239–1251, 2011. View at: Publisher Site  Google Scholar
 Z. Yang, Z. Wu, Z. Yin, T. Quan, and H. Sun, “Hybrid radar emitter recognition based on rough kmeans classifier and relevance vector machine,” Sensors, vol. 13, no. 1, pp. 848–864, 2013. View at: Publisher Site  Google Scholar
 E. E. Azzouz and A. K. Nandi, “Automatic identification of digital modulation types,” Signal Processing, vol. 47, no. 1, pp. 55–69, 1995. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Pingyan Shi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.