WO2015166693A1 - 符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、記録媒体 - Google Patents

符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、記録媒体 Download PDF

Info

Publication number
WO2015166693A1
WO2015166693A1 PCT/JP2015/054717 JP2015054717W WO2015166693A1 WO 2015166693 A1 WO2015166693 A1 WO 2015166693A1 JP 2015054717 W JP2015054717 W JP 2015054717W WO 2015166693 A1 WO2015166693 A1 WO 2015166693A1
Authority
WO
WIPO (PCT)
Prior art keywords
variable length
sequence
decoding
encoding
frequency domain
Prior art date
Application number
PCT/JP2015/054717
Other languages
English (en)
French (fr)
Inventor
守谷　健弘
優 鎌本
登 原田
Original Assignee
日本電信電話株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority to ES15785387T priority Critical patent/ES2732859T3/es
Priority to JP2016515878A priority patent/JP6276845B2/ja
Application filed by 日本電信電話株式会社 filed Critical 日本電信電話株式会社
Priority to CN202110083487.4A priority patent/CN112820305B/zh
Priority to KR1020187006878A priority patent/KR101861787B1/ko
Priority to EP20167431.4A priority patent/EP3703051B1/en
Priority to KR1020187006872A priority patent/KR101861781B1/ko
Priority to EP19157754.3A priority patent/EP3509063B1/en
Priority to PL19157754T priority patent/PL3509063T3/pl
Priority to KR1020167029792A priority patent/KR101839016B1/ko
Priority to US15/304,733 priority patent/US10199046B2/en
Priority to EP20167433.0A priority patent/EP3696812B1/en
Priority to PL20167433T priority patent/PL3696812T3/pl
Priority to PL15785387T priority patent/PL3139380T3/pl
Priority to EP15785387.0A priority patent/EP3139380B1/en
Priority to CN201580022887.7A priority patent/CN106663437B/zh
Priority to PL20167431T priority patent/PL3703051T3/pl
Publication of WO2015166693A1 publication Critical patent/WO2015166693A1/ja
Priority to US16/199,370 priority patent/US10607616B2/en
Priority to US16/199,387 priority patent/US10629214B2/en
Priority to US16/788,539 priority patent/US11164589B2/en

Links

Images

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/02—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/06—Determination or coding of the spectral characteristics, e.g. of the short-term prediction coefficients
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/08—Determination or coding of the excitation function; Determination or coding of the long-term prediction parameters
    • G10L19/12—Determination or coding of the excitation function; Determination or coding of the long-term prediction parameters the excitation function being a code excitation, e.g. in code excited linear prediction [CELP] vocoders
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/90—Pitch determination of speech signals
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
    • H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
    • H03M7/40—Conversion to or from variable length codes, e.g. Shannon-Fano code, Huffman code, Morse code

Definitions

  • the present invention relates to an encoding device, a decoding device, an encoding method, a decoding method, an encoding program, a decoding program, and a recording medium that encode and decode an acoustic signal using a spectral envelope of the acoustic signal.
  • Adaptive coding for orthogonal transform coefficients such as DFT (Discrete Fourier Transform) and MDCT (Modified Discrete Cosine Transform) is known as a coding method for low-bit (for example, about 10 kbit / s to 20 kbit / s) speech and acoustic signals. It has been.
  • TCX transform coded excitation
  • Non-Patent Document 1 a coefficient sequence X [1],..., X [N] that is a frequency domain representation of an input sound signal.
  • a sequence normalized coefficient sequence X N [1],..., X N [N]
  • N in [] is a positive integer.
  • the amplitude spectrum envelope is calculated by the following procedure.
  • Step 1 A linear prediction coefficient ⁇ 1 ,..., ⁇ P is obtained by performing linear prediction analysis on the input time domain acoustic digital signal (hereinafter referred to as input acoustic signal) in units of frames that are predetermined time intervals.
  • P is a positive integer indicating the predicted order.
  • Step2 linear prediction coefficient alpha 1, ..., a alpha P quantizes quantized linear prediction coefficient ⁇ alpha 1, ..., determine the ⁇ alpha P.
  • the quantized linear prediction coefficients ⁇ ⁇ 1 ,..., ⁇ ⁇ P the quantized linear prediction coefficients ⁇ ⁇ 1 ,..., ⁇ ⁇ P .
  • W [n] of the amplitude spectrum envelope series can be obtained by Expression (2).
  • n is an integer of 1 ⁇ n ⁇ N
  • exp (•) is an exponential function with the Napier number as the base
  • j is an imaginary unit
  • ⁇ is the amplitude of the prediction residual signal.
  • a symbol written without a square bracket on the right shoulder represents a power operation. That is, ⁇ 2 represents the square of ⁇ .
  • ⁇ 2 represents the square of ⁇ .
  • the symbols “ ⁇ ”, “ ⁇ ”, etc. used in the sentence should be described immediately above the character that immediately follows, but are described immediately before the character due to restrictions on text notation. In the mathematical expression, these symbols are described in their original positions, that is, directly above the characters.
  • Non-Patent Document 1 the “code corresponding to the spectral envelope” sent to the decoding side is the “code corresponding to the linear prediction coefficient”, and the amount of code can be small.
  • the spectral envelope information obtained by the linear prediction coefficient may deteriorate the approximation accuracy in the vicinity of the peak due to the pitch period of the input acoustic signal, which causes the normalized coefficient sequence to be variable-length encoded. Encoding efficiency may be reduced.
  • the present invention provides an encoding device and a decoding device that can reproduce an envelope sequence in a frequency domain with good approximation accuracy near a peak due to the pitch period of an acoustic signal with a small code amount.
  • the encoding apparatus of the present invention includes a periodic integrated envelope generation unit and a variable length encoding unit.
  • the periodic integrated envelope generation unit includes a spectral envelope sequence that is a frequency domain sequence corresponding to a linear prediction coefficient code obtained from an input acoustic signal in a predetermined time interval, and a frequency domain corresponding to a periodic code obtained from the input acoustic signal. And a periodic integrated envelope sequence that is a frequency domain sequence based on the period.
  • the variable length encoding unit encodes a frequency domain sequence derived from the input acoustic signal on the assumption that the amplitude of the input acoustic signal is larger as the frequency of the periodic integrated envelope sequence is larger.
  • the decoding device of the present invention includes a periodicity integrated envelope generation unit and a variable length decoding unit.
  • the periodic integrated envelope generation unit is a periodic integrated envelope sequence that is a frequency domain sequence based on a spectrum envelope sequence that is a frequency domain sequence corresponding to a linear prediction coefficient code and a frequency domain period that corresponds to a periodic code. Is generated.
  • the variable length decoding unit obtains a frequency domain sequence by decoding the variable length code on the premise that the amplitude of the acoustic signal increases as the frequency of the periodic integrated envelope sequence increases.
  • encoding and decoding that can reproduce an envelope sequence in a frequency domain with good approximation accuracy near the peak due to the pitch period of the input acoustic signal with a small amount of code.
  • the encoding efficiency of the acoustic signal is good.
  • FIG. 6 is a diagram illustrating a functional configuration example of an encoding apparatus according to a second embodiment.
  • FIG. 10 is a diagram illustrating a processing flow of the encoding apparatus according to the second embodiment.
  • FIG. 6 is a diagram illustrating a functional configuration example of a decoding device according to a second embodiment.
  • FIG. 10 is a diagram illustrating a processing flow of the decoding apparatus according to the second embodiment.
  • FIG. 10 is a diagram illustrating a functional configuration example of an encoding apparatus according to a third embodiment.
  • FIG. 10 is a diagram illustrating a processing flow of the encoding apparatus according to the third embodiment.
  • FIG. 10 is a diagram illustrating a functional configuration example of a decoding device according to a third embodiment.
  • FIG. 10 is a diagram illustrating a processing flow of the decoding apparatus according to the third embodiment.
  • FIG. 1 shows a functional configuration example of the periodic integrated envelope sequence generation device of the present invention
  • FIG. 2 shows a processing flow of the periodic integrated envelope sequence generation device of the present invention
  • the periodic integrated envelope sequence generating apparatus 100 includes a spectrum envelope sequence calculating unit 120, a frequency domain converting unit 110, a periodicity analyzing unit 130, a periodic envelope sequence generating unit 140, and a periodic integrated envelope generating unit 150.
  • An acoustic digital signal in the time domain is set as an input acoustic signal x (t), and a periodic integrated envelope sequence is generated by modifying the amplitude spectrum envelope sequence based on the frequency component of the coefficient sequence.
  • the spectrum envelope sequence calculation unit 120 calculates the amplitude spectrum envelope sequence W [1],..., W [N] of the input acoustic signal based on the time domain linear prediction of the input acoustic signal x (t) (S120). However, N is a positive integer.
  • the spectrum envelope sequence calculation unit 120 is the same as that in the prior art, and may be calculated by the following procedure.
  • Step 1 Linear prediction coefficients ⁇ 1 ,..., ⁇ P are obtained by performing linear prediction analysis on the input acoustic signal in units of frames that are predetermined time intervals.
  • P is a positive integer indicating the predicted order.
  • the input acoustic signal x (t) at the time t becomes the past value x (t ⁇ 1),. tP) and the prediction residuals e (t) and the linear prediction coefficients alpha 1, ..., represented by the formula (1) by alpha p.
  • each value W of the amplitude spectral envelope sequence [n] is the linear prediction coefficients alpha 1, ..., alpha quantized linear prediction coefficients corresponding to P ⁇ alpha 1, ..., using ⁇ alpha P Equation (2) Can be obtained.
  • each value W [n] of the amplitude spectrum envelope sequence can be obtained by an equation in which ⁇ ⁇ p in Equation (2) is replaced by ⁇ p using linear prediction coefficients ⁇ 1 ,..., ⁇ P.
  • the frequency domain conversion unit 110 converts the input time domain input acoustic signal into N frequency coefficient sequences X [1],..., X [N] in units of frames that are predetermined time intervals. Output (S110).
  • the conversion to the frequency domain may be performed by a method such as MDCT (Modified Discrete Cosine Transform) or DFT (Discrete Fourier Transform).
  • the periodicity analysis unit 130 receives the coefficient sequence X [1],..., X [N], obtains the cycle T of the coefficient sequence X [1],. S130).
  • the period T is a frequency interval coefficient sequence derived from the input acoustic signal, for example, the interval between the components having periodicity of the coefficient sequence X [1],. Information).
  • the period T may be expressed as the interval T, but only the difference in expression is the same.
  • T is a positive value, and may be an integer or a decimal (for example, 5.0, 5.25, 5.5, 5.75).
  • the periodicity analysis unit 130 may obtain and output an index S indicating the degree of periodicity by inputting the coefficient sequence X [1],..., X [N] as necessary.
  • the index S is an index indicating the degree of periodicity of the frequency domain sample sequence.
  • the periodicity analysis unit 130 may obtain the period T by obtaining the time domain period from the time domain input acoustic signal and converting the obtained time domain period into the frequency domain period. Alternatively, a constant multiple of a time domain period converted to a frequency domain period or a value in the vicinity thereof may be obtained as the period T. Similarly, the periodicity analysis unit 130 is an index indicating the degree of periodicity based on, for example, the magnitude of correlation between signal sequences shifted in time by the period of the time domain from the time domain input acoustic signal. S may be obtained.
  • the periodic envelope sequence generation unit 140 receives the interval T and outputs a periodic envelope sequence P [1],..., P [N] (S140).
  • the periodic envelope sequence P [1],..., P [N] is a discrete sequence in the frequency domain having a peak with a period due to the pitch period, that is, a discrete series corresponding to the harmonic model.
  • FIG. 3 shows an example of the periodic envelope series P [1],..., P [N].
  • the periodic envelope sequence P [1],..., P [N] is divided into an index that is an integer value in the vicinity of an integer multiple of the interval T and a predetermined number of indexes before and after that, as in the waveform shown in FIG.
  • An index that is an integer value in the vicinity of an integer multiple of the interval T periodically has a maximum value (peak), and a value of P [n] corresponding to a predetermined number of indexes before and after the index n corresponds to the peak.
  • On the horizontal axis represent discretized sample point indexes (hereinafter referred to as “frequency indexes”).
  • the peak shape can be expressed by the following function Q (n), where n is a variable representing a frequency index and ⁇ is a frequency index corresponding to a maximum value (peak).
  • Q the number of digits after the decimal point of the interval T is L digits
  • T ′ T ⁇ 2L .
  • h the height of the peak, and the peak height increases as the interval T increases.
  • PD represents the width of the peak portion, and the width increases as the interval T increases.
  • U is a positive integer indicating 1 to the number of peaks (for example, 1 to 10 in the case of FIG. 4), v is an integer of 1 or more (for example, about 1 to 3), and floor ( ⁇ ) is a decimal point
  • the periodic envelope sequence P [n] is, for example, It is sufficient to calculate as follows. However, (U ⁇ T ′) / 2 L ⁇ v ⁇ n ⁇ (U ⁇ T ′) / 2 L + v.
  • the periodic envelope series P [n] is calculated using the function Round (•) that rounds off the first decimal place and returns an integer value. You may ask as follows.
  • the periodic integrated envelope generation unit 150 receives at least the periodic envelope sequences P [1],..., P [N] and the amplitude spectrum envelope sequences W [1],.
  • the series W M [1],..., W M [N] is obtained (S150).
  • the periodic integrated envelope W M [n] is obtained as follows. Note that ⁇ is a value determined such that the shape of the absolute value series of the periodic integrated envelope W M [n] and the coefficient X [n] is close or a predetermined value.
  • the periodic integrated envelope generator 150 determines ⁇ so that the shape of the absolute value series of the periodic integrated envelope W M [n] and the coefficient X [n] is close
  • the periodic integrated envelope generator 150 , coefficient sequence X [1], ..., and X [N] also input, determined ⁇ and its periodicity integrated envelope sequence W M [1] when, ..., may be output W M [N].
  • ⁇ may be determined from several candidates of ⁇ , for example, two of 0.4 and 0.8, to ⁇ that minimizes E defined by the following expression. In other words, it may be determined to be ⁇ that makes the shape of the absolute value series of the periodic integrated envelope W M [n] and the coefficient X [n] close.
  • ⁇ is a value that determines how much the periodic envelope P [n] is considered in the periodic integrated envelope W M [n]. In other words, ⁇ can be said to be a value that determines the mixing ratio of the amplitude spectrum envelope W [n] and the periodic envelope P [n] in the periodic integrated envelope W M [n].
  • G in Equation (9) is the inner product of the absolute value sequence of each coefficient X [n] of the coefficient sequence X [1],..., X [N] and the reciprocal sequence of the periodic integrated envelope sequence.
  • ⁇ W M [n] in Expression (8) is a normalized periodic integrated envelope obtained by normalizing each value W M [n] of the periodic integrated envelope with G.
  • Equation (7) the fourth power of the inner product of the coefficient sequence X [1], ..., X [N] and the normalized periodic integrated envelope sequence ⁇ W M [1], ..., ⁇ W M [N] is calculated.
  • This is intended to reduce the value (distance) obtained by taking the inner product by emphasizing the coefficient X [n] having a particularly large absolute value.
  • ⁇ is determined so that the coefficient X [n] having a particularly large absolute value in the coefficient sequence X [1],..., X [N] is close to the periodic integrated envelope W M [n]. is doing.
  • the periodicity integrated envelope generation unit 150 determines the number of candidates for ⁇ according to the degree of periodicity
  • the periodicity integrated envelope generation unit 150 also receives an index S indicating the degree of periodicity as an input.
  • S indicates that the frame corresponds to a high periodicity
  • a ⁇ that minimizes E defined by Equation (7) is selected from a large number of candidates for ⁇
  • ⁇ may be a predetermined value. That is, when determining the number of candidates for ⁇ according to the degree of periodicity in the periodicity integrated envelope generation unit 150, the higher the periodicity, the larger the number of candidates for ⁇ .
  • FIG. 4 shows an example for explaining the difference between sequences generated for the same acoustic signal.
  • Fig. 4 (A) shows the shape of the curve obtained by interpolating the coefficient sequence X [1], ..., X [N], and
  • Fig. 4 (B) shows the periodic envelope series P [1], ..., P [N].
  • 4C shows the shape of the curve obtained by interpolating the smoothed amplitude spectrum envelope sequence ⁇ W [1],..., W [N]
  • FIG. 4D shows the periodic integrated envelope series. Shows the shape of the curve obtained by interpolating W M [1], ..., W M [N]. As shown in FIG.
  • the periodic integrated envelope sequence W M [1],..., W M [N] is a coefficient compared to the smoothed amplitude spectrum envelope sequence ⁇ W [1],.
  • the shape includes periodic peaks appearing in the rows X [1], ..., X [N].
  • the periodic integrated envelope sequence W M [1],..., W M [N] has an interval T or an interval T in addition to the linear prediction coefficient or the quantized linear prediction coefficient which is information representing the spectrum envelope. And the value ⁇ can be generated.
  • the peak of the amplitude due to the pitch period of the input acoustic signal can be expressed with higher accuracy than the spectral envelope determined by the linear prediction coefficient. it can. That is, the amplitude of the input acoustic signal can be estimated with high accuracy with a small amount of information of the linear prediction coefficient or the quantized linear prediction coefficient and the interval T, or the interval T and the value ⁇ .
  • the smoothed amplitude spectrum envelope to W [n] is an envelope expressed by the following equation, and ⁇ is a positive constant of 1 or less for blunting (smoothing) the amplitude spectrum coefficient.
  • the periodic integrated envelope sequence generation device of the present invention when used in the encoding device and the decoding device, the quantized linear obtained by a processing unit other than the periodic integrated envelope sequence generation device included in the encoding device. Since the code for specifying the prediction coefficient ⁇ ⁇ p (linear prediction coefficient code C L ) and the code for specifying the period T and the time domain period (period code C T ) are input to the decoding apparatus, the periodicity integration of the present invention If the code indicating the information of ⁇ is output from the envelope sequence generation device, the same period as the periodic integrated envelope sequence generated by the encoding-side periodic integrated envelope sequence generation device in the decoding-side periodic integrated envelope sequence generation device A gender integrated envelope sequence can be generated. Therefore, the amount of code that increases when a code is sent from the encoding device to the decoding device is small.
  • the periodic integrated envelope generating unit 150 performs the amplitude spectrum envelope sequence W [1 based on the periodic components of the coefficient sequence X [1],..., X [N]. ], ..., a modification of the W [N], periodicity integrated envelope sequence W M [1], ..., that is set to W M [N] is the most important point.
  • a “neighboring sample” is a sample indicated by an index that is an integer value in the vicinity of an integral multiple of the interval T.
  • the “neighborhood” may be a range determined by a predetermined method such as the equations (3) to (5).
  • the periodic envelope series P [1], represented by the equations (4) and (5) .., P [N] have a large value and a non-zero value in a wide range, that is, an integer multiple of the interval T (period) and many samples in the vicinity thereof. That is, the periodicity integrated envelope generation unit 150 increases the integer multiple of the interval T (period) and the values of samples in the vicinity thereof in the amplitude spectrum envelope sequence as the interval T between the components having periodicity in the coefficient sequence is wider. change.
  • the periodicity integrated envelope generation unit 150 increases the amplitude spectrum envelope sequence with a wider width, that is, an integer multiple of the interval T (period) and the vicinity thereof, as the interval T of the component having periodicity in the coefficient sequence is wider. For many samples, change the sample value. “With many samples in the vicinity” means increasing the number of samples existing in a range corresponding to “neighborhood” (range determined by a predetermined method). That is, the periodic integrated envelope generation unit 150 can easily obtain the above effect by modifying the amplitude spectrum envelope sequence in this way.
  • the first embodiment describes the periodic integrated envelope sequence generation device.
  • Modification 1 Example of periodicity analysis using a normalized coefficient sequence
  • the periodic integrated envelope sequence generation apparatus of Modification 1 is also shown in FIG.
  • the processing flow of the periodic integrated envelope series generation apparatus of the modification 1 is also shown in FIG.
  • the periodic integrated envelope sequence generation apparatus 101 includes a frequency domain sequence normalization unit 111, and the spectrum envelope sequence calculation unit 121 and the periodicity analysis unit 131 differ from the periodic integrated envelope sequence generation apparatus 100, and other configurations are as follows. The same. Only the differences will be described below.
  • the spectrum envelope sequence calculation unit 121 obtains not only the amplitude spectrum envelope sequence W [1],..., W [N] but also the smoothed amplitude spectrum envelope sequence ⁇ W [1],.
  • the spectrum envelope sequence calculation unit 121 performs the following procedure in addition to (step 1) and (step 2) indicated by the spectrum envelope sequence calculation unit 120.
  • Step 3 Multiply each quantized linear prediction coefficient ⁇ ⁇ p by ⁇ p to obtain quantized smoothed linear prediction coefficients ⁇ ⁇ 1 ⁇ , ⁇ ⁇ 2 ⁇ 2 ,..., ⁇ ⁇ P ⁇ P.
  • ⁇ is a positive constant of 1 or less for smoothing.
  • a smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N] is obtained by Expression (10) (S121).
  • the spectral envelope sequence calculation unit 120 may use a linear prediction coefficient alpha p instead of quantized linear prediction coefficient ⁇ alpha p.
  • Periodicity analysis unit 131 Periodicity analysis unit 131, the normalized coefficient sequence X N [1], ..., and enter the X N [N], the normalized coefficient sequence X N [1], ..., a period T of X N [N]
  • the period T is obtained and output (S131). That is, in this modification, the interval between the components having the periodicity of normalized coefficient sequences X N [1],..., X N [N], which is a frequency sequence coefficient sequence derived from the input acoustic signal, is obtained as the period T.
  • the periodicity analysis unit 131 may obtain and output an index S indicating the degree of periodicity by inputting the coefficient sequence X [1],..., X [N] as necessary.
  • periodic integrated envelope generation unit 150 uses smoothed amplitude spectrum envelope sequences ⁇ W [instead of amplitude spectrum envelope sequences W [1], ..., W [N]. 1], ..., ⁇ W [N] may be used. In this case, the following equation is used instead of equation (6).
  • a coefficient other than the periodic integrated envelope sequence generation device included in the encoding device or decoding device Sequence X [1], ..., X [N], normalized coefficient sequence XN [1], ..., XN [N], quantized linear prediction coefficient ⁇ ⁇ p , quantized smoothed linear prediction coefficient ⁇ ⁇ p ⁇ p , amplitude spectrum envelope W [1],..., W [N], smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N], period T, index S, etc. are obtained. There may be.
  • the periodic integrated envelope sequence generation device may be configured not to include at least one of the frequency domain conversion unit, the frequency domain normalization unit, the spectrum envelope sequence calculation unit, and the periodicity analysis unit.
  • a code linear prediction coefficient code C L
  • a code for specifying the period (period code C T ), a code for specifying the index S, and the like are output and input to the decoding device.
  • the periodic integrated envelope sequence generation device in the encoding device receives a code (linear prediction coefficient code C L ) for specifying the quantized linear prediction coefficient ⁇ ⁇ p , the period T and the time domain. There is no need to output a code for specifying the period (periodic code C T ), a code for specifying the index S, and the like.
  • the periodic integrated envelope sequence generation apparatus of the present invention when used in an encoding device or a decoding device, it is necessary to obtain the same periodic integrated envelope sequence in the encoding device and the decoding device. Therefore, it is necessary to obtain a periodic integrated envelope sequence using information that can be specified from codes output from the encoding device and input to the decoding device.
  • the spectral envelope sequence calculator periodicity integrated envelope sequence generator used in the coding device calculates the amplitude spectral envelope sequence using the quantized linear prediction coefficients corresponding to the linear prediction coefficient code C L, in the decoding device
  • the spectrum envelope sequence calculation unit of the periodic integrated envelope sequence generation device to be used uses the decoded linear prediction coefficient corresponding to the linear prediction coefficient code C L output from the encoding device and input to the decoding device to generate an amplitude spectrum envelope sequence. Need to ask.
  • the periodic integrated envelope sequence generating device when using a periodic integrated envelope sequence in an encoding device or a decoding device, is not provided internally as described above, but a necessary requirement in the periodic integrated envelope sequence generating device is required.
  • the processing unit may be provided in the encoding device and the decoding device. Such an encoding device and decoding device will be described in a second embodiment.
  • FIG. 5 shows a functional configuration example of the encoding apparatus according to the second embodiment
  • FIG. 6 shows a processing flow of the encoding apparatus according to the second embodiment.
  • the encoding apparatus 200 includes a spectrum envelope sequence calculation unit 221, a frequency domain conversion unit 110, a frequency domain sequence normalization unit 111, a periodicity analysis unit 230, a periodic envelope sequence generation unit 140, a periodicity integrated envelope generation unit 250, a variable A long coding parameter calculation unit 260 and a variable length coding unit 270 are provided.
  • Encoder 200 an audio digital signal of the input time domain as an input audio signal x (t), at least quantized linear prediction coefficient ⁇ alpha 1, ..., code C L, normalization coefficient indicating a ⁇ alpha P column X N [1], variable length ..., code C T interval T representing a cycle of X N [N], the normalized coefficient sequence X N [1], ..., and variable-length coding the X N [N]
  • the code CX is output.
  • the frequency domain sequence normalization unit 111 is the same as that in the first modification of the first embodiment.
  • the frequency domain transform unit 110 and the periodic envelope sequence generation unit 140 are the same as those in the first embodiment. Hereinafter, different components will be described.
  • the spectrum envelope sequence calculation unit 221 performs the amplitude spectrum envelope sequence W [1],..., W [N] of the input acoustic signal and the smoothed amplitude spectrum envelope sequence based on the time domain linear prediction of the input acoustic signal x (t). ⁇ W [1], ..., and calculates a ⁇ W [N], obtained in the course of calculation quantized linear prediction coefficient ⁇ ⁇ 1, ..., also obtains the code C L indicating the ⁇ ⁇ P (S221). However, N is a positive integer.
  • the spectrum envelope sequence calculation unit 221 may be processed in the following procedure.
  • Step 1 Linear prediction coefficients ⁇ 1 ,..., ⁇ P are obtained by performing linear prediction analysis on the input acoustic signal in units of frames that are predetermined time intervals.
  • P is a positive integer indicating the predicted order.
  • the input acoustic signal x (t) at the time t becomes the past value x (t ⁇ 1),. tP) and the prediction residuals e (t) and the linear prediction coefficients alpha 1, ..., represented by the formula (1) by alpha p.
  • Step2 linear prediction coefficients ⁇ 1, ..., and outputs to obtain a code C L encodes the alpha P, quantized linear prediction coefficients corresponding to the code C L ⁇ ⁇ 1, ..., determine the ⁇ alpha P . Furthermore, quantized linear prediction coefficient ⁇ ⁇ 1, ..., ⁇ ⁇ using P of the input audio signals of N points amplitude spectral envelope sequence W [1], ..., determine the W [N]. For example, each value W [n] of the amplitude spectrum envelope series can be obtained by Expression (2).
  • linear prediction coefficients alpha 1, ..., a method of obtaining a alpha P to be encoded code C L converts the linear prediction coefficients to LSP parameters, such as obtaining a code C L encodes the LSP parameter, linear predictive You may use any method for obtaining the code C L any convertible coefficients in the coefficient is encoded.
  • Step 3 Multiply each quantized linear prediction coefficient ⁇ ⁇ p by ⁇ p to obtain quantized smoothed linear prediction coefficients ⁇ ⁇ 1 ⁇ , ⁇ ⁇ 2 ⁇ 2 ,..., ⁇ ⁇ P ⁇ P.
  • ⁇ is a positive constant of 1 or less for smoothing in advance.
  • the smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N] is obtained by Expression (10).
  • Periodicity analysis unit 230 Periodicity analysis unit 230, the normalized coefficient sequence X N [1], ..., X N as input [N], the normalized coefficient sequence X N [1], ..., X N [N] interval T ( An interval having a periodically large value is obtained, and an interval T and a code CT indicating the interval T are output (S230). Further, the periodicity analysis unit 230 also obtains and outputs an index S indicating the degree of periodicity (that is, an index indicating the degree of periodicity of the sample sequence in the frequency domain) as necessary. The period analysis section 230, if necessary, also obtained an output code C S indicating the index S. The index S and the interval T itself are the same as those of the periodicity analysis unit 131 of the first modification of the first embodiment.
  • the periodic integrated envelope generation unit 250 receives at least the periodic envelope sequence P [1],..., P [N] and the amplitude spectrum envelope sequence W [1],.
  • the sequence W M [1],..., W M [N] is obtained, and the periodic integrated envelope W M [n] is output.
  • the integrated periodic envelope generation unit 150 selects any one of a plurality of predetermined candidate values as the value ⁇ instead of a predetermined value, the coefficient sequence X [1], ..., X [N] is also input, and a candidate value that is close to the shape of the absolute value series of periodic integrated envelope W M [n] and coefficient X [n] among a plurality of predetermined candidate values is obtained as value ⁇ .
  • the code C ⁇ indicating the value ⁇ is also output (S250).
  • the periodic integrated envelope W M [n] and the value ⁇ are the same as those in the first embodiment, and the periodic integrated envelope W M [n] may be obtained as in the equations (6),.
  • the periodicity integrated envelope generation unit 150 determines the number of candidates for ⁇ according to the degree of periodicity, the periodicity integrated envelope generation unit 150 also receives an index S indicating the degree of periodicity, and the index S is In the case of a frame corresponding to high periodicity, ⁇ that minimizes E defined by Equation (7) is selected from among a large number of candidate ⁇ , and the index S is low in periodicity.
  • ⁇ may be set to one predetermined value. When ⁇ is set to a predetermined value, it is not necessary to output the code C ⁇ indicating the value ⁇ .
  • variable length coding parameter calculation unit 260 performs regular integration with the periodic integrated envelope sequence W M [1],..., W M [N] and the smoothed amplitude spectrum envelope sequence ⁇ W [1],. reduction coefficient sequence X n [1], ..., and enter the X n [n], obtaining the variable length coding parameters r n (S260).
  • Variable-length encoding parameter calculating unit 260 periodicity integrated envelope sequence W M [1], ..., characterized in that in dependence on the amplitude value obtained from W M [N] to calculate the variable length coding parameters r n It is said.
  • the variable length coding parameter is a parameter that specifies a possible range of the amplitude of each coefficient of the signal to be coded, that is, the normalized coefficient sequence X N [1],..., X N [N].
  • the Rice parameter corresponds to a variable length coding parameter
  • the range that the amplitude of a signal to be coded can take corresponds to the variable length coding parameter.
  • variable length coding parameter calculating unit 260 for each normalization moiety coefficient sequence which is part of the normalization coefficient sequence, to calculate the variable length coding parameters r n.
  • the normalized partial coefficient sequences include the coefficients of the normalized coefficient sequence without overlapping.
  • Step1 normalized coefficient sequence X N [1], ..., X N the logarithm of the average of the amplitudes of the coefficients of [N], serving as a reference Rice parameter sb equation as (variable length coding parameter as a reference) Calculate as follows.
  • the sb is encoded only once for each frame, and transmitted to the decoding apparatus 400 as a code C sb corresponding to the reference Rice parameter (reference variable-length encoding parameter).
  • the average value of the amplitudes of the normalized coefficient sequences X N [1],..., X N [N] can be estimated from other information transmitted to the decoding apparatus 400, the encoding apparatus 200 and the decoding apparatus 400 are common.
  • a method of approximately determining sb from the estimated value of the average amplitude value may be determined.
  • the average value of the amplitude can be estimated from other information transmitted to the decoding device 400. In this case, it is not necessary to encode sb and output the code C sb corresponding to the reference rice parameter to the decoding device 400.
  • Step 2 The threshold value ⁇ is calculated by the following equation. ⁇ is the logarithm of the average amplitude of values obtained by dividing each value W M [n] of the periodic integrated envelope sequence by each value ⁇ W [n] of the smoothed amplitude spectrum envelope sequence.
  • Step3 W M [n ] / ⁇ W [n]
  • Variable length coding unit 270 the normalization coefficients using variable length coding parameters r n obtained by the variable length coding parameter calculating section 260 columns X N [1], ..., variable-length codes X N [N]
  • the variable length code C X is output (S270).
  • the variable length coding unit 270, the normalized coefficient sequence X N [1] using a Rice parameter r n obtained by the variable length coding parameter calculating section 260, ..., and Rice coding the X N [N] The obtained code is output as a variable length code CX .
  • Rice parameter r n obtained by the variable length coding parameter calculation unit 260, a variable length coding parameter dependent on the amplitude value of the periodicity integration envelope sequence, a larger value as the frequency value greater periodicity integrated envelope sequence It has become.
  • Rice coding is one of the known techniques of variable-length coding which depends on the amplitude value, and performs variable length coding which depends on the amplitude value using the Rice parameter r n.
  • the periodic integrated envelope sequence generated by the periodic integrated envelope generating unit 250 expresses the spectral envelope of the input acoustic signal with high accuracy.
  • variable length coding unit 270 determines that the frequency of the coefficient sequence in the frequency domain of the input acoustic signal X [1],..., X [N] increases as the frequency of the periodic integrated envelope sequence increases. Therefore, the normalized coefficient sequence X N [1], ..., X N [N] is variable-length encoded, in other words, depends on the amplitude value using the variable-length encoding parameter.
  • the normalized coefficient sequence X N [1],..., X N [N] is encoded by variable length coding.
  • the amplitude value here is an average amplitude value of the coefficient sequence to be encoded, an estimated value of the amplitude of each coefficient included in the coefficient sequence, an estimated value of an envelope of the amplitude of the coefficient sequence, or the like.
  • the encoding apparatus 200 includes a code C L indicating the quantized linear prediction coefficients ⁇ ⁇ 1 ,..., ⁇ ⁇ P obtained by such processing, a code C T indicating the interval T, and a normalized coefficient sequence X N [ 1],..., X N [N]
  • a variable length code C X obtained by variable length encoding is output.
  • code C sb showing a variable length coding parameter sb as a code C [delta] and the reference indicating the value [delta], if necessary.
  • the code output from the encoding device 200 is input to the decoding device 400.
  • the encoding device includes only the periodic envelope sequence generation unit 140, the periodic integrated envelope generation unit 250, the variable length encoding parameter calculation unit 260, and the variable length encoding unit 270, and is generated outside the encoding device. Smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N], normalized coefficient sequence XN [1], ..., XN [N], interval T, and amplitude as necessary The spectrum envelope sequence W [1],..., W [N] and the index S as necessary may be input and the variable length code C X may be output.
  • FIG. 7 shows a functional configuration example of the decoding apparatus according to the second embodiment
  • FIG. 8 shows a processing flow of the decoding apparatus according to the second embodiment.
  • Decoding apparatus 400 includes spectrum envelope sequence calculation unit 421, periodic envelope sequence generation unit 440, periodic integrated envelope generation unit 450, variable length coding parameter calculation unit 460, variable length decoding unit 470, frequency domain sequence denormalization unit 411 and a frequency domain inverse transform unit 410.
  • Decoding device 400 quantized linear prediction coefficient ⁇ ⁇ 1, ..., ⁇ ⁇ code indicating the P C L, code C T indicating the interval T, the normalization coefficient sequence X N [1], ..., X N [N ] to receive the variable length code C X which variable length coding, and outputs an acoustic signal.
  • code C S indicating the sign C sb and index S indicating a variable-length coding parameters sb as a code C [delta] and the reference indicating the value [delta], if necessary. Details of each component will be described below.
  • Spectrum envelope series calculation unit 421 inputs the code C L, the amplitude spectral envelope sequence W [1], ..., W [N] and the smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N] Is calculated (S421). More specifically, the processing may be performed according to the following procedure.
  • Step1 decodes the code C L, decodes the linear prediction coefficient ⁇ ⁇ 1, ..., obtaining ⁇ alpha P.
  • Step2 decoded linear prediction coefficient ⁇ ⁇ 1, ..., ⁇ amplitude spectral envelope sequence W of N points using ⁇ P [1], ..., determine the W [N].
  • each value W [n] of the amplitude spectrum envelope series can be obtained by Expression (2).
  • Step 3 Each of the decoded linear prediction coefficients ⁇ ⁇ p is multiplied by ⁇ p to obtain decoded smoothed linear prediction coefficients ⁇ ⁇ 1 ⁇ , ⁇ ⁇ 2 ⁇ 2 ,..., ⁇ ⁇ P ⁇ P.
  • ⁇ is a positive constant of 1 or less for smoothing in advance.
  • the smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N] is obtained by Expression (10).
  • Periodicity envelope sequence generating unit 440 inputs the code C T indicating the interval T, and decodes the code C T, obtaining a spacing T. Then, periodic envelope sequences P [1],..., P [N] are obtained and output by the same method as the periodic envelope sequence generation unit 140 of the encoding device 200 (S440).
  • Periodicity integrated envelope generator 450 decodes the code C [delta], to obtain a value [delta]. However, if the code C [delta] is not input, the decoding of the code C [delta] is not performed, to obtain a pre-stored value [delta] Periodicity integrated envelope generator 450.
  • the periodic integrated envelope generation unit 450 acquires the index S by decoding the code C S , and the acquired index S is a frame corresponding to the high periodicity. In the case of, the code C ⁇ is decoded to obtain the value ⁇ . If the obtained index S is a frame corresponding to low periodicity, the code C ⁇ is not decoded, and the periodic integrated envelope is obtained. The value ⁇ stored in advance in the generation unit 450 is acquired. Then, the periodic integrated envelope generation unit 450 obtains the periodic integrated envelope sequence W M [1],..., W M [N] by Expression (6). (S450)
  • variable length coding parameter calculation unit 460 encodes the periodic integrated envelope sequence W M [1],..., W M [N] and the smoothed amplitude spectrum envelope sequence ⁇ W [1],. as input C sb, to obtain a variable length coding parameters r n (S460).
  • a method of approximately determining sb from the estimated average amplitude value estimated from the other information is determined. Also good. In this case, the code C sb is not input.
  • a variable length coding parameter calculation method will be described by taking as an example the case of performing rice decoding for each sample.
  • Step 1 The code C sb is decoded to obtain a reference rice parameter sb (a reference variable-length encoding parameter).
  • a reference rice parameter sb a reference variable-length encoding parameter
  • Step 2 The threshold value ⁇ is calculated by the equation (14).
  • Step3 W M [n ] / ⁇ W [n]
  • variable length decoding unit 470 obtains variable length coding parameters r n by decoding the variable length codes C X using the decoded normalization coefficient sequence ⁇ X N [1], ... , ⁇ X N [N] is obtained (S470).
  • the variable length decoding unit 470, a variable length coding parameter calculating section 460 obtains the Rice parameter r n by decoding the variable length codes C X using the decoded normalization coefficient sequence ⁇ X N [1], ... , Get ⁇ X N [N].
  • the decoding method of the variable length decoding unit 470 corresponds to the encoding method of the variable length encoding unit 270.
  • the frequency domain sequence inverse normalization unit 411 performs the decoding normalization coefficient sequence ⁇ X N [1],..., ⁇ X N [N] and the smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N].
  • ⁇ X [n] ⁇ X N [n] ⁇ ⁇ W [n]
  • the decoding coefficient sequence ⁇ X [1],..., ⁇ X [N] is obtained and output (S411).
  • the frequency domain inverse transform unit 410 receives the decoded coefficient sequence ⁇ X [1],..., ⁇ X [N] as input, and outputs the decoded coefficient sequence ⁇ X [1], ..., ⁇ X [N] in a predetermined time interval.
  • the sound signal is converted into a certain frame unit (time domain) (S410).
  • the decoding apparatus includes only a periodic envelope sequence generation unit 440, a periodic integrated envelope generation unit 450, a variable length coding parameter calculation unit 460, and a variable length decoding unit 470, and is input to the decoding apparatus as necessary.
  • the smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., W [N] the smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., W [N]
  • the amplitude spectrum envelope series W [1] obtained outside the decoding apparatus.
  • variable-length coding is a coding method that improves coding efficiency by adaptively determining a code in accordance with a possible range of the amplitude of an input value to be coded.
  • normalized coefficient sequences X N [1],..., X N [N] that are frequency sequence coefficient sequences are to be encoded, but the amplitude of each coefficient included in the coefficient sequence to be encoded is If variable-length coding is performed using variable-length coding parameters obtained using information more accurately, the coding efficiency of the variable-length coding itself performed by the coding apparatus is increased.
  • the decoding device in order for the decoding device to obtain the variable-length encoding parameter, it is necessary to send more accurately the amplitude information of each coefficient included in the coefficient sequence to be encoded from the encoding device to the decoding device. Therefore, the amount of code sent from the encoding device to the decoding device increases.
  • a method for obtaining an estimated value of the amplitude of each coefficient included in the coefficient sequence to be encoded from a code with a small code amount is required. Since the periodic integrated envelope sequence W M [1],..., W M [N] of the second embodiment approximates the coefficient sequence X [1],..., X [N] with high accuracy,
  • the envelope including the peak of the amplitude due to the pitch period of the input acoustic signal input to the encoding device is expressed by the code C L , the code C T , and the code C. It can be reproduced by a decoding device with a small amount of information of only ⁇ .
  • the encoding device and the decoding device of the second embodiment are often used in combination with an encoding device and a decoding device that perform encoding and decoding with linear prediction and pitch prediction.
  • the code C L and the code C T are decoded from a coding apparatus that performs coding with linear prediction and pitch prediction outside the coding apparatus 200, and with linear prediction and pitch prediction outside the decoding apparatus 400. It is the code
  • C ⁇ .
  • the code amount of the code C ⁇ is small (each is about 3 bits at most, and an effect can be obtained even with 1 bit), and corresponds to the variable length coding parameter for each partial sequence included in the normalized coefficient sequence to be coded. Less than the total code amount of codes.
  • the encoding efficiency can be improved with a small increase in the code amount.
  • the encoding apparatus 200 includes: A frequency domain based on a spectrum envelope sequence that is a frequency domain sequence corresponding to a linear prediction coefficient code obtained from an input acoustic signal in a predetermined time interval and a frequency domain period corresponding to a periodic code obtained from the input acoustic signal Periodic integrated envelope generator 250 for generating a periodic integrated envelope sequence that is a sequence of A variable length encoding unit 270 that encodes a frequency domain sequence derived from an input acoustic signal on the assumption that the amplitude of the input acoustic signal is larger as the frequency of the periodic integrated envelope sequence is larger.
  • the decoding device 400 has A periodic integrated envelope that generates a periodic integrated envelope sequence that is a frequency domain sequence based on a spectrum envelope sequence that is a frequency domain sequence corresponding to a linear prediction coefficient code and a frequency domain period that corresponds to a periodic code.
  • Generation unit 450 A variable length decoding unit 470 that decodes a variable length code to obtain a frequency domain sequence on the premise that the amplitude of the acoustic signal is larger as the frequency of the periodic integrated envelope sequence is larger, What is necessary is just to make it have.
  • “derived from the input sound signal” means being obtained from the input sound signal or corresponding to the input sound signal. For example, the coefficient sequence X [1],..., X [N] and the normalized coefficient sequence X N [1],..., X N [N] are frequency domain sequences derived from the input acoustic signal.
  • FIG. 9 shows a functional configuration example of the encoding apparatus of the third embodiment
  • FIG. 10 shows a processing flow of the encoding apparatus of the third embodiment.
  • the encoding apparatus 300 includes a spectral envelope sequence calculation unit 221, a frequency domain conversion unit 110, a frequency domain sequence normalization unit 111, a periodicity analysis unit 330, a periodicity envelope sequence generation unit 140, a periodicity integrated envelope generation unit 250, a variable A long coding parameter calculation unit 260, a second variable length coding parameter calculation unit 380, and a variable length coding unit 370 are provided.
  • Encoding apparatus 300 a sound digital signal of the input time domain as an input audio signal x (t), at least quantized linear prediction coefficient ⁇ alpha 1, ..., code C L, normalization coefficient indicating a ⁇ alpha P , X N [1],..., X N [N] representing the period of the code C T , coefficient sequence X [1],..., X [N] or normalized coefficient sequence X N [1],. , X N code C S indicating the predetermined index S and the index S indicating the degree of periodicity of the [N], the normalized coefficient sequence X N [1], variable that ... and variable-length coding the X N [N]
  • the long code CX is output.
  • the frequency domain sequence normalization unit 111 is the same as that in the first modification of the first embodiment.
  • the frequency domain transform unit 110 and the periodic envelope sequence generation unit 140 are the same as those in the first embodiment.
  • the amplitude spectrum envelope sequence calculation unit 221, the periodic integrated envelope generation unit 250, and the variable length coding parameter calculation unit 260 are the same as those in the second embodiment.
  • different components will be described.
  • the periodicity analysis unit 330 receives the normalization coefficient sequence X N [1],..., X N [N] as inputs, An index S indicating the degree of periodicity of the normalization coefficient sequence X N [1],..., X N [N] and an interval T (interval that periodically becomes a large value) are obtained. and it outputs the code C T showing the codes C S and interval T and spacing T shown (S330).
  • the index S and the interval T itself are the same as those of the periodicity analysis unit 131 of the first modification of the first embodiment.
  • variable-length coding parameter calculation unit 260 calculates the variable length coding parameters r n, the index S If it is not in the range that indicates that a large degree of predetermined periodicity, the second variable length coding parameter calculation unit 380 calculates the variable length coding parameters r n (S390).
  • the “range indicating that the predetermined degree of periodicity is large” may be set, for example, when the index S is equal to or greater than a predetermined threshold.
  • the second variable length coding parameter calculation unit 380 normalizes the amplitude spectrum envelope sequence W [1], ..., W [N] and the smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N]. coefficient sequence X n [1], ..., and enter the X n [n], obtaining the variable length coding parameters r n (S380).
  • Variable-length encoding parameter calculating unit 260 periodicity integrated envelope sequence W M [1], ..., characterized in that in dependence on the amplitude value obtained from W M [N] to calculate the variable length coding parameters r n
  • the second variable length coding parameter calculation unit 380 is characterized in that the variable length coding parameter is calculated depending on the amplitude value obtained from the amplitude spectrum envelope sequence.
  • a variable length coding parameter calculation method will be described by taking as an example the case of performing rice coding for each sample.
  • Step 1 The logarithm of the average amplitude of each coefficient of the normalized coefficient sequence X N [1],..., X N [N] is expressed as a reference Rice parameter sb (reference variable-length encoding parameter). 13). This process is the same as that of the variable length coding parameter calculation unit 260.
  • Step 2 The threshold value ⁇ is calculated by the following equation. ⁇ is the logarithm of the average amplitude of values obtained by dividing each value W [n] of the amplitude spectrum envelope series by each value of the smoothed amplitude spectrum envelope series to W [n].
  • Step3 W [n] / ⁇ W [n]
  • Variable length coding unit 370 by using the variable length coding parameters r n normalized coefficient sequence X N [1], ..., the X N [N] and variable length coding, and outputs a variable-length code C X ( S370).
  • variable length coding parameters r n if the range indicates that the degree of periodicity indicator S is predetermined large, a variable length coding parameters r n of the variable length coding parameter calculation unit 260 has calculated There, if the index S is not in the range that indicates that the degree of a predetermined periodicity greater, a variable length coding parameters r n the second variable length coding parameter calculation unit 380 has calculated.
  • the encoding apparatus 300 includes a code C L indicating quantized linear prediction coefficients ⁇ ⁇ 1 ,..., ⁇ ⁇ P obtained by such processing, a code C S indicating an index S indicating the degree of periodicity, and an interval.
  • a variable length code C X obtained by variable length coding the code C T indicating the T and the normalized coefficient sequence X N [1],..., X N [N] is output and transmitted to the decoding side.
  • output code C sb showing a variable length coding parameters sb as a code C [delta] and the reference indicating the value [delta], if necessary, sent to the decoding side.
  • the encoding apparatus includes a periodic envelope sequence generation unit 140, a periodic integrated envelope generation unit 250, a variable length encoding parameter calculation unit 260, a second variable length encoding parameter calculation unit 380, and a variable length encoding unit 370. And a smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N] and a normalized coefficient sequence X N [1], ..., X N [N] generated outside the encoder ., W [N] as necessary, and an index S as necessary, and a variable length code C X may be output.
  • FIG. 11 shows a functional configuration example of the decoding apparatus according to the third embodiment
  • FIG. 12 shows a processing flow of the decoding apparatus according to the third embodiment.
  • the decoding apparatus 500 includes a spectrum envelope sequence calculation unit 421, an index decoding unit 530, a periodic envelope sequence generation unit 440, a periodic integrated envelope generation unit 450, a variable length coding parameter calculation unit 460, and a second variable length coding parameter calculation. 580, variable length decoding section 570, frequency domain sequence inverse normalization section 411, and frequency domain inverse transform section 410.
  • Decoding device 500 quantized linear prediction coefficient ⁇ alpha 1, ..., ⁇ code C L indicating the alpha P, code C S indicating the index S, code C T, normalized coefficient sequence X N [1 showing the interval T ],..., X N [N] variable length code C X is received, and an acoustic signal is output. Even receive code C sb showing a variable length coding parameters sb as a code C [delta] and the reference indicating the value [delta], if necessary.
  • the spectrum envelope sequence calculation unit 421, the periodic envelope sequence generation unit 440, the periodic integrated envelope generation unit 450, the variable length coding parameter calculation unit 460, the frequency domain sequence denormalization unit 411, and the frequency domain inverse conversion unit 410 are examples. Same as 2. Hereinafter, different components will be described.
  • Index decoding unit 530 decodes the code C S, to obtain an index S.
  • the decoding device 500 in the case of a range indicating that the degree of periodicity indicator S is predetermined large, variable-length coding parameter calculation unit 460 calculates the variable length coding parameters r n, defined index S in advance If the degree of periodicity is not the range indicated is greater, the second variable length coding parameter calculation unit 580 calculates the variable length coding parameters r n (S590).
  • the “range indicating that the predetermined degree of periodicity is large” is the same range as that of the encoding apparatus 300.
  • the second variable length coding parameter calculation unit 580 includes the amplitude spectrum envelope sequence W [1],..., W [N] and the smoothed amplitude spectrum envelope sequence ⁇ W [1],. as input sb, obtaining a variable length coding parameters r n (S580).
  • a method of approximately determining sb from the estimated average amplitude value estimated from the other information is determined. Also good. In this case, the code C sb is not input.
  • a variable length coding parameter calculation method will be described by taking as an example the case of performing rice decoding for each sample.
  • Step 1 The code C sb is decoded to obtain a reference rice parameter sb (a reference variable-length encoding parameter).
  • a reference rice parameter sb a reference variable-length encoding parameter
  • Step 2 The threshold value ⁇ is calculated by the equation (16).
  • Step3 W [n] / ⁇ W [n]
  • Variable length decoding unit 570 by using the variable length coding parameters r n variable length code C X decrypted by the decryption normalized coefficient sequence ⁇ X N [1], ... , ⁇ X N Request [N] (S570 ).
  • variable length coding parameters r n if the range indicates that the degree of periodicity indicator S is predetermined large, a variable length coding parameters r n of the variable length coding parameter calculation unit 460 has calculated There, if the index S is not in the range that indicates that the degree of a predetermined periodicity greater, a variable length coding parameters r n the second variable length coding parameter calculation unit 580 has calculated.
  • Modification 1 of Decoding Device (Example in which information is input from the outside)
  • the periodic envelope sequence generation unit 440 the periodic integrated envelope generation unit 450, the variable length coding parameter calculation unit 460, the second variable length coding parameter calculation unit 580, and the variable length decoding unit 570 are included.
  • the smoothed amplitude spectrum envelope sequence ⁇ W [1], ..., ⁇ W [N] obtained outside the decoding apparatus Amplitude spectrum envelope series W [1], ..., W [N], interval T, index S are also input, and normalization coefficient sequence X N [1], ..., X N [N] is output and smoothed externally.
  • the time domain acoustic signal may be converted by multiplying the normalized amplitude spectrum envelope sequence.
  • the encoding device and the decoding device obtain a variable length encoding parameter using the periodic integrated envelope sequence when the degree of periodicity of the acoustic signal to be encoded is large, and perform encoding. If the degree of periodicity of the target acoustic signal is not large, the variable length coding parameter is obtained using the amplitude spectrum envelope sequence, so that the variable length coding parameter can be used for variable length coding. There is an effect that the encoding accuracy can be increased.
  • a power sequence, , W [n], ⁇ W [n], W M [n] may be a power spectrum envelope sequence, a smoothed power spectrum envelope sequence, or a periodic integrated envelope sequence that is a power sequence.
  • the program describing the processing contents can be recorded on a computer-readable recording medium.
  • a computer-readable recording medium any recording medium such as a magnetic recording device, an optical disk, a magneto-optical recording medium, and a semiconductor memory may be used.
  • this 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. Furthermore, 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 a program recorded on a portable recording medium or a program transferred from a server computer in its own storage device.
  • the computer reads a program stored in its own recording medium and executes a process according to the read program.
  • the computer may directly read the program from a portable recording medium and execute processing according to the program, and the program is transferred from the server computer to the computer.
  • the processing according to the received program may be executed sequentially.
  • 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 a processing function only by an execution instruction and result acquisition. It is good.
  • the program in this embodiment includes information that is used for processing by an electronic computer and that conforms to the program (data that is not a direct command to the computer but has a property that defines the processing of the computer).
  • the present apparatus is configured by executing a predetermined program on a computer.
  • a predetermined program on a computer.
  • at least a part of these processing contents may be realized by hardware.
  • Periodic integrated envelope sequence generation device 110 Frequency domain conversion unit 111 Frequency domain coefficient normalization unit 120, 121, 221, 421 Spectral envelope sequence calculation unit 130, 131, 230, 330 Periodicity analysis unit 140, 440 Periodicity Envelope sequence generation unit 150, 250, 450 Periodic integrated envelope generation unit 200, 300 Encoding device 260, 360, 460 Variable length encoding parameter calculation unit 270, 370 Variable length encoding unit 380, 580 Second variable length encoding Parameter calculation unit 400, 500 decoding device 410 frequency domain inverse transformation unit 411 frequency domain sequence inverse normalization unit 470, 570 variable length decoding unit 530 index decoding unit

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

　音響信号のピッチ周期に起因するピークの付近での近似精度がよい周波数領域の包絡系列を少ない符号量で再現できる符号化装置、復号装置を提供する。本発明の符号化装置は、周期性統合包絡生成部と可変長符号化部を備える。周期性統合包絡生成部は、入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する。可変長符号化部は、入力音響信号に由来する周波数領域の系列を符号化する。本発明の復号装置は、周期性統合包絡生成部と可変長復号部を備える。周期性統合包絡生成部は、線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する。可変長復号部は、可変長符号を復号して周波数領域の系列を得る。

Description

符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、記録媒体
　本発明は、音響信号のスペクトル包絡を利用して音響信号を符号化、復号化する符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、および記録媒体に関する。
　低ビット（例えば10kbit/s～20kbit/s程度）の音声信号や音響信号の符号化方法として、ＤＦＴ（離散フーリエ変換）やＭＤＣＴ（変形離散コサイン変換）などの直交変換係数に対する適応符号化が知られている。例えば非特許文献１で用いられているＴＣＸ（transform coded excitation：変換符号化励振）符号化方法では、入力された音信号の周波数領域表現である係数列X[1]，…，X[N]から振幅スペクトル包絡の影響を取り除いた系列（正規化係数列XN[1]，…，XN[N]）を求め、これを可変長符号化する。ただし、[]内のＮは正整数である。
　振幅スペクトル包絡は、以下の手順で算出される。
（ｓｔｅｐ１）所定の時間区間であるフレーム単位で、入力された時間領域の音響ディジタル信号（以下、入力音響信号）に対する線形予測分析を行って線形予測係数α1，…，αPを求める。ただし、Ｐは予測次数を示す正整数である。例えば、全極型モデルであるＰ次自己回帰過程により、時刻ｔでの入力音響信号x(t)は、Ｐ時点まで遡った過去の自分自身の値x(t-1)，…，x(t-P)と予測残差e(t)と線形予測係数α1，…，αpによって式（１）で表される。
　x(t)=α1x(t-1)+…+αp x(t-P)+e(t)   (1)
（ｓｔｅｐ２）線形予測係数α1，…，αPを量子化し、量子化済線形予測係数^α1，…，^αPを求める。量子化済線形予測係数^α1，…，^αPを用いてＮ点の入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]を求める。例えば、振幅スペクトル包絡系列の各値W[n]は、式（２）で求めることができる。ただし、ｎは1≦n≦Nの整数、exp（・）はネイピア数を底とする指数関数、ｊは虚数単位、σは予測残差信号の振幅である。
Figure JPOXMLDOC01-appb-M000001
　なお、本明細書では、右肩に角括弧なしで表記されている記号はべき乗演算を表す。つまり、σ2はσの２乗を表す。また、文中で使用する記号「~」「＾」等は、本来直後の文字の真上に記載されるべきものであるが、テキスト記法の制限により、当該文字の直前に記載する。数式中においてはこれらの記号は本来の位置、すなわち文字の真上に記述している。
Anthony Vetro, "MPEG Unified Speech and Audio Coding", Industry and Standards, IEEE MultiMedia, April-June, 2013.
　音響信号の符号化では、復号側でもスペクトル包絡の情報を得るために、スペクトル包絡に対応する符号を復号側へ送る必要がある。非特許文献１のように線形予測係数によりスペクトル包絡を求める場合には、復号側へ送る「スペクトル包絡に対応する符号」は「線形予測係数に対応する符号」であり、符号量が少なくて済むという利点がある。一方、線形予測係数により求めたスペクトル包絡の情報は、入力音響信号のピッチ周期に起因するピークの付近での近似精度が悪くなることがあり、このことが正規化係数列を可変長符号化する際の符号化効率の低下につながることがある。
　このような問題に鑑み、本発明では、音響信号のピッチ周期に起因するピークの付近での近似精度がよい周波数領域の包絡系列を少ない符号量で再現できる符号化装置、復号装置を提供する。
　本発明の符号化装置は、周期性統合包絡生成部と可変長符号化部を備える。周期性統合包絡生成部は、所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する。可変長符号化部は、周期性統合包絡系列の値が大きい周波数ほど、入力音響信号の振幅が大きいことを前提に、入力音響信号に由来する周波数領域の系列を符号化する。本発明の復号装置は、周期性統合包絡生成部と可変長復号部を備える。周期性統合包絡生成部は、線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する。可変長復号部は、周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、可変長符号を復号して周波数領域の系列を得る。
　本発明の符号化装置と復号装置によれば、入力音響信号のピッチ周期に起因するピーク付近での近似精度がよい周波数領域の包絡系列を少ない符号量で再現できる符号化、復号を行うので、音響信号の符号化効率がよい。
実施例１の周期性統合包絡系列生成装置の機能構成例を示す図。 実施例１の周期性統合包絡系列生成装置の処理フローを示す図。 周期性包絡系列P[1]，…，P[N]の例を示す図。 同じ音響信号に対して生成された系列の違いを説明するための例を示す図。 実施例２の符号化装置の機能構成例を示す図。 実施例２の符号化装置の処理フローを示す図。 実施例２の復号装置の機能構成例を示す図。 実施例２の復号装置の処理フローを示す図。 実施例３の符号化装置の機能構成例を示す図。 実施例３の符号化装置の処理フローを示す図。 実施例３の復号装置の機能構成例を示す図。 実施例３の復号装置の処理フローを示す図。
　以下、本発明の実施の形態について、詳細に説明する。なお、同じ機能を有する構成部には同じ番号を付し、重複説明を省略する。
　図１に本発明の周期性統合包絡系列生成装置の機能構成例を、図２に本発明の周期性統合包絡系列生成装置の処理フローを示す。周期性統合包絡系列生成装置１００は、スペクトル包絡系列計算部１２０、周波数領域変換部１１０、周期性分析部１３０、周期性包絡系列生成部１４０、周期性統合包絡生成部１５０を備え、入力された時間領域の音響ディジタル信号を入力音響信号x(t)とし、係数列の周波数成分に基づいて振幅スペクトル包絡系列を変形した周期性統合包絡系列を生成する。
＜スペクトル包絡系列計算部１２０＞
　スペクトル包絡系列計算部１２０は、入力音響信号x(t)の時間領域の線形予測に基づき、入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]を計算する（Ｓ１２０）。ただし、Ｎは正整数である。スペクトル包絡系列計算部１２０は、従来技術と同じであり、以下の手順で計算すればよい。
　（ｓｔｅｐ１）所定の時間区間であるフレーム単位で、入力音響信号に対する線形予測分析を行って線形予測係数α1，…，αPを求める。ただし、Ｐは予測次数を示す正整数である。例えば、全極型モデルであるＰ次自己回帰過程により、時刻ｔでの入力音響信号x(t)は、Ｐ時点まで遡った過去の自分自身の値x(t-1)，…，x(t-P)と予測残差e(t)と線形予測係数α1，…，αpによって式（１）で表される。
　（ｓｔｅｐ２）線形予測係数α1，…，αPを用いてＮ点の入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]を求める。例えば、振幅スペクトル包絡系列の各値W[n]は、線形予測係数α1，…，αPに対応する量子化済線形予測係数^α1，…，^αPを用いて式（２）で求めることができる。または、振幅スペクトル包絡系列の各値W[n]は、線形予測係数α1，…，αPを用いて、式（２）の^αpをαpに置き換えた式で求めることができる。
＜周波数領域変換部１１０＞
　周波数領域変換部１１０は、所定の時間区間であるフレーム単位で、入力された時間領域の入力音響信号を周波数領域のＮ点の係数列X[1]，…，X[N]に変換して出力する（Ｓ１１０）。周波数領域への変換は、ＭＤＣＴ（変形離散コサイン変換）やＤＦＴ（離散フーリエ変換）などの方法で行えばよい。
＜周期性分析部１３０＞
　周期性分析部１３０は、係数列X[1]，…，X[N]を入力とし、当該係数列X[1]，…，X[N]の周期Ｔを求め、周期Ｔを出力する（Ｓ１３０）。
　周期Ｔは、入力音響信号に由来する周波数領域の係数列、例えば、係数列X[1]，…，X[N]、の周期性を有する成分の間隔（係数列が周期的に大きな値となる間隔）に対応する情報である。以下では周期Ｔを間隔Ｔと表現する場合もあるが、表現上の違いだけであり、同じものである。Ｔは正値であり、整数であってもよいし、小数（例えば、5.0、5.25、5.5、5.75）であってもよい。
　また、周期性分析部１３０は、必要に応じて、係数列X[1]，…，X[N]を入力とし、周期性の程度を示す指標Ｓも求めて出力してもよい。この場合、例えば、係数列X[1]，…，X[N]の周期性を有する成分の部分のエネルギーとそれ以外の部分のエネルギーとの比など基づいて周期性の程度を示す指標Ｓを求める。この場合は、指標Ｓは周波数領域のサンプル列の周期性の程度を示す指標となる。なお、周期性を有する成分の大きさが大きいほど、すなわち、周期Ｔの整数倍のサンプルやその近傍にあるサンプルの振幅（サンプル値の絶対値）が大きいほど、周波数領域のサンプル列の「周期性の程度」は大きい。
　なお、周期性分析部１３０は、時間領域の入力音響信号から時間領域の周期を求め、求めた時間領域の周期を周波数領域の周期に変換することで、周期Ｔを求めてもよい。また、時間領域の周期を周波数領域の周期に変換したものの定数倍やその近傍の値を周期Ｔとして求めてもよい。同様に、周期性分析部１３０は、時間領域の入力音響信号から、例えば、時間領域の周期分だけ時間がずれた信号列間の相関の大きさ等に基づいて、周期性の程度を示す指標Ｓを求めてもよい。
　要は、時間領域の入力音響信号やそれに由来する周波数領域係数列から周期Ｔや指標Ｓを求める方法は、従来より様々な方法が存在するので、その何れの方法を選択して利用してもよい。
＜周期性包絡系列生成部１４０＞
　周期性包絡系列生成部１４０は、間隔Ｔを入力とし、周期性包絡系列P[1]，…，P[N]を出力する（Ｓ１４０）。周期性包絡系列P[1]，…，P[N]は、ピッチ周期に起因する周期でピークを持つ周波数領域の離散系列、すなわち調波モデルに対応する離散系列である。図３に周期性包絡系列P[1]，…，P[N]の例を示す。周期性包絡系列P[1]，…，P[N]は、図３に示された波形のように、間隔Ｔの整数倍の近傍の整数値であるインデックスと、その前後所定数のインデックスに対応する周期性包絡の値のみ正の値を持ち、それ以外は０であるような系列である。間隔Ｔの整数倍の近傍の整数値であるインデックスが周期的に最大値（ピーク）をとり、その前後所定数のインデックスに対応するP[n]の値は、そのインデックスｎがピークに対応するインデックスから離れるにつれて単調減少する関係にある。図３の横軸の1,2,…,は離散化サンプル点のインデックス（以下、「周波数インデックス」）を表す。
　例えば、ｎを周波数インデックスを表す変数とし、τを極大値（ピーク）に対応する周波数インデックスとして、ピークの形状は以下の関数Q(n)で表せる。ただし、間隔Ｔの小数点以下の桁数がＬ桁であり、間隔Ｔ’をＴ’＝Ｔ×２Ｌとする。
Figure JPOXMLDOC01-appb-M000002
 ｈはピークの高さを表し、間隔Ｔが大きいほどピークの高さが高くなる。また、ＰＤはピーク部分の幅を表し、間隔Ｔが大きいほど幅が広くなる。
　Ｕを１からピークの数までを示す正整数（例えば、図４の場合は１～１０）とし、ｖを１以上の整数（例えば、１から３程度）とし、floor(・)を小数点以下を切り捨てて整数値を返す関数とすると、周期性包絡系列P[n]は、例えば、
Figure JPOXMLDOC01-appb-M000003
 のように計算すればよい。ただし、（Ｕ×Ｔ’）／２Ｌ－ｖ≦ｎ≦（Ｕ×Ｔ’）／２Ｌ＋ｖである。例えば、L=2の場合、T=20.00であればT’=80、T=20.25であればT’=81、T=20.50であればT’=82、T=20.75であればT’=83である。なお、周期性包絡系列P[n]は、小数点第一位を四捨五入して整数値を返す関数Round(・)を用いて、
Figure JPOXMLDOC01-appb-M000004
 のように求めてもよい。
＜周期性統合包絡生成部１５０＞
　周期性統合包絡生成部１５０は、少なくとも、周期性包絡系列P[1]，…，P[N]、振幅スペクトル包絡系列W[1]，…，W[N]を入力とし、周期性統合包絡系列WM[1]，…，WM[N]を求める（Ｓ１５０）。具体的には、周期性統合包絡WM[n]を次式のように求める。
Figure JPOXMLDOC01-appb-M000005
 なお、δは、周期性統合包絡WM[n]と係数X[n]の絶対値系列の形状が近くなるように決定される値または予め定めた値である。
　周期性統合包絡生成部１５０において周期性統合包絡WM[n]と係数X[n]の絶対値系列の形状が近くなるようにδを決定する場合には、周期性統合包絡生成部１５０は、係数列X[1], …, X[N]も入力とし、決定されたδとそのときの周期性統合包絡系列WM[1]，…，WM[N]を出力すればよい。例えば、δは、いくつかのδの候補、例えば、０．４と０．８の２つをδの候補、の中から以下の式により定義されるＥが最小となるδに決めればよい。言い換えると、周期性統合包絡WM[n]と係数X[n]の絶対値系列の形状が近くなるδに決めればよい。
Figure JPOXMLDOC01-appb-M000006
　δは、周期性統合包絡WM[n]において周期性包絡P[n]をどの程度考慮するかを決める値である。言い換えれば、δは周期性統合包絡WM[n]における振幅スペクトル包絡W[n]と周期性包絡P[n]の混合比率を決める値といえる。また、式（９）のＧは係数列X[1]，…，X[N]の各係数X[n]の絶対値の系列と周期性統合包絡系列の逆数の系列との内積である。式（８）の~WM[n]は、周期性統合包絡の各値WM[n]をＧで正規化した正規化周期性統合包絡である。式（７）において、係数列X[1]，…，X[N]と正規化周期性統合包絡系列~WM[1]，…，~WM[N]の内積の４乗を計算しているのは、特に絶対値の大きい係数X[n]を強調して内積をとった値（距離）を小さくすることを意図している。つまり、係数列X[1]，…，X[N]の中で特に絶対値の大きい係数X[n]と周期性統合包絡WM[n]が近くなるようにδを決定することを意味している。
　また、周期性統合包絡生成部１５０において周期性の程度に応じてδの候補数を決定する場合には、周期性統合包絡生成部１５０は、周期性の程度を示す指標Ｓも入力とし、指標Ｓが、周期性が高いことに対応するフレームであることを示している場合には多くの候補数のδの候補の中から式（７）で定義されるＥが最小となるδを選び、指標Ｓが、周期性が低いことに対応するフレームであることを示している場合にはδを予め定めた値としてもよい。すなわち、周期性統合包絡生成部１５０において周期性の程度に応じてδの候補数を決定する場合には、周期性が高いほどδの候補の数を多くすればよい。
＜実施例１の発明の効果＞
　図４に同じ音響信号に対して生成された系列の違いを説明するための例を示す。図４（Ａ）に係数列X[1]，…，X[N]を補間した曲線の形状を、図４（Ｂ）に周期性包絡系列P[1]，…，P[N]を補間した曲線の形状を、図４（Ｃ）に平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を補間した曲線の形状を、図４（Ｄ）に周期性統合包絡系列WM[1]，…，WM[N]を補間した曲線の形状を示す。図４に示すとおり、周期性統合包絡系列WM[1]，…，WM[N]は、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]に比べて、係数列X[1]，…，X[N]に現れる周期的なピークを含んだ形状となっている。また、周期性統合包絡系列WM[1]，…，WM[N]は、スペクトル包絡を表す情報である線形予測係数または量子化済線形予測係数の他に、間隔Ｔ、または、間隔Ｔと値δの情報があれば生成できる。したがって、入力音響信号のスペクトル包絡を表す情報に少ない情報量を追加するだけで、入力音響信号のピッチ周期に起因する振幅のピークを、線形予測係数により求まるスペクトル包絡より高精度に表現することができる。すなわち、線形予測係数または量子化済線形予測係数と、間隔Ｔ、または、間隔Ｔと値δと、の少ない情報量で入力音響信号の振幅を高精度に推定することができることになる。なお、平滑化振幅スペクトル包絡~W[n]は次式で表現される包絡であり、γは振幅スペクトル係数を鈍らせる（平滑化する）ための１以下の正の定数である。
Figure JPOXMLDOC01-appb-M000007
　また、本発明の周期性統合包絡系列生成装置を符号化装置と復号装置で用いる場合には、符号化装置に含まれる周期性統合包絡系列生成装置以外の処理部で得られた量子化済線形予測係数^αｐを特定する符号（線形予測係数符号CL）と周期Ｔや時間領域の周期を特定する符号（周期符号CT）が復号装置に入力されるので、本発明の周期性統合包絡系列生成装置からはδの情報を示す符号を出力すれば、復号側の周期性統合包絡系列生成装置でも符号化側の周期性統合包絡系列生成装置で生成した周期性統合包絡系列と同じ周期性統合包絡系列を生成できる。したがって、符号化装置から復号装置に符号を送る際に増加する符号量は少ない。
＜実施例１の発明のポイント＞
　実施例１の周期性統合包絡系列生成装置１００では、周期性統合包絡生成部１５０が係数列X[1]，…，X[N]の周期性成分に基づいて、振幅スペクトル包絡系列W[1]，…，W[N]を変形し、周期性統合包絡系列WM[1]，…，WM[N]としている点が最も重要なポイントである。特に、係数列X[1]，…，X[N]の周期性の程度が大きいほど、すなわち、周期性を有する成分の大きさが大きいほど、振幅スペクトル包絡系列W[1]，…，W[N]のうち間隔Ｔ（周期）の整数倍およびそれらの近傍のサンプルの値を大きく変更すれば、上記の効果を得やすい。「近傍のサンプル」とは、間隔Ｔの整数倍の近傍の整数値であるインデックスで示されるサンプルである。また、「近傍」とは、例えば、式（３）～（５）などのあらかじめ定めた方法で決まる範囲とすればよい。
　また、係数列X[1]，…，X[N]の周期性を有する成分の間隔Ｔが広いほど、式（４）と式（５）に示された周期性包絡系列P[1]，…，P[N]は、大きい値を持ち、広い幅で、すなわち、間隔Ｔ（周期）の整数倍およびそれらの近傍の多くのサンプルで、０以外の値を持つ。つまり、周期性統合包絡生成部１５０は、係数列の周期性を有する成分の間隔Ｔが広いほど、振幅スペクトル包絡系列のうち間隔Ｔ（周期）の整数倍およびそれらの近傍のサンプルの値を大きく変更する。また、周期性統合包絡生成部１５０は、係数列の周期性を有する成分の間隔Ｔが広いほど、振幅スペクトル包絡系列を広い幅で、すなわち、間隔Ｔ（周期）の整数倍およびそれらの近傍の多くのサンプルで、サンプル値を変更する。「近傍の多くのサンプルで」とは、「近傍」に該当する範囲（あらかじめ定めた方法で決まる範囲）に存在するサンプルを多くすることを意味している。つまり、周期性統合包絡生成部１５０は、このように振幅スペクトル包絡系列を変形すれば、上記の効果を得やすい。
　なお、周期性統合包絡系列が持つ「入力音響信号のピッチ周期に起因する振幅のピークをより高精度に表現することができる。」という特徴を効果的に利用する例としては、符号化装置と復号装置があり、この例を実施例２，３に示している。ただし、周期性統合包絡系列の特徴の利用例は、符号化装置と復号装置以外にも、雑音除去装置やポストフィルタなどがあり得る。したがって、実施例１では周期性統合包絡系列生成装置を説明している。
［変形例１］（正規化係数列で周期性分析する例）
　変形例１の周期性統合包絡系列生成装置も図１に示す。また、変形例１の周期性統合包絡系列生成装置の処理フローも図２に示す。周期性統合包絡系列生成装置１０１は、周波数領域系列正規化部１１１も備える点と、スペクトル包絡系列計算部１２１、周期性分析部１３１が周期性統合包絡系列生成装置１００と異なり、その他の構成は同じである。以下では相違点についてのみ説明する。
＜スペクトル包絡系列計算部１２１＞
　スペクトル包絡系列計算部１２１は、振幅スペクトル包絡系列W[1]，…，W[N]だけではなく、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]も求める。
　具体的には、スペクトル包絡系列計算部１２１は、スペクトル包絡系列計算部１２０で示した（ｓｔｅｐ１），（ｓｔｅｐ２）に加えて、以下の手順の処理を行う。
　（ｓｔｅｐ３）量子化済線形予測係数^αｐのそれぞれにγｐを乗算し、量子化済平滑化線形予測係数^α1γ，^α2γ2，…，^αPγPを求める。γは平滑化するための１以下の正の定数である。そして、式（１０）によって、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を求める（Ｓ１２１）。もちろん、スペクトル包絡系列計算部１２０と同様に、量子化済線形予測係数^αｐに代えて線形予測係数αｐを用いてもよい。
＜周波数領域系列正規化部１１１＞
　周波数領域系列正規化部１１１は、係数列X[1]，…，X[N]の各係数を平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]の各係数で除算して正規化係数列XN[1]，…，XN[N]を得る。すなわち、ｎ＝１，…，Ｎに対して
　XN[n]=X[n]/~W[n]　　　（１１）
の計算を行い、正規化係数列XN[1]，…，XN[N]を求める（Ｓ１１１）。
＜周期性分析部１３１＞
　周期性分析部１３１は、正規化係数列XN[1]，…，XN[N]を入力とし、当該正規化係数列XN[1]，…，XN[N]の周期Ｔを求め、周期Ｔを出力する（Ｓ１３１）。すなわち、本変形例では、入力音響信号に由来する周波数領域の係数列である正規化係数列XN[1]，…，XN[N]の周期性を有する成分の間隔を周期Ｔとして求める。また、周期性分析部１３１は、必要に応じて、係数列X[1]，…，X[N]を入力とし、周期性の程度を示す指標Ｓも求めて出力してもよい。
　その他の処理は周期性統合包絡系列生成装置１００と同じである。したがって、実施例１と同様の効果が得られる。なお、周期性統合包絡系列生成装置１０１の場合は、周期性統合包絡生成部１５０は、振幅スペクトル包絡系列W[1]，…，W[N]の代わりに平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を用いてもよい。この場合は、式（６）の代わりに次式の計算となる。
Figure JPOXMLDOC01-appb-M000008
［変形例２］（外部から情報が入力される例）
　本発明の周期性統合包絡系列生成装置を符号化装置や復号装置が内部に備えている場合には、符号化装置や復号装置に含まれる周期性統合包絡系列生成装置以外の処理部で、係数列X[1]，…，X[N]、正規化係数列XN[1]，…，XN[N]、量子化済線形予測係数^αｐ、量子化済平滑化線形予測係数^αｐγｐ、振幅スペクトル包絡W[1]，…，W[N]、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]、周期Ｔ、指標Ｓなどが求められていることがある。このような場合は、周期性統合包絡系列生成装置に、周波数領域変換部、周波数領域正規化部、スペクトル包絡系列計算部、周期性分析部の少なくとも何れかを備えない構成としてもよい。この場合には、符号化装置内の周期性統合包絡系列生成装置以外の処理部から、量子化済線形予測係数^αｐを特定する符号（線形予測係数符号CL）、周期Ｔや時間領域の周期を特定する符号（周期符号CT）、指標Ｓを特定する符号、などが出力され、復号装置に入力される。したがって、この場合には、符号化装置内の周期性統合包絡系列生成装置からは、量子化済線形予測係数^αｐを特定する符号（線形予測係数符号CL）、周期Ｔや時間領域の周期を特定する符号（周期符号CT）、指標Ｓを特定する符号、などを出力する必要がない。
　また、本発明の周期性統合包絡系列生成装置を符号化装置や復号装置で用いる場合には、符号化装置と復号装置とで同一の周期性統合包絡系列を得られるようにする必要がある。したがって、符号化装置が出力し復号装置に入力される符号から特定可能な情報を用いて周期性統合包絡系列を得る必要がある。たとえば、符号化装置で用いる周期性統合包絡系列生成装置のスペクトル包絡系列計算部では、線形予測係数符号CLに対応する量子化済線形予測係数を用いて振幅スペクトル包絡系列を求め、復号装置で用いる周期性統合包絡系列生成装置のスペクトル包絡系列計算部では、符号化装置から出力されて復号装置に入力される線形予測係数符号CLに対応する復号線形予測係数を用いて振幅スペクトル包絡系列を求める必要がある。
　なお、符号化装置や復号装置で周期性統合包絡系列を用いる場合には、上述のように周期性統合包絡系列生成装置を内部に備えるのではなく、周期性統合包絡系列生成装置内の必要な処理部を符号化装置と復号装置に備えるようにすればよい。このような符号化装置や復号装置は実施例２で説明する。
≪符号化装置≫
　図５に実施例２の符号化装置の機能構成例を、図６に実施例２の符号化装置の処理フローを示す。符号化装置２００は、スペクトル包絡系列計算部２２１、周波数領域変換部１１０、周波数領域系列正規化部１１１、周期性分析部２３０、周期性包絡系列生成部１４０、周期性統合包絡生成部２５０、可変長符号化パラメータ計算部２６０、可変長符号化部２７０を備える。符号化装置２００は、入力された時間領域の音響ディジタル信号を入力音響信号x(t)とし、少なくとも量子化済線形予測係数^α1，…，^αPを示す符号ＣＬ、正規化係数列XN[1]，…，XN[N]の周期を表す間隔Ｔの符号ＣＴ、正規化係数列XN[1]，…，XN[N]を可変長符号化した可変長符号ＣＸを出力する。周波数領域系列正規化部１１１は実施例１変形例１と同じである。周波数領域変換部１１０と周期性包絡系列生成部１４０は実施例１と同じである。以下では異なる構成部について説明する。
＜スペクトル包絡系列計算部２２１＞
　スペクトル包絡系列計算部２２１は、入力音響信号x(t)の時間領域の線形予測に基づき、入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を計算し、計算の過程で得た量子化済線形予測係数^α1，…，^αPを示す符号ＣＬも求める（Ｓ２２１）。ただし、Ｎは正整数である。スペクトル包絡系列計算部２２１は、以下の手順で処理すればよい。
　（ｓｔｅｐ１）所定の時間区間であるフレーム単位で、入力音響信号に対する線形予測分析を行って線形予測係数α1，…，αPを求める。ただし、Ｐは予測次数を示す正整数である。例えば、全極型モデルであるＰ次自己回帰過程により、時刻ｔでの入力音響信号x(t)は、Ｐ時点まで遡った過去の自分自身の値x(t-1)，…，x(t-P)と予測残差e(t)と線形予測係数α1，…，αpによって式（１）で表される。
　（ｓｔｅｐ２）線形予測係数α1，…，αPを符号化して符号ＣＬを得て出力するとともに、符号ＣＬに対応する量子化済線形予測係数^α1，…，^αPを求める。また、量子化済線形予測係数^α1，…，^αPを用いてＮ点の入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]を求める。例えば、振幅スペクトル包絡系列の各値W[n]は、式（２）で求めることができる。なお、線形予測係数α1，…，αPを符号化して符号ＣＬを得る方法は、線形予測係数をLSPパラメータに変換して、LSPパラメータを符号化して符号ＣＬを得るなど、線形予測係数に変換可能な係数の何れを符号化して符号ＣＬを得る何れの方法を用いてもよい。
　（ｓｔｅｐ３）量子化済線形予測係数^αｐのそれぞれにγｐを乗算し、量子化済平滑化線形予測係数^α1γ，^α2γ2，…，^αPγPを求める。γはあらかじめ定めた平滑化するための１以下の正の定数である。そして、式（１０）によって、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を求める。
＜周期性分析部２３０＞
　周期性分析部２３０は、正規化係数列XN[1]，…，XN[N]を入力とし、当該正規化係数列XN[1]，…，XN[N]の間隔Ｔ（周期的に大きな値となる間隔）を求め、間隔Ｔと間隔Ｔを示す符号ＣＴを出力する（Ｓ２３０）。また、周期性分析部２３０は、必要に応じて、周期性の程度を示す指標Ｓ（すなわち、周波数領域のサンプル列の周期性の程度を示す指標）、も求めて出力する。また、周期性分析部２３０は、必要に応じて、指標Ｓを示す符号ＣＳも得て出力する。なお、指標Ｓと間隔Ｔ自体は実施例１変形例１の周期性分析部１３１と同じである。
＜周期性統合包絡生成部２５０＞
　周期性統合包絡生成部２５０は、少なくとも、周期性包絡系列P[1]，…，P[N]、振幅スペクトル包絡系列W[1]，…，W[N]を入力とし、周期性統合包絡系列WM[1]，…，WM[N]を求めて周期性統合包絡WM[n]を出力する。また、周期性統合包絡生成部１５０は、値δとして、予め定めた１つの値ではなく、予め定めた複数の候補値のうちの何れかを選択する場合には、係数列X[1], …, X[N]も入力とし、予め定めた複数の候補値のうち周期性統合包絡WM[n]と係数X[n]の絶対値系列の形状が近くなる候補値を値δとして求め、値δを示す符号Ｃδも出力する（Ｓ２５０）。
　周期性統合包絡WM[n]と値δは実施例１と同じであり、周期性統合包絡WM[n]は式（６），…，（９）のように求めればよい。周期性統合包絡生成部１５０において周期性の程度に応じてδの候補数を決定する場合には、周期性統合包絡生成部１５０は、周期性の程度を示す指標Ｓも入力とし、指標Ｓが周期性が高いことに対応するフレームの場合には多くの候補数のδの候補の中から式（７）で定義されるＥが最小となるδを選び、指標Ｓが周期性が低いことに対応するフレームである場合にはδを１つの予め定めた値としてもよい。なお、δを予め定めた値にする場合は、値δを示す符号Ｃδを出力する必要はない。
＜可変長符号化パラメータ計算部２６０＞
　可変長符号化パラメータ計算部２６０は、周期性統合包絡系列WM[1]，…，WM[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と正規化係数列XN[1]，…，XN[N]を入力とし、可変長符号化パラメータｒｎを求める（Ｓ２６０）。可変長符号化パラメータ計算部２６０は、周期性統合包絡系列WM[1]，…，WM[N]から求めた振幅値に依存して可変長符号化パラメータｒｎを計算することを特徴としている。
　可変長符号化パラメータは、符号化対象の信号、すなわち、正規化係数列XN[1]，…，XN[N]の各係数の振幅の取り得る範囲を特定するパラメータである。例えば、ライス符号化の場合にはライスパラメータが可変長符号化パラメータに相当し、算術符号化の場合は符号化対象の信号の振幅の取り得る範囲が可変長符号化パラメータに相当する。
　１サンプルごとに可変長符号化を行う場合には、正規化係数列の各係数XN[n]について可変長符号化パラメータが計算される。複数のサンプルからなるサンプル群ごとに（例えば２サンプルずつ）まとめて可変長符号化を行う場合には、サンプル群ごとに可変長符号化パラメータが計算される。つまり、可変長符号化パラメータ計算部２６０は、正規化係数列の一部である正規化部分係数列ごとに、可変長符号化パラメータｒｎを計算する。ここで、正規化部分係数列は複数個あり、複数個の正規化部分係数列には正規化係数列の係数が重複されずに含まれるものとする。以下に、１サンプルごとにライス符号化を行う場合を例に、可変長符号化パラメータの計算方法を説明する。
　（ｓｔｅｐ１）正規化係数列XN[1],…,XN[N]の各係数の振幅の平均の対数を、基準となるライスパラメータｓｂ（基準となる可変長符号化パラメータ）として次式のように算出する。
Figure JPOXMLDOC01-appb-M000009
 ｓｂはフレームごとに１度だけ符号化されて、基準となるライスパラメータ（基準となる可変長符号化パラメータ）に対応する符号Ｃｓｂとして復号装置４００に伝送される。あるいは復号装置４００に伝送される別の情報から正規化係数列XN[1],…,XN[N]の振幅の平均値を推定できる場合は、符号化装置２００と復号装置４００で共通に振幅の平均値の推定値からｓｂを近似的に決定する方法を決めておいてもよい。例えば、包絡の傾きを表すパラメータ、区分帯域ごとの平均包絡の大きさを表すパラメータを別途使う符号化の場合には、復号装置４００に伝送される別の情報から振幅の平均値を推定できる。この場合は、ｓｂを符号化し、基準となるライスパラメータに対応する符号Ｃｓｂを復号装置４００へ出力しなくてもよい。
　（ｓｔｅｐ２）下記式により閾値θを算出する。
Figure JPOXMLDOC01-appb-M000010
 θは、周期性統合包絡系列の各値WM[n]を平滑化振幅スペクトル包絡系列の各値~W[n]で除算した値の振幅の平均の対数である。
　（ｓｔｅｐ３） |WM[n]/~W[n]|がθより大きいほど、正規化係数XN[n]をライス符号化するためのライスパラメータｒｎをｓｂよりも大きな値として決定する。|WM[n]/~W[n]|がθより小さいほど、正規化係数XN[n]をライス符号化するためのライスパラメータｒｎをｓｂよりも小さな値として決定する。
　（ｓｔｅｐ４）ｓｔｅｐ３の処理を全てのn=1,2,…,Nについて繰り返して、各XN[n]についてのライスパラメータｒｎを求める。
＜可変長符号化部２７０＞
　可変長符号化部２７０は、可変長符号化パラメータ計算部２６０で求めた可変長符号化パラメータｒｎを用いて正規化係数列XN[1]，…，XN[N]を可変長符号化し、可変長符号ＣＸを出力する（Ｓ２７０）。例えば、可変長符号化部２７０は、可変長符号化パラメータ計算部２６０で求めたライスパラメータｒｎを用いて正規化係数列XN[1]，…，XN[N]をライス符号化し、得られた符号を可変長符号ＣＸとして出力する。可変長符号化パラメータ計算部２６０で求めたライスパラメータｒｎは、周期性統合包絡系列の振幅値に依存する可変長符号化パラメータであり、周期性統合包絡系列の値が大きい周波数ほど大きな値となっている。ライス符号化は、振幅値に依存する可変長符号化の公知技術のうちの１つであり、ライスパラメータｒｎを用いて振幅値に依存する可変長符号化を行うものである。また、周期性統合包絡生成部２５０で生成した周期性統合包絡系列は、入力音響信号のスペクトル包絡を高精度に表現するものである。すなわち、可変長符号化部２７０は、周期性統合包絡系列の値が大きい周波数ほど、前記入力音響信号の周波数領域の係数列であるX[1]，…，X[N]の振幅が大きいとことを前提に、正規化係数列XN[1]，…，XN[N]を可変長符号化していることになり、言い換えれば、可変長符号化パラメータを用いて、振幅値に依存する可変長符号化により、正規化係数列XN[1]，…，XN[N]を符号化していることになる。ここでいう振幅値とは、符号化対象の係数列の平均振幅値、係数列に含まれる各係数の振幅の推定値、係数列の振幅の包絡の推定値などである。
　符号化装置２００は、このような処理によって得られた量子化済線形予測係数^α1，…，^αPを示す符号ＣＬ、間隔Ｔを示す符号ＣＴ、正規化係数列XN[1]，…，XN[N]を可変長符号化した可変長符号ＣＸを出力する。また、必要に応じて値δを示す符号Ｃδと基準となる可変長符号化パラメータｓｂを示す符号Ｃｓｂも出力する。符号化装置２００から出力された符号は、復号装置４００に入力される。
［符号化装置の変形例１］（外部から情報が入力される例）
　なお、符号化装置としては、周期性包絡系列生成部１４０と周期性統合包絡生成部２５０と可変長符号化パラメータ計算部２６０と可変長符号化部２７０だけを備え、符号化装置の外部で生成された平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と、正規化係数列XN[1]，…，XN[N]、間隔Ｔと、必要に応じて振幅スペクトル包絡系列W[1]，…， W[N]と、必要に応じて指標Ｓとを入力とし、可変長符号ＣＸを出力してもよい。
［符号化装置の変形例２]（係数列X[n]から間隔Ｔを求める例）
　上述の周期性分析部２３０では正規化係数列XN[1]，…，XN[N]を入力として間隔Ｔを求めているが、周期性分析部２３０では周波数領域変換部１１０が出力した係数列X[1]，…，X[N]を入力として間隔Ｔを求めてもよい。この場合は、実施例１の周期性分析部１３０と同じ方法で間隔Ｔを求める。
≪復号装置≫
　図７に実施例２の復号装置の機能構成例を、図８に実施例２の復号装置の処理フローを示す。復号装置４００は、スペクトル包絡系列計算部４２１、周期性包絡系列生成部４４０、周期性統合包絡生成部４５０、可変長符号化パラメータ計算部４６０、可変長復号部４７０、周波数領域系列逆正規化部４１１、周波数領域逆変換部４１０を備える。復号装置４００は、量子化済線形予測係数^α1，…，^αPを示す符号ＣＬ、間隔Ｔを示す符号ＣＴ、正規化係数列XN[1]，…，XN[N]を可変長符号化した可変長符号ＣＸを受け取り、音響信号を出力する。なお、必要に応じて値δを示す符号Ｃδと基準となる可変長符号化パラメータｓｂを示す符号Ｃｓｂと指標Ｓを示す符号ＣＳも受け取る。以下に、各構成部の詳細を示す。
＜スペクトル包絡系列計算部４２１＞
　スペクトル包絡系列計算部４２１は、符号ＣＬを入力とし、振幅スペクトル包絡系列W[1]，…，W[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を計算する（Ｓ４２１）。より具体的には、以下の手順で処理すればよい。
　（ｓｔｅｐ１）符号ＣＬを復号し、復号線形予測係数^α1，…，^αPを得る。
　（ｓｔｅｐ２）復号線形予測係数^α1，…，^αPを用いてＮ点の振幅スペクトル包絡系列W[1]，…，W[N]を求める。例えば、振幅スペクトル包絡系列の各値W[n]は、式（２）で求めることができる。
　（ｓｔｅｐ３）復号線形予測係数^αｐのそれぞれにγｐを乗算し、復号平滑化線形予測係数^α1γ，^α2γ2，…，^αPγPを求める。γはあらかじめ定めた平滑化するための１以下の正の定数である。そして、式（１０）によって、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を求める。
＜周期性包絡系列生成部４４０＞
　周期性包絡系列生成部４４０は、間隔Ｔを示す符号ＣＴを入力とし、符号ＣＴを復号し、間隔Ｔを得る。そして、符号化装置２００の周期性包絡系列生成部１４０と同じ方法で周期性包絡系列P[1]，…，P[N]を求め、出力する（Ｓ４４０）。
＜周期性統合包絡生成部４５０＞
　周期性統合包絡生成部４５０には、周期性包絡系列P[1]，…，P[N]、振幅スペクトル包絡系列W[1]，…，W[N]、符号Ｃδ、符号ＣＳが入力される。ただし、符号Ｃδ、符号ＣＳは入力されない場合もある。周期性統合包絡生成部４５０は、符号Ｃδを復号し、値δを取得する。ただし、符号Ｃδが入力されない場合は、符号Ｃδの復号は行わず、周期性統合包絡生成部４５０に予め記憶された値δを取得する。なお、周期性統合包絡生成部４５０は、符号ＣＳが入力された場合には、符号ＣＳを復号して指標Ｓを取得し、取得した指標Ｓが、周期性が高いことに対応するフレームの場合には符号Ｃδを復号して値δを取得し、取得した指標Ｓが、周期性が低いことに対応するフレームである場合には符号Ｃδの復号は行わず、周期性統合包絡生成部４５０に予め記憶された値δを取得する。そして、周期性統合包絡生成部４５０は、式（６）によって、周期性統合包絡系列WM[1]，…，WM[N]を求める。（Ｓ４５０）
＜可変長符号化パラメータ計算部４６０＞
　可変長符号化パラメータ計算部４６０は、周期性統合包絡系列WM[1]，…，WM[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と符号Ｃｓｂを入力とし、可変長符号化パラメータｒｎを得る（Ｓ４６０）。ただし、復号装置４００に伝送される別の情報から振幅の平均値を推定できる場合は、別の情報から推定した振幅の平均値の推定値からｓｂを近似的に決定する方法を決めておいてもよい。この場合は、符号Ｃｓｂは入力されない。以下に、１サンプルごとにライス復号を行う場合を例に、可変長符号化パラメータの計算方法を説明する。
　（ｓｔｅｐ１）符号Ｃｓｂを復号して、基準となるライスパラメータｓｂ（基準となる可変長符号化パラメータ）を得る。なお、符号化装置２００と復号装置４００で共通に振幅の平均値の推定値からｓｂを近似的に決定する方法を決めている場合は、その方法で求める。
　（ｓｔｅｐ２）閾値θを式（１４）で算出する。
　（ｓｔｅｐ３） |WM[n]/~W[n]|がθより大きいほど、ライスパラメータｒｎをｓｂよりも大きな値として、符号化装置２００の可変長符号化パラメータ計算部２６０と同じ方法で決定する。|WM[n]/~W[n]|がθより小さいほど、ライスパラメータｒｎをｓｂよりも小さな値として、符号化装置２００の可変長符号化パラメータ計算部２６０と同じ方法で決定する。
　（ｓｔｅｐ４）ｓｔｅｐ３の処理を全てのn=1,2,…,Nについて繰り返して、各XN[n]についてのライスパラメータｒｎを求める。
＜可変長復号部４７０＞
　可変長復号部４７０は、可変長符号化パラメータ計算部４６０で求めた可変長符号化パラメータｒｎを用いて可変長符号ＣＸを復号して復号正規化係数列^XN[1]，…，^XN[N]を得る（Ｓ４７０）。例えば、可変長復号部４７０は、可変長符号化パラメータ計算部４６０で求めたライスパラメータｒｎを用いて可変長符号ＣＸを復号して復号正規化係数列^XN[1]，…，^XN[N]を得る。可変長復号部４７０の復号方法は、可変長符号化部２７０の符号化方法に対応するものである。
＜周波数領域系列逆正規化部４１１＞
　周波数領域系列逆正規化部４１１は、復号正規化係数列^XN[1]，…，^XN[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を入力とし、
　^X[n]=^XN[n]・~W[n]　　　（１５）
のように、復号係数列^X[1]，…，^X[N]を求めて出力する（Ｓ４１１）。
＜周波数領域逆変換部４１０＞
周波数領域逆変換部４１０は、復号係数列^X[1]，…，^X[N]を入力とし、復号係数列^X[1]，…，^X[N]を所定の時間区間であるフレーム単位の音響信号（時間領域）に変換する（Ｓ４１０）。
［復号装置の変形例１］（外部から情報が入力される例）
　なお、復号装置としては、周期性包絡系列生成部４４０と周期性統合包絡生成部４５０と可変長符号化パラメータ計算部４６０と可変長復号部４７０だけを備え、復号装置に必要に応じて入力される符号Ｃδと符号Ｃｓｂに加えて、復号装置の外部で得られた平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]、振幅スペクトル包絡系列W[1]，…，W[N]、間隔Ｔ、必要に応じて指標Ｓも入力とし、正規化係数列XN[1]，…，XN[N]を出力し、外部で平滑化振幅スペクトル包絡系列を乗算して時間領域の音響信号に変換してもよい。
＜実施例２の発明の効果＞
　可変長符号化は、符号化対象の入力値の振幅の取りうる範囲に合わせて適応的に符号を決定することで符号化効率を向上させる符号化方法である。実施例２では周波数領域の係数列である正規化係数列XN[1]，…，XN[N]を符号化対象としているが、符号化対象の係数列に含まれる各係数の振幅の情報をより正確に用いて求めた可変長符号化パラメータを用いて可変長符号化をすれば符号化装置が行う可変長符号化自体の符号化効率は高くなる。しかし、復号装置が可変長符号化パラメータを求めるために、符号化装置から復号装置に対して符号化対象の係数列に含まれる各係数の振幅の情報をより正確に送る必要があり、その分だけ符号化装置から復号装置に送る符号量が増大してしまう。
　符号量の増大を抑えるためには、少ない符号量の符号から符号化対象の係数列に含まれる各係数の振幅の推定値を得る方法が必要である。実施例２の周期性統合包絡系列WM[1]，…，WM[N]は係数列X[1]，…，X[N]を高精度に近似するので、|WM[1]/~W[1]|，…，|WM[N]/~W[N]|は可変長符号化対象の係数であるXN[1], XN[2],…, XN[N]の振幅包絡を高精度に近似できる。つまり、|WM[1]/~W[1]|，…，|WM[N]/~W[N]|は、符号化対象の各係数の振幅と正の相関を持つ系列となっている。
　また、|WM[1]/~W[1]|, |WM[2]/~W[2]|,…,|WM[N]/~W[N]|を復号装置側で復元するために必要な情報は、
・量子化済線形予測係数^α1，…，^αPの情報（符号ＣＬ）
・間隔Ｔを示す情報（符号ＣＴ）
・値δを示す情報（符号Ｃδ）
である。すなわち、実施例２の符号化装置と復号装置によれば、符号化装置に入力された入力音響信号のピッチ周期に起因する振幅のピークを含む包絡を、符号ＣＬ、符号ＣＴ、符号Ｃδのみの少ない情報量で、復号装置で再現することが可能となる。
　なお、実施例２の符号化装置と復号装置は、線形予測やピッチ予測を伴う符号化及び復号を行う符号化装置及び復号装置と併用して用いられることが多い。この場合は、符号ＣＬと符号ＣＴは、符号化装置２００外にある線形予測やピッチ予測を伴う符号化を行う符号化装置から、復号装置４００外にある線形予測やピッチ予測を伴う復号を行う復号装置に送られている符号である。したがって、符号化装置側に入力された入力音響信号のピッチ周期に起因する振幅のピークを含む包絡を復号装置側で復元するために符号化装置２００から復号装置４００に送る必要があるのは符号Ｃδである。符号Ｃδの符号量は小さく（それぞれ、せいぜい３ビット程度であり、１ビットでも効果が得られる）、符号化対象の正規化係数列に含まれる部分系列ごとの可変長符号化パラメータに対応する符号の総符号量よりも少ない。
　よって、実施例２の符号化装置、復号装置によれば、少ない符号量の増加で、符号化効率を向上させることができる。
＜実施例２の発明のポイント＞
　上述の効果を得るというポイントで実施例２の符号化装置、復号装置を考えると、符号化装置２００が、
　・所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部２５０
　・周期性統合包絡系列の値が大きい周波数ほど、入力音響信号の振幅が大きいことを前提に、入力音響信号に由来する周波数領域の系列を符号化する可変長符号化部２７０
を有し、復号装置４００が、
　・線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部４５０
　・周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、可変長符号を復号して周波数領域の系列を得る可変長復号部４７０、
を有することを特徴とすればよい。なお、「周期性統合包絡系列の値が大きい周波数ほど、入力音響信号の振幅が大きいことを前提に」と「周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に」とは、周期性統合包絡系列が、入力音響信号または音響信号の振幅の大きい周波数において大きい値になることを特徴としていることを示している。また、「入力音響信号に由来する」とは、入力音響信号から求められることや入力音響信号に対応していることを意味している。例えば、係数列X[1]，…，X[N]や正規化係数列XN[1]，…，XN[N]は、入力音響信号に由来する周波数領域の系列である。
≪符号化装置≫
　図９に実施例３の符号化装置の機能構成例を、図１０に実施例３の符号化装置の処理フローを示す。符号化装置３００は、スペクトル包絡系列計算部２２１、周波数領域変換部１１０、周波数領域系列正規化部１１１、周期性分析部３３０、周期性包絡系列生成部１４０、周期性統合包絡生成部２５０、可変長符号化パラメータ計算部２６０、第２可変長符号化パラメータ計算部３８０、可変長符号化部３７０を備える。符号化装置３００は、入力された時間領域の音響ディジタル信号を入力音響信号x(t)とし、少なくとも量子化済線形予測係数^α1，…，^αPを示す符号ＣＬ、正規化係数列XN[1]，…，XN[N]の周期を表す間隔Ｔの符号ＣＴ、係数列X[1]，…，X[N]または正規化係数列XN[1]，…，XN[N]の周期性の程度を示す所定の指標Ｓと指標Ｓを示す符号ＣＳ、正規化係数列XN[1]，…，XN[N]を可変長符号化した可変長符号ＣＸを出力する。周波数領域系列正規化部１１１は実施例１変形例１と同じである。周波数領域変換部１１０と周期性包絡系列生成部１４０は実施例１と同じである。振幅スペクトル包絡系列計算部２２１、周期性統合包絡生成部２５０、可変長符号化パラメータ計算部２６０は、実施例２と同じである。以下では異なる構成部について説明する。
＜周期性分析部３３０＞
　周期性分析部３３０は、正規化係数列XN[1]，…，XN[N]を入力とし、
当該正規化係数列XN[1]，…，XN[N]の周期性の程度を示す指標Ｓと間隔Ｔ（周期的に大きな値となる間隔）とを求め、指標Ｓと指標Ｓを示す符号ＣＳと間隔Ｔと間隔Ｔを示す符号ＣＴを出力する（Ｓ３３０）。なお、指標Ｓと間隔Ｔ自体は実施例１変形例１の周期性分析部１３１と同じである。
　そして、符号化装置３００では、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、可変長符号化パラメータ計算部２６０が可変長符号化パラメータｒｎを計算し、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、第２可変長符号化パラメータ計算部３８０が可変長符号化パラメータｒｎを計算する（Ｓ３９０）。「あらかじめ定めた周期性の程度が大きいことを示す範囲」は、例えば、指標Ｓが所定の閾値以上のときとすればよい。
＜第２可変長符号化パラメータ計算部３８０＞
　第２可変長符号化パラメータ計算部３８０は、振幅スペクトル包絡系列W[1]，…，W[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と正規化係数列XN[1]，…，XN[N]を入力とし、可変長符号化パラメータｒｎを求める（Ｓ３８０）。可変長符号化パラメータ計算部２６０は、周期性統合包絡系列WM[1]，…，WM[N]から求めた振幅値に依存して可変長符号化パラメータｒｎを計算することを特徴としているのに対して、第２可変長符号化パラメータ計算部３８０は、振幅スペクトル包絡系列から求めた振幅値に依存して可変長符号化パラメータを計算することを特徴としている。以下に、１サンプルごとにライス符号化を行う場合を例に、可変長符号化パラメータの計算方法を説明する。
　（ｓｔｅｐ１）正規化係数列XN[1],…,XN[N]の各係数の振幅の平均の対数を、基準となるライスパラメータｓｂ（基準となる可変長符号化パラメータ）として式（１３）のように算出する。この処理は、可変長符号化パラメータ計算部２６０と同じである。
　（ｓｔｅｐ２）下記式により閾値θを算出する。
Figure JPOXMLDOC01-appb-M000011
 θは、振幅スペクトル包絡系列の各値W[n]を平滑化振幅スペクトル包絡系列の各値~W[n]で除算した値の振幅の平均の対数である。
　（ｓｔｅｐ３） |W[n]/~W[n]|がθより大きいほど、正規化係数XN[n]をライス符号化するためのライスパラメータｒｎをｓｂよりも大きな値として決定する。|W[n]/~W[n]|がθより小さいほど、正規化係数XN[n]をライス符号化するためのライスパラメータｒｎをｓｂよりも小さな値として決定する。
　（ｓｔｅｐ４）ｓｔｅｐ３の処理を全てのn=1,2,…,Nについて繰り返して、各XN[n]についてのライスパラメータｒｎを求める。
＜可変長符号化部３７０＞
　可変長符号化部３７０は、可変長符号化パラメータｒｎを用いて正規化係数列XN[1],…,XN[N]を可変長符号化し、可変長符号ＣＸを出力する（Ｓ３７０）。ただし、可変長符号化パラメータｒｎは、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、可変長符号化パラメータ計算部２６０が計算した可変長符号化パラメータｒｎであり、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、第２可変長符号化パラメータ計算部３８０が計算した可変長符号化パラメータｒｎである。
　符号化装置３００は、このような処理によって得られた量子化済線形予測係数^α1，…，^αPを示す符号ＣＬ、周期性の程度を示す指標Ｓを示す符号ＣＳ、間隔Ｔを示す符号ＣＴ、正規化係数列XN[1],…,XN[N]を可変長符号化した可変長符号ＣＸを出力し、復号側に送信する。また、必要に応じて値δを示す符号Ｃδと基準となる可変長符号化パラメータｓｂを示す符号Ｃｓｂも出力し、復号側に送信する。
［符号化装置の変形例１］（外部から情報が入力される例）
　なお、符号化装置としては、周期性包絡系列生成部１４０と周期性統合包絡生成部２５０と可変長符号化パラメータ計算部２６０と第２可変長符号化パラメータ計算部３８０と可変長符号化部３７０だけを備え、符号化装置の外部で生成された平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と正規化係数列XN[1],…,XN[N]、間隔Ｔと、必要に応じて振幅スペクトル包絡系列W[1]，…， W[N]と、必要に応じて指標Ｓとを入力とし、可変長符号ＣＸを出力してもよい。
［符号化装置の変形例２]（係数列X[n]から間隔Ｔを求める例）
　上述の周期性分析部３３０では正規化係数列XN[1]，…，XN[N]を入力として間隔Ｔを求めているが、周期性分析部３３０では周波数領域変換部１１０が出力した係数列X [1]，…，X [N]を入力として間隔Ｔを求めてもよい。この場合は、実施例１の周期性分析部１３０と同じ方法で間隔Ｔを求める。
≪復号装置≫
　図１１に実施例３の復号装置の機能構成例を、図１２に実施例３の復号装置の処理フローを示す。復号装置５００は、スペクトル包絡系列計算部４２１、指標復号部５３０、周期性包絡系列生成部４４０、周期性統合包絡生成部４５０、可変長符号化パラメータ計算部４６０、第２可変長符号化パラメータ計算部５８０、可変長復号部５７０、周波数領域系列逆正規化部４１１、周波数領域逆変換部４１０を備える。復号装置５００は、量子化済線形予測係数^α1，…，^αPを示す符号ＣＬ、指標Ｓを示す符号ＣＳ、間隔Ｔを示す符号ＣＴ、正規化係数列XN[1]，…，XN[N]を可変長符号化した可変長符号ＣＸを受け取り、音響信号を出力する。なお、必要に応じて値δを示す符号Ｃδと基準となる可変長符号化パラメータｓｂを示す符号Ｃｓｂも受け取る。スペクトル包絡系列計算部４２１、周期性包絡系列生成部４４０、周期性統合包絡生成部４５０、可変長符号化パラメータ計算部４６０、周波数領域系列逆正規化部４１１、周波数領域逆変換部４１０は実施例２と同じである。以下では異なる構成部について説明する。
＜指標復号部５３０＞
　指標復号部５３０は、符号ＣＳを復号し、指標Ｓを得る。復号装置５００では、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、可変長符号化パラメータ計算部４６０が可変長符号化パラメータｒｎを計算し、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、第２可変長符号化パラメータ計算部５８０が可変長符号化パラメータｒｎを計算する（Ｓ５９０）。なお、「あらかじめ定めた周期性の程度が大きいことを示す範囲」は、符号化装置３００と同じ範囲である。
＜第２可変長符号化パラメータ計算部５８０＞
　第２可変長符号化パラメータ計算部５８０は、振幅スペクトル包絡系列W[1]，…，W[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と符号Ｃｓｂを入力とし、可変長符号化パラメータｒｎを求める（Ｓ５８０）。ただし、復号装置５００に伝送される別の情報から振幅の平均値を推定できる場合は、別の情報から推定した振幅の平均値の推定値からｓｂを近似的に決定する方法を決めておいてもよい。この場合は、符号Ｃｓｂは入力されない。以下に、１サンプルごとにライス復号を行う場合を例に、可変長符号化パラメータの計算方法を説明する。
　（ｓｔｅｐ１）符号Ｃｓｂを復号して、基準となるライスパラメータｓｂ（基準となる可変長符号化パラメータ）を得る。なお、符号化装置３００と復号装置５００で共通に振幅の推定値からｓｂを近似的に決定する方法を決めている場合は、その方法で求める。
　（ｓｔｅｐ２）閾値θを式（１６）で算出する。
　（ｓｔｅｐ３） |W[n]/~W[n]|がθより大きいほど、ライスパラメータｒｎをｓｂよりも大きな値として、符号化装置３００の第２可変長符号化パラメータ計算部３８０と同じ方法で決定する。|W[n]/~W[n]|がθより小さいほど、ライスパラメータｒｎをｓｂよりも小さな値として、符号化装置３００の第２可変長符号化パラメータ計算部３８０と同じ方法で決定する。
　（ｓｔｅｐ４）ｓｔｅｐ３の処理を全てのn=1,2,…,Nについて繰り返して、各XN[n]についてのライスパラメータｒｎを求める。
＜可変長復号部５７０＞
　可変長復号部５７０は、可変長符号化パラメータｒｎを用いて可変長符号ＣＸを復号して復号正規化係数列^XN[1]，…，^XN[N]を求める（Ｓ５７０）。ただし、可変長符号化パラメータｒｎは、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、可変長符号化パラメータ計算部４６０が計算した可変長符号化パラメータｒｎであり、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、第２可変長符号化パラメータ計算部５８０が計算した可変長符号化パラメータｒｎである。
［復号装置の変形例１］（外部から情報が入力される例）
　なお、復号装置としては、周期性包絡系列生成部４４０と周期性統合包絡生成部４５０と可変長符号化パラメータ計算部４６０と第２可変長符号化パラメータ計算部５８０と可変長復号部５７０だけを備え、復号装置に必要に応じて入力される符号Ｃδと符号Ｃｓｂに加えて、復号装置の外部で得られた平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]、振幅スペクトル包絡系列W[1]，…，W[N]、間隔Ｔ、指標Ｓも入力とし、正規化係数列XN[1]，…，XN[N]を出力し、外部で平滑化振幅スペクトル包絡系列を乗算して時間領域の音響信号に変換してもよい。
＜実施例３の発明の効果＞
　入力音響信号の周期性の程度が小さい場合には、入力音響信号のピッチ周期に起因する振幅のピークは小さい。そこで、実施例３の符号化装置、復号装置は、符号化の対象となる音響信号の周期性の程度が大きい場合には周期性統合包絡系列を用いて可変長符号化パラメータを求め、符号化の対象となる音響信号の周期性の程度が大きくない場合には振幅スペクトル包絡系列を用いて可変長符号化パラメータを求めるため、より適した可変長符号化パラメータを用いて可変長符号化でき、符号化精度を上げることができるという効果がある。
　上述の実施例１～３では、振幅スペクトル包絡系列、平滑化振幅スペクトル包絡系列、周期性統合包絡系列等について振幅の系列を用いる例を説明したが、振幅の系列に代えてパワーの系列、すなわち、W[n]、~W[n]、WM[n]としてパワースペクトル包絡系列、平滑化パワースペクトル包絡系列、パワーの系列である周期性統合包絡系列を用いてもよい。
［プログラム、記録媒体］
　上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。
　また、上述の構成をコンピュータによって実現する場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。
　この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。
　また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ－ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。
　このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。
　また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、本装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。
１００、１０１　周期性統合包絡系列生成装置
１１０　周波数領域変換部　　　　　　１１１　周波数領域係数正規化部
１２０、１２１、２２１、４２１　スペクトル包絡系列計算部
１３０，１３１、２３０、３３０　周期性分析部
１４０、４４０　周期性包絡系列生成部
１５０、２５０、４５０　周期性統合包絡生成部
２００、３００　符号化装置
２６０、３６０、４６０　可変長符号化パラメータ計算部
２７０、３７０　可変長符号化部
３８０、５８０　第２可変長符号化パラメータ計算部
４００、５００　復号装置
４１０　周波数領域逆変換部　　　　　４１１　周波数領域系列逆正規化部
４７０、５７０　可変長復号部　　　　５３０　指標復号部

Claims (20)

  1. 　所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、前記入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部と、
    　前記周期性統合包絡系列から振幅値に依存する可変長符号化パラメータを計算する可変長符号化パラメータ計算部と、
    　可変長符号化パラメータを用いて、振幅値に依存する可変長符号化により、前記入力音響信号に由来する周波数領域の系列を符号化する可変長符号化部、
    　を備える符号化装置。
  2. 　所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、前記入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部と、
    　前記周期性統合包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化する可変長符号化部、
    　を備える符号化装置。
  3. 　請求項１記載の符号化装置であって、
    　前記スペクトル包絡系列から振幅値に依存する可変長符号化パラメータを計算する第２可変長符号化パラメータ計算部も備え、
    　前記可変長符号化部は、
    　前記入力音響信号の周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記可変長符号化パラメータ計算部が計算した可変長符号化パラメータを用いて符号化し、
    　前記周期性の程度を示す指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記第２可変長符号化パラメータ計算部が計算した可変長符号化パラメータを用いて符号化する
    　符号化装置。
  4. 　請求項２記載の符号化装置であって、
    　前記可変長符号化部は、
    　前記入力音響信号の周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記周期性統合包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化し、
    　前記指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記スペクトル包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化する
    　符号化装置。
  5. 　線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部と、
    　前記周期性統合包絡系列から振幅値に依存する可変長符号化パラメータを計算する可変長符号化パラメータ計算部と、
    　可変長符号化パラメータを用いて復号する可変長復号部、
    　を備える復号装置。
  6. 　線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部と、
    　前記周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、可変長符号を復号して周波数領域の系列を得る可変長復号部、
    　を備える復号装置。
  7. 　請求項５記載の復号装置であって、
    　入力された指標符号を復号して周期性の程度を示す指標を得る指標復号部と、
    　前記スペクトル包絡系列から振幅値に依存する可変長符号化パラメータを計算する第２可変長符号化パラメータ計算部
    　も備え、
    　前記可変長復号部は、
    　前記指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記可変長符号化パラメータ計算部が計算した可変長符号化パラメータを用いて復号し、
    　前記指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記第２可変長符号化パラメータ計算部が計算した可変長符号化パラメータを用いて復号する
    　復号装置。
  8. 　請求項６記載の復号装置であって、
    　前記可変長復号部は、
    　入力された周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、前記可変長符号を復号して周波数領域の系列を得て、
    　前記指標があらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、前記スペクトル包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、前記可変長符号を復号して周波数領域の系列を得る
    　復号装置。
  9. 　所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、前記入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成ステップと、
    　前記周期性統合包絡系列から振幅値に依存する可変長符号化パラメータを計算する可変長符号化パラメータ計算ステップと、
    　可変長符号化パラメータを用いて、振幅値に依存する可変長符号化により、前記入力音響信号に由来する周波数領域の系列を符号化する可変長符号化ステップ、
    　を実行する符号化方法。
  10. 　所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、前記入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成ステップと、
    　前記周期性統合包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化する可変長符号化ステップ、
    　を実行する符号化方法。
  11. 　請求項９記載の符号化方法であって、
    　前記スペクトル包絡系列から振幅値に依存する可変長符号化パラメータを計算する第２可変長符号化パラメータ計算ステップも有し、
    　前記可変長符号化ステップは、
    　前記入力音響信号の周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記可変長符号化パラメータ計算ステップで計算した可変長符号化パラメータを用いて符号化し、
    　前記周期性の程度を示す指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記第２可変長符号化パラメータ計算ステップで計算した可変長符号化パラメータを用いて符号化する
    　符号化方法。
  12. 　請求項１０記載の符号化方法であって、
    　前記可変長符号化ステップは、
    　前記入力音響信号の周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記周期性統合包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化し、
    　前記指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記スペクトル包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化する
    　符号化方法。
  13. 　線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成ステップと、
    　前記周期性統合包絡系列から振幅値に依存する可変長符号化パラメータを計算する可変長符号化パラメータ計算ステップと、
    　可変長符号化パラメータを用いて復号する可変長復号ステップ、
    　を実行する復号方法。
  14. 　線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成ステップと、
    　前記周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、可変長符号を復号して周波数領域の系列を得る可変長復号ステップ、
    　を実行する復号方法。
  15. 　請求項１３記載の復号方法であって、
    　入力された指標符号を復号して周期性の程度を示す指標を得る指標復号ステップと、
    　前記スペクトル包絡系列から振幅値に依存する可変長符号化パラメータを計算する第２可変長符号化パラメータ計算ステップ
    　も有し、
    　前記可変長復号ステップは、
    　前記指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記可変長符号化パラメータ計算ステップで計算した可変長符号化パラメータを用いて復号し、
    　前記指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記第２可変長符号化パラメータ計算ステップで計算した可変長符号化パラメータを用いて復号する
    　復号方法。
  16. 　請求項１４記載の復号方法であって、
    　前記可変長復号ステップは、
    　入力された周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、前記可変長符号を復号して周波数領域の系列を得て、
    　前記指標があらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、前記スペクトル包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、前記可変長符号を復号して周波数領域の系列を得る
    　復号方法。
  17. 　請求項１から４のいずれかの符号化装置としてコンピュータを機能させるための符号化プログラム。
  18. 　請求項５から８のいずれかの復号装置としてコンピュータを機能させるための復号プログラム。
  19. 　請求項１から４のいずれかの符号化装置としてコンピュータを機能させるための符号化プログラムを記録したコンピュータが読み取り可能な記録媒体。
  20. 　請求項５から８のいずれかの復号装置としてコンピュータを機能させるための復号プログラムを記録したコンピュータが読み取り可能な記録媒体。
PCT/JP2015/054717 2014-05-01 2015-02-20 符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、記録媒体 WO2015166693A1 (ja)

Priority Applications (19)

Application Number Priority Date Filing Date Title
EP20167433.0A EP3696812B1 (en) 2014-05-01 2015-02-20 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium
US15/304,733 US10199046B2 (en) 2014-05-01 2015-02-20 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium
PL20167433T PL3696812T3 (pl) 2014-05-01 2015-02-20 Koder, dekoder, sposób kodowania, sposób dekodowania, program kodujący, program dekodujący i nośnik rejestrujący
JP2016515878A JP6276845B2 (ja) 2014-05-01 2015-02-20 符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、記録媒体
EP20167431.4A EP3703051B1 (en) 2014-05-01 2015-02-20 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium
KR1020187006872A KR101861781B1 (ko) 2014-05-01 2015-02-20 부호화 장치, 복호 장치, 부호화 방법, 복호 방법, 부호화 프로그램, 복호 프로그램, 기록매체
EP19157754.3A EP3509063B1 (en) 2014-05-01 2015-02-20 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium
PL19157754T PL3509063T3 (pl) 2014-05-01 2015-02-20 Koder, dekoder, sposób kodowania, sposób dekodowania, program kodujący, program dekodujący i nośnik rejestrujący
KR1020167029792A KR101839016B1 (ko) 2014-05-01 2015-02-20 부호화 장치, 복호 장치, 부호화 방법, 복호 방법, 부호화 프로그램, 복호 프로그램, 기록매체
ES15785387T ES2732859T3 (es) 2014-05-01 2015-02-20 Codificador, descodificador, método de codificación, método de descodificación, programa de codificación, programa de descodificación y medio de grabación
KR1020187006878A KR101861787B1 (ko) 2014-05-01 2015-02-20 부호화 장치, 복호 장치, 부호화 방법, 복호 방법, 부호화 프로그램, 복호 프로그램, 기록매체
CN202110083487.4A CN112820305B (zh) 2014-05-01 2015-02-20 编码装置、编码方法、编码程序、记录介质
PL15785387T PL3139380T3 (pl) 2014-05-01 2015-02-20 Koder, dekoder, sposób kodowania, sposób dekodowania, program kodujący, program dekodujący i nośnik rejestrujący
EP15785387.0A EP3139380B1 (en) 2014-05-01 2015-02-20 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium
CN201580022887.7A CN106663437B (zh) 2014-05-01 2015-02-20 编码装置、解码装置、编码方法、解码方法、记录介质
PL20167431T PL3703051T3 (pl) 2014-05-01 2015-02-20 Koder, dekoder, sposób kodowania, sposób dekodowania, program kodujący, program dekodujący i nośnik rejestrujący
US16/199,370 US10607616B2 (en) 2014-05-01 2018-11-26 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium
US16/199,387 US10629214B2 (en) 2014-05-01 2018-11-26 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium
US16/788,539 US11164589B2 (en) 2014-05-01 2020-02-12 Periodic-combined-envelope-sequence generating device, encoder, periodic-combined-envelope-sequence generating method, coding method, and recording medium

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2014094881 2014-05-01
JP2014-094881 2014-05-01

Related Child Applications (3)

Application Number Title Priority Date Filing Date
US15/304,733 A-371-Of-International US10199046B2 (en) 2014-05-01 2015-02-20 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium
US16/199,370 Continuation US10607616B2 (en) 2014-05-01 2018-11-26 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium
US16/199,387 Continuation US10629214B2 (en) 2014-05-01 2018-11-26 Encoder, decoder, coding method, decoding method, coding program, decoding program and recording medium

Publications (1)

Publication Number Publication Date
WO2015166693A1 true WO2015166693A1 (ja) 2015-11-05

Family

ID=54358434

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/054717 WO2015166693A1 (ja) 2014-05-01 2015-02-20 符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、記録媒体

Country Status (9)

Country Link
US (4) US10199046B2 (ja)
EP (4) EP3703051B1 (ja)
JP (4) JP6276845B2 (ja)
KR (3) KR101861781B1 (ja)
CN (3) CN112820304A (ja)
ES (4) ES2790278T3 (ja)
PL (4) PL3139380T3 (ja)
TR (1) TR201909906T4 (ja)
WO (1) WO2015166693A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ES2768090T3 (es) 2014-03-24 2020-06-19 Nippon Telegraph & Telephone Método de codificación, codificador, programa y soporte de registro
KR102529747B1 (ko) * 2016-05-23 2023-05-08 삼성에스디아이 주식회사 무선 충전 장치 및 방법
EP3644515B1 (en) * 2017-06-22 2022-08-24 Nippon Telegraph And Telephone Corporation Encoding device, decoding device, encoding method, decoding method and program
WO2019028883A1 (en) * 2017-08-11 2019-02-14 Lenovo (Beijing) Limited ENCODING POWER RECEIVED FROM REFERENCE SIGNAL

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58168094A (ja) * 1982-03-29 1983-10-04 藤崎　博也 音声分析処理方式
JPS5994795A (ja) * 1982-11-22 1984-05-31 藤崎　博也 音声分析処理方式
WO2012137617A1 (ja) * 2011-04-05 2012-10-11 日本電信電話株式会社 符号化方法、復号方法、符号化装置、復号装置、プログラム、記録媒体

Family Cites Families (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2615862B2 (ja) * 1988-06-14 1997-06-04 日本電気株式会社 音声符号化復号化方法とその装置
KR940002854B1 (ko) * 1991-11-06 1994-04-04 한국전기통신공사 음성 합성시스팀의 음성단편 코딩 및 그의 피치조절 방법과 그의 유성음 합성장치
IT1277194B1 (it) * 1995-06-28 1997-11-05 Alcatel Italia Metodo e relativi apparati di codifica e di decodifica di un segnale vocale campionato
US6345145B1 (en) * 1995-08-25 2002-02-05 Sony Corporation Signal recording/reproducing method and apparatus, signal record medium and signal transmission/reception method and apparatus
JP3259759B2 (ja) * 1996-07-22 2002-02-25 日本電気株式会社 音声信号伝送方法及び音声符号復号化システム
FI964975A (fi) * 1996-12-12 1998-06-13 Nokia Mobile Phones Ltd Menetelmä ja laite puheen koodaamiseksi
EP0993674B1 (en) * 1998-05-11 2006-08-16 Philips Electronics N.V. Pitch detection
US6377915B1 (en) * 1999-03-17 2002-04-23 Yrp Advanced Mobile Communication Systems Research Laboratories Co., Ltd. Speech decoding using mix ratio table
US20030187663A1 (en) * 2002-03-28 2003-10-02 Truman Michael Mead Broadband frequency translation for high frequency regeneration
JP4296753B2 (ja) * 2002-05-20 2009-07-15 ソニー株式会社 音響信号符号化方法及び装置、音響信号復号方法及び装置、並びにプログラム及び記録媒体
US20040002859A1 (en) * 2002-06-26 2004-01-01 Chi-Min Liu Method and architecture of digital conding for transmitting and packing audio signals
CN101552007B (zh) * 2004-03-01 2013-06-05 杜比实验室特许公司 用于对编码音频信道和空间参数进行解码的方法和设备
US8396717B2 (en) * 2005-09-30 2013-03-12 Panasonic Corporation Speech encoding apparatus and speech encoding method
MX2008010836A (es) * 2006-02-24 2008-11-26 France Telecom Un metodo para codificacion binaria de indices de cuantificacion de una envoltura de señal, un metodo para descodificar una envoltura de señal, y modulos de codificacion y descodificacion correspondiente.
US7761900B2 (en) * 2006-08-02 2010-07-20 Clarendon Foundation, Inc. Distribution of content and advertisement
US8239190B2 (en) * 2006-08-22 2012-08-07 Qualcomm Incorporated Time-warping frames of wideband vocoder
US7953595B2 (en) * 2006-10-18 2011-05-31 Polycom, Inc. Dual-transform coding of audio signals
CN101192408A (zh) * 2006-11-24 2008-06-04 华为技术有限公司 选择导谱频率系数矢量量化的方法及装置
KR100868763B1 (ko) * 2006-12-04 2008-11-13 삼성전자주식회사 오디오 신호의 중요 주파수 성분 추출 방법 및 장치와 이를이용한 오디오 신호의 부호화/복호화 방법 및 장치
US8688437B2 (en) * 2006-12-26 2014-04-01 Huawei Technologies Co., Ltd. Packet loss concealment for speech coding
JP5275612B2 (ja) * 2007-07-18 2013-08-28 国立大学法人 和歌山大学 周期信号処理方法、周期信号変換方法および周期信号処理装置ならびに周期信号の分析方法
US8515767B2 (en) * 2007-11-04 2013-08-20 Qualcomm Incorporated Technique for encoding/decoding of codebook indices for quantized MDCT spectrum in scalable speech and audio codecs
EP2214165A3 (en) * 2009-01-30 2010-09-15 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Apparatus, method and computer program for manipulating an audio signal comprising a transient event
US8463599B2 (en) * 2009-02-04 2013-06-11 Motorola Mobility Llc Bandwidth extension method and apparatus for a modified discrete cosine transform audio coder
US8665945B2 (en) * 2009-03-10 2014-03-04 Nippon Telegraph And Telephone Corporation Encoding method, decoding method, encoding device, decoding device, program, and recording medium
JP5223786B2 (ja) * 2009-06-10 2013-06-26 富士通株式会社 音声帯域拡張装置、音声帯域拡張方法及び音声帯域拡張用コンピュータプログラムならびに電話機
US8670990B2 (en) * 2009-08-03 2014-03-11 Broadcom Corporation Dynamic time scale modification for reduced bit rate audio coding
WO2011086923A1 (ja) * 2010-01-14 2011-07-21 パナソニック株式会社 符号化装置、復号装置、スペクトル変動量算出方法及びスペクトル振幅調整方法
CN102194458B (zh) * 2010-03-02 2013-02-27 中兴通讯股份有限公司 频带复制方法、装置及音频解码方法、系统
JP5694751B2 (ja) * 2010-12-13 2015-04-01 日本電信電話株式会社 符号化方法、復号方法、符号化装置、復号装置、プログラム、記録媒体
US9711158B2 (en) * 2011-01-25 2017-07-18 Nippon Telegraph And Telephone Corporation Encoding method, encoder, periodic feature amount determination method, periodic feature amount determination apparatus, program and recording medium
WO2012111512A1 (ja) * 2011-02-16 2012-08-23 日本電信電話株式会社 符号化方法、復号方法、符号化装置、復号装置、プログラム及び記録媒体
US20130006644A1 (en) * 2011-06-30 2013-01-03 Zte Corporation Method and device for spectral band replication, and method and system for audio decoding
ES2657039T3 (es) * 2012-10-01 2018-03-01 Nippon Telegraph And Telephone Corporation Método de codificación, dispositivo de codificación, programa, y medio de grabación

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58168094A (ja) * 1982-03-29 1983-10-04 藤崎　博也 音声分析処理方式
JPS5994795A (ja) * 1982-11-22 1984-05-31 藤崎　博也 音声分析処理方式
WO2012137617A1 (ja) * 2011-04-05 2012-10-11 日本電信電話株式会社 符号化方法、復号方法、符号化装置、復号装置、プログラム、記録媒体

Also Published As

Publication number Publication date
US11164589B2 (en) 2021-11-02
EP3509063A2 (en) 2019-07-10
JP6744471B2 (ja) 2020-08-19
CN106663437A (zh) 2017-05-10
ES2732859T3 (es) 2019-11-26
US10629214B2 (en) 2020-04-21
EP3509063A3 (en) 2019-08-07
ES2790278T3 (es) 2020-10-27
JP6412993B2 (ja) 2018-10-24
CN112820305B (zh) 2023-12-15
US20170040023A1 (en) 2017-02-09
US20190096414A1 (en) 2019-03-28
ES2884626T3 (es) 2021-12-10
EP3696812B1 (en) 2021-06-09
EP3139380A1 (en) 2017-03-08
KR101839016B1 (ko) 2018-03-16
EP3139380A4 (en) 2017-11-01
US10199046B2 (en) 2019-02-05
EP3703051A1 (en) 2020-09-02
EP3139380B1 (en) 2019-04-10
PL3139380T3 (pl) 2019-09-30
PL3703051T3 (pl) 2021-11-22
PL3696812T3 (pl) 2021-09-27
CN112820304A (zh) 2021-05-18
CN106663437B (zh) 2021-02-02
US10607616B2 (en) 2020-03-31
KR20160137615A (ko) 2016-11-30
ES2883848T3 (es) 2021-12-09
US20190096413A1 (en) 2019-03-28
KR20180030243A (ko) 2018-03-21
EP3696812A1 (en) 2020-08-19
KR101861781B1 (ko) 2018-05-28
JPWO2015166693A1 (ja) 2017-04-20
TR201909906T4 (tr) 2019-07-22
JP2018013795A (ja) 2018-01-25
JP2019215587A (ja) 2019-12-19
EP3703051B1 (en) 2021-06-09
US20200202876A1 (en) 2020-06-25
JP2019003219A (ja) 2019-01-10
JP6595684B2 (ja) 2019-10-23
KR101861787B1 (ko) 2018-05-28
KR20180030240A (ko) 2018-03-21
EP3509063B1 (en) 2020-04-22
CN112820305A (zh) 2021-05-18
PL3509063T3 (pl) 2020-08-24
JP6276845B2 (ja) 2018-02-07

Similar Documents

Publication Publication Date Title
JP6412994B2 (ja) 周期性統合包絡系列生成装置、周期性統合包絡系列生成方法、周期性統合包絡系列生成プログラム、記録媒体
JP6595684B2 (ja) 符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、記録媒体

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15785387

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2016515878

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 15304733

Country of ref document: US

ENP Entry into the national phase

Ref document number: 20167029792

Country of ref document: KR

Kind code of ref document: A

REEP Request for entry into the european phase

Ref document number: 2015785387

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2015785387

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE