Next Article in Journal
Adaptive Probabilistic Optimization Approach for Vibration-Based Structural Health Monitoring
Next Article in Special Issue
Mixture Density Conditional Generative Adversarial Network Models (MD-CGAN)
Previous Article in Journal
Robustness and Sensitivity Tuning of the Kalman Filter for Speech Enhancement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Voice Transformation Using Two-Level Dynamic Warping and Neural Networks

Department of Electrical and Computer Engineering, Utah State University (USU), 4120 Old Main Hill, Logan, UT 84322-4120, USA
*
Author to whom correspondence should be addressed.
Signals 2021, 2(3), 456-474; https://doi.org/10.3390/signals2030028
Submission received: 12 April 2021 / Revised: 24 May 2021 / Accepted: 13 July 2021 / Published: 14 July 2021
(This article belongs to the Special Issue Machine Learning and Signal Processing)

Abstract

:
Voice transformation, for example, from a male speaker to a female speaker, is achieved here using a two-level dynamic warping algorithm in conjunction with an artificial neural network. An outer warping process which temporally aligns blocks of speech (dynamic time warp, DTW) invokes an inner warping process, which spectrally aligns based on magnitude spectra (dynamic frequency warp, DFW). The mapping function produced by inner dynamic frequency warp is used to move spectral information from a source speaker to a target speaker. Artifacts arising from this amplitude spectral mapping are reduced by reconstructing phase information. Information obtained by this process is used to train an artificial neural network to produce spectral warping information based on spectral input data. The performance of the speech mapping compared using Mel-Cepstral Distortion (MCD) with previous voice transformation research, and it is shown to perform better than other methods, based on their reported MCD scores.

1. Introduction

Voice transformation (VT) refers to the process of changing speech so that speech uttered by one speaker (the source speaker) sounds as if another speaker (the target speaker) had spoken it, for example, transforming from a male voice to a female voice [1]. VT has applications such as text-to-speech synthesis (TTS), international dubbing, health-care, multi-media, language education, music, security-related usage, vocal restoration, speech-to-speech translation, and preprocessing for speech recognition, etc. [2,3,4].
The general framework for VT is shown in Figure 1 [5]. In a training stage, performed offline, features are extracted from both source and target speech sources. VT, which maps information from the source speech to target speech, is trained to form conversion rules. In the actual operation, these conversion rules are used to convert features from a source to a target.
In this work, the spectral features computed using the DFT are transformed using a neural network (NN). Speech signals are aligned using dynamic warping (DW) which aligns both temporally and spectrally to produce training data which is used to train a NN. The NN is trained to produce the spectral mapping functions determined by the frequency-domain warping. Several NN structures (differing numbers of layers, and differing numbers of neurons per layer) were examined to determine which had best performance, as measured by mel-cepstral distortion (MCD). On the basis of MCD, this method was compared with previously developed methods, and is shown to perform well. The test set was based on transformation of a male speaker to a female speaker. Future work will test this architecture transformations between different speakers.

2. Previous Voice Transformation Research

Early work on VT dates back to 1990, when the author in [6] proposed using a vector quantization approach, where the vector-quantized speech features are mapped spectrally via codebook mapping. The features transformed include acoustic features such as pitch, formant frequencies and bandwidth, spectral tilt, etc. [7,8]. These features are extracted and represented using vector quantized codebooks. Generally, codebooks are generated for both source speakers and target speaker. DTW is used to find the correspondence between these vectors. These correspondences are then accumulated to find a histogram, which is used as a weighting function to compute the mapping codebook from source speaker to a target speaker. In the synthesis step, the speech iwas quantized using the speaker’s codebook, then all feature parameters are transformed with the mapping codebook from the training step, and finally the speech is synthesized with LPC vocoder.
To improve this general approach, the authors in [9] propose a method of doing VT based on Pitch Synchronous Overlap and Add (PSOLA) to improve the quality of the converted speech signal. In this work, the authors introduced DFW to achieve VT.
Gaussian Mixture Models (GMMs) have been used for VT [10,11,12,13,14,15,16,17]. GMMs model the source and target data with a GMM to produce transformation functions for each Gaussian component [11]. Two approaches have been used for GMM-based VT techniques: first modeling the source and target information separately using a GMM [18], and second modeling the joint density GMM (JDGMM) between source and target speakers [10]. The conversion function is typical assumed to be probabilistic piecewise linear mapping function for each Gaussian. The unknown parameters are calculated by solving the normal linear conversion equations for a least squares solution [18]. The combination of the spectral source aligned vectors and the corresponding target aligned spectral vectors is used to estimate GMM parameters for the joint density model by using expectation maximization algorithm [11]. Although this algorithm works effectively, its conversion performance is still insufficient [19]. The performance of the standard GMM based VT method have been improved by implementing sophisticated techniques, such as Gaussian process regression [20,21], also using WaveNet vocoder in GMM-based voice conversion to improve the naturalness of the synthetic voice [15,16]. An improvement can be achieved by incorporating a trajectory-based conversion algorithm to achieve temporal correlation in the conversion process [18]. Additional features were modeled in order to reduce the over-smoothing effect of the transformed signal parameters, such us transformation based on matrix variate GMMs, in which joint features from the source and target speakers are represented as matrices using multiple frame features [14], global variance [22] and modulation spectrum [23]. Furthermore, a real-time conversion process has also been successfully implemented for state-of-the-art GMM-based VT [24].
As an alternative to the GMM, which relies on probabilistic mappings between the source and target parameters, DFW seeks to warp the source spectral information in frequency, in order to compensate for the differences between the spectral features of the source and the target speakers, to most closely match that of the target speaker [25]. Cepstral information was extracted from both source and target speakers to do DFW [26,27]. In [15] a new framework was proposed by combining frequency warping and exemplar-based method to achieve VT, where the target was represented by applying frequency warping on the high-resolution spectrum. The warping function is generated by a sparse interpolation from a dictionary of exemplar warping functions.
Neural networks have also been used for VT, where the NNs map the source speech feature vectors into the feature vectors of the target speaker [28,29]. The work [2] is an early effort using NN to transform the source speaker formants to target speaker formants using a feedforward NN. In [30], the authors proposed a method for VT using NN with three layers. This method is based on LPC spectral features using radial basis function NN. In [31] the authors compared using a NN and GMM for VT. GMMs capture the joint distribution of the source features and the target feature, while the work in [31] directly maps the spectral source feature information onto the spectral target feature information. Additionally, with GMM they use maximum likelihood parameters generation (MLPG) to obtain a smooth trajectory of spectral features, while the mapping with NN provide best transformation results without using MLPG.
Long short-term memory (LSTM) NNs have been used for VT [32]. The joint use of deep bidirectional (DB) LSTM and identify vectors (i-vectors) [33], Kullback–Leibler divergence and deep NN-based approach [34], and DBLSTM based Recurrent NNs [35,36], have been recently used for VT. More recently, Generative Adversarial Networks (GAN) such as VAW-GAN [37], CycleGAN [38,39], and many-to-many mapping with StarGAN [40] further advance the state-of-the-art.
Previous methods generally work by extracting a parametric acoustic parameter (such as pitch, MFCC, etc.) and transforming this parameter from one speaker to another based on estimated joint densities. This introduces complexities associated with this initial parameter estimation. The approach taken here avoids the need to find parametric acoustic parameters, or to model the joint density, and instead deals directly with spectral information. (That is, this is analogous to spectral estimation, where a parametric model might be used, such as AR or ARMA, or the spectrum may be estimated directly using the DFT.) NNs are used to learn the transformation which maps the spectral information of one speaker to another speaker. We have found that this basic idea yields improvements when the spectral information from the preceding and following speech segments are used to provide context. The transformation is achieved using a two-level DW. Based on the two-level DW it is straightforward to map the source speech to target speech when both are available.

3. Voice Transformation Process

3.1. Two-Level Dynamic Warping

The transformation in our VT is accomplished using a two-level DW transformation which aligns signals both temporally (DTW) and spectrally (DFW).
Let i = 1 , 2 denote the index of a speaker, with i = 1 representing the source speaker and i = 2 representing the target speaker. Let be the index of a segment of speech, where successive segments may be overlapping. Let r i ( , : ) denote the time-domain signal associated with the th segment of speech. This is a vector of length 2 K = 512 . This was transformed using a 512-point FFT. This was normalized as follows. Let D i ( , j ) , i = 1 , 2 , = 1 , 2 , , L i , j = 1 , 2 , , K be the magnitude of the positive frequency coefficients in the DFT of r i ( , : ) . The feature vector is produced by
s i ( , j ) = D i ( , j ) 1 K L i = 1 L i j = 1 K D i ( , j ) 2 , = 1 , 2 , , L i , j = 1 , 2 , , K .
A sequence of these feature vectors is denoted as
S i = { s i ( 1 , : ) , s i ( 2 , : ) , , s i ( L i , : ) } , i = 1 , 2 .
Here, L 1 and L 2 denote the number of source and target feature vectors. Let d T ( 1 , 2 ) denote the minimum cost between the sequences S 1 (source) and S 2 (target), up to segment indices 1 and 2 , computed using DTW to align the segments in time. (The subscript T refers to time warping.) As usual for DTW [41], d T ( 1 , 2 ) is recursively computed as
d T ( 1 , 2 ) = d ( s 1 ( 1 , : ) , s 2 ( 2 , : ) ) + min { d T ( 1 1 , 2 ) , d T ( 1 , 2 1 ) , d T ( 1 1 , 2 1 ) } .
Here, d ( s 1 ( 1 , : ) , s 2 ( 2 , : ) ) is the distance between individual feature vectors at segments 1 and 2 . This is computed using DFW, as described below. Conceptually, the dynamic warping aligns the peaks and valleys between the source and target data, as suggested by Figure 2. When 1 = L 1 and 2 = L 2 , the temporal warping process determines an overall warped metric distance d T ( L 1 , L 2 ) between the sequences S 1 and S 2 . DTW computes a sequence of indices a T = ( a T ( 1 ) , a T ( 2 ) , , a T ( N T ) ) and b T = ( b T ( 1 ) , b T ( 2 ) , , b T ( N T ) ) , which are called the temporal warping function paths. The temporal warping function paths define a time-aligned (TA) function s 1 , T A such that the sequence
s 1 , T A ( a T ( j ) , : ) = s 1 ( b T ( j ) , : ) , j = 1 , 2 , , N T
is as similar to S 2 as possible.
The distance between two feature vectors in (2), d ( s 1 ( 1 , : ) , s 2 ( 2 , : ) ) , can itself be computed using warping between the components of these individual vectors. Because these vectors represent frequency information, this is DFW. DFW is computed in a way similar to the DTW. Let s i = s i ( i , : ) , i = 1 , 2 denote the spectral information vectors at segment i that is passed to the DFW function from the (outer) DTW function. The elements of s i are denoted by s i ( k ) . DFW is applied to calculate the distance between s 1 and s 2 as
d F ( k 1 , k 2 ) = dist ( s 1 ( k 1 ) , s 2 ( k 2 ) ) + min { d F ( k 1 1 , k 2 ) , d F ( k 1 , k 2 1 ) , d F ( k 1 1 , k 2 1 ) } ,
where the subscript F refers to a frequency (or spectral) warping process, and dist ( s 1 ( k 1 ) , s 2 ( k 2 ) ) represents the metric distance between elements of the feature spectral vectors. The metric distance, dist ( s 1 ( k 1 ) , s 2 ( k 2 ) ) , used looks only at the magnitude of the spectral information, so
dist ( s 1 ( k 1 ) , s 2 ( k 2 ) ) = ( | | s 1 ( k 1 ) | | s 2 ( k 2 ) | | ) .
When there are low energy speech segments (such as from unvoiced speech) which have no particular spectral information to match, the distance d ( s 1 , s 2 ) is generally noninformative. In order to downplay this effect the small energy segments, the minimimum energy E min = min ( s 1 2 , s 2 2 ) , is used to scale the distance, as follows. At the end of frequency warping process, the metric distance between spectral vectors to be used in (2) is computed as d ( s 1 , s 2 ) = E m i n d F ( K , K ) .
The DFW on each feature vector produces a sequence of indices a F = ( a F ( 1 ) , a F ( 2 ) , , a F ( N F ) ) and b F = ( b F ( 1 ) , b F ( 2 ) , , b F ( N F ) ) . a F and b F are called the spectral warping function paths. Here N F is the length of the spectral warping paths, which may be different from segment to segment.
The temporally aligned feature vector of speaker s 1 T A is transformed to spectrally match the target spectrum information of speaker s 2 by creating a modified, spectrally-aligned (SA), feature vector according to
s ^ 1 , T F ( a F ( i ) ) = s 1 , T A ( b F ( i ) ) , i = 1 , 2 , , N F .
This SA map drags spectral components of source vectors to produce transformed voice spectra. Each time-aligned and spectrally-aligned segment is inverse Fourier transformed to produce the time signal for segment , r ^ 2 ( , : ) . This time signal is windowed (we used a Hamming window of length 512 with 16-ms overlapping when the sample rate was 16,000 samples/sec) and the overlapping segments are added together.
The combination of the inner and outer warping processes applied to speech signal is portrayed in Figure 3. Starting at the bottom of the diagram, a speech signal is split into different overlapping segments. Spectral features are computed for each segment. These spectral feature vectors are passed through DTW for time alignment, and at every stage of DTW, the metric between vectors is computed using DFW. The method for doing two-level DW to achieve VT is given in Algorithm 1.
Algorithm 1: Dynamic Warping ( D W ) for Voice Transformation.
Input:
First spectral sequence, S 1 = { s 1 ( 1 , 1 : n T ) , s 1 ( 2 , 1 : n T ) , , s 1 ( n F , 1 : n T ) }
Second spectral sequence, S 2 = { s 2 ( 1 , 1 : m T ) , s 1 ( 2 , 1 : m T ) , , s 1 ( m F , 1 : m T ) }
window size w T for DTW process
window size w F for DFW process
Output:
Distance between S 1 and S 2
Indices a ( N ) and b ( N ) for temporal alignment
Indices b ( M ) and b ( M ) for spectral alignment
Begin
Initialize DTW array, D T W = a r r a y [ 0 n T , 0 m T ]
Initialize DFW array, D F W = a r r a y [ 0 n S , 0 m F ]
Adapt window size, w T = m a x ( w T , a b s ( n T m T ) )
For i T = 1 t o n T
For j T = 1 t o m T
   D T W [ i T , j T ] =
 End For j T
End For i T
Set D T W [ 0 , 0 ] = 0
For i F = 1 t o n F
 For j F = 1 t o m F
   D T W [ i F , j F ] =
 End For j F
End For i F
Set D F W [ 0 , 0 ] = 0
For i T = 1 t o n T
 For j T = m a x ( 1 , i T w T ) t o m i n ( m T , i T + w T )
   Adapt window size, w F = m a x ( w F , a b s ( n S m F ) )
   For i F = 1 t o n F
    For j F = m a x ( 1 , i F w F ) t o m i n ( m F , i F + w F )
     c o s t = a b s ( a b s ( S 1 ( i T , i F ) ) ) ( a b s ( S 2 ( j T , j F ) ) )
     D F W [ i F , j F ] : = c o s t + m i n i m u m ( D T W [ i F 1 , j F ] ,
                       D F W [ i F , j F 1 ] ,
                       D F W [ i F 1 , j F 1 ] )
  End For j F
 End For i F
 Searching minimum path through D F W [ i F , j F ] , save a F , b F
c o s t n e w = D F W [ n F , m F ]
 Minimum energy [ s 1 and s 2 ]
d = E m i n c o s t n e w
D T W [ i T , j T ] : = d + m i n i m u m ( D T W [ i T 1 , j T ] ,
                     D T W [ i T , j T 1 ] ,
                     D T W [ i T 1 , j T 1 ] )
End For j T
End For i T
Searching minimum path through D T W [ i T , j T ] , save a T , b T

3.2. Experiment Description and a Proof of Concept Result

Experiments were carried out on CMU ARCTIC database [42]. The database consisted of seven primary sets of recordings, recorded by two US males, one Canadian male with English accent, one Scottish male with English accent, one Indian male, and two US females. Each speaker recorded a set of 1132 English utterances, most being between one and four seconds long and between 5 to 15 words long. The prompt list was designed to be phonetically balanced American English and have diphone coverage representative of the source material. The sound files were recorded in a soundproof booth at 16,000 samples/sec with simultaneous EGG measurements (the speech needed to be clean, or else background disturbances affected the transformed speech). The lexical and phonetic descriptions derived from the US English front-end module. In experiments presented here, one US male and one US female were chosen to do VT from male to female. While this was a limited experiment, since the transformation process was general (being based on DW of spectral features), this was considered a valid proof of concept for this initial experiment.
After identifying the starting point and ending point of the first 720 sentences of CMU ARCTIC database, each sentence was temporally segmented into 32-ms segments, using a Hamming window with 16-ms overlapping, zero-padded, then transformed using a 512-point FFT then normalized as described above. DW was applied on these features and the features of the first 600 sentences were used to train the NN and a separate set of 120 sentences were used for testing the trained NN.
Figure 4a shows a typical spectrogram for a male speaker, using the phrase “Author of the danger trail, Philip Steels, etc.” from the CMU ARCTIC database. Figure 4b shows the spectrogram for a female speaker for the same phrase. Figure 4c shows the VT from male speaker to female speaker. After transformation, the overall energy distribution looked similar to the spectrogram of the female speaker. The fine-scale variations (the ripples in frequency) evident in Figure 4a,b are due to pitch. These lines are also in Figure 4c, but the pitch details are fuzzier. The resulting speech had very severe processing artifacts. (It did not sound very good.) Since the transformation only involved shifting magnitude spectral information, we were led to suspect that the problem might result from a mismatch between the magnitude and phase information of the warped signal.
Figure 5a shows the spectral feature information for one segment of speech for the female speaker. Figure 5b shows the spectral feature information for one time-aligned segment of speech for the male speaker. By applying the DFW, the peaks and valleys of the male spectrum were generally aligned to the peaks and valleys of the female spectrum, as portrayed in Figure 5c. In warping the spectral information, artifacts were introduced into the signal which were not typical of speech. For example, Figure 5c shows that the warping smeared some spectral values across several frequency bins (indicated with green circles) in the plot, resulting in regions of the spectrum that had constant spectral magnitude.
We found that the quality of the speech was significantly improved by using a phase reconstruction algorithm on the warped speech. In this work, the Griffin–Lim algorithm [43] was applied to do phase reconstruction after the spectral transformation has been accomplished [44].
Figure 6 shows the spectrogram for the warped speech after phase reconstruction. The pitch lines were stronger (less fuzzy) after the GLA phase reconstruction. Figure 7 shows the spectral feature information for one segment of speech for the female speaker with the time aligned spectral warped segment for male from the phrase mentioned above. The peaks and valleys were still aligned, but the regions where there were spectral flats in Figure 5c, denoted by the green circles, were replaced with more speech-like spectra.
Table 1 shows the MCD score obtained using warping with and without the GLA phase reconstruction. The MCD was nearly 3 dB better when GLA phase reconstruction was used. The audio file also sounded significantly better. This verified that the MCD was valid as a measure of audio quality, and that the phase reconstruction was an important part of the VT. Phase reconstruction was used in all transformation experiments described below.

4. Spectral Warping Using NN

The proof-of-concept experiments above demonstrated that spectral mapping followed by phase reconstruction produces effective VT. For a complete, viable VT system, it is necessary to establish conversion rules for transforming speech when the target speaker saying the desired phrase is not available. These conversion rules are represented using a NN which is trained to accept spectral inputs from the source speech and to produce the spectral warping function paths, a F and b F . The process is outlined in Figure 8. The input data for training the NN is the magnitude of the spectrum of a segment of speech that has been time-aligned with the target speech. The spectral information for segments of speech from source and target are used in DW to find the ( a F , b F ) spectral warping information, as suggested by Figure 8a. The NN is trained to produce this ( a F ( i ) , b F ( i ) ) information when the magnitude spectrum for the segment s 1 T A is the input. This is suggested by the NN diagram in Figure 8b.
It was found that better VT was produced when the NN uses the spectral magnitude from the preceding segment s 1 T A ( i ) and the following segment s 1 T A ( i + 1 ) as inputs, to provide context for the current segment s 1 T A ( i + 2 ) . This configuration is suggested in Figure 8c.

Interpolating the Data to Achieve Constant Length

As noted above, the length of the warping path ( a F , b F ) may vary among the different segments of speech. This varying length poses a problem for a NN which has a fixed number of outputs. Rather than zero-pad the outputs of the NN to produce a desired length (which would create artifacts) the a F and b F are interpolated to produce a modified a F and b F which is the same (maximum) length for all segments of the speech vector. In our experiments the maximum length of the warping paths was 482 (as determined from training data). The interpolation is described using Matlab-like notation and functions. Let max index denote the maximum length of the interpolated vectors. The interpolation procedes in two steps. First, interpolated indices are computed according to
n ˜ = [ 1 : max index ] × length ( a ) max index ;
Then interpolation among the a values using these interpolated indices is computed by
a ˜ = interpolate ( a , [ 1 : length ( a ) ] , n ˜ )
The sequence b is similarly interpolated. Figure 9 illustrates un-interpolated function paths for a and b with interpolated paths ( a ˜ and b ˜ respectively). The interpolated a ˜ and b ˜ have the same general shape as the original a and b , but have length 482, as desired. The NN is trained to produce the interpolated warping sequences at its output.

5. Mel-Cepstral Distortion as an Objective Measure

We use Mel-Cepstral Distortion (MCD) to evaluate the quality of the transformed speech, and to be able to compare our work to other work in the literature. MCD has been used as an objective error measure for evaluating the quality of synthetic voice [10]. It is a measure of the difference between two sequences of mel-cepstra. It is computed as
M C D = ( 10 / ln 10 ) 2 i = 1 25 ( m c i ( t ) m c i ( w ) ) 2
where m c i ( t ) and m c i ( w ) represent the mel-cepstrum of the target speech and warped speech, respectively.
We used MCD both as a way of evaluating the performance of different NN architectures, as well as a tool for comparing with previous VT research which reports MCD scores.

6. NN Architecture Experiments

A series of experiments was conducted to determine an effective NN structure. Two different input models were considered. In the first model, a single input block of time-aligned spectral magnitude data, s 1 , T A ( : , i ) , is presented as the input of the NN. There are 256 inputs in the block (corresponding to the positive frequencies of the DFT of a 512-sample segment of speech). The NN has 482 outputs for each element of a and b for a total of 964 outputs. (482 is the largest length of the a , b data that was found among the training data.) This architecture is referred to as the 1-IP/1-OP model.
In the second model, three consecutive time-aligned vectors ( s 1 , T A ( : , i ) , s 1 , T A ( : , i + 1 ) , s 1 , T A ( : , i + 2 ) ) were used as the input, for a total of 768 inputs. The three inputs provided temporal/frequency context for the central block. There were still 964 outputs. This architecture is called the 3-IP/1-OP model.
Training was performed by Tensorflow and Keras to estimate the warping paths. The neural architecture was developed in a series of phases, as described below. Initialize weights were selected at random. Each optimizer was trained for different learning rates (1, 0.1, 0.01, 0.001, 0.0001 and 0.00001); we found that a learning rate of 0.001 achieved the best results. Our NN models were trained with different optimizers (Gradient Descent, Adam and RMS Prop). The best optimizer for this work was the Gradient Descent. The NN architecture was developed in a series of phases, as described below.

6.1. Phase One

This phase was designed as a proof of concept: can a NN learn the spectral warping function paths? For this phase, a single phrase was selected to be trained using the two input models mentioned above. This phase also provided a basis for exploring different NN architectures. For all these tests, the phrase, “Author of the danger trail, Philip Steels, etc.” was selected from the CMU-ARCTIC database.
We experimented with both input models on different number of layers in order determine a good architecture. Table 2 shows the number of nodes in each layer and the output function used for that layer for both architecture. The notation 256L 500R 900R 1000R 964L means that this NN had 5-layers with 256 inputs, 964 output nodes with 500, 900, and 1000 neuron nodes in the hidden layers. L and R represent “linear” and “ReLU” activation output functions, respectively. This table also shows results for architectures described in later phases of development.
Figure 10 shows how the NN learned warping paths at various segments of the selected speech signal for different training iterations. In this figure, each subplot represents the mapping for one spectral segment of speech. In each of these, the blue plots shows the “truth” from DW the paths; the NN-learned paths are in orange. In this figure, the horizontal axes are the sequence of indices (1–468) and the vertical axes are the values of the a F , truth and learned. Figure 10a uses the 19-layers-1-IP/1-OP architecture. There was fairly good matching between the training data and the NN output, but Figure 10b, produced using the 19-layers-3-IP/1-OP architecture showed better results, with the training data being virtually indistinguishable from the NN outputs. Figure 10c,d shows how the NN learned a b warping path at various segments of the selected speech signal after 5000 training iterations. Again an excellent match was achieved with 3-IP/1-Op architecture. Not only was the NN better trained as demonstrated by these figures, but the quality of the audible transformed signal was significantly better when the three-input NN was used. Figure 11 shows the MSE between the true and estimated warping paths for the selected architecture (3-IP/1-OP) with 5000 training iterations, demonstrating that the learning essentially converged by that point.
The use of three input blocks provided results superior to the single input block model. On the basis of the MCD values from Table 2 and Figure 10b,d, we concluded that the NN used in the phases of this experiment should be based on the 19-layers-3-IP/1-OP architecture.

6.2. Phase Two

After validating from Phase One that the NN was able to produce good transformed speech for a single phrase, the NN was trained to deal with multiple phrases. In this phase, training was conducted using 600 phrases from a single speaker for training, and a separate set of 120 phrases from a single speaker was used to test. The selected NN architecture was trained for 10,000 iterations.
Figure 12 shows the warping paths for some of the segments, where the DW path (truth) is shown in blue and the learned path is in orange. The NN learned some of the warping paths well. Figure 12b, which was produced using the selected architecture after 10,000 training iterations, showed more matching between the a F warping path and the learned path than Figure 12a, which was produced at 2500 training iterations. Figure 12d, which was produced at 10,000 training iterations, showed more matching between the b F warping path and the learned path than Figure 12c, which produced at 2500 training iterations. By 10,000 training iterations, many of the warping functions were learned well, but there was still a significant number that were not learned well. Figure 13 shows the MSE between the true values of the warping paths and the estimated one.

6.3. Phase Three

While the results of Phase Two were promising, there was still significant error in the estimates of the warping functions. To further improve the performance, it was determined to assist the NN by decomposing the problem into smaller pieces by clustering. To this end, several NNs were trained, each to be responsive to data pertaining to a particular cluster, where clustering was accomplished using k-means clustering. Different numbers of clusters were tried: k = 4 , 5 , 6 , 10 , 15 , 20 . It was found that k = 6 clusters gave good performance. Figure 14 shows the k = 6 clusters for the selected spectral data from the 600 training phrases. The six clusters are shown using different colors.
In the training stage, data from each of the k clusters were used to train k NNs. In the synthesis stage, for each segment of speech data it was determined which cluster the associated feature vector best matched, and the ( a F , b F ) warping functions were obtained from the corresponding NN. These were used to provide the spectral warping.
Figure 15a,b show how the clustering NN learned the paths a and b , respectively, at 5000 training iterations. The DW path (truth) is shown in blue and the NN-learned path is in orange. At this point, the final sound at the 5000 training iterations is much better than the final sound produced from Phase Two, with little signal processing artifacts. Figure 16 shows the MSE between the true values of the warping paths ( a and b ) and the estimated paths for each of the six clusters. Interestingly, the learning demonstrated by the MSE varied from cluster to cluster. The MCD value for this phase is reported in Phase Three row of Table 2.

6.4. Other NN Structures

The fully connected NN structures described above were shown to produce the warping functions with good fidelity. However, they took a lot of training and a lot of memory. In the interest of reducing the NN complexity, we examined convolutional NN structures (CNNs). However, CNNs failed to produce good acoustic results. The MCD scores for CNN architectures are reported in Table 2 as Phase Four. The CNN structures performed significantly worse than the fully connected NN structures.

7. Comparison

From Table 2 the average MCD across the test data was MCD = 2.8 dB. By contrast, in [32], the MCD for male to female conversion was 5.5 dB. The performance of our method was better than reported in [5], which used NN and reported an average MCD of 6.55. The work of VT using NN in [31] reported MCD of 6.1 dB. In all comparisons made, our method was significantly better (as measured by MCD) than previous work. Table 3 shows more comparison between this work and other works, mentioned in the references part, that used MCD scores to evaluate their works. The method presented here worked better than all previous work that reported performance using MCD.

8. Summary and Conclusions

VT has been achieved by employing two-level DW in combination with NN structures. The DW does both temporal and spectral alignment. The path warping functions from DFW are used to train a NN to produce the warping functions. These warping functions were interpolated to provide vectors of constant length. The warped amplitude spectra are shaped by phase reconstruction, which was shown to dramatically improve the speech quality.
To provide sufficient representation capability on larger training data sets, we found it necessary to cluster the training data, training each segment separately then combining the results. This produced low error in the representation of the warping functions.
Based on comparisons with mel-cepstral scores, the VT approach provides performance competitive with previous works.
The primary focus of this work has been on the development of the two-level DW and design of a NN architecture to do VT. Future work will consider transformation among different classes of speakers, including dealing with speakers with different accents.

Author Contributions

Conceptualization, A.-W.A.-D., T.K.M., J.H.G.; methodology, A.-W.A.-D., T.K.M.; software, A.-W.A.-D.; validation, A.-W.A.-D.; formal analysis, A.-W.A.-D.; investigation, A.-W.A.-D.; resources, T.K.M., J.H.G.; writing—original draft preparation, A.-W.A.-D.; writing—review and editing, T.K.M., J.H.G.; supervision, T.K.M.; project administration, T.K.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The code, data, and transformed sounds can be found here: https://digitalcommons.usu.edu/etd/8030, https://github.com/AlwaledAldulaimi/Voice-Transformation-.git (accessed on 23 May 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Helander, E.; Schwarz, J.; Nurminen, J.; Silen, H.; Gabbouj, M. On the impact of alignment on voice conversion performance. In Proceedings of the 2008 9th Annual Conference of the International Speech Communication Association, Brisbane, Australia, 22–26 September 2008; Volume 1, pp. 1453–1456. [Google Scholar]
  2. Narendranath, M.; Murthy, H.A.; Rajendran, S.; Yegnanarayana, B. Transformation of formants for voice conversion using artificial neural networks. Speech Commun. 1995, 16, 207–216. [Google Scholar] [CrossRef] [Green Version]
  3. Turk, O.; Buyuk, O.; Haznedaroglu, A.; Arslan, L.M. Application of voice conversion for cross-language rap singing transformation. In Proceedings of the IEEE 2009 International Conference on Acoustics, Speech and Signal Processing, (ICASSP), Taipei, Taiwan, 19–24 April 2009; pp. 3597–3600. [Google Scholar]
  4. Stylianou, Y. Voice transformation: A survey. In Proceedings of the IEEE 2009 International Conference on Acoustics, Speech and Signal Processing (ICASSP), Taipei, Taiwan, 19–24 April 2009; pp. 3585–3588. [Google Scholar]
  5. Ramos, M.V. Voice Conversion with Deep Learning. Masters of Science Thesis, Universidade do Lisboa, Lisboa, Portugal, 2016. [Google Scholar]
  6. Abe, M.; Nakamura, S.; Shikano, K.; Kuwabara, H. Voice conversion through vector quantization. J. Acoust. Soc. Jpn. 1990, 11, 71–76. [Google Scholar] [CrossRef] [Green Version]
  7. Kuwabara, H. Quality control of speech by modifying formant frequencies and bandwidth. In Proceedings of the 11th International Congress of Phonetic Sciences (ICPhS), Tallinn, Estonia, 1–7 August 1987; pp. 281–284. [Google Scholar]
  8. Childers, D.; Yegnanarayana, B.; Wu, K. Voice conversion: Factors responsible for quality. In Proceedings of the 1985 IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP), Tampa, FL, USA, 26 March 1985; Volume 10, pp. 748–751. [Google Scholar]
  9. Valbret, H.; Moulines, E.; Tubach, J.P. Voice transformation using PSOLA technique. Speech Commun. 1992, 11, 175–187. [Google Scholar] [CrossRef]
  10. Toth, A.R.; Black, A.W. Using articulatory position data in voice transformation. In Proceedings of the 6th International Speech and Communication Association (ISCA) Workshop on Speech Synthesis, Antwerp, Belgium, 27–31 August 2007; pp. 182–187. [Google Scholar]
  11. Helander, E.; Virtanen, T.; Nurminen, J.; Gabbouj, M. Voice conversion using partial least squares regression. IEEE Trans. Audio Speech Lang. Process. 2010, 18, 912–921. [Google Scholar] [CrossRef] [Green Version]
  12. Benisty, H.; Malah, D. Voice conversion using GMM with enhanced global variance. In Proceedings of the 2011 12th Annual Conference of the International Speech Communication Association (ISCA), Florence, Italy, 28–31 August 2011. [Google Scholar]
  13. Sathiarekha, K.; Kumaresan, S. A survey on the evolution of various voice conversion techniques. In Proceedings of the 2016 IEEE 3rd International Conference on Advanced Computing and Communication Systems (ICACCS), Coimbatore, India, 22–23 January 2016; Volume 1, pp. 1–5. [Google Scholar]
  14. Yang, Y.; Uchida, H.; Saito, D.; Minematsu, N. Voice Conversion Based on Matrix Variate Gaussian Mixture Model Using Multiple Frame Features. In Proceedings of the 2016 17th Annual Conference of the International Speech Communication Association (ISCA), San Francisco, CA, USA, 8–12 September 2016. [Google Scholar]
  15. Tian, X.; Lee, S.W.; Wu, Z.; Chng, E.S.; Li, H. An exemplar-based approach to frequency warping for voice conversion. IEEE Trans. Audio Speech Lang. Process. 2017, 25, 1863–1876. [Google Scholar] [CrossRef]
  16. Sisman, B.; Zhang, M.; Li, H. Group sparse representation with wavenet vocoder adaptation for spectrum and prosody conversion. IEEE Trans. Audio Speech Lang. Process. 2019, 27, 1085–1097. [Google Scholar] [CrossRef]
  17. Sisman, B.; Yamagishi, J.; King, S.; Li, H. An overview of voice conversion and its challenges: From statistical modeling to deep learning. IEEE Trans. Audio Speech Lang. Process. 2020, 29, 132–157. [Google Scholar] [CrossRef]
  18. Toda, T.; Black, A.W.; Tokuda, K. Voice conversion based on maximum-likelihood estimation of spectral parameter trajectory. IEEE Trans. Audio Speech Lang. Process. 2007, 15, 2222–2235. [Google Scholar] [CrossRef]
  19. Zen, H.; Nankaku, Y.; Tokuda, K. Probabilistic feature mapping based on trajectory HMMs. In Proceedings of the 2008 9th Annual Conference of the International Speech Communication Association, Brisbane, Australia, 22–26 September 2008. [Google Scholar]
  20. Pilkington, N.C.; Zen, H.; Gales, M.J. Gaussian process experts for voice conversion. In Proceedings of the 2011 12th Annual Conference of the International Speech Communication Association, (ISCA), Florence, Italy, 28–31 August 2011. [Google Scholar]
  21. Xu, N.; Tang, Y.; Bao, J.; Jiang, A.; Liu, X.; Yang, Z. Voice conversion based on Gaussian processes by coherent and asymmetric training with limited training data. Speech Commun. 2014, 58, 124–138. [Google Scholar] [CrossRef]
  22. Takamichi, S.; Toda, T.; Black, A.W.; Nakamura, S. Modulation spectrum-based post-filter for GMM-based voice conversion. In Proceedings of the IEEE Signal and Information Processing Association Annual Summit and Conference (APSIPA), Chiang Mai, Thailand, 9–12 December 2014; pp. 1–4. [Google Scholar]
  23. Takamichi, S.; Toda, T.; Black, A.W.; Neubig, G.; Sakti, S.; Nakamura, S. Postfilters to modify the modulation spectrum for statistical parametric speech synthesis. IEEE/ACM Trans. Audio Speech Lang. Process. 2016, 24, 755–767. [Google Scholar] [CrossRef]
  24. Toda, T.; Muramatsu, T.; Banno, H. Implementation of computationally efficient real-time voice conversion. In Proceedings of the 2012 13th Annual Conference of the International Speech Communication Association, Portland, OR, USA, 9–13 September 2012. [Google Scholar]
  25. Godoy, E.; Rosec, O.; Chonavel, T. Voice conversion using dynamic frequency warping with amplitude scaling, for parallel or nonparallel corpora. IEEE Trans. Audio Speech Lang. Process. 2012, 20, 1313–1323. [Google Scholar] [CrossRef]
  26. Mohammadi, S.H.; Kain, A. Transmutative voice conversion. In Proceedings of the 38th 2013 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Vancouver, BC, Canada, 26–31 May 2013; pp. 6920–6924. [Google Scholar]
  27. Erro, D.; Navas, E.; Hernaez, I. Parametric voice conversion based on bilinear frequency warping plus amplitude scaling. IEEE Trans. Audio Speech Lang. Process. 2012, 21, 556–566. [Google Scholar] [CrossRef]
  28. Ayodeji, A.O.; Oyetunji, S.A. Voice conversion using coefficient mapping and neural network. In Proceedings of the 2016 IEEE International Conference for Students on Applied Engineering (ICSAE), Newcastle, UK, 20–21 October 2016; pp. 479–483. [Google Scholar]
  29. Zhang, M.; Sisman, B.; Zhao, L.; Li, H. Deepconversion: Voice conversion with limited parallel training data. Speech Commun. 2020, 122, 31–43. [Google Scholar] [CrossRef]
  30. Watanabe, T.; Murakami, T.; Namba, M.; Hoya, T.; Ishida, Y. Transformation of spectral envelope for voice conversion based on radial basis function networks. In Proceedings of the 2002 7th International Conference on Spoken Language Processing, Denver, CO, USA, 16–20 September 2002. [Google Scholar]
  31. Desai, S.; Raghavendra, E.V.; Yegnanarayana, B.; Black, A.W.; Prahallad, K. Voice conversion using artificial neural networks. In Proceedings of the 2009 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Taipei, Taiwan, 19–24 April 2009; pp. 3893–3896. [Google Scholar]
  32. Sun, L.; Kang, S.; Li, K.; Meng, H. Voice conversion using deep bidirectional long short-term memory based recurrent neural networks. In Proceedings of the 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Brisbane, Australia, 19–24 April 2015; pp. 4869–4873. [Google Scholar]
  33. Wu, J.; Wu, Z.; Xie, L. On the use of i-vectors and average voice model for voice conversion without parallel data. In Proceedings of the 2016 IEEE Asia-Pacific Signal and Information Processing Association Annual Summit and Conference (APSIPA), Jeju, Korea, 13–16 December 2016; pp. 1–6. [Google Scholar]
  34. Xie, F.L.; Soong, F.K.; Li, H. A KL Divergence and DNN-Based Approach to Voice Conversion without Parallel Training Sentences. In Proceedings of the 2016 17th Annual Conference of the International Speech Communication Association (ISCA), San Francisco, CA, USA, 8–12 September 2016; pp. 287–291. [Google Scholar]
  35. Sun, L.; Li, K.; Wang, H.; Kang, S.; Meng, H. Phonetic posteriorgrams for many-to-one voice conversion without parallel data training. In Proceedings of the 2016 IEEE International Conference on Multimedia and Expo (ICME), Seattle, WA, USA, 11–15 July 2016; pp. 1–6. [Google Scholar]
  36. Sun, L.; Wang, H.; Kang, S.; Li, K.; Meng, H.M. Personalized, Cross-Lingual TTS Using Phonetic Posteriorgrams. In Proceedings of the 2016 17th Annual Conference of the International Speech Communication Association (ISCA), San Francisco, CA, USA, 8–12 September 2016; pp. 322–326. [Google Scholar]
  37. Hsu, C.C.; Hwang, H.T.; Wu, Y.C.; Tsao, Y.; Wang, H.M. Voice conversion from unaligned corpora using variational autoencoding wasserstein generative adversarial networks. arXiv 2017, arXiv:1704.00849. [Google Scholar]
  38. Fang, F.; Yamagishi, J.; Echizen, I.; Lorenzo-Trueba, J. High-quality nonparallel voice conversion based on cycle-consistent adversarial network. In Proceedings of the 2018 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Calgary, AB, Canada, 15–20 April 2018; pp. 5279–5283. [Google Scholar]
  39. Lorenzo-Trueba, J.; Fang, F.; Wang, X.; Echizen, I.; Yamagishi, J.; Kinnunen, T. Can we steal your vocal identity from the Internet?: Initial investigation of cloning Obama’s voice using GAN, WaveNet and low-quality found data. arXiv 2018, arXiv:1803.00860. [Google Scholar]
  40. Kameoka, H.; Kaneko, T.; Tanaka, K.; Hojo, N. Stargan-vc: Non-parallel many-to-many voice conversion using star generative adversarial networks. In Proceedings of the 2018 IEEE Spoken Language Technology Workshop (SLT), Athens, Greece, 18–21 December 2018; pp. 266–273. [Google Scholar]
  41. Parsons, T. Voice and Speech Processing; McGraw Hill: New York, NY, USA, 1987. [Google Scholar]
  42. Kominek, J.; Black, A.W. The CMU Arctic speech databases. In Proceedings of the 5th International Speech and Communication Association (ISCA) Workshop on Speech Synthesis, Pittsburgh, PA, USA, 14–16 June 2004. [Google Scholar]
  43. Griffin, D.; Lim, J. Signal estimation from modified short-time Fourier transform. IEEE Trans. Acoust. Speech Signal Process. 1984, 32, 236–243. [Google Scholar] [CrossRef]
  44. Al-Dulaimi, A.W.; Moon, T.K.; Gunther, J.H. Phase Effects on Speech and Its Influence on Warped Speech. In Proceedings of the 2020 IEEE Intermountain Engineering, Technology and Computing (IETC), Orem, MD, USA, 4–5 May 2020; pp. 1–5. [Google Scholar]
Figure 1. Typical voice transformation system.
Figure 1. Typical voice transformation system.
Signals 02 00028 g001
Figure 2. Distance measure for outer dynamic warping for VT.
Figure 2. Distance measure for outer dynamic warping for VT.
Signals 02 00028 g002
Figure 3. Inner/outer dynamic warping.
Figure 3. Inner/outer dynamic warping.
Signals 02 00028 g003
Figure 4. Spectrogram for male, female and warped male to female.
Figure 4. Spectrogram for male, female and warped male to female.
Signals 02 00028 g004
Figure 5. One spectral segment feature for male, female and warped male to female.
Figure 5. One spectral segment feature for male, female and warped male to female.
Signals 02 00028 g005
Figure 6. Spectrogram information for warped male to female with using GLA.
Figure 6. Spectrogram information for warped male to female with using GLA.
Signals 02 00028 g006
Figure 7. One spectral segment feature for warped male to female with GLA.
Figure 7. One spectral segment feature for warped male to female with GLA.
Signals 02 00028 g007
Figure 8. Spectrogram for male, female and warped male to female. (a) Making DFW training data. (b) Training a neural network—one input spectral vector. (c) Training a neural network—three input spectral vectors.
Figure 8. Spectrogram for male, female and warped male to female. (a) Making DFW training data. (b) Training a neural network—one input spectral vector. (c) Training a neural network—three input spectral vectors.
Signals 02 00028 g008
Figure 9. Un-interpolated and interpolated paths. (a) a and a ˜ . (b) b and b ˜ .
Figure 9. Un-interpolated and interpolated paths. (a) a and a ˜ . (b) b and b ˜ .
Signals 02 00028 g009
Figure 10. Learned warping paths, phase one, 5000 training iterations. (a) Warping path a produced using 1-IP/1-OP architecture. (b) Warping path a produced using 3-IP/1-OP architecture. (c) Warping path b produced using 1-IP/-1-OP architecture. (d) Warping path b produced using 3-IP/1-OP architecture.
Figure 10. Learned warping paths, phase one, 5000 training iterations. (a) Warping path a produced using 1-IP/1-OP architecture. (b) Warping path a produced using 3-IP/1-OP architecture. (c) Warping path b produced using 1-IP/-1-OP architecture. (d) Warping path b produced using 3-IP/1-OP architecture.
Signals 02 00028 g010
Figure 11. Mean squared error (l curve), phase one.
Figure 11. Mean squared error (l curve), phase one.
Signals 02 00028 g011
Figure 12. Learned warping path a , Phase Two. (a) Warping path a , 2500 iterations. (b) Warping path a , 10,000 iterations. (c) Warping path b , 2500 iterations (d) Warping path b , 10,000 iterations.
Figure 12. Learned warping path a , Phase Two. (a) Warping path a , 2500 iterations. (b) Warping path a , 10,000 iterations. (c) Warping path b , 2500 iterations (d) Warping path b , 10,000 iterations.
Signals 02 00028 g012
Figure 13. Mean squared error learning curve, Phase Two.
Figure 13. Mean squared error learning curve, Phase Two.
Signals 02 00028 g013
Figure 14. Clustering of spectral information from 600 phrases ( k = 6 ).
Figure 14. Clustering of spectral information from 600 phrases ( k = 6 ).
Signals 02 00028 g014
Figure 15. Learned warping paths a and b , Phase Three. (a) Warping path a , 5000 iterations. (b) Warping path b , 5000 iterations.
Figure 15. Learned warping paths a and b , Phase Three. (a) Warping path a , 5000 iterations. (b) Warping path b , 5000 iterations.
Signals 02 00028 g015
Figure 16. Mean squared error (learning curve) using cluster method.
Figure 16. Mean squared error (learning curve) using cluster method.
Signals 02 00028 g016
Table 1. MCDs obtained for 2-Level DW with and without GLA.
Table 1. MCDs obtained for 2-Level DW with and without GLA.
No.Two-Level DW MethodMCD [dB]
1DW, with phase reconstruction2.4
2DW, without phase reconstruction6.3
Table 2. MCDs obtained for different architectures.
Table 2. MCDs obtained for different architectures.
Phase No.Input ModelArchitecture TypeNN ArchitectureMCD [dB]
5256L 600R 1200R 1000R 964L20.5
1-IP-1-OP6256L 600R 1000R 1500R 1000R 964L18.01
7256L 600R 1000R 1500R 2000R 1500R 964L18.002
19500L 1500R 2500R 3500R 4500R 5040R 6459R 7459R3.3
8459R 9459R 10459R 8459R 7939R 6939R 5500R 4500R
3500R 2500R 964L
Phase One 20500L 1500R 2500R 3500R 4500R 5040R 6459R 7459R5.5
8459R 9459R 10459R 9459R 8459R 7939R 6939R 5500R
4500R 3500R 2500R 964L
5768L 4500R 6500R 6939R 964L16.85
6768L 4500R 6500R 7040R 6939R 964L15.3
3-IP-1-OP7768L 4500R 6500R 8040R 8939R 6939R 964L12.92
9768L 4500R 6500R 7040R 9459R 8939R 7939R 6939R 964L7.5
10768L 4500R 6500R 7040R 9459R 10459R 8939R 7939R 6939R 964L4.7
191500L 2500R 3500R 4500R 6500R 7040R 8459R2.25
9459R 10459R 10459R 9459R 8459R 7939R 6939R 5500R 4500R
3500R 2500R 964L
201500L 2500R 3500R 4500R 6500R 7040R 8459R 9459R3.5
10459R 11459R 10459R 9459R 8459R 7939R 6939R 5500R 4500R
3500R 2500R 964L
Phase Two3-IP-1-OP191500L 2500R 3500R 4500R 6500R 7040R 8459R 9459R4.8
10459R 9459R 8459R 7939R 6939R 5500R 4500R 3500R
2500R 964L
Phase Three3-IP-1-OPSix-Clusters.Each Cluster trained for 19 layers1500L 2500R 3500R 4500R 6500R 7040R 8459R 9459R 10459R2.8
10459R 9459R 8459R 7939R 6939R 5500R 4500R 3500R 2500R
964L
Phase Four1-IP-1-OP7256L 64R (k = 3) 100R (k = 3) 150R (k = 3) 100R (k = 3) FL D100R 964L13.33
3-IP-1-OP5768L 64R (k = 3) 100R (k = 3) FL D100R 964L13.5
Table 3. Comparative MCDs from previous work.
Table 3. Comparative MCDs from previous work.
Reference NumberReference YearMCD [dB]
This Works20202.8
[30]20025.6
[32]20155.9
[33]20165.1
[16]20195.09
[29]20204.11
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Al-Dulaimi, A.-W.; Moon, T.K.; Gunther, J.H. Voice Transformation Using Two-Level Dynamic Warping and Neural Networks. Signals 2021, 2, 456-474. https://doi.org/10.3390/signals2030028

AMA Style

Al-Dulaimi A-W, Moon TK, Gunther JH. Voice Transformation Using Two-Level Dynamic Warping and Neural Networks. Signals. 2021; 2(3):456-474. https://doi.org/10.3390/signals2030028

Chicago/Turabian Style

Al-Dulaimi, Al-Waled, Todd K. Moon, and Jacob H. Gunther. 2021. "Voice Transformation Using Two-Level Dynamic Warping and Neural Networks" Signals 2, no. 3: 456-474. https://doi.org/10.3390/signals2030028

APA Style

Al-Dulaimi, A. -W., Moon, T. K., & Gunther, J. H. (2021). Voice Transformation Using Two-Level Dynamic Warping and Neural Networks. Signals, 2(3), 456-474. https://doi.org/10.3390/signals2030028

Article Metrics

Back to TopTop