Modulation Classification Based on Extensible Neural Networks

Qing Yang, Guan

doi:https://doi.org/10.1155/2017/6416019

Mathematical Problems in Engineering

On this page

Abstract Introduction System Model Analysis Conclusion Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2017 | Article ID 6416019 | https://doi.org/10.1155/2017/6416019

Modulation Classification Based on Extensible Neural Networks

Guan Qing Yang¹

Academic Editor: Paolo Crippa

Received13 Apr 2017

Revised03 Aug 2017

Accepted29 Aug 2017

Published09 Oct 2017

Abstract

A deep learning architecture based on Extensible Neural Networks is proposed for modulation classification in multipath fading channel. Expanded Neural Networks (ENN) are established based on energy natural logarithm model. The model is set up using hidden layers. Modulation classification based on ENN is implemented through the amplitude, phase, and frequency hidden network, respectively. In order to improve Probability of Correct classification (PCC), one or more communication signal features are extracted using hidden networks. Through theoretical proof, ENN learning network is demonstrated to be effective in improving PCC using amplitude, phase, and the frequency weight subnetwork, respectively. Compared with the traditional algorithms, the simulation results show that the proposed ENN has higher PCC than traditional algorithm for modulation classification within the same training sequence and Signal to Noise Ratio (SNR).

1. Introduction

In the communication transmission technology, modulation classification for communication signals is an important research direction, especially in the field of space communication, satellite communication, underwater communication, and so forth. Modulation classification is constrained by the transmission environmental parameters, such as signal amplitude, signal frequency, and signal phase, especially fading channel.

In the related research literature, there are two kinds of research trends for communication signal modulation classification. One trend is based on the Maximum Likelihood function method based on feature parameter estimation [1]. The other trend is signal statistical features [2]. In order to improve PCC for modulation classification, [3–6] established different Maximum Likelihood (ML) functions. This method was optimized to estimate parameters, such as modulation mode, carrier frequency, phase, and channel response. In [3], a generalized ML function was established. To solve this likelihood function, the QAM signal was verified to obtain lower complexity. In [4], a hybrid ML function was proposed. Instantaneous signal feature was obtained to reduce the complexity. In [5, 6], the optimization method based on different Maximum Likelihood was proposed for MFSK signal modulation. Thus a higher PCC had been obtained. In [7], a hybrid ML algorithm was proposed to classify AMC communication signals.

ML function for modulation classification is generally assisted with Bayesian framework. When the numbers of unknown parameters are too large, the complexity for established likelihood functions would be increasing. Different ML functions and Bayesian decision networks must be evaluated carefully. According to different modulation types, literature [3–7] proposed different likelihood ratio functions, which could obtain the theoretical optimal result under the Gaussian white noise condition. However, multipath fading channel would seriously impact on modulation classification.

Another common method is signal statistical features, which is an easy method to achieve suboptimal result by obtaining statistical characteristics features, such as high-order cumulants, wavelet features, and multidimensional signal features. In [8], an instantaneous signal statistical characteristic was proposed, which included instantaneous amplitude, instantaneous phase, and instantaneous carrier frequency. Reference [9] proposed wavelet feature, which was decomposed into wavelet coefficients by Haar transform. Haar transformed wavelet coefficients have better antinoise features than the higher order cumulants.

Literature [10, 11] acquired high-order cumulants features. However, higher computation improved the complexity. In [12], HOS statistics feature was obtained, which included time accumulation statistics feature. HOS statistics features had excellent anti-Gaussian white noise. Similarly, the literature [13] adopted HOS features to combat multipath fading. In [14], a statistic classifier within a multipath channel had been proposed to estimate signal parameters. Within multipath fading channel, literature [15] proposed a blind channel estimation method based on higher order statistics. However, this blind channel estimation is only effective in fading channels within two-path channel. In [14], within multipath fading channel, the signal parameters were firstly estimated. In [16], modulation classification in multiuser system was researched. Most of modulation classification algorithms are within Gaussian white noise channel. In fact, due to multipath fading channel, the application for modulation classification algorithms was greatly reduced.

In the past few years, due to rapid development of artificial intelligence, researchers introduced machine learning algorithms into modulation classification and gave the relevant theoretical derivation proof. In [17], the -nearest neighbor clustering classifier was used for AMC modulation classification. Under ideal conditions, better PCC based on -nearest neighbor clustering classifier was only limited to ideal conditions.

In [18], a multilayer neural network was used for the modulation classification. On this basis, literature [19] used ANN neural network for modulation classification. However, the ANN network was only limited to the shallow network. Literature [18, 20] implemented modulation classification combined with GA algorithm. On this basis, [17] proposed a KNN algorithm, which was used to improve PCC.

In this paper, a deep learning architecture based on Extensible Neural Networks is proposed for modulation classification in multipath fading channel. Expanded Neural Networks are established through the energy natural logarithm model, which works through the amplitude, phase, and frequency subnetwork, respectively, in order to improve the PCC.

2. System Model and Problem Formulation

The Gaussian white noise is according to the Gaussian distribution with mean 0 and variance , which represents the mean and variance using A/D sampling additive Gaussian white noise (AWGN) signal. Specifically, Gaussian white noise signal is passing through correlated Receive or Matched filter and forming discrete random variable. So we could obtain that the mean is 0 and variance is . The block diagram end-to-end signal flow chart of the system for ENN is shown in Figure 1.

The received signal passes through ideal matched filter. After A/D sampling, intersymbol interference is cancelled by the equalization. Similarly, the intersymbol interference is also eliminated by ideal symbol equalizer.

For discussion convenience, suppose that the symbol period is less than the channel coherence time. The system adopts ideal channel equalizer to reduce the intersymbol interference.

The received signal at the tth time gap can be expressed as

is defined as the transmitting signal data set for modulation mode , , is defined as amplitude gain, is defined as carrier offset, is defined as phase deviation, and is defined as Gauss white noise for power .

The received signal is the input for ENN; suppose that is the output estimator for the ENN, is the output value estimated for the lth order network, and is defined as the ideal output value for the lth order network.

Here, a nonlinear logarithmic error energy function has been proposed as the objective function. Suppose that the signal is in the form of polar coordinates, including the amplitude features, phase features, and frequency features of the signal. Define the system’s energy loss as :

Similarly, the th order network output is for estimating the output value, which can be expressed as the exponential energy loss function,

could be defined as the weight between the kth hidden layer network and the lth output network. The weight amplitude can be expressed as

The corresponding weight for the phase and frequency can be expressed, respectively, as , . can be defined as a hidden layer activation function; the hidden layer needs to be obtained by the following activation function, which can be expressed as

Here, the input for ENN is the sequence , which is the received multipath fading signal defined as . Thus, the kth hidden layer unit can be expressed as

is defined as the error attenuation factor for the kth neural network, which could be expressed as . is the kth hidden layer average center processing unit, which can be expressed as .

In this paper, the hidden layer needs to be obtained by the following activation function, which can be expressed as . The activation function needs to be differential everywhere. The condition is as below: , while satisfying and .

Currently, the available function is considered as activation function, including the sigmod function and the tanh function. As learning network construction, each layer requires faster convergence. So the ENN select tanh function as an activation function.

Define , , as the system network training library, in which is defined as the input of the th training unit. is defined as the target output of the nth training unit. Define signal error as , after th layer network training, which can be defined asDefine amplitude error as , which can be defined asDefine phase error as , which can be defined asDefine phase error as , which can be defined asThe framework for deep learning ENN is shown as in Figure 2.

The output feature for the th training sample is defined as . It can be expressed in polar coordinates to obtain amplitude, phase, and frequency information based on energy function model. The signal can be expressed as

is defined as signal amplitude, which also can be expressed as , is defined as signal phase, and is defined as frequency. Substituting (11) into (2), we could obtain that

For training sequences, the target energy function can be defined aswhere is the weight of the network training system. In order to obtain the optimal training objectives, the objective function is minimized, which can be expressed asSimilarly, the connection weight can be expressed asIt could be obtained asExpanding formula (16), we could obtainSubstituting (2) into (17), which is available, we could obtainWe could obtain the gradient of (18) for , and we could obtainAnd then we simplified (19) asThe received multipath fading signal enters the different training network through amplitude, phase, and frequency subnetwork, respectively. Among them, the amplitude weight network can be expressed asAmong them, the phase weight network can be expressed asAmong them, the frequency weight network can be expressed asThe above formulas could be represented by matrix, which can be simplified aswhereIn accordance with different training network, the matrix form for amplitude can be expressed asThe matrix form for phase can be expressed asThe matrix form for frequency can be expressed asAccording to formulas (24)–(28), we could obtain the learning network as shown in Figure 3.

3. ENN Algorithm

According to Figure 2, the network of ENN algorithm can be implemented as follows.

Algorithm 1 (network learning based on amplitude information). Given a training set, , one has the following.
Select the layer hidden layer network for training, in which the hidden layer parameters are .
Calculating the th hidden layer response, ,Calculate the front matrix ,Calculate the posterior matrix ,Calculate the network weight matrix,

Algorithm 2 (network learning based on phase information). Given a training set, , one has the following.
Select the layer hidden layer network for training, in which the hidden layer parameters are .
Calculating the th hidden layer response, ,Calculate the front matrix ,Calculate the posterior matrix ,Calculate the network weight matrix,

Algorithm 3 (network learning based on frequency information). Given a training set, , one has the following.
Select the layer hidden layer network for training, in which the hidden layer parameters are .
Calculating the th hidden layer response, ,Calculate the front matrix ,Calculate the posterior matrix ,Calculate the network weight matrix,

4. Experimental Classification Results and Analysis

The network of deep learning ENN is established for modulation classification. The modulation consists of BPSK, QPSK, 8PSK, and QAM. The signal modulation mode is generated by a random generator. Signal carrier frequency is defined as KHz. Sampling rate is defined as KHz.

The experiment is simulated by Monte Carlo method. SNR ranges from −5 dB to 20 dB. Suppose that 512 signal samples are constructed into one data frame, and there are 1000 sets. The first 100 sets are used as training, and 900 sets left are used for testing.

Channel model has been used corresponding to a typical multipath model shown in Table 1, which is modeled as Rayleigh fading channels. According to the three-diameter model, the main channel is adopted as the Rice channel model, and the second-diameter model is adopted as the Rayleigh channel model.

In order to simplify the simulation, we set ENN learning network using three subnetworks, each subnetwork including 30 layers. The number of hidden nodes is defined as 100. In this paper, we use the HML algorithm [7], comparison with the ANN algorithm based on the constant modulus statistics algorithm [14], and literature [19]. We demonstrate the performance based on the PCC under different SNR, different training length, and the convergence performance.

4.1. Classification for Different Modulation Modes with ENN Deep Learning Network

Define that the SNR varies from 5 dB to 15 dB, and each of the 512 simulation signal samples forms a data frame. Figure 4 shows the PCC based on BPSK, QPSK, 8PSK, and 16QAM at different SNR. When the SNR is greater than 0 dB, the PCC is above 93%. When the SNR is greater than 10 dB, PCC for all modulation methods is more than 80%.

Figure 5 shows the PCC based on BPSK, QPSK, 8PSK, and 16QAM modulation signals at different SNR. Define different modulation modes occurring at equal probability. Suppose that SNR ranges from −5 dB to 20 dB, and 1024 signal sampling points constitute a data frame. The number of data frames is set as 1000. Different algorithms are carried out according to different training lengths. Define the length of the training symbols, respectively, , 20, 30. Enough training lengths can effectively improve the PCC. Under the same PCC, the distance for different training lengths can be nearly 3 dB.

4.2. PCC Simulation Based on Different Signal to Noise Ratio

Figure 6 shows the PCC performance of ENN, and algorithms in literature [7, 14, 19] for comparison at SNR ranged from the −5 dB to 20 dB. Suppose that all the algorithms have adopted the same training sequence, ; when the SNR = 5 dB, the PCC is below 90%. The performance of the algorithm is greatly affected by noise. When SNR = 10 dB, the PCC is increasing. In the multipath fading channel, with the increasing of SNR, the overall performance of different algorithms is closer to the ideal value. The proposed ENN algorithm uses hierarchical identification, and the overall PCC is higher.

HML algorithm has established the likelihood function model in the Gaussian white noise, but the PCC has greatly impacted on the conditions of fading channel. The cumulants of amplitude are affected by the fading multipath channel, especially under low SNR condition. The proposed deep learning ENN algorithm improves the fitness for the amplitude in the fading channel, which is implemented by the learning of the amplitude weight subnetwork, especially within the slow fading channel.

In addition, PCC of HML increases rapidly with SNR increasing. HML has better classification ability. However, HML algorithm performance is constrained by multipath fading channels. Thus the performance has not shown advantages. It can be seen from the simulation results that the proposed algorithm has high accuracy and stability under different SNR. Compared with the traditional neural network ANN algorithm, the proposed algorithm has higher PCC and reliability. This is because the deep learning ENN algorithm is based on signal features, which is based on different dimensions of subnetwork identification. The overall performance has shown more advantages. In the SNR at more than 5 dB, the PCC can reach more than 90%.

4.3. Probability of Correct Classification Simulation Based on Different Lengths of Training Sequences

In order to demonstrate the performance of different modulation classification algorithms with different training samples, suppose that the 1024 signal sampling points constitute a data frame, and the total number for classified data frames is identified as 1000 according to the different training length.

Figure 7 shows the PCC performance under the condition of training samples from to at SNR = 5 dB. We can obtain from the simulation figure that, as the training number increases, the PCC is near optimal. When the number of samples is , the PCC is higher than 90%. When the training number is small, the neural network ANN structure weight for training is not sufficient, so the performance is poor. The cumulants algorithm has low performance, mainly because of the feature extraction based on amplitude feature, which reduces the performance in multipath channel.

The ML is based on the likelihood criterion and makes full use of the estimation about unknown parameters. Therefore, under the same training number condition, the performance is higher than ANN algorithm. The proposed ENN algorithm has higher accuracy and stability, with high reliability and better performance advantages. In the case of small sample group, the correct classification rate can reach more than 90%. Thus, the proposed ENN also has excellent classification performance with small sample. The PCC of the signal modulation mode is higher than the literature algorithm [7, 14, 19].

4.4. Convergence Performance Simulation Based on Different Algorithm

Figure 8 shows the convergence performance of the proposed deep learning network ENN and ANN classifier algorithm proposed in [19] within different training sequences. With the length of the training sequence increasing, the ENN algorithm has faster convergence performance. This is mainly because the ENN algorithm has faster convergence using multidimensional features and quickly enters fast convergence direction. Therefore, the convergence performance of the ENN algorithm can be improved rapidly. In addition, when the number of training sequences is small, ENN classifier algorithm performance is better than ANN algorithm according to PCC. So the proposed ENN algorithm has better adaptability with small training sequences. Compared with the existing literature [19], the ENN algorithm proposed in this paper has better classification performance and has a better training sequence dynamic range.

4.5. Probability of Correct Classification Simulation Based on MIMO-OFDM

In order to analyze the modulation classification for OFDM, the communication system based on OFDM modulation is established. In the simulation, establishing a 2-channel channel, the main path has a direct component. The second path has a fading rate of 0.4 with a maximum delay 180 ns. Amplitude fading is supposed as Rician distribution. The link transmission bit rate is 40 Mbit/s. The frame length is 16384 bits.

Suppose that the cyclic prefix length of OFDM is larger than the maximum spreading delay; the guard interval should be 250 ns multiplied by 4, which is 1 us. Then the symbol period of OFDM can be 6 us. So each OFDM symbol needs to be transmitted with (40 Mbit/s)/(1/6 us) = 240 bits. OFDM symbols within each subcarrier can be transmitted with 2 bits, which should need 240/2 = 120 subcarriers to meet with transmission rate requirements. The zero-padded subcarriers should be added to eight digit number to facilitate the realization of 128-point FFT/IFFT. For the OFDM modulation, the normalized frequency offset factor is used to measure frequency offset. The normalized frequency offset factor is defined as the ratio of the carrier frequency deviation to the subcarrier spacing. Normally, if the normalized frequency offset is an integral, the orthogonality for subcarrier will not be destroyed.

The orthogonality for subcarrier will be destroyed. When the frequency deviation of the subcarrier is 0.5 times for the carrier interval, the orthogonality of the subcarrier will be minimized, resulting in serious loss of demodulation signal performance. The OFDM system can obtain a carrier frequency interval of approximately 0.5 times the frequency deviation, it will seriously damage the subcarrier orthogonality, and the demodulation of the signal performance has a great impact.

The paper is mainly for the normalized frequency offset factor within 0.5 times. For the deviation analysis, the normalized frequency offset factor is set as , , , .

It can be obtained from the simulation in Figure 9, when the normalized frequency offset factor , which is lower, the PCC for ENN is lower due to the different normalized frequency offsets. This is because the proposed ENN has correction capability for the small frequency deviation. When the normalized frequency offset , the frequency offset has introduced a larger orthogonal subcarrier failure. For the signal modulation classification, the error has exceeded the relevant decision limits, and PCC is low.

It can be obtained from the simulation in Figure 9, when the normalized frequency offset factor , which is lower, the PCC for ENN is lower due to the different normalized frequency offsets. This is because the proposed ENN has some error correction capability for the small frequency deviation. When the normalized frequency offset , the frequency offset has introduced a larger orthogonal subcarrier failure. For the signal modulation classification, the error has exceeded the relevant decision limits, and the PCC is low.

4.6. Probability of Correct Classification Simulation Based on MIMO

Based on OFDM modulation mode link communication, the antenna model is extended from SISO to MIMO. MIMO adopts and antenna group, respectively, to establish STBC coding mode. The PCC for OFDM-MIMO has been analyzed. In the simulation, a 2-channel channel is used. The main path has a direct component, and the amplitude fading follows the Rician distribution. The fading amplitude for the second diameter is 0.4 for the first diameter, the maximum delay is 180 ns, the channel bandwidth is 10 MHz, and the frame length is 128-symbol cycle.

From the simulation in Figure 10, we can obtain that MIMO-STBC system can combat with a certain channel fading characteristic. The ENN with MIMO-STBC can be used to improve the Probability of Correct classification. Therefore, combined with the STBC coding of MIMO system, the proposed ENN depth modulation pattern recognition network can obtain high PCC.

4.7. Complexity Analysis of Different Algorithms

The proposed ENN deep learning network performs the best performance because of sacrificing the complexity, and the multiple identifiers are trained by different signal features including amplitude, phase, and frequency and are judged by the combined features. So the PCC results are available in both small sample and low SNR conditions.

Viewing the classification time with the same experimental environment, the constant modulus statistics algorithm has the faster velocity because of only extracting the amplitude features information. Followed by the ANN network, ANN network needs to establish the neural networks through the amplitude characteristics.

5. Conclusion

In this paper, an adaptive Extensible Neural Network has been proposed for modulation classification in the multipath fading channel. By establishing the natural logarithmic energy model of the desired model, the deep learning ENN has been divided into the amplitude, phase, and frequency weight subnetwork. Compared with the traditional algorithm, the simulation results show that the proposed ENN algorithm has efficient classification performance and small sample classification ability and could be adaptive in multipath fading channels with the same length training sequence and SNR.

Conflicts of Interest

The author declares that they have no conflicts of interest.

Acknowledgments

The work reported here is supported by National Natural Science Foundation of China (no. 61501306), The Natural Science Foundation of Liaoning Province of China (no. 2015020026), and The Natural Science Foundation of Liaoning Provincial Department of Education (no. L2015402).

References

J. L. Xu, W. Su, and M. Zhou, “Likelihood-ratio approaches to automatic modulation classification,” IEEE Transactions on Systems, Man and Cybernetics Part C: Applications and Reviews, vol. 41, no. 4, pp. 455–469, 2011.
View at: Publisher Site | Google Scholar
W. Su, “Feature space analysis of modulation classification using very high-order statistics,” IEEE Communications Letters, vol. 17, no. 9, pp. 1688–1691, 2013.
View at: Publisher Site | Google Scholar
O. A. Dobre, A. Abdi, Y. Bar-Ness, and W. Su, “Blind modulation classification: a concept whose time has come,” in Proceedings of the IEEE/Sarnoff Symposium on Advances in Wired and Wireless Communication, pp. 223–228, April 2005.
View at: Publisher Site | Google Scholar
F. Hameed, O. A. Dobre, and D. C. Popescu, “On the likelihood-based approach to modulation classification,” IEEE Transactions on Wireless Communications, vol. 8, no. 12, pp. 5884–5892, 2009.
View at: Publisher Site | Google Scholar
B. F. Beidas and C. L. Weber, “Asynchronous classification of mfsk signals using the higher order correlation domain,” IEEE Transactions on Communications, vol. 46, no. 4, pp. 480–493, 1998.
View at: Publisher Site | Google Scholar
A. E. El-Mahdy and N. M. Namazi, “Classification of multiple M-ary frequency-shift keying signals over a Rayleigh fading channel,” IEEE Transactions on Communications, vol. 50, no. 6, pp. 967–974, 2002.
View at: Publisher Site | Google Scholar
O. Ozdemir, P. K. Varshney, and W. Su, “Asynchronous hybrid maximum likelihood classification of linear modulations,” in Proceedings of the 2013 IEEE Global Communications Conference, GLOBECOM 2013, pp. 3235–3240, December 2013.
View at: Publisher Site | Google Scholar
A. K. Nandi and E. E. Azzouz, “Algorithms for automatic modulation recognition of communication signals,” IEEE Transactions on Communications, vol. 46, no. 4, pp. 431–436, 1998.
View at: Publisher Site | Google Scholar
Z. Fucai, H. Yihua, and H. Shiqi, “Classification using wavelet packet decomposition and support vector machine for digital modulations,” Journal of Systems Engineering and Electronics, vol. 19, no. 5, pp. 914–918, 2008.
View at: Publisher Site | Google Scholar
P. Li, F. Wang, and Z. Wang, “Algorithm for modulation recognition based on high-order cumulants and subspace decomposition,” in Proceedings of the 8th International Conference on Signal Processing, ICSP 2006, November 2006.
View at: Publisher Site | Google Scholar
M. R. Mirarab and M. A. Sobhani, “Robust modulation classification for PSK/QAM/ASK using higher-order cumulants,” in Proceedings of the 6th International Conference on Information, Communications and Signal Processing, ICICS, December 2007.
View at: Publisher Site | Google Scholar
M. Pedzisz and A. Mansour, “Automatic modulation recognition of MPSK signals using constellation rotation and its 4th order cumulant,” Digital Signal Processing: A Review Journal, vol. 15, no. 3, pp. 295–304, 2005.
View at: Publisher Site | Google Scholar
B. Ramkumar, T. Bose, and M. S. Radenkovic, “Robust multiuser automatic modulation classifier for multipath fading channels,” in Proceedings of the IEEE Symposium on New Frontiers in Dynamic Spectrum, DySPAN 2010, April 2010.
View at: Publisher Site | Google Scholar
D. C. Chang and P. K. Shih, “Cumulants-based modulation classification technique in multipath fading channels,” IET Communications, vol. 9, no. 6, pp. 828–835, 2015.
View at: Publisher Site | Google Scholar
V. Orlić and M. L. Dukić, “Algorithm for automatic modulation classification in multipath channel based on sixth-order cumulants,” in Proceedings of the 9th International Conference on Telecommunication in Modern Satellite, Cable, and Broadcasting Services (TELSIKS '09), pp. 423–426, IEEE, Niš, Serbia, October 2009.
View at: Publisher Site | Google Scholar
M. S. Mühlhaus, M. Öner, O. A. Dobre, and F. K. Jondral, “A low complexity modulation classification algorithm for MIMO systems,” IEEE Communications Letters, vol. 17, no. 10, pp. 1881–1884, 2013.
View at: Publisher Site | Google Scholar
M. W. Aslam, Z. Zhu, and A. K. Nandi, “Automatic modulation classification using combination of genetic programming and KNN,” IEEE Transactions on Wireless Communications, vol. 11, no. 8, pp. 2742–2750, 2012.
View at: Publisher Site | Google Scholar
A. K. Nandi and E. E. Azzouz, “Modulation recognition using artificial neural networks,” Signal Processing, vol. 56, no. 2, pp. 165–175, 1997.
View at: Publisher Site | Google Scholar
M. L. D. Wong and A. K. Nandi, “Automatic digital modulation recognition using artificial neural network and genetic algorithm,” Signal Processing, vol. 84, no. 2, pp. 351–365, 2004.
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 © 2017 Guan Qing Yang. 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.

PDF Download Citation

Download other formats

Order printed copies

Views

2381

Downloads

1111

Citations