JP2016200750A - Encoding device, decoding device and method and program therefor - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide an encoding or decoding technique capable of performing more efficient encoding or decoding than before.SOLUTION: An encoding device comprises: a linear predictive analysis part 22 obtaining a coefficient which can be converted into a linear prediction coefficient and energy of a predictive residue by taking a linear predictive analysis of a pseudo-correlation function signal sequence obtained by performing inverse Fourier transformation in which the absolute value, raised to ηth power, of a frequency region sample sequence corresponding to a time-series signal is regarded as a power spectrum; an encoding part 26 performing encoding, in which bit assignment substantially changes, on each coefficient of the frequency region sample sequence corresponding to the time-series signal based upon information corresponding to a spectrum envelope and a model parameter estimated based upon a coefficient which can be transformed into a linear predictive coefficient while using the square root of the energy of the predictive residue as the model parameter; and a model parameter information encoding part 29 obtaining a model parameter information code representing the information corresponding to the model parameter.SELECTED DRAWING: Figure 4

Description

  The present invention relates to a technique for encoding or decoding a time series signal such as a sound signal.

  As an encoding method for sound signals of low bits (for example, about 10 kbit / s to 20 kbit / s), adaptive encoding for orthogonal transform coefficients in the frequency domain such as DFT (Discrete Fourier Transform) and MDCT (Modified Discrete Cosine Transform) is available. Are known. For example, MEPG USAC (Unified Speech and Audio Coding), a standard technology, has a TCX (transform coded excitation) coding mode, in which MDCT coefficients are normalized for each frame and variable after quantization. Long encoding is performed (for example, see Non-Patent Document 1).

  A configuration example of a conventional TCX-based encoding device is shown in FIG. Hereinafter, each part of FIG. 1 will be described.

<Frequency domain converter 11>
A sound signal which is a time-series signal in the time domain is input to the frequency domain converter 11. The sound signal is, for example, an audio signal or an acoustic signal.

  The frequency domain transform unit 11 converts an input time domain sound signal into N frequency MDCT coefficient sequences X (0), X (1),..., X (N− Convert to 1). N is a positive integer.

  The converted MDCT coefficient sequences X (0), X (1),..., X (N−1) are output to the envelope normalization unit 14.

<Linear prediction analysis unit 12>
The linear prediction analysis unit 12 receives a sound signal that is a time-series signal in the time domain.

The linear prediction analysis unit 12 generates linear prediction coefficients α ₁ , α ₂ ,..., Α _p by performing linear prediction analysis on the sound signal input in units of frames. Further, the linear prediction analysis unit 12 encodes the generated linear prediction coefficients α ₁ , α ₂ ,..., Α _p to generate a linear prediction coefficient code. Examples of the linear prediction coefficient code is the linear prediction coefficients α _1, α _2, ..., a LSP code is a code corresponding to the column of the quantized value of the LSP (Line Spectrum Pairs) parameter sequence corresponding to alpha _p. p is an integer of 2 or more.

Further, the linear prediction analysis unit 12 generates quantized linear prediction coefficients ^ α ₁ , ^ α ₂ ,..., ^ Α _p that are linear prediction coefficients corresponding to the generated linear prediction coefficient code.

The generated quantized linear prediction coefficients ^ α ₁ , ^ α ₂ ,..., ^ Α _p are output to the smoothed amplitude spectrum envelope sequence generation unit 14 and the non-smoothed amplitude spectrum envelope sequence generation unit 13. The generated linear prediction coefficient code is output to the decoding device.

  For the linear prediction analysis, for example, a method of obtaining a linear prediction coefficient by obtaining an autocorrelation for a sound signal input in units of frames and performing a Levinson-Durbin algorithm using the obtained autocorrelation is used. Alternatively, the MDCT coefficient sequence obtained by the frequency domain conversion unit 11 is input to the linear prediction analysis unit 12, and the Levinson-Durbin algorithm is performed on the inverse Fourier transform of the square value series of each coefficient of the MDCT coefficient sequence. A method of obtaining a linear prediction coefficient may be used.

<Smoothing Amplitude Spectrum Envelope Sequence Generation Unit 14>
Quantized linear prediction coefficients ^ α ₁ , ^ α ₂ ,..., ^ Α _p generated by the linear prediction analysis unit 12 are input to the smoothed amplitude spectrum envelope sequence generation unit 14.

The smoothed amplitude spectrum envelope sequence generation unit 14 uses the quantized linear prediction coefficients ^ α ₁ , ^ α ₂ ,..., ^ Α _p to smooth the smoothed amplitude spectrum envelope sequence defined by the following equation (B1) ^ W _γ (0), ^ W _γ (1), ..., ^ W _γ (N-1) are generated. Exp (·) is an exponential function with the Napier number as the base, and j is an imaginary unit. γ is a positive constant of 1 or less, and the amplitude spectrum envelope sequence ^ W (0), ^ W (1),…, ^ W (N-1) defined by the following formula (B2) It is a coefficient for smoothing the unevenness, in other words, a coefficient for smoothing the amplitude spectrum envelope series.

The generated smoothed amplitude spectrum envelope sequences ^ _Wγ (0), ^ _Wγ (1),..., ^ _Wγ (N-1) are the envelope normalization unit 15 and the dispersion parameter determination unit of the encoding unit 16. It is output to 163.

<Non-smoothed amplitude spectrum envelope sequence generation unit 13>
The textured amplitude spectral envelope sequence generating unit 13, quantized linear prediction coefficients the linear prediction analyzer 12 generates _{^ α 1, ^ α 2,} ..., ^ α p is input.

The non-smoothed amplitude spectrum envelope sequence generation unit 13 uses the quantized linear prediction coefficients ^ α ₁ , ^ α ₂ ,..., ^ Α _p and uses the unsmoothed amplitude spectrum envelope defined by the above equation (B2). Generate the sequence ^ W (0), ^ W (1), ..., ^ W (N-1).

  The generated non-smoothed amplitude spectrum envelope sequences ^ W (0), ^ W (1),..., ^ W (N-1) are output to the dispersion parameter determination unit 163 of the encoding unit 16.

<Envelope normalization unit 15>
The envelope normalization unit 15 outputs the MDCT coefficient sequence X (0), X (1),..., X (N-1) generated by the frequency domain conversion unit 11 and the smoothed amplitude spectrum envelope sequence generation unit 14. The smoothed amplitude spectrum envelope sequence ^ _Wγ (0), ^ _Wγ (1), ..., ^ _Wγ (N-1) is input.

Envelope normalization unit 15, by normalizing the respective values ^ W γ _(k) of each coefficient X (k) the smoothed amplitude spectrum envelope sequences of MDCT coefficients, normalized MDCT coefficients X _N (0) , X _N (1), ..., X _N (N-1) are generated. That is, X _N (k) = X (k) / ^ W _γ (k) [k = 0, 1,..., N−1].

The generated normalized MDCT coefficient sequences X _N (0), X _N (1),..., X _N (N−1) are output to the encoding unit 16.

Here, in order to realize quantization that audibly reduces distortion, the envelope normalization unit 15 performs a smoothed amplitude spectrum envelope sequence ^ W _γ (0), which is a sequence in which the amplitude spectrum envelope is blunted. ^ W _γ (1), ..., ^ W _γ (N-1) is used to normalize MDCT coefficient sequence X (0), X (1), ..., X (N-1) in units of frames .

<Encoding unit 16>
The encoding unit 16 includes normalized MDCT coefficient sequences X _N (0), X _N (1),..., X _N (N−1) generated by the envelope normalization unit 15, a smoothed amplitude spectrum envelope sequence generation unit 14, the smoothed amplitude spectrum envelope sequence ^ _Wγ (0), ^ _Wγ (1),..., ^ _Wγ (N-1), the non-smoothed amplitude spectrum envelope sequence generator 13 outputs ^ W (0), ^ W (1), ..., ^ W (N-1) is input.

The encoding unit 16 generates a code corresponding to the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1).

Codes corresponding to the generated normalized MDCT coefficient sequences X _N (0), X _N (1),..., X _N (N−1) are output to the decoding device.

A sequence of integer values obtained by dividing the coefficients of normalized MDCT coefficient sequences X _N (0), X _N (1), ..., X _N (N-1) by gain (global gain) g and quantizing the result A code obtained by encoding the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1) is an integer signal code. In the technique of Non-Patent Document 1, the encoding unit 16 determines a gain g such that the number of bits of the integer signal code is equal to or less than the allocated bit number B, which is the number of bits allocated in advance, and as large as possible. To do. Then, the encoding unit 16 generates a gain code corresponding to the determined gain g and an integer signal code corresponding to the determined gain g.

The generated gain code and integer signal code are output to the decoding apparatus as codes corresponding to the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1).

[Specific Example of Encoding Process Performed by Encoding Unit 16]
A specific example of the encoding process performed by the encoding unit 16 will be described.

  A configuration example of a specific example of the encoding unit 16 is shown in FIG. As shown in FIG. 2, the encoding unit 16 includes a gain acquisition unit 161, a quantization unit 162, a dispersion parameter determination unit 168, an arithmetic encoding unit 169, a gain encoding unit 165, and a determination unit 166. For example, a gain updating unit 167 is provided. Hereinafter, each part of FIG. 2 will be described.

<Gain acquisition unit 161>
The gain acquisition unit 161 is a bit in which the number of bits of the integer signal code is allocated in advance from the input normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1). A global gain g that is equal to or smaller than the distribution bit number B, which is a number, and is as large as possible is determined and output. The global gain g obtained by the gain acquisition unit 161 is an initial value of the global gain used in the quantization unit 162.

<Quantization unit 162>
The quantization unit 162 obtains each coefficient of the input normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) by the gain acquisition unit 161 or the gain update unit 167. Then, quantized normalized coefficient sequences X _Q (0), X _Q (1),..., X _Q (N−1), which are sequences based on the integer part of the result of division by the global gain g, are obtained and output.

  Here, the global gain g used when the quantization unit 162 is executed for the first time is the global gain g obtained by the gain acquisition unit 161, that is, the initial value of the global gain. The global gain g used when the quantization unit 162 is executed for the second time or later is the global gain g obtained by the gain update unit 167, that is, the updated value of the global gain.

<Dispersion parameter determination unit 163>
The dispersion parameter determination unit 163 inputs the input unsmoothed amplitude spectrum envelope sequence ^ W (0), ^ W (1),..., ^ W (N-1) and the input smoothed amplitude spectrum envelope sequence ^ From W _γ (0), ^ W _γ (1), ..., ^ W _γ (N-1), the dispersion parameters φ (0), φ (1), ..., Obtain φ (N-1) and output.

<Arithmetic coding unit 164>
The arithmetic encoding unit 164 uses the dispersion parameters φ (0), φ (1),..., Φ (N−1) obtained by the dispersion parameter determination unit 163 to perform quantization normalization obtained by the quantization unit 162. , X _Q (0), X _Q (1), ..., X _Q (N-1) are arithmetically encoded to obtain an integer signal code, and the integer signal code and the number of bits of the integer signal code are consumed The number of bits C is output. This arithmetic code is optimal when the quantized normalized coefficient sequence at each frequency k (= 0,..., N-1) follows a Laplace distribution represented by the following equation for the following random variable X, for example: This bit allocation is performed.

<Determining unit 166>
The determination unit 166 outputs an integer signal code when the number of gain updates is a predetermined number, and also instructs the gain encoding unit 165 to encode the global gain g obtained by the gain updating unit 167. And the number of consumed bits C measured by the arithmetic encoding unit 164 is output to the gain updating unit 167.

<Gain Updater 167>
When the number of consumed bits C measured by the arithmetic coding unit 164 is larger than the allocated bit number B, the gain updating unit 167 updates the global gain g value to a larger value and outputs the updated value. When the number is smaller than the number B, the value of the global gain g is updated to a small value, and the updated value of the global gain g is output.

<Gain Encoding Unit 165>
The gain encoder 165 encodes the global gain g obtained by the gain updater 167 in accordance with the instruction signal output from the determination unit 166, obtains a gain code, and outputs the gain code.

  The integer signal code output from the determination unit 166 and the gain code output from the gain encoding unit 165 are output to the decoding apparatus as codes corresponding to the normalized MDCT coefficient sequence.

  As described above, in the coding based on the conventional TCX, after normalizing the MDCT coefficient sequence using the smoothed amplitude spectrum envelope sequence in which the non-smoothed amplitude spectrum envelope is blunted, the normalized MDCT coefficient sequence is encoded. ing. This encoding method is employed in the above MPEG-4 USAC and the like.

M. Neuendorf, et al., “MPEG Unified Speech and Audio Coding- The ISO / MPEG Standard for High-Efficiency Audio Coding of all Content Types”, AES 132nd Convention, Budapest, Hungary, 2012.

  In the conventional coding apparatus, optimal bit allocation is performed for the Laplace distribution by arithmetic codes. Then, in order to use the information on the unevenness of the spectrum envelope at the time of arithmetic coding, a dispersion parameter corresponding to the dispersion of the Laplace distribution is generated from the envelope value. However, there is diversity in the probability distribution to which the encoding target belongs, and it does not generally follow the Laplace distribution. As described above, if the same bit allocation is performed on the encoding target belonging to the distribution deviating from the assumption, the compression efficiency may be lowered. Also, when other distributions are introduced, it is difficult to improve efficiency unless a dispersion parameter is generated for the distributions and information on the unevenness of the spectrum envelope is correctly incorporated, as in the case of a conventional encoding device.

  It is an object of the present invention to provide an encoding device, a decoding device, and a method and program thereof that can perform encoding or decoding more efficiently than conventional methods.

  According to the encoding apparatus according to one aspect of the present invention, an encoding apparatus that encodes a time-series signal in the frequency domain, where η is a positive number, and an absolute value of a frequency-domain sample sequence corresponding to the time-series signal Linear prediction that obtains the energy of coefficients and prediction residuals that can be converted into linear prediction coefficients by performing linear prediction analysis on the pseudo correlation function signal sequence obtained by performing inverse Fourier transform with the power of η as the power spectrum A non-smoothed spectral envelope that estimates a spectral envelope by obtaining a non-smoothed spectral envelope sequence that is a series obtained by raising an amplitude spectral envelope sequence corresponding to a coefficient that can be converted into a linear prediction coefficient to the 1 / η power Based on the estimated spectral envelope and the information corresponding to the model parameter, the series generator and the square root of the energy of the prediction residual are used as model parameters. A coding unit that performs coding with different allocations on each coefficient of a frequency domain sample sequence corresponding to a time-series signal, a model parameter information coding unit that obtains a model parameter information code representing information according to model parameters, It is equipped with.

  According to the decoding device according to one aspect of the present invention, a decoding device that obtains a frequency domain sample sequence corresponding to a time-series signal by decoding in the frequency domain, and decodes an input linear prediction coefficient code to generate a linear prediction coefficient A linear prediction coefficient decoding unit that obtains a coefficient that can be converted into a model parameter, a model parameter information decoding unit that decodes an input model parameter code to obtain information according to the model parameter, and linear prediction using η as a predetermined positive number A non-smoothed spectrum envelope sequence generation unit for obtaining a non-smoothed spectrum envelope sequence that is a sequence obtained by raising the amplitude spectrum envelope sequence corresponding to a coefficient that can be converted into a coefficient to the 1 / η power, and a non-smoothed spectrum envelope sequence and model parameters Corresponding to the time series signal by decoding the input integer signal code according to the bit allocation which changes substantially based on the information according to the And a decoding unit for obtaining a frequency domain sample sequence.

  Encoding or decoding can be performed more efficiently than in the past.

The block diagram for demonstrating the example of the conventional encoding apparatus. The block diagram for demonstrating the example of the conventional encoding part. The figure for demonstrating generalized Gaussian distribution. The block diagram for demonstrating the example of an encoding apparatus. The flowchart for demonstrating the example of the encoding method. The block diagram for demonstrating the example of an encoding part. The block diagram for demonstrating the example of an encoding part. The flowchart for demonstrating the example of a process of an encoding part. The block diagram for demonstrating the example of a decoding apparatus. The flowchart for demonstrating the example of a decoding method. The flowchart for demonstrating the example of a process of a decoding part. The block diagram for demonstrating the example of an encoding apparatus. The flowchart for demonstrating the example of the encoding method. The block diagram for demonstrating the example of a parameter determination apparatus. The flowchart for demonstrating the example of the parameter determination method. Histogram to explain the technical background.

[Technical background]
Normalization of MDCT sequence X (0), X (1), ..., X (N-1) with smoothed amplitude spectrum envelope is more accurate than MDCT sequence X (0), X (1), ..., X (N-1) is not whitened. Specifically, MDCT coefficient sequences X (0), X (1), ..., X (N-1) are smoothed amplitude spectrum envelope sequences ^ W _γ (0), ^ W _γ (1), ..., ^ Normalized MDCT coefficient sequence X _N (0) = X (0) / ^ W _γ (0), X _N (1) = X (1) / ^ W _γ obtained by normalizing with W _γ (N-1) (1), ..., X _N (N-1) = X (N-1) / ^ W _γ (N-1) is the MDCT coefficient sequence X (0), X (1), ..., X (N- Normalized sequence X (0) / ^ W () obtained by normalizing 1) with unsmoothed amplitude spectrum envelope sequence ^ W (0), ^ W (1), ..., ^ W (N-1) Than 0), X (1) / ^ W (1), ..., X (N-1) / ^ W (N-1), ^ W (0) / ^ W _γ (0), ^ W (1 ) / ^ W _γ (1),…, ^ W (N-1) / ^ W _γ (N-1) is large. Therefore, the MDCT coefficient sequence X (0), X (1), ..., X (N-1) is transformed into the unsmoothed amplitude spectrum envelope sequence ^ W (0), ^ W (1), ..., ^ W (N- Normalized sequence X (0) / ^ W (0), X (1) / ^ W (1),…, X (N-1) / ^ W (N-1 ) Is an envelope irregularity flattened to an extent suitable for encoding in the encoding unit 16, normalized MDCT coefficient sequences X _N (0), X _N (1), …, X _N (N-1) has ^ W (0) / ^ W _γ (0), ^ W (1) / ^ W _γ (1),…, ^ W (N-1) / ^ W _The envelope irregularities represented by the sequence of _γ (N-1) (hereinafter referred to as normalized amplitude spectrum envelope sequence ^ W _N (0), ^ W _N (1),…, W _N (N-1)) It is left.

Fig. 16 shows the envelope irregularities of the normalized MDCT sequence ^ W (0) / ^ _Wγ (0), ^ W (1) / ^ _Wγ (1),…, ^ W (N-1) / ^ _Wγ The frequency of appearance of the value of each coefficient included in the normalized MDCT coefficient sequence when (N-1) takes each value is shown. envelope: The curve of 0.2-0.3 is the normalized MDCT coefficient X _N (k) corresponding to sample k where the irregularity of the normalized MDCT sequence ^ W (k) / ^ _Wγ (k) is 0.2 or more and less than 0.3 Represents the frequency of the values. envelope: The curve of 0.3-0.4 is the normalized MDCT coefficient X _N (k) corresponding to the sample k whose envelope irregularities ^ W (k) / ^ _Wγ (k) of the normalized MDCT sequence is 0.3 or more and less than 0.4 Represents the frequency of the values. envelope: The curve of 0.4-0.5 is the normalized MDCT coefficient X _N (k) corresponding to the sample k whose envelope irregularity ^ W (k) / ^ _Wγ (k) is 0.4 or more and less than 0.5 Represents the frequency of the values.

  Referring to FIG. 16, it can be seen that the average value of each coefficient included in the normalized MDCT coefficient sequence is almost 0, but the variance is related to the envelope value. That is, the larger the envelope unevenness of the normalized MDCT sequence, the wider the base of the curve representing the frequency. It is clear that there is a relationship. In order to realize more efficient compression, encoding using this relationship is performed. Specifically, encoding in which the bit allocation is substantially changed based on the spectrum envelope is performed on each coefficient of the frequency domain coefficient sequence to be encoded.

For this reason, for example, when the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N-1) is arithmetically encoded, the variance determined based on the spectrum envelope is used. Use parameters.

  In addition, since there is diversity in the probability distribution to which the encoding target belongs, an optimal bit allocation assuming the encoding target belonging to a certain probability distribution (for example, Laplace distribution) is assigned to a code belonging to the probability distribution that deviates from the assumption. If it is performed on the conversion target, the compression efficiency may decrease.

  Therefore, a generalized Gaussian distribution represented by the following equation, which is a distribution that can express various probability distributions, is used as the probability distribution to which the encoding target belongs.

In the generalized Gaussian distribution, by changing the parameter η (> 0) which is a shape parameter, various distributions such as a Laplace distribution when η = 1 and a Gaussian distribution when η = 2 as shown in FIG. Can be expressed. η is a predetermined number greater than zero. The value of η may be determined in advance, or may be selected or varied for each frame that is a predetermined time interval. Also, φ in the above equation is a value corresponding to the dispersion of the distribution, and information on the unevenness of the spectrum envelope is incorporated using this value as a dispersion parameter. That is, the dispersion parameters φ (0), φ (1),..., Φ (N−1) are generated from the spectrum envelope, and for the quantized normalized coefficient X _Q (k) at each frequency k, f An arithmetic code that is optimal when _GG (X | φ (k), η) is followed is configured, and encoding is performed using the arithmetic code based on this configuration.

For example, in addition to the information of the prediction residual energy σ ² and the global gain g, information on the distribution to be used is further incorporated, and the quantized normalized coefficient series X _Q (0), X _Q (1) _,. For example, the dispersion parameter for each coefficient of (N-1) is calculated by the following equation (A1).

However, the model parameter σ is the square root of the energy σ ² of the prediction residual.

Specifically, the Levinson-Durbin algorithm is performed on the inverse Fourier transform of a series of values obtained by raising the absolute value of the MDCT coefficient to the power of η, and the resulting linear prediction coefficient is quantized β ₁ , ^ β ₂ ,…, ^ β _p is used in place of the quantized linear prediction coefficients ^ α ₁ , ^ α ₂ ,…, ^ α _{p and} the unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1 ), ..., ^ H (N-1) and smoothed amplitude spectrum envelope sequence ^ H _γ (0), ^ H _γ (1),…, ^ H _γ (N-1) And formula (A3)

Unsmoothed amplitude spectrum envelope series ^ H (0), ^ H (1),…, ^ H (N-1) coefficients corresponding to the smoothed amplitude spectrum envelope series ^ _Hγ (0 ), ^ H _γ (1),…, ^ H _γ (N-1) divided by the coefficients, normalized amplitude spectrum envelope sequence ^ H _N (0) = ^ H (0) / ^ H _γ (0 ), ^ H _N (1) = ^ H (1) / ^ H _γ (1),…, ^ H _N (N-1) = ^ H (N-1) / ^ H _γ (N-1) Then, the dispersion parameter is calculated from the normalized amplitude spectrum envelope sequence and the global gain g by the above equation (A1).

Here, since the value of σ ^{2 / η} / g in the equation (A1) obtained from actual observation varies, in the present invention, σ ^{2 / η} / g is variable for each frame in a predetermined time interval. And By considering σ ^{2 / η 2 /} g, it is possible to assign a more accurate dispersion parameter, and it is possible to realize more efficient encoding or decoding. Therefore, as will be described later, the encoding apparatus transmits information corresponding to the model parameter σ to the decoding apparatus.

The above technique is based on the minimization problem based on the code length when the quantized normalized coefficient sequence X _Q (0), X _Q (1), ..., X _Q (N-1) is arithmetically coded. Is. The derivation of the above technique is described below.

The code length when the quantized normalized coefficient X _Q (k) is encoded by the arithmetic code using the generalized Gaussian distribution of the shape parameter η by the dispersion parameter φ (k) is sufficiently finely quantized. Assuming

Is proportional to In order to reduce the code length, it is assumed that the dispersion parameter sequence φ (0), φ (1),..., Φ (N−1) is obtained based on the linear prediction coefficients that have already been quantized and encoded. The above equation (A4) is transformed into

Can be rewritten. Where ln is the logarithm based on the number of Napiers, C is a constant for the dispersion parameter, and D _IS (X | Y) is the distance from Y to Itakura Saito

Suppose that That is, the minimization problem of code length L for distributed parameter sequence is ^{φ η (k) / (ηB} η (η)) and | is reduced to minimize problems of the sum of Itakura Saito distance between ^η | X _Q ^(k) The Here, the dispersion parameter sequence φ (0), φ (1 ), ..., φ (N-1) and the linear prediction coefficients _{β 1, β 2, ...,} β p, correspondence between the energy sigma ² prediction residual If one is determined, an optimization problem for obtaining a linear prediction coefficient that minimizes the code length can be established. However, in order to use a conventional fast solution, the following correspondence is made.

Quantized normalized coefficient series X _Q (0), X _Q (1), ..., X _Q (N-1) ignore the influence of quantization, MDCT series X (0), X (1), ... , X (N-1) and smoothed amplitude spectrum envelope ^ _Hγ (0), ^ _Hγ (1),…, ^ _Hγ (N-1), respectively, using global gain g, X _Q (k) = X (k) / (g ^ H _γ (k)), the term that depends on the dispersion parameter in equation (A5) is

As shown, it is expressed as the Itakura Saito distance between the absolute value of the MDCT coefficient series and the all-pole spectral envelope. The conventional linear prediction analysis, that is, applying the Levinson-Durbin algorithm to the inverse Fourier transform of the power spectrum, has a linear prediction coefficient that minimizes the Itakura Saito distance between the power spectrum and the all-pole spectral envelope. It is known that this is a desired operation. Therefore, the above code length minimization problem is the same as in the conventional method by applying the Levinson-Durbin algorithm to the ηth power of the amplitude spectrum, that is, the ηth power of the absolute value of the MDCT coefficient series. An optimal solution can be obtained.

[First embodiment]
(Coding)
A configuration example of the encoding apparatus according to the first embodiment is shown in FIG. As shown in FIG. 4, the encoding apparatus of the third embodiment includes a frequency domain transform unit 21, a linear prediction analysis unit 22, a non-smoothed amplitude spectrum envelope sequence generation unit 23, and a smoothed amplitude spectrum envelope sequence generation. For example, a unit 24, an envelope normalization unit 25, an encoding unit 26, a parameter determination unit 27, and a model parameter information encoding unit 29 are provided. An example of each process of the encoding method according to the first embodiment realized by this encoding apparatus is shown in FIG.

  Hereinafter, each part of FIG. 4 will be described.

<Parameter determining unit 27>
In the first embodiment, any one of a plurality of parameters η can be selected by the parameter determination unit 27 for each predetermined time interval.

  It is assumed that the parameter determination unit 27 stores a plurality of parameters η as parameter η candidates. The parameter determination unit 27 sequentially reads one parameter η among the plurality of parameters, and sends it to the linear prediction analysis unit 22, the unsmoothed amplitude spectrum envelope sequence generation unit 23, the decoding unit 26, and the model parameter information encoding unit 29. Output (step A0).

  Frequency domain transform unit 21, linear prediction analysis unit 22, unsmoothed amplitude spectrum envelope sequence generation unit 23, smoothed amplitude spectrum envelope sequence generation unit 24, envelope normalization unit 25, encoding unit 26, and model parameter information encoding unit 29 is a frequency region corresponding to time-series signals in the same predetermined time interval by performing the processing of step A1 to step A6 and step A9 described below, for example, based on each parameter η sequentially read by the parameter determination unit 27 A code is generated for the sample sequence. In general, given the parameter η, two or more codes may be obtained for frequency domain sample sequences corresponding to time-series signals in the same predetermined time interval. In this case, the codes for the frequency domain sample sequences corresponding to the time-series signals in the same predetermined time section are a combination of these two or more obtained codes. In this example, the code is a combination of a linear prediction coefficient code, a gain code, an integer signal code, and a model parameter information code. Thereby, the code | symbol for each parameter (eta) with respect to the frequency domain sample sequence corresponding to the time series signal of the same predetermined time interval is obtained.

  After the process of step A9, the parameter determination unit 27 selects one code from the codes obtained for each parameter η with respect to the frequency domain sample sequence corresponding to the time-series signal in the same predetermined time interval. Then, the parameter η corresponding to the selected code is determined (step A7). This determined parameter η becomes the parameter η for the frequency domain sample sequence corresponding to the time-series signal in the same predetermined time interval. Then, the parameter determining unit 27 outputs the selected code and the code representing the determined parameter η to the decoding device. Details of the process of step A7 by the parameter determination unit 27 will be described later.

  In the following, it is assumed that one parameter η is read by the parameter determination unit 27, and processing is performed on the read one parameter η.

<Frequency domain converter 21>
The frequency domain converter 21 receives a sound signal that is a time-series signal in the time domain. Examples of sound signals are voice digital signals or acoustic digital signals.

  The frequency domain transform unit 21 converts the input time domain sound signal into N frequency MDCT coefficient sequences X (0), X (1),..., X (N− 1) (step A1). N is a positive integer.

  The obtained MDCT coefficient sequences X (0), X (1),..., X (N-1) are output to the linear prediction analysis unit 22 and the envelope normalization unit 25.

  Unless otherwise specified, the subsequent processing is performed in units of frames.

  In this way, the frequency domain conversion unit 21 obtains a frequency domain sample sequence corresponding to the sound signal, for example, an MDCT coefficient sequence.

<Linear prediction analysis unit 22>
The linear prediction analysis unit 22 receives the MDCT coefficient sequence X (0), X (1),..., X (N-1) obtained by the frequency domain conversion unit 21.

The linear prediction analysis unit 22 uses the MDCT coefficient sequence X (0), X (1),..., X (N-1) to define ~ R (0), ~ R defined by the following equation (A7): (1),..., ~ R (N-1) are subjected to linear prediction analysis to generate linear prediction coefficients β ₁ , β ₂ ,..., Β _p, and the generated linear prediction coefficients β ₁ , β ₂ ,. Encode β _p to generate linear prediction coefficient code and quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,…, ^ β _p , which are quantized linear prediction coefficients corresponding to the linear prediction coefficient code ( Step A2).

The generated quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p are output to the non-smoothed spectrum envelope sequence generation unit 23 and the smoothed amplitude spectrum envelope sequence generation unit 24. Note that the energy σ ² of the prediction residual is calculated in the course of the linear prediction analysis process. In this case, the calculated energy σ ² of the prediction residual is output to the model parameter information encoding unit 29 and the variance parameter determining unit 268 of the encoding unit 26.

  Further, the generated linear prediction coefficient code is transmitted to the parameter determination unit 27.

Specifically, the linear prediction analysis unit 22 first converts the absolute value of the MDCT coefficient sequence X (0), X (1),..., X (N-1) to the inverse FFT considering the power spectrum as a power spectrum. By performing the corresponding operation, that is, the operation of equation (A7), the time domain signal corresponding to the absolute value of the MDCT coefficient sequence X (0), X (1),. The pseudo correlation function signal sequence ~ R (0), ~ R (1), ..., ~ R (N-1), which is a sequence, is obtained. Then, the linear prediction analysis unit 22 performs linear prediction analysis using the obtained pseudo correlation function signal sequence ~ R (0), ~ R (1), ..., ~ R (N-1) to obtain a linear prediction coefficient. β ₁ , β ₂ ,..., β _p are generated. Then, the linear prediction analysis unit 22 encodes the generated linear prediction coefficients β ₁ , β ₂ ,..., Β _{p so as} to encode the linear prediction coefficient code and the quantized linear prediction coefficient corresponding to the linear prediction coefficient code. ^ β ₁ , ^ β ₂ ,…, ^ β _p are obtained.

The linear prediction coefficients β ₁ , β ₂ , ..., β _p are obtained when the absolute value of the MDCT coefficient sequence X (0), X (1), ..., X (N-1) is considered as the power spectrum Is a linear prediction coefficient corresponding to a signal in the time domain.

  The generation of the linear prediction coefficient code by the linear prediction analysis unit 22 is performed by, for example, a conventional encoding technique. The conventional encoding technique is, for example, an encoding technique in which a code corresponding to the linear prediction coefficient itself is a linear prediction coefficient code, and a code corresponding to the LSP parameter by converting the linear prediction coefficient into an LSP parameter. For example, an encoding technique for converting a linear prediction coefficient into a PARCOR coefficient and a code corresponding to the PARCOR coefficient as a linear prediction coefficient code. For example, in a coding technique in which a code corresponding to a linear prediction coefficient itself is a linear prediction coefficient code, a plurality of quantized linear prediction coefficient candidates are determined in advance, and each candidate is stored in association with a linear prediction coefficient code in advance. In this technique, any one of candidates is determined as a quantized linear prediction coefficient for the generated linear prediction coefficient, and a quantized linear prediction coefficient and a linear prediction coefficient code are obtained.

  In this way, the linear prediction analysis unit 22 converts the absolute value of the frequency domain sample sequence, which is an MDCT coefficient sequence, for example, to the power spectrum and performs an inverse Fourier transform to obtain a pseudo correlation function signal sequence. On the other hand, linear prediction analysis is performed to generate coefficients that can be converted into linear prediction coefficients.

<Non-smoothed Amplitude Spectrum Envelope Sequence Generation Unit 23>
The unsmoothed amplitude spectrum envelope sequence generation unit 23 receives the quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p generated by the linear prediction analysis unit 22.

Textured amplitude spectral envelope sequence generating unit 23, the quantized linear prediction coefficient _{^ β 1, ^ β 2,} ..., ^ β is the sequence of the amplitude spectrum envelope corresponding to _p textured amplitude spectral envelope sequence ^ H ( 0), ^ H (1),..., ^ H (N-1) are generated (step A3).

  The generated non-smoothed amplitude spectrum envelope sequences ^ H (0), ^ H (1),..., ^ H (N-1) are output to the encoding unit 26.

The unsmoothed amplitude spectrum envelope sequence generation unit 23 uses the quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p to generate the unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H ( 1), ..., ^ H (N-1), the unsmoothed amplitude spectrum envelope sequence defined by equation (A2) ^ H (0), ^ H (1), ..., ^ H (N-1) Is generated.

  In this way, the unsmoothed amplitude spectrum envelope sequence generation unit 23 is a sequence obtained by raising the amplitude spectrum envelope sequence corresponding to the coefficient that can be converted into the linear prediction coefficient generated by the linear prediction analysis unit 22 to the 1 / η power. The spectral envelope is estimated by obtaining a non-smoothed spectral envelope sequence. Here, the sequence obtained by raising c to a power of a sequence composed of a plurality of values, where c is an arbitrary number, is a sequence composed of values obtained by raising each of the plurality of values to the c-th power. For example, a series obtained by raising the amplitude spectrum envelope series to the power of 1 / η is a series constituted by values obtained by raising each coefficient of the amplitude spectrum envelope to the power of 1 / η.

  The 1 / η power processing by the non-smoothed amplitude spectrum envelope sequence generation unit 23 is caused by processing in which the absolute value η power of the frequency domain sample sequence is regarded as a power spectrum performed by the linear prediction analysis unit 22. It is. That is, the process of the 1 / η power by the non-smoothed amplitude spectrum envelope sequence generation unit 23 is performed by the process in which the absolute value of the frequency domain sample sequence performed by the linear prediction analysis unit 22 is regarded as the power spectrum as η. This is done to return the raised value to its original value.

<Smoothing Amplitude Spectrum Envelope Sequence Generation Unit 24>
Quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p generated by the linear prediction analysis unit 22 are input to the smoothed amplitude spectrum envelope sequence generation unit 24.

Smoothing the amplitude spectral envelope sequence generating unit 24, the quantized linear prediction coefficient _{^ β 1, ^ β 2,} ..., smoothing a series blunted amplitude of irregularities of the amplitude spectral envelope of the sequence corresponding to the ^ beta _p Amplitude spectrum envelope sequences ^ _Hγ (0), ^ _Hγ (1),..., ^ _Hγ (N-1) are generated (step A4).

The generated smoothed amplitude spectrum envelope sequences ^ _Hγ (0), ^ _Hγ (1),..., ^ _Hγ (N−1) are output to the envelope normalization unit 25 and the encoding unit 26.

The smoothed amplitude spectrum envelope sequence generation unit 24 uses the quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p and the correction coefficient γ to smooth the smoothed amplitude spectrum envelope sequence ^ H _γ (0), ^ H _γ (1),…, ^ H _γ (N-1), the smoothed amplitude spectrum envelope sequence defined by equation (A3) ^ H _γ (0), ^ H _γ (1),…, ^ H _γ (N-1) is generated.

  Here, the correction coefficient γ is a predetermined constant less than 1, and the amplitude unevenness of the unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1),…, ^ H (N-1) The coefficient for blunting, in other words, the coefficient for smoothing the unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1), ..., ^ H (N-1).

<Envelope normalization unit 25>
The envelope normalization unit 25 includes the MDCT coefficient sequence X (0), X (1),..., X (N-1) obtained by the frequency domain conversion unit 21 and the smoothed amplitude spectrum envelope generation unit 24. ^ H _γ (0), ^ H _γ (1), ..., ^ H _γ (N-1) are input.

The envelope normalization unit 25 converts each coefficient of the MDCT coefficient sequence X (0), X (1),..., X (N-1) into a corresponding smoothed amplitude spectrum envelope sequence ^ H _γ (0), ^ H. Normalized MDCT coefficient sequence X _N (0), X _N (1), ..., X _N (N-1 by normalizing with each value of _γ (1), ..., ^ H _γ (N-1) ) Is generated (step A5).

  The generated normalized MDCT coefficient sequence is output to the encoding unit 26.

For example, the envelope normalization unit 25 sets each coefficient X (k) of the MDCT coefficient sequence X (0), X (1),..., X (N-1) as k = 0, 1,. Is divided by the smoothed amplitude spectrum envelope series ^ H _γ (0), ^ H _γ (1),…, ^ H _γ (N-1) to obtain the normalized MDCT coefficient sequence X _N (0), X _N Each coefficient X _N (k) of (1),..., X _N (N−1) is generated. That is, X _N (k) = X (k) / ^ H _γ (k) where k = 0, 1,..., N−1.

<Encoding unit 26>
The encoding unit 26 includes normalized MDCT coefficient sequences X _N (0), X _N (1),..., X _N (N−1) generated by the envelope normalization unit 25, an unsmoothed amplitude spectrum envelope generation unit. 23, the non-smoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1),..., ^ H (N-1), and the smoothed amplitude spectrum envelope sequence generated by the smoothed amplitude spectrum envelope generation unit 24 ^ _Hγ (0), ^ _Hγ (1),..., ^ _Hγ (N−1) and the prediction residual energy σ ² calculated by the linear prediction analysis unit 22 are input.

  The encoding unit 26 performs encoding, for example, by performing the processing from step A61 to step A65 shown in FIG. 8 (step A6).

The encoding unit 26 obtains a global gain g corresponding to the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) (step A61), and the normalized MDCT coefficient sequence Quantized normalized coefficient series X, which is a series of integer values obtained by quantizing the result of dividing each coefficient of X _N (0), X _N (1), ..., X _N (N-1) by global gain g _Q (0), X _Q (1), ..., X _Q (N-1) is obtained (step A62), and the quantized normalized coefficient series X _Q (0), X _Q (1), ..., X _Q Dispersion parameters φ (0), φ (1), ..., φ (N-1) corresponding to each coefficient of (N-1) are set to global gain g and unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1),…, ^ H (N-1) and smoothed amplitude spectrum envelope series ^ H _γ (0), ^ H _γ (1),…, ^ H _γ (N-1) and the average residual It is obtained from the energy σ ^{2 by} the equation (A1) (step A63), and the quantized normalized coefficient series X _Q (0) using the dispersion parameters φ (0), φ (1),. ), X _Q (1), ..., X _Q (N-1) are arithmetically encoded to obtain an integer signal code ( Step A64), a gain code corresponding to the global gain g is obtained (Step A65).

Here, the normalized amplitude spectrum envelope sequence ^ H _N (0), ^ H _N (1), ..., ^ H _N in the above equation (A1) is the unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1),…, ^ H (N-1) values are converted into corresponding smoothed amplitude spectrum envelope sequences ^ H _γ (0), ^ H _γ (1),…, ^ H _γ (N- Divided by each value of 1), that is, obtained by the following equation (A8).

  The generated integer signal code and gain code are output to the parameter determination unit 27 as codes corresponding to the normalized MDCT coefficient sequence.

  In step A61 to step A65, the encoding unit 26 determines a global gain g such that the number of bits of the integer signal code is equal to or smaller than the allocated bit number B, which is the number of bits allocated in advance, and as large as possible. A function of generating a gain code corresponding to the determined global gain g and an integer signal code corresponding to the determined global gain g is realized.

Of the steps A61 to A65 performed by the encoding unit 26, the characteristic processing is included in step A63, where the global gain g and the quantized normalized coefficient series X _Q (0), X _Q (1 ),..., X _Q (N-1) are encoded to obtain a code corresponding to the normalized MDCT coefficient sequence. The encoding process itself includes various techniques including those described in Non-Patent Document 1. Known techniques exist. Two specific examples of the encoding process performed by the encoding unit 26 will be described below.

[Specific Example 1 of Encoding Process Performed by Encoder 26]
As a specific example 1 of the encoding process performed by the encoding unit 26, an example not including a loop process will be described.

  A configuration example of the encoding unit 26 of the first specific example is shown in FIG. As shown in FIG. 6, the encoding unit 26 of the first specific example includes a gain acquisition unit 261, a quantization unit 262, a dispersion parameter determination unit 268, an arithmetic encoding unit 269, and a gain encoding unit 265. For example. Hereinafter, each part of FIG. 6 will be described.

<Gain acquisition unit 261>
The gain acquisition unit 261 receives the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) generated by the envelope normalization unit 25.

The gain acquisition unit 261 is the number of bits allocated in advance from the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1). A global gain g that is equal to or less than the number of allocated bits B and that is as large as possible is determined and output (step S261). The gain acquisition unit 261 has, for example, a negative correlation between the square root of the total energy of the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) and the allocated bit number B. The multiplication value with a certain constant is obtained as the global gain g and output. Alternatively, the gain acquisition unit 261 calculates the total energy of the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1), the number of allocated bits B, and the global gain g. , And a global gain g may be obtained and output by referring to the table.

  In this way, the gain acquisition unit 261 obtains a gain for dividing all samples of the normalized frequency domain sample sequence that is a normalized MDCT coefficient sequence, for example.

  The obtained global gain g is output to the quantization unit 262, the dispersion parameter determination unit 268, and the model parameter information encoding unit 29.

<Quantization unit 262>
The quantization unit 262 includes the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) generated by the envelope normalization unit 25 and the global obtained by the gain acquisition unit 261. Gain g is input.

The quantization unit 262 is a series of integer parts as a result of dividing each coefficient of the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) by the global gain g. Quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1) are obtained and output (step S262).

  In this way, the quantization unit 262 divides each sample of the normalized frequency domain sample sequence, which is a normalized MDCT coefficient sequence, for example, by the gain and quantizes it to obtain a quantized normalized coefficient sequence.

The obtained quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1) are output to the arithmetic coding unit 269.

<Dispersion parameter determination unit 268>
The variance parameter determination unit 268 includes the parameter η read by the parameter determination unit 27, the global gain g obtained by the gain acquisition unit 261, and the unsmoothed amplitude spectrum envelope sequence ^ H generated by the unsmoothed amplitude spectrum envelope generation unit 23. (0), ^ H (1), ..., ^ H (N-1), the smoothed amplitude spectrum envelope sequence generated by the smoothed amplitude spectrum envelope generator 24 ^ _Hγ (0), ^ _Hγ (1) ,..., ^ H _γ (N−1) and the prediction residual energy σ ² obtained by the linear prediction analysis unit 22 are input.

The dispersion parameter determination unit 268 calculates the global gain g, the unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1), ..., ^ H (N-1), and the smoothed amplitude spectrum envelope sequence ^ H. _{From γ} (0), ^ H _γ (1),…, ^ H _γ (N-1) and the prediction residual energy σ ² , the dispersion parameter sequence φ is obtained by the above formulas (A1) and (A8). Each of the dispersion parameters (0), φ (1),..., Φ (N−1) is obtained and output (step S268).

  The obtained dispersion parameter series φ (0), φ (1),..., Φ (N−1) are output to arithmetic coding section 269.

<Arithmetic Coding Unit 269>
The arithmetic coding unit 269 includes the parameter η read by the parameter determining unit 27, the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N -1) and the dispersion parameter series φ (0), φ (1),..., Φ (N−1) obtained by the dispersion parameter determination unit 268 are input.

The arithmetic coding unit 269 uses a dispersion parameter sequence φ (0) as a dispersion parameter corresponding to each coefficient of the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1). ), φ (1), ..., φ (N-1) using the respective dispersion parameters, the quantized normalized coefficient series X _Q (0), X _Q (1), ..., X _Q (N-1 ) Is arithmetically encoded to obtain and output an integer signal code (step S269).

The arithmetic coding unit 269 performs generalized Gaussian distribution on each coefficient of the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1) during arithmetic coding. An arithmetic code that is optimal when following f _GG (X | φ (k), η) is configured, and encoding is performed using the arithmetic code based on this configuration. As a result, the expected value of the bit allocation to each coefficient of the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N-1) is expressed as the dispersion parameter series φ (0), φ (1),..., φ (N−1).

  The obtained integer signal code is output to the parameter determination unit 27.

Arithmetic coding may be performed across a plurality of coefficients in the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1). In this case, the dispersion parameters of the dispersion parameter series φ (0), φ (1),..., Φ (N-1) are unsmoothed amplitude spectrum envelopes as can be seen from equations (A1) and (A8). Since it is based on the sequence ^ H (0), ^ H (1), ..., ^ H (N-1), the arithmetic coding unit 269 is based on the estimated spectral envelope (unsmoothed amplitude spectral envelope). Thus, it can be said that encoding is performed in which the bit allocation is substantially changed.

<Gain Encoding Unit 265>
The gain encoder 265 receives the global gain g obtained by the gain acquisition unit 261.

  The gain encoder 265 encodes the global gain g to obtain and output a gain code (step S265).

  The generated integer signal code and gain code are output to the parameter determination unit 27 as codes corresponding to the normalized MDCT coefficient sequence.

  Steps S261, S262, S268, S269, and S265 of the first specific example correspond to the above steps A61, A62, A63, A64, and A65, respectively.

[Specific Example 2 of Encoding Process Performed by Encoder 26]
As a specific example 2 of the encoding process performed by the encoding unit 26, an example including a loop process will be described.

  A configuration example of the encoding unit 26 of the specific example 2 is shown in FIG. As illustrated in FIG. 7, the encoding unit 26 of the specific example 2 includes a gain acquisition unit 261, a quantization unit 262, a dispersion parameter determination unit 268, an arithmetic encoding unit 269, a gain encoding unit 265, For example, a determination unit 266 and a gain update unit 267 are provided. Hereinafter, each part of FIG. 7 will be described.

<Gain acquisition unit 261>
The gain unit 261 receives the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) generated by the envelope normalization unit 25.

The gain acquisition unit 261 is the number of bits allocated in advance from the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1). A global gain g that is equal to or less than the number of allocated bits B and that is as large as possible is determined and output (step S261). The gain acquisition unit 261 has, for example, a negative correlation between the square root of the total energy of the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) and the allocated bit number B. The multiplication value with a certain constant is obtained as the global gain g and output.

  The obtained global gain g is output to the quantization unit 262 and the dispersion parameter determination unit 268.

  The global gain g obtained by the gain acquisition unit 261 is an initial value of the global gain used by the quantization unit 262 and the dispersion parameter determination unit 268.

<Quantization unit 262>
The quantization unit 262 includes a normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) generated by the envelope normalization unit 25 and a gain acquisition unit 261 or a gain update unit. The global gain g obtained by 267 is input.

The quantization unit 262 is a series of integer parts as a result of dividing each coefficient of the normalized MDCT coefficient sequence X _N (0), X _N (1),..., X _N (N−1) by the global gain g. Quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1) are obtained and output (step S262).

  Here, the global gain g used when the quantization unit 262 is executed for the first time is the global gain g obtained by the gain acquisition unit 261, that is, the initial value of the global gain. The global gain g used when the quantizing unit 262 is executed for the second time or later is the global gain g obtained by the gain updating unit 267, that is, the updated value of the global gain.

The obtained quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1) are output to the arithmetic coding unit 269.

<Dispersion parameter determination unit 268>
The dispersion parameter determination unit 268 includes the parameter η read by the parameter determination unit 27, the global gain g obtained by the gain acquisition unit 261 or the gain update unit 267, and the non-smoothed amplitude generated by the non-smoothed amplitude spectrum envelope generation unit 23. Spectral envelope sequence ^ H (0), ^ H (1), ..., ^ H (N-1), smoothed amplitude spectrum envelope sequence generated by the smoothed amplitude spectrum envelope generator 24 ^ _Hγ (0), ^ H _γ (1),..., ^ H _γ (N−1) and the prediction residual energy σ ² obtained by the linear prediction analysis unit 22 are input.

The dispersion parameter determination unit 268 calculates the global gain g, the unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1), ..., ^ H (N-1), and the smoothed amplitude spectrum envelope sequence ^ H. _{From γ} (0), ^ H _γ (1),…, ^ H _γ (N-1) and the prediction residual energy σ ² , the dispersion parameter sequence φ is obtained by the above formulas (A1) and (A8). Each of the dispersion parameters (0), φ (1),..., Φ (N−1) is obtained and output (step S268).

  Here, the global gain g used when the dispersion parameter determination unit 268 is executed for the first time is the global gain g obtained by the gain acquisition unit 261, that is, the initial value of the global gain. The global gain g used when the dispersion parameter determination unit 268 is executed for the second time or later is the global gain g obtained by the gain update unit 267, that is, the updated value of the global gain.

  The obtained dispersion parameter series φ (0), φ (1),..., Φ (N−1) are output to arithmetic coding section 269.

<Arithmetic Coding Unit 269>
The arithmetic coding unit 269 includes the parameter η read by the parameter determining unit 27, the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N -1) and the dispersion parameter series φ (0), φ (1),..., Φ (N−1) obtained by the dispersion parameter determination unit 268 are input.

The arithmetic coding unit 269 uses a dispersion parameter sequence φ (0) as a dispersion parameter corresponding to each coefficient of the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1). ), φ (1), ..., φ (N-1) using the respective dispersion parameters, the quantized normalized coefficient series X _Q (0), X _Q (1), ..., X _Q (N-1 ) Are arithmetically encoded to obtain and output an integer signal code and a consumed bit number C that is the number of bits of the integer signal code (step S269).

The arithmetic coding unit 269 performs generalized Gaussian distribution on each coefficient of the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1) during arithmetic coding. An arithmetic code that is optimal when following f _GG (X | φ (k), η) is configured, and encoding is performed using the arithmetic code based on this configuration. As a result, the expected value of the bit allocation to each coefficient of the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N-1) is expressed as the dispersion parameter series φ (0), φ (1),..., φ (N−1).

  The obtained integer signal code and consumed bit number C are output to the determination unit 266.

Arithmetic coding may be performed across a plurality of coefficients in the quantized normalized coefficient series X _Q (0), X _Q (1),..., X _Q (N−1). In this case, the dispersion parameters of the dispersion parameter series φ (0), φ (1),..., Φ (N-1) are unsmoothed amplitude spectrum envelopes as can be seen from equations (A1) and (A8). Since it is based on the sequence ^ H (0), ^ H (1), ..., ^ H (N-1), the arithmetic coding unit 269 is based on the estimated spectral envelope (unsmoothed amplitude spectral envelope). Thus, it can be said that encoding is performed in which the bit allocation is substantially changed.

<Determining unit 266>
The integer signal code obtained by the arithmetic coding unit 269 is input to the determination unit 266.

  The determination unit 266 outputs an integer signal code when the number of gain updates is a predetermined number, and also instructs the gain encoding unit 265 to encode the global gain g obtained by the gain updating unit 267. When the gain update count is less than the predetermined count, the consumed bit count C measured by the arithmetic encoding section 264 is output to the gain update section 267 (step S266).

<Gain Update Unit 267>
The gain updating unit 267 receives the number of consumed bits C measured by the arithmetic coding unit 264.

  The gain updating unit 267 updates the global gain g value when the consumed bit number C is larger than the allocated bit number B, and outputs the updated value. If the consumed bit number C is smaller than the allocated bit number B, the gain update unit 267 The gain g is updated to a smaller value, and the updated global gain g is output (step S267).

  The updated global gain g obtained by the gain updating unit 267 is output to the quantization unit 262, the gain encoding unit 265, and the model parameter information encoding unit 29.

<Gain Encoding Unit 265>
The gain encoding unit 265 receives the output instruction from the determination unit 266 and the global gain g obtained by the gain update unit 267.

  The gain encoding unit 265 encodes the global gain g according to the instruction signal to obtain and output a gain code (step 265).

  The integer signal code output from the determination unit 266 and the gain code output from the gain encoding unit 265 are output to the parameter determination unit 27 as codes corresponding to the normalized MDCT coefficient sequence.

  That is, in this specific example 2, step S267 performed last corresponds to step A61, and steps S262, S263, S264, and S265 correspond to steps A62, A63, A64, and A65, respectively.

  Specific example 2 of the encoding process performed by the encoding unit 26 is described in more detail in International Publication No. WO2014 / 054556 and the like.

<Model Parameter Information Encoding Unit 29>
The model parameter information encoding unit 29 includes the prediction residual energy σ ² obtained by the linear prediction analysis unit 22, the global gain g obtained by the gain acquisition unit 261 or the gain update unit 267, and the parameters determined by the parameter determination unit 27. A parameter η corresponding to the code is input.

The model parameter information encoding unit 29 obtains and outputs a model parameter information code representing information according to the model parameter σ using at least the prediction residual energy σ ² (step A9). The model parameter information code is output to the parameter determination unit 27.

The information corresponding to the model parameter σ is, for example, σ ^{2 / η} / g in the formula (A1). That is, information according to the model parameter σ, where m is a positive number, is, for example, the value of the ratio between σ ^m and the global gain g. For example, m is 2 / η. The model parameter σ and the global gain g are in a mutually dependent relationship. More specifically, since the dynamic range of the spectrum increases as the model parameter σ increases, there is a mutually dependent relationship that the quantization width, that is, the global gain g increases. For this reason, the possible range of the value of the ratio between σ ^m and the global gain g is often limited. Therefore, if the value of the ratio between σ ^m and the global gain g is information corresponding to the model parameter σ, the information corresponding to the model parameter σ can be efficiently quantized (encoded).

  Of course, the information corresponding to the model parameter σ may be the model parameter σ itself, or may be information depending on the other model parameter σ.

  For example, the model parameter information encoding unit 29 stores a plurality of candidates for information according to model parameters. The model parameter information decoding unit 39 selects the model parameter closest to the information corresponding to the model parameter determined by the inputted σ (and further g, η as necessary) from the plurality of information candidates corresponding to the model parameter. A corresponding information candidate is selected, and information for identifying a candidate information corresponding to the selected model parameter is set as a model parameter information code.

<Parameter determining unit 27>
Through the processing from step A1 to step A6 and step A9, a code (in this example, a linear prediction coefficient code) generated for each parameter η with respect to a frequency domain sample sequence corresponding to a time-series signal in the same predetermined time interval. , Gain code, integer signal code, and model parameter information code) are input to the parameter determination unit 27.

  The parameter determination unit 27 selects one code from the codes obtained for each parameter η with respect to the frequency domain sample sequence corresponding to the time-series signal in the same predetermined time interval, and sets the selected code as the selected code. The corresponding parameter η is determined (step A7). This determined parameter η becomes the parameter η for the frequency domain sample sequence corresponding to the time-series signal in the same predetermined time interval. Then, the parameter determining unit 27 outputs the selected code and the parameter code representing the determined parameter η to the decoding device. The selection of the code is performed based on at least one of the code amount of the code and the coding distortion corresponding to the code. For example, the code with the smallest code amount or the code with the smallest coding distortion is selected.

  Here, the coding distortion is an error between a frequency domain sample sequence obtained from an input signal and a frequency domain sample sequence obtained by local decoding of a generated code. The encoding apparatus may include an encoding distortion calculation unit for calculating encoding distortion. The encoding distortion calculation unit includes a decoding unit that performs processing similar to that of the decoding device described below, and locally decodes the code generated by the decoding unit. Thereafter, the coding distortion calculation unit calculates an error between the frequency domain sample sequence obtained from the input signal and the frequency domain sample sequence obtained by local decoding, and obtains the coding distortion.

(Decryption)
A configuration example of a decoding device corresponding to the encoding device is shown in FIG. As shown in FIG. 9, the decoding device of the first embodiment includes a linear prediction coefficient decoding unit 31, a non-smoothed amplitude spectrum envelope sequence generating unit 32, a smoothed amplitude spectrum envelope sequence generating unit 33, and a decoding unit 34. And an envelope denormalization unit 35, a time domain conversion unit 36, a parameter decoding unit 37, and a model parameter information decoding unit 39, for example. An example of each process of the decoding method according to the first embodiment realized by this decoding apparatus is shown in FIG.

  The decoding apparatus receives at least the parameter code, the code corresponding to the normalized MDCT coefficient sequence, and the linear prediction coefficient code output from the encoding apparatus.

  Hereinafter, each part of FIG. 9 will be described.

<Parameter decoding unit 37>
The parameter code output from the encoding device is input to the parameter decoding unit 37.

  The parameter decoding unit 37 obtains a decoding parameter η by decoding the parameter code. The obtained decoding parameter η is output to the non-smoothed amplitude spectrum envelope sequence generation unit 32, the smoothed amplitude spectrum envelope sequence generation unit 33, and the decoding unit 34. The parameter decoding unit 37 stores a plurality of decoding parameters η as candidates. The parameter decoding unit 37 obtains a decoding parameter η candidate corresponding to the parameter code as a decoding parameter η. The plurality of decoding parameters η stored in the parameter decoding unit 37 are the same as the plurality of parameters η stored in the parameter determining unit 27 of the encoding device.

<Model Parameter Information Decoding Unit 39>
The model parameter information decoding unit 39 receives the model parameter information code output from the encoding device.

  The model parameter information decoding unit 39 obtains information corresponding to the model parameter by decoding the model parameter information code. Information corresponding to the obtained model parameter is output to the dispersion parameter determination unit 342 of the decoding unit 34.

  For example, the model parameter information decoding unit 39 stores a plurality of information candidates corresponding to the model parameters. The model parameter information decoding unit 39 selects information candidates corresponding to the model parameters specified by the model parameter information code from among a plurality of information candidates corresponding to the model parameters, thereby obtaining information corresponding to the model parameters. obtain. The plurality of decoding parameters η stored in the model parameter information decoding unit 39 are the same as the plurality of parameters η stored in the model parameter information encoding unit 29 of the encoding device.

<Linear prediction coefficient decoding unit 31>
The linear prediction coefficient decoding unit 31 receives the linear prediction coefficient code output from the encoding device.

Linear prediction coefficient decoding unit 31, for each frame, the linear prediction coefficient code that has been entered for example by decoding by conventional decoding technique decodes the linear prediction coefficient _{^ β 1, ^ β 2,} ..., obtaining ^ beta _p ( Step B1).

The obtained decoded linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p are output to the non-smoothed amplitude spectrum envelope sequence generation unit 32 and the non-smoothed amplitude spectrum envelope sequence generation unit 33.

  Here, the conventional decoding technique is the same as the linear prediction coefficient quantized by decoding the linear prediction coefficient code when the linear prediction coefficient code is a code corresponding to the quantized linear prediction coefficient, for example. A technique for obtaining a decoded linear prediction coefficient, a technique for decoding the linear prediction coefficient code and obtaining the same decoded LSP parameter as the quantized LSP parameter when the linear prediction coefficient code is a code corresponding to the quantized LSP parameter, etc. It is. In addition, linear prediction coefficients and LSP parameters can be converted to each other, and conversion between decoded linear prediction coefficients and decoded LSP parameters is performed according to the input linear prediction coefficient code and information necessary for subsequent processing. What is necessary is just to perform a process. From the above, what includes the decoding process of the linear prediction coefficient code and the conversion process performed as necessary is “decoding by a conventional decoding technique”.

  In this way, the linear prediction coefficient decoding unit 31 decodes the input linear prediction coefficient code, thereby reversing the absolute value of the frequency domain sample sequence corresponding to the time-series signal as a power spectrum. A coefficient that can be converted into a linear prediction coefficient corresponding to a pseudo correlation function signal sequence obtained by performing Fourier transform is generated.

<Non-smoothed Amplitude Spectrum Envelope Sequence Generation Unit 32>
The unsmoothed amplitude spectrum envelope sequence generation unit 32 includes the decoding parameter η obtained by the parameter decoding unit 37 and the decoded linear prediction coefficients ^ β ₁ , ^ β _{2 ,.} Is entered.

Textured amplitude spectral envelope sequence generating unit 32, decodes the linear prediction coefficient _{^ β 1, ^ β 2,} ..., ^ β unsmoothed amplitude spectrum is a series of amplitude spectrum envelope corresponding to _p envelope sequence ^ H (0 ), ^ H (1),..., ^ H (N-1) are generated by the above equation (A2) (step B2).

  The generated non-smoothed amplitude spectrum envelope sequences ^ H (0), ^ H (1),..., ^ H (N-1) are output to the decoding unit 34.

  In this way, the unsmoothed amplitude spectrum envelope sequence generation unit 32 converts the amplitude spectrum envelope sequence corresponding to the coefficient that can be converted into the linear prediction coefficient generated by the linear prediction coefficient decoding unit 31 to 1 / η. A non-smoothed spectral envelope sequence which is a raised sequence is obtained.

<Smoothing Amplitude Spectrum Envelope Sequence Generation Unit 33>
The smoothed amplitude spectrum envelope sequence generation unit 33 receives the decoding parameter η obtained by the parameter decoding unit 37 and the decoded linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p obtained by the linear prediction coefficient decoding unit 31. Entered.

Smoothing the amplitude spectral envelope sequence generating unit 33 decodes the linear prediction coefficient _{^ β 1, ^ β 2,} ..., smoothing the amplitude is a sequence blunted amplitude of irregularities of the amplitude spectral envelope of the sequence corresponding to the ^ beta _p Spectral envelope sequences ^ H _γ (0), ^ H _γ (1),..., ^ H _γ (N-1) are generated by the above equation A (3) (step B3).

The generated smoothed amplitude spectrum envelope sequences ^ _Hγ (0), ^ _Hγ (1),..., ^ _Hγ (N-1) are output to the decoding unit 34 and the envelope denormalization unit 35.

<Decoding unit 34>
The decoding unit 34 includes a decoding parameter η obtained by the parameter decoding unit 37, a code corresponding to the normalized MDCT coefficient sequence output by the encoding device, and a non-smoothed amplitude spectrum generated by the non-smoothed amplitude spectrum envelope generating unit 32. Envelope sequence ^ H (0), ^ H (1),…, ^ H (N-1), smoothed amplitude spectrum envelope sequence generated by the smoothed amplitude spectrum envelope generator 33 ^ H _γ (0), ^ H _γ (1),..., ^ H _γ (N−1) and information corresponding to the model parameters obtained by the model parameter information decoding unit 39 are input.

  The decoding unit 34 includes a dispersion parameter determination unit 342.

The decoding unit 34 performs decoding by performing, for example, the processing from step B41 to step B44 shown in FIG. 11 (step B4). That is, the decoding unit 34 decodes the gain code included in the code corresponding to the input normalized MDCT coefficient sequence for each frame to obtain the global gain g (step B41). The dispersion parameter determination unit 342 of the decoding unit 34 includes a global gain g, a non-smoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1),..., ^ H (N-1) and a smoothed amplitude spectrum envelope sequence. ^ H _γ (0), ^ H _γ (1),…, ^ H _γ (N-1), parameter η, and information corresponding to the model parameter, the dispersion parameter sequence φ (0) by the above equation (A1) , φ (1),..., φ (N−1) are obtained (step B42). The decoding unit 34 converts the integer signal code included in the code corresponding to the normalized MDCT coefficient sequence to arithmetic corresponding to each dispersion parameter of the dispersion parameter sequence φ (0), φ (1),..., Φ (N−1). According to the coding configuration, arithmetic decoding is performed to obtain a decoded normalized coefficient sequence ^ X _Q (0), ^ X _Q (1), ..., ^ X _Q (N-1) (step B43), and decoding normalization Coefficient series ^ X _Q (0), ^ X _Q (1), ..., ^ X _Q (N-1) is multiplied by global gain g and decoded normalized MDCT coefficient sequence ^ X _N (0) , ^ X _N (1),..., ^ X _N (N-1) are generated (step B44). As described above, the decoding unit 34 may perform decoding of the input integer signal code according to bit allocation that substantially changes based on the non-smoothed spectrum envelope sequence.

The generated decoded normalized MDCT coefficient sequence ^ X _N (0), ^ X _N (1),..., ^ X _N (N−1) is output to the envelope denormalization unit 35.

<Envelope inverse normalization unit 35>
The envelope denormalization unit 35 includes a smoothed amplitude spectrum envelope sequence ^ _Hγ (0), ^ _Hγ (1), ..., ^ _Hγ (N-1) generated by the smoothed amplitude spectrum envelope generation unit 33. The decoding normalization MDCT coefficient sequence ^ X _N (0), ^ X _N (1),..., ^ X _N (N-1) generated by the decoding unit 34 is input.

The envelope denormalization unit 35 uses the smoothed amplitude spectrum envelope sequence ^ _Hγ (0), ^ _Hγ (1),…, ^ _Hγ (N-1) to decode the normalized MDCT coefficient sequence ^ X By denormalizing _N (0), ^ X _N (1), ..., ^ X _N (N-1), the decoded MDCT coefficient sequence ^ X (0), ^ X (1), ..., ^ X (N-1) is generated (step B5).

  The generated decoded MDCT coefficient sequence ^ X (0), ^ X (1),..., ^ X (N-1) is output to the time domain conversion unit 36.

For example, the envelope inverse normalization unit 35, k = 0, 1, ..., a N-1, decoding the normalized MDCT coefficients _{^ X N (0), ^} X N (1), ..., ^ X N (N -1) for each coefficient ^ X _N (k), the smoothed amplitude spectrum envelope series ^ H _γ (0), ^ H _γ (1),…, ^ H _γ (N-1) envelope values ^ H _The decoded MDCT coefficient sequence ^ X (0), ^ X (1), ..., ^ X (N-1) is generated by multiplying by _γ (k). That is, ^ X (k) = ^ X _N (k) × ^ H _γ (k) where k = 0, 1,..., N−1.

<Time domain conversion unit 36>
The time domain transform unit 36 receives the decoded MDCT coefficient sequence ^ X (0), ^ X (1),..., ^ X (N-1) generated by the envelope denormalization unit 35.

  The time domain transform unit 36 transforms the decoded MDCT coefficient sequence ^ X (0), ^ X (1), ..., ^ X (N-1) obtained by the envelope denormalization unit 35 into the time domain for each frame. Thus, a sound signal (decoded sound signal) in units of frames is obtained (step B6).

  In this way, the decoding device obtains a time series signal by decoding in the frequency domain.

[Second Embodiment]
The encoding apparatus and method according to the first embodiment generate a code by performing encoding for each of a plurality of parameters η, select an optimal code from the codes generated for each parameter η, and select the selected code. And a parameter code corresponding to the selected code.

  On the other hand, in the encoding apparatus and method of the second embodiment, the parameter determination unit 27 first determines the parameter η, performs encoding based on the determined parameter η, generates a code, and outputs it. . In the second embodiment, the parameter η is made variable by the parameter determination unit 27 for each predetermined time interval. Here, the parameter η being variable for each predetermined time interval means that the parameter η can be changed if the predetermined time interval is changed, and the value of the parameter η is not changed in the same time interval.

  Hereinafter, a description will be given centering on differences from the first embodiment. A duplicate description of the same parts as in the first embodiment is omitted.

(Coding)
A configuration example of the encoding apparatus according to the second embodiment is shown in FIG. As shown in FIG. 12, the encoding device includes a frequency domain transform unit 21, a linear prediction analysis unit 22, a non-smoothed amplitude spectrum envelope sequence generation unit 23, a smoothed amplitude spectrum envelope sequence generation unit 24, and an envelope. For example, a normalization unit 25, an encoding unit 26, a parameter determination unit 27 ′, and a model parameter information encoding unit 29 are provided. An example of each process of the encoding method realized by this encoding apparatus is shown in FIG.

  Hereinafter, each part of FIG. 12 will be described.

<Parameter determining unit 27 '>
A time domain sound signal, which is a time-series signal, is input to the parameter determination unit 27 ′. Examples of sound signals are voice digital signals or acoustic digital signals.

The parameter determining unit 27 ′ determines the parameter η by a process described later based on the input time series signal (step A7 ′).
Η determined by the parameter determination unit 27 ′ is output to the linear prediction analysis unit 22, the non-smoothed amplitude spectrum envelope estimation unit 23, the smoothed amplitude spectrum envelope estimation unit 24, and the encoding unit 26.

  The parameter determination unit 27 ′ generates a parameter code by encoding the determined η. The generated parameter code is transmitted to the decoding device. Also, the parameter η corresponding to the generated parameter code is output to the model parameter information encoding unit 29.

  Details of the parameter determination unit 27 'will be described later.

  Frequency domain transform unit 21, linear prediction analysis unit 22, unsmoothed amplitude spectrum envelope sequence generation unit 23, smoothed amplitude spectrum envelope sequence generation unit 24, envelope normalization unit 25 encoding unit 26, and model parameter information encoding unit 29 Generates a code by the same process as in the first embodiment based on the parameter η determined by the parameter determination unit 27 (step A1 to step A6 and step A9). In this example, the code is a combination of a linear prediction coefficient code, a gain code, an integer signal code, and a model parameter information code. The generated code is transmitted to the decoding device.

  A configuration example of the parameter determination unit 27 'is shown in FIG. As illustrated in FIG. 14, the parameter determination unit 27 ′ includes, for example, a frequency domain conversion unit 41, a spectrum envelope estimation unit 42, a whitened spectrum sequence generation unit 43, and a parameter acquisition unit 44. The spectrum envelope estimation unit 42 includes, for example, a linear prediction analysis unit 421 and a non-smoothed amplitude spectrum envelope sequence generation unit 422. For example, FIG. 2 shows an example of each process of the parameter determination method realized by the parameter determination unit 27 '.

  Hereinafter, each part of FIG. 14 will be described.

<Frequency domain conversion unit 41>
The time domain sound signal, which is a time series signal, is input to the frequency domain transform unit 41. Examples of sound signals are voice digital signals or acoustic digital signals.

  The frequency domain conversion unit 41 converts the input time domain sound signal into N frequency MDCT coefficient sequences X (0), X (1),..., X (N− Convert to 1). N is a positive integer.

  The obtained MDCT coefficient sequences X (0), X (1),..., X (N−1) are output to the spectrum envelope estimation unit 42 and the whitened spectrum sequence generation unit 43.

  Unless otherwise specified, the subsequent processing is performed in units of frames.

  In this way, the frequency domain conversion unit 41 obtains a frequency domain sample sequence corresponding to the sound signal, for example, an MDCT coefficient sequence (step C41).

<Spectrum envelope estimation unit 42>
The spectrum envelope estimation unit 42 receives the MDCT coefficient sequence X (0), X (1),..., X (N−1) obtained by the frequency domain conversion unit 21.

Based on the parameter η ₀ determined by a predetermined method, the spectrum envelope estimation unit 42 performs spectrum envelope estimation using the absolute value η ₀ of the frequency domain sample sequence corresponding to the time-series signal as a power spectrum ( Step C42).

  The estimated spectrum envelope is output to the whitened spectrum sequence generation unit 43.

  The spectrum envelope estimation unit 42 estimates the spectrum envelope by generating a non-smoothed amplitude spectrum envelope sequence, for example, by processing of a linear prediction analysis unit 421 and a non-smoothed amplitude spectrum envelope sequence generation unit 422 described below. .

It is assumed that the parameter η ₀ is determined by a predetermined method. For example, η ₀ is a predetermined number greater than zero. For example, η ₀ = 1. Moreover, you may use (eta) calculated | required by the flame | frame before the frame which is calculating | requiring the present parameter (eta). The frame before the frame for which the current parameter η is to be obtained (hereinafter referred to as the current frame) is, for example, a frame before the current frame and in the vicinity of the current frame. The frame in the vicinity of the current frame is, for example, a frame immediately before the current frame.

<Linear prediction analysis unit 421>
MDCT coefficient sequences X (0), X (1),..., X (N−1) obtained by the frequency domain transform unit 41 are input to the linear prediction analysis unit 421.

The linear prediction analysis unit 421 uses the MDCT coefficient sequence X (0), X (1),..., X (N-1) to define ~ R (0), ~ R defined by the following equation (C1). (1),..., ~ R (N-1) are subjected to linear prediction analysis to generate linear prediction coefficients β ₁ , β ₂ ,..., Β _p, and the generated linear prediction coefficients β ₁ , β ₂ ,. β _p is encoded to generate a linear prediction coefficient code and quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ β _p which are quantized linear prediction coefficients corresponding to the linear prediction coefficient code.

The generated quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p are output to the non-smoothed spectrum envelope sequence generation unit 422.

Specifically, the linear prediction analysis unit 421 first performs an inverse FFT considering the absolute value of the MDCT coefficient sequence X (0), X (1),..., X (N-1) to the power η ₀ as the power spectrum. , That is, in the time domain corresponding to the absolute value of MDCT coefficient sequence X (0), X (1), ..., X (N-1) to the ηth power A pseudo-correlation function signal sequence ~ R (0), ~ R (1), ..., ~ R (N-1) which is a signal string is obtained. Then, the linear prediction analysis unit 421 performs linear prediction analysis using the obtained pseudo correlation function signal sequence ~ R (0), ~ R (1), ..., ~ R (N-1) to obtain a linear prediction coefficient. β ₁ , β ₂ ,..., β _p are generated. Then, the linear prediction analysis unit 421 encodes the generated linear prediction coefficients β ₁ , β ₂ ,..., Β _{p so as} to encode a linear prediction coefficient code and a quantized linear prediction coefficient corresponding to the linear prediction coefficient code. ^ β ₁ , ^ β ₂ ,…, ^ β _p are obtained.

Linear prediction coefficients _{β 1, β 2, ...,} β p is, MDCT coefficient sequence X (0), X (1 ), ..., and the eta ₀ square of the absolute value of X (N-1) was regarded as a power spectrum It is a linear prediction coefficient corresponding to the time domain signal.

  The generation of the linear prediction coefficient code by the linear prediction analysis unit 421 is performed by, for example, a conventional encoding technique. The conventional encoding technique is, for example, an encoding technique in which a code corresponding to the linear prediction coefficient itself is a linear prediction coefficient code, and a code corresponding to the LSP parameter by converting the linear prediction coefficient into an LSP parameter. For example, an encoding technique for converting a linear prediction coefficient into a PARCOR coefficient and a code corresponding to the PARCOR coefficient as a linear prediction coefficient code.

  In this way, the linear prediction analysis unit 421 generates a pseudo correlation function signal sequence obtained by performing an inverse Fourier transform assuming that the absolute value of the absolute value of the frequency domain sample sequence, which is an MDCT coefficient sequence, is a power spectrum, for example. On the other hand, linear prediction analysis is performed to generate a coefficient that can be converted into a linear prediction coefficient (step C421).

<Non-smoothed Amplitude Spectrum Envelope Sequence Generation Unit 422>
Quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p generated by the linear prediction analysis unit 421 are input to the unsmoothed amplitude spectrum envelope sequence generation unit 422.

Textured amplitude spectral envelope sequence generation unit 422, the quantized linear prediction coefficient _{^ β 1, ^ β 2,} ..., ^ β is the sequence of the amplitude spectrum envelope corresponding to _p textured amplitude spectral envelope sequence ^ H ( 0), ^ H (1), ..., ^ H (N-1) are generated.

  The generated non-smoothed amplitude spectrum envelope sequences ^ H (0), ^ H (1),..., ^ H (N-1) are output to the whitened spectrum sequence generation unit 43.

The unsmoothed amplitude spectrum envelope sequence generation unit 422 uses the quantized linear prediction coefficients ^ β ₁ , ^ β ₂ ,..., ^ Β _p to generate the unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H ( 1),…, ^ H (N-1) as unsmoothed amplitude spectrum envelope sequence defined by equation (C2) ^ H (0), ^ H (1),…, ^ H (N-1) Is generated.

In this way, the unsmoothed amplitude spectrum envelope sequence generation unit 422 performs linear prediction analysis on the unsmoothed spectrum envelope sequence that is a sequence obtained by raising the amplitude spectrum envelope sequence corresponding to the pseudo correlation function signal sequence to the 1 / η ₀ power. The spectral envelope is estimated by obtaining the coefficient based on the coefficient that can be converted into the linear prediction coefficient generated by the unit 421 (step C422).

<Whitening spectrum series generation unit 43>
The whitened spectrum sequence generation unit 43 includes an MDCT coefficient sequence X (0), X (1),..., X (N-1) obtained by the frequency domain conversion unit 41 and a non-smoothed amplitude spectrum envelope generation unit 422. The generated non-smoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1), ..., ^ H (N-1) is input.

The whitened spectrum sequence generation unit 43 converts each coefficient of the MDCT coefficient sequence X (0), X (1),..., X (N-1) into a corresponding non-smoothed amplitude spectrum envelope sequence ^ H (0), By dividing each value of ^ H (1), ..., ^ H (N-1), the whitened spectrum series X _W (0), X _W (1), ..., X _W (N-1) Generate.

The generated whitening spectrum series X _W (0), X _W (1),..., X _W (N−1) are output to the parameter acquisition unit 44.

For example, the whitening spectrum sequence generation unit 43 sets k = 0, 1,..., N−1 as the coefficients X (()) of the MDCT coefficient sequence X (0), X (1),. By dividing k) by the unsmoothed amplitude spectrum envelope sequence ^ H (0), ^ H (1),…, ^ H (N-1) values ^ H (k), the whitened spectrum sequence X Each value X _W (k) of _W (0), X _W (1),..., X _W (N−1) is generated. That is, X _W (k) = X (k) / ^ H (k) where k = 0, 1,..., N−1.

  In this way, the whitened spectrum sequence generation unit 43 obtains a whitened spectrum sequence that is a sequence obtained by dividing a frequency domain sample sequence that is an MDCT coefficient sequence, for example, by a spectrum envelope that is an unsmoothed amplitude spectrum envelope sequence, for example ( Step C43).

<Parameter acquisition unit 44>
The parameter acquisition unit 44 receives the whitened spectrum series X _W (0), X _W (1),..., X _W (N−1) generated by the whitened spectrum series generating unit 43.

The parameter acquisition unit 44 approximates the histogram of the whitened spectrum series X _W (0), X _W (1),..., X _W (N−1) with the generalized Gaussian distribution having the parameter η as a shape parameter. Is obtained (step C44). In other words, the parameter acquisition unit 44 is a distribution of histograms in which the generalized Gaussian distribution having the parameter η as a shape parameter is a whitened spectrum series X _W (0), X _W (1), ..., X _W (N-1). The parameter η that is close to is determined.

  A generalized Gaussian distribution with the parameter η as a shape parameter is defined as follows, for example. Γ is a gamma function.

  By changing the shape parameter η, the generalized Gaussian distribution can represent various distributions such as a Laplace distribution when η = 1 and a Gaussian distribution when η = 2 as shown in FIG. It can be done. φ is a parameter corresponding to the variance.

Here, η obtained by the parameter acquisition unit 44 is defined by the following equation. F ⁻¹ is an inverse function of the function F. This equation is derived by the so-called moment method.

When the inverse function F ⁻¹ is formulated, the parameter acquisition unit 44 inputs the value of m ₁ / ((m ₂ ) ^1/2 ) into the formulated inverse function F ⁻¹ . The parameter η can be obtained by calculating the output value.

If the inverse function F ⁻¹ is not formulated, the parameter acquisition unit 44 calculates, for example, the first method or the second method described below in order to calculate the value of η defined by the equation (C3). The parameter η may be obtained by

A first method for obtaining the parameter η will be described. In the first method, the parameter acquisition unit 44 calculates m ₁ / ((m ₂ ) ^1/2 ) based on the whitened spectrum sequence, and a plurality of different F prepared in advance corresponding to η. Η corresponding to F (η) closest to the calculated m ₁ / ((m ₂ ) ^1/2 ) is obtained with reference to the pair of (η).

A plurality of different pairs of F (η) corresponding to η prepared in advance are stored in advance in the storage unit 441 of the parameter acquisition unit 44. The parameter acquisition unit 44 refers to the storage unit 441, finds F (η) closest to the calculated m ₁ / ((m ₂ ) ^1/2 ), and stores η corresponding to the found F (η). Read from the unit 441 and output.

The calculated m ₁ / closest to ^{_{((m 2) 1/2) F}} (η) , the absolute value of the difference between the calculated _{m 1 / ((m 2)} 1/2) is smallest F (η).

A second method for obtaining the parameter η will be described. In the second method, the approximate curve function of the inverse function F ⁻¹ is set to ~ F ⁻¹ , and the parameter acquisition unit 44 calculates m ₁ / ((m ₂ ) ^1/2 ) based on the whitened spectrum series. Then, η is obtained by calculating an output value when m ₁ / ((m ₂ ) ^1/2 ) calculated in the approximate curve function ~ F ⁻¹ is input.

The parameter determination unit 27 ′ may obtain the parameter η by loop processing. That is, the parameter determination unit 27 ′ sets the parameter η obtained by the parameter acquisition unit 44 as the parameter η ₀ determined by a predetermined method, and performs processing by the spectrum envelope estimation unit 42, the whitened spectrum sequence generation unit 43, and the parameter acquisition unit 44. May be performed once more.

In this case, for example, as indicated by a broken line in FIG. 14, the parameter η obtained by the parameter acquisition unit 44 is output to the spectrum envelope estimation unit 42. The spectrum envelope estimation unit 42 estimates the spectrum envelope by performing the same process as described above using η obtained by the parameter acquisition unit 44 as the parameter η ₀ . Based on the newly estimated spectrum envelope, the whitened spectrum sequence generation unit 43 generates a whitened spectrum sequence by performing the same process as described above. The parameter acquisition unit 44 performs a process similar to the process described above based on the newly generated whitened spectrum sequence to obtain the parameter η.

  For example, the processing of the spectrum envelope estimation unit 42, the whitened spectrum series generation unit 43, and the parameter acquisition unit 44 may be further performed by τ times that is a predetermined number of times. τ is a predetermined positive integer, for example, τ = 1 or τ = 2.

  Further, the spectrum envelope estimation unit 42 performs the spectrum envelope estimation unit 42, the whitened spectrum sequence generation unit 43, and the parameter until the absolute value of the difference between the parameter η obtained this time and the parameter η obtained last time is equal to or less than a predetermined threshold. You may repeat the process of the acquisition part 44. FIG.

(Decryption)
Since the decoding apparatus and method of the second embodiment are the same as those of the first embodiment, redundant description is omitted.

[Modifications, etc.]
When the linear prediction analysis unit 22 and the unsmoothed amplitude spectrum envelope sequence generation unit 23 are regarded as one spectrum envelope estimation unit 2A, the spectrum envelope estimation unit 2A is a frequency domain that is, for example, an MDCT coefficient sequence corresponding to a time series signal. It can be said that the spectrum envelope (unsmoothed amplitude spectrum envelope sequence) is estimated by regarding the absolute value of the sample string to the power of η as the power spectrum. Here, “considered as a power spectrum” means to use a power of η where a power spectrum is normally used.

  In this case, the linear prediction analysis unit 22 of the spectrum envelope estimation unit 2A performs, for example, a pseudo correlation obtained by performing an inverse Fourier transform in which the absolute value of η power of a frequency domain sample sequence that is an MDCT coefficient sequence is regarded as a power spectrum. It can be said that a coefficient that can be converted into a linear prediction coefficient is obtained by performing linear prediction analysis on the function signal sequence. Further, the unsmoothed amplitude spectrum envelope sequence generation unit 23 of the spectrum envelope estimation unit 2A converts the amplitude spectrum envelope sequence corresponding to the coefficient that can be converted into the linear prediction coefficient obtained by the linear prediction analysis unit 22 to the 1 / ηth power. It can be said that the spectral envelope is estimated by obtaining the non-smoothed spectral envelope sequence which is the obtained sequence.

  Further, if the smoothed amplitude spectrum envelope sequence generation unit 24, the envelope normalization unit 25, and the encoding unit 26 are regarded as one encoding unit 2B, the encoding unit 2B is a spectrum estimated by the spectrum envelope estimation unit 2A. It can be said that the coding in which the bit allocation is substantially changed based on the envelope (non-smoothed amplitude spectrum envelope sequence) is performed for each coefficient of the frequency domain sample sequence corresponding to the time series signal, for example, an MDCT coefficient sequence. .

  When the decoding unit 34 and the envelope denormalization unit 35 are regarded as a single decoding unit 3A, the decoding unit 3A is configured to input the input integer signal code according to bit allocation that substantially changes based on the non-smoothed spectrum envelope sequence. It can be said that the frequency domain sample sequence corresponding to the time-series signal is obtained by performing the decoding.

  The encoding unit 2B may perform encoding processing other than the arithmetic encoding described above as long as it performs encoding in which the bit allocation substantially changes based on the spectrum envelope (unsmoothed amplitude spectrum envelope sequence). Good. In this case, the decoding unit 3A performs a decoding process corresponding to the encoding process performed by the encoding unit 2B.

  For example, the encoding unit 2B may perform Golomb-Rice encoding on the frequency domain sample sequence using the determined Rice parameter determined based on the spectrum envelope (unsmoothed amplitude spectrum envelope sequence). In this case, the decoding unit 3A may perform Golomb-Rice decoding using the determined Rice parameter determined based on the spectrum envelope (unsmoothed amplitude spectrum envelope sequence).

  In the first embodiment, the encoding device may not perform the encoding process to the end when determining the parameter η. In other words, the parameter determination unit 27 may determine the parameter η based on the estimated code amount. In this case, the encoding unit 2B uses each of the plurality of parameters η to estimate the code obtained by the same encoding process as described above for the frequency domain sample sequence corresponding to the time-series signal in the same predetermined time interval. Get quantity. The parameter determination unit 27 selects one of a plurality of parameters η based on the obtained estimated code amount. For example, the parameter η having the smallest estimated code amount is selected. The encoding unit 2B obtains and outputs a code by performing the same encoding process as described above using the selected parameter η.

  The same predetermined parameter η may be predetermined for the encoding device and the decoding device. In this case, the parameter determination units 27 and 27 'may not be provided in the encoding device and the decoding device. Each unit other than the parameter determination units 27 and 27 ′ of the encoding device and the decoding device performs, for example, the same encoding process and decoding process as in the second embodiment based on the predetermined parameter η.

  The encoding apparatus may further include a dividing unit 28 indicated by a broken line in FIG. 4 or FIG. Based on the frequency domain sample sequence that is, for example, the MDCT coefficient sequence generated by the frequency domain transform unit 21, the dividing unit 28 includes a first frequency domain sample sequence that includes samples corresponding to periodic components of the frequency domain sample sequence, Generating a second frequency domain sample sequence composed of samples other than the sample corresponding to the periodic component of the frequency domain sample sequence, and outputting information representing the sample corresponding to the periodic component to the decoding apparatus as auxiliary information .

  In other words, the first frequency domain sample sequence is a sample sequence composed of samples corresponding to the peaks of the frequency domain sample sequence, and the second frequency domain sample sequence is a sample corresponding to the valleys of the frequency domain sample sequence. It is a sample sequence composed of

  For example, one or a plurality of consecutive samples including a periodicity of a time-series signal corresponding to a frequency domain sample sequence in a frequency domain sample sequence or a sample corresponding to a fundamental frequency, and a frequency domain in a frequency domain sample sequence The first frequency domain includes a sample sequence composed of all or part of one or a plurality of consecutive samples including samples corresponding to the periodicity of the time series signal corresponding to the sample sequence or an integer multiple of the fundamental frequency A second frequency domain sample sequence is generated from the sample sequence and a sample sequence composed of samples not included in the first frequency domain sample sequence of the frequency domain sample sequences. The generation of the first frequency domain sample sequence and the second frequency domain sample sequence can be performed using a method described in International Publication WO2012 / 046685.

  The linear prediction analysis unit 22, the unsmoothed amplitude spectrum envelope sequence generation unit 23, the smoothed amplitude spectrum envelope sequence generation unit 24, the envelope normalization unit 25, the encoding unit 26, and the parameter determination unit 27 include a first frequency domain sample sequence And about each of a 2nd frequency domain sample sequence, the encoding process demonstrated in 1st embodiment or 2nd embodiment is performed, and a code | symbol is produced | generated. That is, for example, when arithmetic coding is performed, a parameter code, a linear prediction coefficient code, an integer signal code, and a gain code corresponding to the first frequency domain sample sequence are generated, and parameters corresponding to the second frequency domain sample sequence are generated. A code, a linear prediction coefficient code, an integer signal code, and a gain code are generated.

  Thus, encoding can be performed more efficiently by encoding each of the first frequency domain sample sequence and the second frequency domain sample sequence.

  In this case, the decoding apparatus may further include a combining unit 38 indicated by a broken line in FIG. The decoding apparatus performs the decoding process described in the first embodiment or the second embodiment based on a code (for example, a parameter code, a linear prediction coefficient code, an integer signal code, and a gain code) corresponding to the first frequency domain sample sequence. To obtain a decoded first frequency domain sample sequence. Further, the decoding apparatus has been described in the first embodiment or the second embodiment based on a code (for example, a parameter code, a linear prediction coefficient code, an integer signal code, and a gain code) corresponding to the second frequency domain sample sequence. A decoding process is performed to obtain a decoded second frequency domain sample sequence. The combining unit 38 appropriately combines the decoded first frequency domain sample sequence and the decoded second frequency domain sample sequence using the input auxiliary information, for example, by decoding MDCT coefficient sequences ^ X (0), ^ X (1 ), ..., ^ X (N-1) is obtained as a decoded frequency domain sample sequence. The time domain transform unit obtains a time series signal by transforming the decoded frequency domain sample sequence into the time domain. The combination using the auxiliary information can be performed using the method described in International Publication No. WO2012 / 046685.

  When the bit rate is low or when it is desired to further reduce the code amount, the encoding device encodes only the first frequency domain sample sequence, generates only the code corresponding to the first frequency domain sample sequence, The code corresponding to the two frequency domain sample sequences is not generated, and the decoding apparatus uses the first frequency domain sample sequence obtained from the code and the second frequency domain sample sequence with the sample value of 0 as a decoded frequency domain sample. A column may be obtained.

  Further, the linear prediction analysis unit 22, the non-smoothed amplitude spectrum envelope sequence generation unit 23, the smoothed amplitude spectrum envelope sequence generation unit 24, the envelope normalization unit 25, the encoding unit 26, and the parameter determination unit 27 are included in the first frequency domain. Even if the sample sequence after the rearrangement is a sample sequence obtained by combining the sample sequence and the second frequency domain sample sequence, the encoding process described in the first embodiment or the second embodiment is performed to generate a code. Good. For example, when arithmetic coding is performed, a parameter code, a linear prediction coefficient code, an integer signal code, and a gain code corresponding to the rearranged sample sequence are generated.

  In this way, encoding can be performed more efficiently by encoding the rearranged sample sequence.

  In this case, the decoding device performs the decoding process described in the first embodiment or the second embodiment, obtains a sample string after decoding rearrangement, and uses the input auxiliary information to obtain the sample string after decoding rearrangement, The first frequency domain sample sequence and the second frequency domain sample sequence are rearranged according to the rule corresponding to the rule generated by the encoder, for example, the decoded MDCT coefficient sequence ^ X (0), ^ X (1),. A decoded frequency domain sample sequence which is ^ X (N-1) is obtained. The time domain transform unit 36 transforms the decoded frequency domain sample sequence into the time domain to obtain a time series signal. Rearrangement using auxiliary information can be performed using the method described in International Publication No. WO2012 / 046685.

  Further, the encoding device (1) performs a coding process on the frequency domain sample sequence to generate a code, and (2) performs a coding process on each of the first frequency domain sample sequence and the second frequency domain sample sequence. A method of generating a code, (3) a method of generating a code by encoding only the first frequency domain sample sequence, and (4) obtained by combining the first frequency domain sample sequence and the second frequency domain sample sequence. Any method may be selected for each frame from among the methods for generating a code by performing an encoding process on the rearranged sample sequence that is a sample sequence. In this case, the encoding apparatus also outputs a code indicating which method (1) to (4) is selected, and the decoding apparatus corresponds to any of the above methods according to the code input for each frame. Perform decryption.

  The parameter determination unit 27 of the encoding device and the parameter decoding unit 37 of the decoding device may store parameter η candidates corresponding to the above methods (1) to (4). . Similarly, the linear prediction analysis unit 22 of the encoding device and the linear prediction coefficient decoding unit 31 of the decoding device include quantized linear prediction coefficient candidates corresponding to the methods (1) to (4), and Decoded linear prediction coefficient candidates may be stored.

  The unsmoothed amplitude spectrum envelope sequence generation unit 23 and the unsmoothed amplitude spectrum envelope sequence generation unit 422 are, for example, MDCT coefficient sequences ^ X (0), ^ X (1),..., ^ X (N-1). The periodic integrated envelope sequence may be generated by modifying the spectrum envelope sequence (unsmoothed amplitude spectrum envelope sequence) based on the periodic component of the frequency domain sample sequence. Similarly, the non-smoothed amplitude spectrum envelope sequence generator 32 generates a period of a decoded frequency domain sample sequence that is, for example, a decoded MDCT coefficient sequence ^ X (0), ^ X (1), ..., ^ X (N-1). The periodic integrated envelope sequence may be generated by modifying the spectral envelope sequence (non-smoothed amplitude spectral envelope sequence) based on the sex component. In this case, the dispersion parameter determination unit 268, the decoding unit 34, and the whitened spectrum sequence generation unit 43 of the encoding unit 26 use the periodic integrated envelope sequence instead of the spectrum envelope sequence (unsmoothed amplitude spectrum envelope sequence). The same processing as above is performed. Since the periodic integrated envelope sequence has good approximation accuracy near the peak due to the pitch period of the time series signal, the use of the periodic integrated envelope sequence can increase the coding efficiency.

  For example, the greater the period of the frequency domain sample sequence, the greater the periodicity of at least the integer multiple of the frequency domain sample sequence in the spectrum envelope sequence and the value of samples in the vicinity of the integer multiple of the period. The integrated envelope series. In addition, as the degree of periodicity of the time-series signal is larger, a sequence obtained by greatly changing the value of a sample in the vicinity of at least an integer multiple of the period of the frequency domain sample sequence and the integer multiple of the period of the spectrum envelope sequence. It may be a periodic integrated envelope sequence. Further, as the period of the frequency domain sample sequence is larger, a sequence obtained by changing the values of many samples in the vicinity of an integer multiple of the period of the frequency domain sample sequence in the spectrum envelope sequence may be used as the periodic integrated envelope sequence. .

Furthermore, N and U are positive integers, T is the interval of the periodic component of the frequency domain sample sequence, L is the number of digits after the decimal point of interval T, v is an integer of 1 or more, and floor (·) is the decimal point A function that returns an integer value by rounding down, Round (·) rounds off the first decimal place, and returns an integer value, T ′ = T × 2 ^L , ^ H [0],…, ^ H [N-1 ] Is a spectral envelope sequence, δ is a value that determines the mixing ratio of spectral envelope ^ H [n] and periodic envelope P [k],
(U × T ′) / 2 ^L −v−1 ≦ k ≦ (U × T ′) / 2 ^L + v−1
For an integer k in the range

P [1], ..., P [N] is obtained as shown below, and the periodicity defined by the following equation using the obtained periodic envelope series P [1], ..., P [N] The integrated envelope sequence ^ H _M [1], ..., ^ H _M [N] may be obtained. h and PD may be predetermined values other than the above example.

  Δ that is a value that determines the mixing ratio of the spectral envelope ^ H [n] and the periodic envelope P [k] may be determined in advance by the encoding device and the decoding device, or may be a value of δ determined by the encoding device. A code indicating information may be generated and output to the decoding device. In the latter case, the decoding apparatus obtains δ by decoding the code indicating the input information of δ. The non-smoothed amplitude spectrum envelope sequence generation unit 32 of the decoding device can obtain the same periodic integrated envelope sequence as the periodic integrated envelope sequence generated by the encoding device by using the obtained δ.

  The processes described above are not only executed in chronological order according to the order of description, but may be executed in parallel or individually as required by the processing capability of the apparatus that executes the processes.

  Various processes in each method may be realized by a computer. In that case, the processing content of each method is described by a program. Then, by executing this program on a computer, various processes in each method are realized on the computer.

  The program describing the processing contents can be recorded on a computer-readable recording medium. As the computer-readable recording medium, for example, any recording medium such as a magnetic recording device, an optical disk, a magneto-optical recording medium, and a semiconductor memory may be used.

  The program is distributed by selling, transferring, or lending a portable recording medium such as a DVD or CD-ROM in which the program is recorded. Further, the program may be distributed by storing the program in a storage device of the server computer and transferring the program from the server computer to another computer via a network.

  A computer that executes such a program first stores, for example, a program recorded on a portable recording medium or a program transferred from a server computer in its storage unit. When executing the process, this computer reads the program stored in its own storage unit and executes the process according to the read program. As another embodiment of this program, a computer may read a program directly from a portable recording medium and execute processing according to the program. Further, each time a program is transferred from the server computer to the computer, processing according to the received program may be executed sequentially. Also, the program is not transferred from the server computer to the computer, and the above-described processing is executed by a so-called ASP (Application Service Provider) type service that realizes the processing function only by the execution instruction and result acquisition. It is good. Note that the program includes information provided for processing by the electronic computer and equivalent to the program (data that is not a direct command to the computer but has a property that defines the processing of the computer).

  In addition, although each device is configured by executing a predetermined program on a computer, at least a part of these processing contents may be realized by hardware.

Claims (7)

An encoding device that encodes a time-series signal in a frequency domain,
η is a positive number, and linear prediction is performed for a pseudo correlation function signal sequence obtained by performing an inverse Fourier transform in which the absolute value of the frequency domain sample sequence corresponding to the time series signal is assumed to be a power spectrum as a power spectrum. A linear prediction analysis unit for performing analysis and obtaining energy of coefficients and prediction residuals that can be converted into linear prediction coefficients;
A non-smoothed spectrum envelope sequence generation unit that estimates a spectrum envelope by obtaining a non-smoothed spectrum envelope sequence that is a sequence obtained by raising the amplitude spectrum envelope sequence corresponding to a coefficient that can be converted into the linear prediction coefficient to the power of 1 / η. When,
A frequency domain sample corresponding to the time-series signal encoding with a bit allocation substantially changed based on the estimated spectral envelope and information according to the model parameter, using the square root of the energy of the prediction residual as a model parameter An encoding unit for each coefficient of the sequence;
A model parameter information encoding unit for obtaining a model parameter information code representing information according to the model parameter;
An encoding device including: The encoding device according to claim 1, wherein
The encoding unit performs encoding on a sequence by an integer part as a result of dividing each coefficient of a frequency domain sample sequence corresponding to the time series signal by a global gain,
The model parameter is σ, m is a positive number, and the information according to the model parameter is a value of a ratio between σ ^m and the global gain.
Encoding device. A decoding device that obtains a frequency domain sample sequence corresponding to a time-series signal by decoding in the frequency domain,
A linear prediction coefficient decoding unit that obtains a coefficient that can be converted into a linear prediction coefficient by decoding the input linear prediction coefficient code;
A model parameter information decoding unit that decodes the input model parameter code to obtain information according to the model parameter;
Generation of a non-smoothed spectral envelope sequence that obtains a non-smoothed spectral envelope sequence that is a sequence obtained by raising the amplitude spectral envelope sequence corresponding to the coefficient that can be converted to the linear prediction coefficient to the 1 / η power, where η is a predetermined positive number And
The frequency domain sample sequence corresponding to the time-series signal is obtained by decoding the input integer signal code according to the bit allocation that substantially changes based on the non-smoothed spectral envelope sequence and the information corresponding to the model parameter. A decoding unit to obtain;
A decoding device. The decoding device according to claim 3,
The decoding unit obtains a frequency domain sample sequence corresponding to the time series signal by multiplying each coefficient of the coefficient series obtained by decoding the input integer signal code by a global gain,
The model parameter is σ, m is a positive number, and the information according to the model parameter is a value of a ratio between σ ^m and the global gain.
Decoding device. An encoding method for encoding a time series signal in a frequency domain,
η is a positive number, and linear prediction is performed for a pseudo correlation function signal sequence obtained by performing an inverse Fourier transform in which the absolute value of the frequency domain sample sequence corresponding to the time series signal is assumed to be a power spectrum as a power spectrum. A linear prediction analysis step for obtaining an energy of a coefficient and a prediction residual that can be analyzed and converted into a linear prediction coefficient;
A non-smoothed spectrum envelope sequence generating step for estimating a spectrum envelope by obtaining a non-smoothed spectrum envelope sequence that is a sequence obtained by raising the amplitude spectrum envelope sequence corresponding to a coefficient that can be converted into the linear prediction coefficient to the 1 / η power. When,
A frequency domain sample corresponding to the time-series signal encoding with a bit allocation substantially changed based on the estimated spectral envelope and information according to the model parameter, using the square root of the energy of the prediction residual as a model parameter An encoding step for each coefficient of the sequence;
A model parameter information encoding step for obtaining a model parameter information code representing information according to the model parameter;
An encoding method including: A decoding method for obtaining a frequency domain sample sequence corresponding to a time-series signal by decoding in the frequency domain,
A linear prediction coefficient decoding step for obtaining a coefficient that can be converted into a linear prediction coefficient by decoding the input linear prediction coefficient code;
A model parameter information decoding step of decoding the input model parameter code to obtain information according to the model parameter;
Generation of a non-smoothed spectral envelope sequence that obtains a non-smoothed spectral envelope sequence that is a sequence obtained by raising the amplitude spectral envelope sequence corresponding to the coefficient that can be converted to the linear prediction coefficient to the 1 / η power, where η is a predetermined positive number Steps,
The frequency domain sample sequence corresponding to the time-series signal is obtained by decoding the input integer signal code according to the bit allocation that substantially changes based on the non-smoothed spectral envelope sequence and the information corresponding to the model parameter. Obtaining a decryption step;
A decoding method including:   A program for causing a computer to function as each unit of the encoding device according to claim 1 or 2 or the decoding device according to claim 3 or 4.

Images