Research Article  Open Access
Falin Zeng, Sunmin Sun, "Study on the Sound Quality of Steady and Unsteady Exhaust Noise", Mathematical Problems in Engineering, vol. 2018, Article ID 6205140, 11 pages, 2018. https://doi.org/10.1155/2018/6205140
Study on the Sound Quality of Steady and Unsteady Exhaust Noise
Abstract
In order to predict and study the sound quality of automobile exhaust noise, Zwicker steadystate and timevarying method were applied to calculate the psychoacoustic objective parameter values in terms of the exhaust noise of sample cars at uniform velocity and accelerated velocity; Thereby, a prediction model of GABP sound quality based on psychoacoustic objective parameters was established. At the same time, wavelet analysis was used to decompose the accelerated signal; in order to overcome the shortcomings such as Heisenberg uncertainty, the RNR (regularization nonstationary regression technique) was applied to compute the WVD distribution (RNRWVD), therefrom obtaining the coefficient matrices of differentband signals after wavelet decomposition, and then A weighting was carried out on the coefficient matrices, so as to establish a new sound quality parameter SQPWRW (sound quality parameter base on wavelet and then proceed to RNRWVD) as the input of GABP model, and therefrom a sound quality prediction model was established. The results indicate that the model based on SQPWRW has higher precision for predicting the sound quality of acceleration signal, and it can better reflect the characteristics of acceleration signal and sound quality.
1. Introduction
The research of automobile NVH has worked its way from noise control to the new stage emphasizing the design of noise and sound quality; the traditional research of vehicle noise aiming for sound pressure level can no longer satisfy the demand of the contemporary consumers.
The sound quality of automobiles reflects the subjective feeling of people towards the noise; the present research on sound quality is mostly based on subjective evaluation test, which can accurately and directly reflect the quality of voice, but it is timeconsuming and laborconsuming. On the ground of that, domestic and foreign scholars put forward prediction models of automobile sound quality based on psychoacoustics parameters. Zhang et al. [1] optimized the genetic algorithm of support vector machines and established the prediction model of diesel engine sound quality based on psychoacoustic objective parameters. Based on the psychoacoustic objective parameters, Zuo Shuguang et al. [2] proposed the multiple linear regressions, neural network, and support vector machine, 3 prediction models for vehicle interior noise, and carried out comparison study. Bi Fengrong et al. [3] established the least squares support vector machine model on the basis of psychoacoustic objective parameters and EEMD signal features, thereby conducting the research of the acoustic quality of diesel engine radiated noise. Lee et al. [4] proposed the wavelet transformbased evaluation parameters of impact sound quality HFEC and the objective parameters with both roughness and volatility as multivariate linear regression model, which were applied to predict the acoustic quality of suspension system components.
To sum up, the present research on the sound quality of automobile mainly stays in the prediction model based on psychoacoustic objective parameters; however, the automobile noise mainly belongs to unsteady signal, and the single use of the analysis of time domain or frequency domain cannot accurately reflect the characteristic of vehicle noise. We should study and extract the signal features from time and frequency domains, while the sound quality prediction model based on advanced time frequency signal processing (such as wavelets, EMD, and WVD) has not yet been widely used [5].
The exhaust noise of vehicles is one of the most important noise sources, the research on noise quality of which is of certain significance to noise pollution control. Firstly, this paper obtained relative psychoacoustic objective parameters of steadystate and acceleration signals in Artermis based on Zwicker steadystate and timevarying algorithm [6], and the noise quality GABP prediction model based on psychoacoustic parameters was established. Later on the methods of wavelet analysis and regularization unsteady regression were introduced to calculate the WVD distribution, and the characteristic parameters SQPWRW of the acceleration exhaust signal were set up as the input of GABP model. The result suggests that the latter can predict the sound quality of nonstationary signal more accurately, which can provide a reference for the research of unsteady exhaust sound quality.
2. Subjective Evaluation Model of Exhaust Noise
Based on GB149679 Measuring Method for Noise of PowerDriven Vehicles, this paper adopted LMS to test and collect the steadystate exhaust noise signals as well as the acceleration noise signals of 10 domestic vehicles. The signals with engine speed of 1000rpm, 2000rpm, 3000rpm, 4000rpm, and 5000rpm were collected, and signals in the entire process in which engine speed changes from 1000rpm to 5000rpm were collected. According to the research experience [7], loudness, sharpness, roughness, volatility, kurtosis, and A sound level were used to be the input parameters of the model. The steadystate noise test results of sample car 1 is shown in Figure 1; it can be concluded that the sound pressure level of noise generally increased with the increase of rotational speed.
2.1. Subjective Evaluation Test
Pair comparison method was utilized in this subjective test. 5s of steadystate signals was intercepted, and 15s of characteristic signals in acceleration was extracted. In order to eliminate the effect of time duration on the subjective evaluation, the steadystate signal samples were treated with delayed time processing. Samples were compared in pair, in which the sample with decent evaluation will obtain 1 point, and the ones with unsatisfying evaluation will gain nothing, so that each sample will get a definite value representing its sound quality, which is more intuitive and conducive for the modelling later on.
A total of 42 people participated in the subjective evaluation; they are postgraduates majoring in vehiclerelated field from a university and related workers in research institutes. There are 24 males and 18 females aged roughly within 24–40. According to sample coincidence degree and consistency coefficient [8], the data of 4 reviewers were eliminated, the average coincidence degree of the samples was 0.783, and the average consistency coefficient was 0.928. Thereof the consistency coefficients were computed by Kendall method [9], and the psychoacoustic objective parameters of the samples and the results of the subjective tests are shown in Table 1. The normalization formula of satisfaction is showed as formula (2).

2.2. Correlation Analysis
In order to study the relationship between the subjective satisfaction degree of exhaust noise and the psychoacoustic objective parameters, a comparison study was carried out between the subjective satisfaction score and psychoacoustic objective parameters. SPSS.19 software was used to do the correlation analysis. The results are shown in Table 2. Because the sample contains unsteady noise, it will produce extreme values with the change of time. Therefore, we applied twotailed Spearman rank correlation to carry on the correlation analysis, wherein the spearman rank correlation formula is as follows:
 
(Note: indicates that when the bilateral confidence level is 0.01, the correlation is significant; indicates that when the bilateral confidence level is 0.05, the correlation is significant.) 
In the formula: and are the ranks of the two variables, the role of which is to change the fixeddistance variable to nonfixeddistance, thereby reducing the effect of extreme values on the results; is the sample number, and is the spearman rank correlation coefficient.
According to the results of correlation analysis, it can be concluded that the correlations between loudness, sharpness, and satisfaction degree are relatively great. The coefficients have reached 0.875 and 0.682, respectively, the influence of these two parameters has a great weight on the noise satisfaction, and, except for the positive correlations between roughness, kurtosis, and subjective satisfaction degree, other parameters have inverse correlations with satisfaction degree.
3. Establishment of GABP Prediction Model for Sound Quality
3.1. Establishment of Steady GABP Model
In this paper, BP (backpropagation) network was used to establish a complex nonlinear mapping relationship between psychoacoustic parameters and subjective satisfaction degree, and GA (genetic algorithms) was applied to optimize the weights and thresholds of neural networks. This not only can solve the problem that the nonlinear model of BP algorithm is easy to fall into local minimum but also can improve the efficiency of the calculation and the accuracy of the model. The model established is shown in Figure 2. The number of hidden layer nodes was calculated according to the formula and the number of nodes was selected 7 according to the training results, wherein and are the number of input and output layer nodes, respectively. In order to improve the accuracy of acoustic quality prediction model, this paper used steadystate noises as the model’s training samples; 46 steadystate samples were chosen as training samples, and the 4 remaining samples were used to validate the model’s accuracy.
Tansig was chosen as the transfer function of the hidden layer, and purelin was chosen as the transfer function of the output layer for the model. Gradient descent algorithm traingd was chosen as network learning algorithm; the learning efficiency was set as 0.1, and the momentum coefficient was set as 0.9; the common mean square error (MSE) was used as the network training object function; the training target was set as 0.001; the maximum number of generations was 200, and the population size was 40, the generation gap (GGAP) was set as 0.85, the crossover rate was set as 0.7, and the mutation rate was set as 0.01. The input and output samples were all normalized before training. The normalized formula is as follows:
The training results of the model are shown in Figure 3. It can be seen from the training results that, with the increase of the iterative times of the genetic algorithm, the target value of the population is decreasing; that is, the adaptability is increasing, and it tends to be stable in the 50th iteration. After the optimization of the genetic algorithm, the object function value of the network training becomes smaller with the increase of training times. At the 100th iteration, the target value of the model tends to stabilize and the error reaches the set target. The value of fitting check of training result(R^{2}) is 0.994; the sample expectation and the training value almost perfectly coincide. The psychoacoustic objective parameter values of the remaining four steady noises were used as the input of the training GABP model, and the validation errors’ percentages are 2.3%, 1.8%, 1.6%, and 3.5%, respectively, and the average verification error is only 2.3%, which proves that the GABP network model has higher precision and can satisfy the requirements of the sound quality research and prediction.
The GABP model established was used to predict the sound quality of the acceleration noise of the 10 vehicle samples, and the predicted results are shown in Table 3. It can be found that the prediction errors of the acceleration noise signals of the 10 vehicles are basically greater than 5%, and the error RMS value is 7.13%.

3.2. Establishment of Unsteady GABP Model
The sound quality prediction model based on steadystate sample training is not ideal when it comes to predicting the sound quality of the acceleration exhaust noise, and, in order to study whether the training samples result in the difference in the model established and lead to insufficient accuracy, 8 groups of the acceleration unsteady noise samples were used as the training samples of GABP model, and the remaining two groups were used in the prediction of the new model.
The network structure of the model is still 671; the training parameters, transfer functions, network training methods, and genetic algorithm parameters are all consistent with model 1. After the training of the model, the objective parameter values of the predicting samples were introduced into the predicting model, and the comparison results of the two predictions are shown in Figure 4.
From the comparison results, it can be perceived that the application of GABP sound quality model trained by acceleration noise samples improved the precision compared with the previous model, but the prediction accuracy was still not satisfying.
4. Based on WVTRNRWVD Sound Quality Model
4.1. Based on WVTRNRWVD Signal Analysis
Wavelet analysis, as a nearly mature analysis method for signal timefrequency, can effectively process acceleration signals, but the classical wavelet analysis is rooted in the superposition of specific basis functions, and it is difficult to approximate the local signal characteristics of the small wave functions derived by a single basis function on different levels. This may result in the loss of nonstationary signal characteristics [10]. In this regard, we can introduce WVD (WignerVille) which is a bilinear time frequency distribution with a higher time frequency resolution and a certain noise suppression capability, to improve the accuracy of signal analysis. However, the WVD spectral decomposition has the disadvantage of cross noise, so the regularization unsteadystate regression (RNR) [11, 12] was proposed in this paper to compute WVD spectral decomposition. Based on the above methods, a new acoustic quality evaluation parameter SQPWRW (sound quality parameter based on wavelet and then proceed) was proposed to evaluate the sound quality of the acceleration exhaust noise; the process chart is shown in Figure 5.
4.1.1. Wavelet Transform
Wavelet is a function with oscillatory attenuation. Given a small wave function ψ(t), the wavelet transform sequence function is a function family [13] transferred from a single primary image wavelet through expansion and translation, as shown in the following form:
wherein, is an expansion factor, reflecting the scale or width of the function, and is a translation factor, reflecting the translation position of the function along the time axis. (t) is the function family obtained by the expansion and translation of wavelet function ψ(t) under the continuous change of and b. Given a squareintegrable signal x(t), then the wavelet transform is as follows:
(a,b) contains the information of x(t) and (t), so the selection of the wavelet function is very important.
Based on the experience, the Daubechies (dbN) wavelets with orthogonality, compactness, and approximate symmetry was selected to analyze the acceleration signals. Firstly, the frequency below 20hz of signals was filtered by Pasteurized highpass filter, so as to eliminate the influence of infrasound, and then highfrequency denoising process was carried out to the filtered noise signals based on Shannon Criterion. Figure 6 shows the contrast spectrum before and after processing.
From the frequency spectrum of the acceleration signals, it can be observed that its energy is basically concentrated in lowfrequency range below 500Hz, but the sharpness was calculated based on Zwicker which focuses on the proportion of mediumhighfrequency components in the entire signals. Roughness and fluctuation reflect the feelings of the human ear on the modulation frequency and amplitude of the sound. Kurtosis reflects the signal abrupt or flat degree; according to Table 2, the roughness, fluctuation, and kurtosis are less related to subjective satisfaction. However, the unsteadystate acceleration signal is very volatile and unordered, so these parameters calculated based on Zwicker could not express the fluctuation and steepness of the unstable signal accurately.
Carrying out the wavelet transform to the processed noise signals based on the Daubechies in MATLAB is then followed by reconstructing the obtained wavelet transform coefficients by wrcoef function. Thereby the result of wavelet decomposition can be obtained, as shown in Figure 7. We can see that the original signal is decomposed into 8layer detail wavelet di (i=1,.8) and approximate wavelet a8 by wavelet transform. Through the wavelet transform, the original signal was processed with multiscale decomposition, and the relationship between the original signal and decomposition wavelets can be expressed in the following formula:
4.1.2. WVD Distribution and Regularization Theories
WVD is proposed by Wigner in the study of quantum mechanics in 1932, which was later applied by Ville in signal analysis; WVD is a quadratic distribution, which can meet many mathematical properties expected by time frequency analysis, and its transformation form is relatively simple, which belongs to the time frequency analysis of strict sense [14]. Given a signal x(t), its WignerVille distribution is defined as
wherein z(t) is the analytic signal; in the time domain, the analytic signal z(t) is defined as
wherein x(t) is the real signal, is the time, is the delayed time, is the frequency, is the transpose of z, and H[x(t)] is the Hilbert transformation [15] of the real signal x(t). From the above formula, it can be found that product term appears in the WVD distribution, which leads to cross interference in the analysis of multifrequency signals, and this cross interference will affect the readability of the signals.
In order to eliminate the cross interference caused by the WVD distribution, regularization regression technique was introduced to compute WVD, the core of which is the shaping regularization theory; firstly, the theory of shaping regularization is briefly expounded.
Regularization technique is a practical mathematical method, which aims to strengthen the limit of estimated model, so as to solve the problem of illposed inverse problem. The most commonly used regularization method is Tikhonov regularization [16]. Fomel [17] proposes shaping regularization theory by considering the function of shaping operator. The method can be used to select the regularization operator in a simple way, such as Gaussian smooth operator and bandpass filter operator. Subsequently, Fomel extends the theory of shaping regularization to nonlinear inverse problem, thereby establishing the theoretical basis of shaping regularization in inverse problem. One assumes that the vector represents data, represents model parameters, and the relationship between the data and the model is defined by forwardmodelling operator , which can be expressed as
With the least square method, the optimization problem can be solved as follows:
wherein indicates the norm of . The aim of least squares optimization method is to estimate the optimal solution under the condition of known data . When the condition number of the operator is large, the direct solution of the inverse problem to calculate is unstable; considering Tikhonov regularization method to be used to strengthen the constraint of the model , there are
In the formula, is the regularization term of Tikhonov, and then the optimization problem of the above formula has the following theoretical solution:
Smooth operator was taken into account in the shaping regularization; in general, smooth operator can be regarded as the mapping of the constraint model in an acceptable space, which is called shaping by Fomel [18]. The shaping operator can be written as
wherein is shaping operator, from which we can derive that
Introducing the above formula into formula (11), we can obtain the theoretical solution under shaping regularization:
Set the discrete WVD distribution as
Its crosscorrelation function is defined as
wherein , and the inverse transformation of WVD can be deduced as follows:
The least squares optimal solution of the formula above is
The correlation functions R(n,m) and (n,k) in the above formula are complex and real, respectively, the above optimization problem is an illposed problem in mathematics, because the unknown quantity is more than the constraint equation, and therefore the shaping regularization algorithm is introduced to solve this problem.
In this paper, the Gaussian smooth operator with adjustable size was used as the shaping operator. What needs to be emphasized is that although WVD expressions are mathematical real functions, the WV(n,k) are complex after using RNR iterations. Therefore, the absolute values of WV(n,k) were taken as time frequency distribution characteristic quantity of WVD. Figure 8 shows the analytic wavelet d8 analysis results of sample car 2 by WVD and RNRWVD. The results show that there are a lot of “burr” in WVD analysis, which is mainly due to the cross noise caused by WVD decomposition algorithm. Through contrasting contour graphs, it can be clearly observed that the interference of cross noises can be effectively overcome by RNRWVD method which is introduced with smooth operator. And this can bring the signal a better smoothness and clearer time frequency resolution. Some false characteristic interference in the signal analysis was also eliminated, and this can make the signal characteristic extraction more accurate.
4.2. Establishment of Sound Quality Model Based on SQPWRW
Based on the above analysis, a new sound quality parameter SQPWRW was established by using wavelet analysis and RNRWVD method, and it was taken as input parameter to establish sound quality prediction model. And a comparison study was carried out with the GABP model based on psychoacoustic parameter of the acceleration signals. The main establishment steps of SQPWRW model are as follows.
Step 1. Wavelet Decomposition. After filtering and denoising the collected acceleration exhaust signals, the wavelet transform was used to decompose the processed signals, and 9 analytic wavelets were obtained, which are an approximate signal and 8 detail signals.
Step 2. RNRWVD Transform. Processing the analytic wavelets obtained by wavelet decomposition based on RNRWVD and then k matrices R (k=9) with m rows and n columns were obtained, in which K is the number of the decomposition wavelets.
Step 3 (weighting). In order to simulate the acoustic filtering characteristic of human ear, the matrices R were treated with A weighting, and the weighting coefficient matrices can be obtained.
Step 4 (SQPWRW calculation). Calculating the acceleration noise quality parameters SQPWRW based on WVT and RNRWVD, the formula is as follows:wherein K is set within 1–9, which is the number of analytic wavelets. RMS] is used to obtain the valid value of each coefficient matrix; the SQPWRW values of each acceleration sample are shown in Table 4.
Unsteadystate signal is different from steadystate signal. Its distribution parameter and distribution rule will change dramatically with time. The SQPWRW parameter presented in this paper can reflect the signal intensity, disorder degree, and jitter degree in time and frequency domain. A coefficient matrix will be obtained after the RNRWVD procession, which records the characteristics of the signal. Taking onedimensional time and onedimensional frequency as an example, but the actual coefficient matrix is m rows and n columns, the schematic diagram is shown in Figure 9. Different curve has different max and min value, different distribution rule and parameters, and different SQPWRW followed. The parameter can almost describe the characteristics of nonstationary signals.

Step 5 (modelling). Normalizing the SQPWR (k=1,2...,9), the results of which were taken as the inputs of the GABP model, and a 9101 sound quality prediction model was established. The algorithm parameters and the objective function were consistent with the model based on psychoacoustic parameters. The 18 acceleration signals were used as the training samples of model, and the left samples were taken as the test data. After training, the mean square error of the prediction values and the target values is 0.0006, as shown in Figure 10, and the value of correlation coefficient (R^{2}) reached 0.996 after the model training.
Table 5 shows a comparison of the prediction results of the three built models, from which it can be found that the model trained by the parameter SQPWRW values is more accurate in predicting the sound quality of acceleration exhaust noise signals, and it can be concluded that the parameter SQPWRW built in this paper can be used to extract the characteristics of the unsteady noise signal, and it is suitable to be used as the training parameter of the prediction model of sound quality.

5. Conclusions
(1) Based on Zwicker steadystate and unsteadystate algorithm, the psychoacoustic objective parameters of steadystate exhaust noise and acceleration exhaust noise were calculated, and the GABP sound quality prediction models were established by steadystate and unsteadystate samples, respectively, which were applied to predict the unsteady acceleration exhaust noise.
(2) By combining WVT and RNRWVD methods, the coefficient matrices featuring the information of time frequency signal characteristics in different frequency range can be obtained, and the obtained coefficient matrices were processed with A weighting, thereby obtaining a new sound quality parameter SQPWRW, which was used as GABP input to train the sound quality model. And a comparison study was conducted with the prediction model based on the psychoacoustic objective parameters. The results suggest that the model trained by SQPWRW has higher precision in terms of the sound quality prediction of the unsteady signal, and the parameter SQPWRW built in this paper can be used to extract the characteristics of the unsteady noise signal, and it is suitable to be used as the training parameter of the prediction model of sound quality, which can provide certain reference for the future research on the sound quality of unsteady noise signal.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
References
 H. Liu, J. Zhang, P. Guo, F. Bi, H. Yu, and G. Ni, “Sound quality prediction for engineradiated noise,” Mechanical Systems and Signal Processing, vol. 56, pp. 277–287, 2015. View at: Publisher Site  Google Scholar
 X. Shen M, S. Zuo G, and L. Lin I, “Interiror sound quality forecast for veheicles based on support vector machine,” Journal of Vibration & Shock, 2010. View at: Google Scholar
 F. Bi, L. Li, J. Zhang, and T. Ma, “Sound quality prediction for diesel engine radiated noise based on EEMDHT and LSSVM,” Tianjin Daxue Xuebao (Ziran Kexue yu Gongcheng Jishu Ban)/Journal of Tianjin University Science and Technology, vol. 50, no. 1, pp. 28–34, 2017. View at: Google Scholar
 S.K. Lee, H.W. Kim, and E.W. Na, “Improvement of impact noise in a passenger car utilizing sound metric based on wavelet transform,” Journal of Sound and Vibration, vol. 329, no. 17, pp. 3606–3619, 2010. View at: Publisher Site  Google Scholar
 U. K. Chandrika and J. H. Kim, “Development of an algorithm for automatic detection and rating of squeak and rattle events,” Journal of Sound and Vibration, vol. 329, no. 21, pp. 4567–4577, 2010. View at: Publisher Site  Google Scholar
 C. Liu, Y. He, and H. Yu, “Loudness characteristics for vehicle interior timevarying noise under braking condition,” Chinese Journal of Automotive Engineering [CJAE], 2013. View at: Google Scholar
 S.K. Lee, “Objective evaluation of interior sound quality in passenger cars during acceleration,” Journal of Sound and Vibration, vol. 310, no. 12, pp. 149–168, 2008. View at: Publisher Site  Google Scholar
 S. Amman, N. Otto, and C. Jones, “Sound quality analysis of vehicle windshield wiper systems,” in Proceedings of the Noise & Vibration Conference & Exposition, 1993. View at: Publisher Site  Google Scholar
 M. Hussain, J. Gölles, A. Ronacher, and H. Schiffbänker, “Statistical evaluation of an annoyance index for engine noise recordings,” in Proceedings of the Noise & Vibration Conference & Exposition, vol. 100, pp. 527–532, 1991. View at: Publisher Site  Google Scholar
 P. Jiang, Q. Shi, W. Chen et al., “A research on the construction of city road driving cycle based on wavelet analysis,” Automotive Engineering, vol. 33, no. 1, pp. 6970, 2011. View at: Google Scholar
 X. Wu and T. Liu, “Spectral decomposition of seismic data with reassigned smoothed pseudo WignerVille distribution,” Journal of Applied Geophysics, vol. 68, no. 3, pp. 386–393, 2009. View at: Publisher Site  Google Scholar
 H. Bai, F. Yingxiong U, and I. Weifu L, “Optimal rate for least squaress regularized regression with Markov chain samples,” Journal of Hubei University, 2016. View at: Google Scholar
 Y. He, C. Liu, Z. Xu, and Z. Zhang, “Combined wavelet denoising for vehicle braking noise signal based on both soft threshold and genetic algorithmbased adaptive threshold,” Qiche Gongcheng/Automotive Engineering, vol. 36, no. 6, pp. 703–708, 2014. View at: Google Scholar
 H. Wang, T. Qiu, and Z. Chen, Non Stationary Random Signal Analysis and Processing, vol. 2, National Defense Industry Press, Beijing, China, 2008.
 X. Wang and Q. Cao W, “The Hilbert transform and its characters,” Journal of Hubei University, 2008. View at: Google Scholar
 Z. Zhang, S. Chen, and Z. Xu, “Iterative regularization method in generalized inverse beamforming,” Journal of Sound Vibration, p. 396, 2017. View at: Google Scholar
 S. Fomel, “Adaptive multiple subtraction using regularized nonstationary regression,” Geophysics, vol. 74, no. 1, pp. V25–V33, 2009. View at: Publisher Site  Google Scholar
 S. Fomel, “Shaping regularization in geophysicalestimation problems,” Geophysics, vol. 72, no. 2, pp. R29–R36, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Falin Zeng and Sunmin Sun. 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.