JPH0715640B2 - Sound analyzer synthesizer - Google Patents

Sound analyzer synthesizer

Info

Publication number
JPH0715640B2
JPH0715640B2 JP58199816A JP19981683A JPH0715640B2 JP H0715640 B2 JPH0715640 B2 JP H0715640B2 JP 58199816 A JP58199816 A JP 58199816A JP 19981683 A JP19981683 A JP 19981683A JP H0715640 B2 JPH0715640 B2 JP H0715640B2
Authority
JP
Japan
Prior art keywords
time
main
spectrum
frequency
value
Prior art date
Legal status (The legal status 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 status listed.)
Expired - Lifetime
Application number
JP58199816A
Other languages
Japanese (ja)
Other versions
JPS6091227A (en
Inventor
弘志 北川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kawai Musical Instrument Manufacturing Co Ltd
Original Assignee
Kawai Musical Instrument Manufacturing Co Ltd
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
Application filed by Kawai Musical Instrument Manufacturing Co Ltd filed Critical Kawai Musical Instrument Manufacturing Co Ltd
Priority to JP58199816A priority Critical patent/JPH0715640B2/en
Publication of JPS6091227A publication Critical patent/JPS6091227A/en
Publication of JPH0715640B2 publication Critical patent/JPH0715640B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Description

【発明の詳細な説明】 (1)発明の技術分野 本発明は音の分析,抽出,合成を行なうため、分析装
置,主スペクトル抽出装置,合成装置を具え、とくに非
調和性の音に対し高速に解析できる音響解析装置にける
合成装置の改良に関するものである。
DETAILED DESCRIPTION OF THE INVENTION (1) Technical Field of the Invention The present invention comprises an analyzing device, a main spectrum extracting device, and a synthesizing device for analyzing, extracting, and synthesizing sounds, and is particularly fast for anharmonic sounds. The present invention relates to an improvement of a synthesizer in an acoustic analyzer that can analyze the above.

(2)従来技術と問題点 一般に自然音に対しては、解析は余りなされていない
が、音声波形はそのスペクトルエンベロープ(声道の周
波数特性)と音源に分離して考え、情報成分としてのス
ペクトルエンベロープを帯域フイルタ群によつて分析抽
出し、エネルギー成分としての音源をパネル系列とラン
ダム雑音に分離する。この2つの情報すなわちスペクト
ル情報と音源情報を記録しまたは伝送する。合成側では
音源情報から音源をパルス系列とランダム雑音で近似
し、その各周波数成分をスペクトル情報によつて振幅変
調し、それを分析に対応した帯域フイルタ群を通して加
え合せて音声波形を再生する。このようなボコーダ方式
や、ターミナルアナログ,声道アナログ方式等のアナロ
グ方式や、RARCOR等のデジタル方式がすでに提案または
実用されている。
(2) Conventional technology and problems Generally, natural sound is not analyzed so much, but a speech waveform is considered as a spectral envelope (frequency characteristic of vocal tract) and a sound source, and a spectrum as an information component is considered. The envelope is analyzed and extracted by a band filter group, and the sound source as an energy component is separated into a panel series and random noise. These two pieces of information, namely spectrum information and sound source information, are recorded or transmitted. On the synthesis side, the sound source is approximated by a pulse sequence and random noise from the sound source information, each frequency component is amplitude-modulated by the spectrum information, and the result is added through the band filter group corresponding to the analysis to reproduce the voice waveform. Such vocoder methods, analog methods such as terminal analog and vocal tract analog methods, and digital methods such as RARCOR have already been proposed or put into practical use.

しかし、これらは全て音声という調和性の音,すなわち
倍音関係が整数倍となる音に限られており、またそれに
対する解析方法として開発されてきたものである。ここ
で前述の方式をトランペツトやトロンボン等調和性の高
い音の解析に適用した場合には有効であると考えられ
る。しかし、ピアノやチヤイムで代表されるような非調
和性の音に対する解析には適しないし、また非常に困難
なことである。このため、自然音全体に対しとくに非調
和性の音に対して解析可能に方式が望まれている。この
場合、従来の調和性の音と同程度に情報量を減少させる
ことも重要でありかつ高速処理が望まれる。
However, these are all limited to harmonic sounds such as voices, that is, sounds whose harmonic relationship is an integral multiple, and have been developed as an analysis method therefor. Here, it is considered to be effective when the above-mentioned method is applied to the analysis of a highly harmonic sound such as a trumpet or thrombone. However, it is not suitable for analysis of anharmonic sounds represented by pianos and chimes, and it is very difficult. Therefore, there is a demand for a method that can analyze all natural sounds, especially nonharmonic sounds. In this case, it is important to reduce the amount of information to the same degree as that of a conventional harmonic sound, and high-speed processing is desired.

これに対し、本出願人は別提案により、詳しくは後述す
るような音響解析装置を提案した。その概略を述べる
と、要部は分析装置と主スペクトル抽出装置と合成装置
より成り、分析装置では音響信号のデイジタルサンプル
値を所定の時間窓で切取り、これをずらしながら時系列
的にスペクトルを算出し、主スペクトル抽出装置ではこ
の算出された周波数スペクトルからスペクトルエンベロ
ープを求め、これにより最も近似するN個以内の正弦波
成分を表わす主スペクトル成分の周波数値と振幅値と位
相値を時系列的に抽出し、次に合成装置ではこれらの主
スペクトル成分値に基づき、N個分のデジタル正弦波発
生器により周波数,位相,振幅を設定し音響波形を合成
するものである。この提案発明により、非調和性の音に
対しても調和性の音の場合と同様に少ない情報量で高速
処理することができるようになつた。しかし、ここで提
案された合成装置は主スペクトル成分値をN個とした結
果、抽出時間間隔がたとえば10msec毎に合成されること
になり、音色変化が大きい場合には忠実度が悪くノイズ
を発生することが起る。
On the other hand, the present applicant has proposed a sound analysis device, which will be described in detail later, by another proposal. In summary, the main part consists of an analysis device, a main spectrum extraction device, and a synthesis device.In the analysis device, the digital sample value of the acoustic signal is cut out in a predetermined time window, and the spectrum is calculated in time series while shifting it. Then, the main spectrum extracting device obtains a spectrum envelope from the calculated frequency spectrum, and thereby obtains the frequency value, the amplitude value, and the phase value of the main spectrum component representing the most approximate N or less sine wave components in time series. After extraction, the synthesizer sets the frequency, phase and amplitude by N digital sine wave generators based on these main spectral component values to synthesize an acoustic waveform. According to the proposed invention, it is possible to process anharmonic sounds at high speed with a small amount of information as in the case of harmonic sounds. However, as a result of setting the number of main spectrum component values to N, the synthesizer proposed here synthesizes at intervals of, for example, 10 msec, and when the timbre change is large, the fidelity is poor and noise is generated. Things will happen.

(3)発明の目的 本発明の目的は自然音全体とくに非調和性の音に対し分
析,抽出,合成を高速に行なうことのできる音響解析装
置において、主スペクトル成分をN個としたまま大きな
音色変化に対しても忠実度が良くノイズの少ない合成装
置を提供することである。
(3) Object of the Invention An object of the present invention is to provide a sound analysis apparatus capable of analyzing, extracting, and synthesizing a whole natural sound, in particular, anharmonic sound at high speed. It is an object of the present invention to provide a synthesizer having high fidelity against changes and having less noise.

(4)発明の構成 前記目的を達成するため、本発明の音響解析装置の合成
装置は、時間的に変化する楽音信号のディジタルサンプ
ル値を所定の時間窓で切取り、これをずらして1ブロッ
ク宛の周波数スペクトルをフーリェ変換手法により所定
の時間間隔で時系列的に周波数スペクトルを算出する分
析装置と、 「該周波数スペクトルの数より十分少ないN個以内の正
弦波成分を表わし周波数,位相,振幅データより成る主
スペクトル成分値を、該周波数スペクトルのエンベロー
プよりピークレベルの順に抽出し、前記時系列に従い抽
出する主スペクトル抽出装置と、 周波数,位相,振幅を個別に制御可能なN個分のディジ
タル正弦波発生器を有し、前記主スペクトル抽出装置よ
り抽出された主スペクトル成分値の各値に基づき前記時
系列に従い時間的に変化する楽音信号を合成する合成装
置と、 からなり調和又は非調和な楽音を分析,合成する音響解
析装置において、 前記合成装置は、 ある時刻においてN個の周波数,位相,振幅より成る1
組の主スペクトル成分値を、時系列的に各値を複数組記
憶する主データメモリと、 時刻Tにおける1組の主スペクトル成分値を記憶し、前
記合成装置の合成速度に合せて該主スペクトル成分値の
各値を読み出す第1メモリと、 単位時刻前のT−1における1組の主スペクトル成分値
を記憶し、前記合成装置の合成速度に合せて該主スペク
トル成分値の各値を読み出す第2メモリと、 前記主データメモリをアドレス時刻Tと時刻T−1にお
けるN個の主スペクトル成分値を合成に用いる速度より
高速度で順次読みだし第1および第2メモリに転送する
アドレス制御手段と、 時刻T−1から時刻Tまでの分割した時刻を出力する補
間カウンタと、 該出力を受け前記第1および第2メモリより合成速度に
合せて読み出された主スペクトル成分値の各値を時刻T
−1から時刻Tまで直線的に補間するための補間手段
と、を具え、該補間手段により補間された主スペクトル
成分値データを前記正弦波発生器に入力することにより
ノイズの少ない調和又は非調和な時間的に変化する楽音
を合成することを特徴とするものである。
(4) Configuration of the Invention In order to achieve the above object, the synthesizer of the acoustic analysis device of the present invention cuts a digital sample value of a musical tone signal which changes with time in a predetermined time window and shifts the digital sample value to address one block. An analyzer for calculating the frequency spectrum of the frequency spectrum of the above in a time series at a predetermined time interval by the Fourier transform method, and "the frequency, phase and amplitude data representing N or less sine wave components sufficiently smaller than the number of the frequency spectrum. A main spectrum extracting device for extracting the main spectrum component value consisting of the frequency spectrum from the envelope in the order of peak levels, and extracting the main spectrum component values according to the time series, and N digital sine capable of individually controlling frequency, phase, and amplitude. It has a wave generator and follows the time series based on each value of the main spectrum component values extracted by the main spectrum extraction device. A synthesizer for synthesizing a musical sound signal that changes with time, and an acoustic analyzer for analyzing and synthesizing harmonic or inharmonic musical tones, wherein the synthesizer comprises N frequencies, phases and amplitudes at a certain time. 1
A main data memory that stores a plurality of sets of main spectrum component values in time series, and a set of main spectrum component values at time T, and stores the main spectra in accordance with the synthesis speed of the synthesizer. A first memory for reading each value of the component value, and a set of main spectral component values at T-1 before the unit time is stored, and each value of the main spectral component value is read in accordance with the combining speed of the combining device. Address control means for sequentially reading the second memory and the main data memory at a speed higher than the speed used for combining the N main spectrum component values at the address times T and T-1 and transferring them to the first and second memories. And an interpolation counter for outputting divided times from time T-1 to time T, and a main spectrum component value read out from the first and second memories in accordance with the combination speed, the interpolation counter outputting the times. Each value time T
-1 to time T for linearly interpolating means, and by inputting the main spectrum component value data interpolated by the interpolating means to the sine wave generator, harmonic or anharmonic with less noise It is characterized by synthesizing various musical sounds that change with time.

(5)発明の実施例 第1図は本発明の基本となる別提案の音響解析装置全体
の概略説明図である。音響信号は分析装置100に入力
し、ここで入力信号の周波数スペクトル成分が算出され
る。入力信号はデイジタル変換された後一定時間窓のサ
ンプル区間を1ブロツクとして分析される。続いて、数
サンプルずつずらして次々分析を行なう。これにより、
時系列的な周波数スペクトルが得られる。この分析結果
は、いくつもの山(フオルマント)を持つスペクトルエ
ンベロープで表わされる。この山は、分析区間に対して
波形周期が非整数倍であることやフーリエ変換で用いら
れる窓関数による影響で生じるものである。このため、
山のピークに相当するスペクトル成分を抽出し、さらに
有効な数のスペクトラムにデータ圧縮を行なう。
(5) Embodiment of the Invention FIG. 1 is a schematic explanatory view of the whole of another proposed acoustic analysis apparatus which is the basis of the present invention. The acoustic signal is input to the analyzer 100, where the frequency spectrum component of the input signal is calculated. The input signal is digitally converted and then analyzed with one block as a sample section of a constant time window. Subsequently, the samples are shifted by several samples and analyzed one after another. This allows
A time series frequency spectrum is obtained. The result of this analysis is represented by a spectral envelope with many peaks (formants). This peak is caused by the fact that the waveform period is a non-integer multiple with respect to the analysis section and the influence of the window function used in the Fourier transform. For this reason,
The spectral component corresponding to the peak of the mountain is extracted, and the data is compressed into an effective number of spectra.

ここでまとめてこの明細書で用いる主要な波形用語を第
2図(a)〜(c)により説明しておく。同図(a)は
横軸の周波数に対して縦軸に音響信号のデイジタルサン
プル値をとつたいわゆる周波数スペクトルを示す。
The main waveform terms used in this specification will be collectively described here with reference to FIGS. FIG. 5A shows a so-called frequency spectrum in which the vertical axis represents the frequency and the vertical axis represents the digital sample value of the acoustic signal.

単に「スペクトル」というときは個々の周波数スペクト
ルと同意であり、各線スペクトルを指している。これに
対し「周波数スペクトル」というときはこのスペクトル
の集合した状態を示し、この状態中には周波数,振幅,
位相の各情報を含んでいる。次に「スペクトルエンベロ
ープ」は同図(b)に示すように周波数スペクトルの包
絡線を指すものである。さらに、以下に示す「主スペク
トル」は本発明で新用語として用いるもので同図(c)
に示すように、前記スペクトルエンベロープより抽出し
たN個以内の任意の周波数の正弦波成分に対応する線ス
ペクトルを意味し、合成のための周波数,位相,振幅デ
ータの1線分を示すものである。
The term "spectrum" is synonymous with the individual frequency spectrum and refers to each line spectrum. On the other hand, the term "frequency spectrum" refers to the state in which this spectrum is aggregated, and the frequency, amplitude, and
It contains each phase information. Next, the "spectral envelope" refers to the envelope of the frequency spectrum as shown in FIG. Further, the "main spectrum" shown below is used as a new term in the present invention and is shown in FIG.
As shown in FIG. 5, it means a line spectrum corresponding to sine wave components of arbitrary frequencies within N extracted from the spectrum envelope, and shows one line segment of frequency, phase, and amplitude data for synthesis. .

主スペクトル抽出装置200は、分析装置100より出力され
るスペクトルエンベロープを表わすデータのそのデータ
数より相当に少ない数N個の主スペクトル抽出する。N
個の主スペクトルは周波数,振幅,位相データより成
り、合成装置300に転送される。合成装置300は周波数,
振幅,位相を任意に制御できるN個の正弦波発振器を有
し、前述のN個の抽出された主スペクトルデータがN個
の発振器にそれぞれ割当てられ、分析と同様の時間間隔
で音色の変化する波形を時々刻々合成し、次にサウンド
システム400で音響が再生される。
The main spectrum extraction device 200 extracts a number N of main spectra, which is considerably smaller than the number of data of the data representing the spectrum envelope output from the analysis device 100. N
Each main spectrum consists of frequency, amplitude, and phase data, and is transferred to the synthesizer 300. The synthesizer 300 has a frequency
It has N sine wave oscillators whose amplitude and phase can be controlled arbitrarily, and the above-mentioned N extracted main spectrum data are assigned to the N oscillators respectively, and the timbre changes at the same time intervals as in the analysis. The waveforms are synthesized moment by moment, and then the sound is reproduced by the sound system 400.

このように音響解析装置において、分析から合成までの
過程はフーリエ変換した音響信号の周波数スペクトルの
中から、重要でかつ相当に少ない数N個の主スペクトル
を限定的に抽出し、これを合成再生することによりなさ
れる。
As described above, in the acoustic analysis device, in the process from analysis to synthesis, a number N of important and considerably small number of main spectra are limitedly extracted from the frequency spectrum of the Fourier-transformed acoustic signal, and this is synthesized and reproduced. It is done by doing.

たとえば、1回の分析区間を1024サンプルとすると、分
析結果の周波数スペクトラムは最大512個算出される。
For example, if one analysis interval is 1024 samples, 512 frequency spectra at the maximum will be calculated.

上記Xkで示されるスペクトルエンベロープ(周波数スペ
クトラムの全体の包絡曲線)よりN個のたとえば32個の
主スペクトルを抽出する。この時、主スペクトルの値は
補間されたスペクトルエンベロープやその形状から推理
する方法によつて算出した値でもよい。続いて抽出した
周波数,位相,振幅を包む主スペクトルは、32個の独立
した正弦波発振器にそれぞれ割当てられ、 (αq:振幅,ωq:周波数,θq:位相) により示される波形x(t)が合成される。また、32個
一組のαq,ωq,θqは所定の時間間隔たとえば10msec
毎に、分析結果に応じて変化してゆく。
N, for example, 32 main spectra are extracted from the spectrum envelope (envelope curve of the entire frequency spectrum) indicated by X k . At this time, the value of the main spectrum may be a value calculated by a method of inferring from the interpolated spectrum envelope or its shape. Subsequently, the main spectrum including the extracted frequency, phase, and amplitude is assigned to 32 independent sine wave oscillators, respectively. A waveform x (t) represented by (α q : amplitude, ω q : frequency, θ q : phase) is synthesized. In addition, a set of 32 α q , ω q , θ q is a predetermined time interval, for example, 10 msec.
Each time, it changes according to the analysis result.

このようにして、単に32倍音の調和性倍音の合成のみな
らず、非調和性音の合成や有効なデータのみを利用する
ため、たとえば方形波においては偶数倍音を含まないた
め、実質的には32個の合成において、64倍音まで再生可
能となり高品質な合成が行なえる。
In this way, not only the synthesis of the harmonic overtones of the 32nd harmonic, but also the synthesis of the nonharmonic tones and the use of only valid data, for example, since the square wave does not include even harmonics, With 32 composites, up to 64 overtones can be played, and high quality composites can be performed.

第3図は分析装置100の実施例の詳細説明図である。FIG. 3 is a detailed explanatory diagram of an embodiment of the analyzer 100.

音響入力信号はADC101でアナログ−デイジタル変換され
た後一旦波形メモリ102に記憶される。波形メモリ102は
たとえばサンプリング周波数51.2KHzで5秒間記憶でき
るものであれば256Kワードの波形メモリである。記憶さ
れた波形は時間窓の1ブロツクとしてたとえば1024サン
プルずつ高速フーリエ変換(Fast Furier Transform)
するFFT演算回路103に入力し周波数スペクトルが算出さ
れる。このブロツクを順次数ワードずつずらして同様に
演算することにより、時間的に変化するいくつものスペ
クトルエンベロープが算出される。このスペクトルデー
タは主スペクトル抽出装置200に転送される。
The acoustic input signal is analog-digital converted by the ADC 101 and then temporarily stored in the waveform memory 102. The waveform memory 102 is a waveform memory of 256 K words if it can store the sampling frequency of 51.2 KHz for 5 seconds. The stored waveform is, for example, 1024 samples of Fast Fourier Transform as one block of the time window.
The frequency spectrum is calculated by inputting it to the FFT calculation circuit 103. By sequentially shifting this block by several words and performing the same calculation, several spectral envelopes that change with time are calculated. This spectrum data is transferred to the main spectrum extraction device 200.

第4図は主スペクトル抽出装置200の実施例の詳細説明
図である。分析装置100により算出されたスペクトルデ
ータは、各時刻毎に1組のスペクトルエンベロープをス
ペクトルエンベロープメモリ201に転送する。
FIG. 4 is a detailed explanatory diagram of an embodiment of the main spectrum extracting apparatus 200. The spectrum data calculated by the analyzer 100 transfers a set of spectrum envelopes to the spectrum envelope memory 201 at each time.

次にn次ピークレベル検出回路202により、スペクトル
エンベロープメモリ201にあるスペクトルエンベロープ
の山のピーク値を検出し、最大のものから順にN個たと
えば32個のピークレベルが検出される。これにより、そ
の検出される32個の点での振幅値と周波数値と位相値が
それぞれ振幅抽出回路203,周波数抽出回路204,位相抽出
回路205により抽出される。
Next, the peak value of the peak of the spectrum envelope in the spectrum envelope memory 201 is detected by the n-th peak level detection circuit 202, and N peak levels, for example 32 peak levels, are detected in order from the largest peak value. As a result, the amplitude value, the frequency value, and the phase value at the detected 32 points are extracted by the amplitude extraction circuit 203, the frequency extraction circuit 204, and the phase extraction circuit 205, respectively.

ここで32個の主スペクトル成分は、その波形も最も良く
代表する成分が選択されるであろう。次に抽出されたN
個の主スペクトルデータは合成装置300にあるN個たと
えば32個の正弦波発振器に割当てられる。
Here, the 32 main spectral components will be selected such that their waveforms are also best represented. Next extracted N
The main spectrum data are assigned to N, for example, 32 sinusoidal oscillators in the synthesizer 300.

第5図はスペクトルエンベロープより8個の主スペクト
ルを抽出する方法を示す波形図である。多数のポイント
(たとえば512ポイント)で示されるなめらかなスペク
トルエンベロープより山の最大レベルを示す点が〜
にその順位で示される。ここで、抽出周波数範囲をfc
でにすると、7番目のスペクトルは不要となり、かわり
に9番目のf9が抽出される。
FIG. 5 is a waveform diagram showing a method of extracting eight main spectra from the spectrum envelope. The point showing the highest level of the mountain is higher than the smooth spectral envelope shown by many points (for example, 512 points).
Are shown in that order. Here, when the extraction frequency range up to f c, 7 th spectrum is unnecessary, 9 th f 9 are extracted instead.

このように抽出周波数帯域を制限することによつて、小
さなレベルのスペクトルも抽出できるようになり利用可
能性が向上する。このようにレベル順で検出された周波
数に基づいて、位相,振幅も同時に抽出される。取りこ
ぼされたスペクトルは相当にレベルの低いものであり、
全体の音の特徴にはそれ程重要なものではない。このた
め少ないスペクトルで合成しても十分な再現性を得るこ
とができる。
By limiting the extraction frequency band in this manner, it becomes possible to extract a spectrum of a small level and the usability is improved. In this way, the phase and the amplitude are simultaneously extracted based on the frequencies detected in order of level. The missed spectrum is of a much lower level,
It is not so important to the overall sound characteristics. Therefore, sufficient reproducibility can be obtained even if synthesis is performed with a small number of spectra.

第6図は合成装置300の実施例の詳細説明図である。主
スペクトル抽出装置200により各時刻毎に抽出された主
スペクトルデータ(振幅,位相,周波数)は全て抽出デ
ータメモリ301に一旦記憶される。これは各時刻毎の分
析および抽出が実時間で行なわれずその時間間隔には間
に合わず次の処理にはいり込むためである。また、何度
も再生するために、全抽出データは一時記憶しておくこ
とが必要となる。
FIG. 6 is a detailed explanatory diagram of an embodiment of the synthesizer 300. All the main spectrum data (amplitude, phase, frequency) extracted at each time by the main spectrum extraction device 200 is temporarily stored in the extracted data memory 301. This is because the analysis and extraction for each time is not performed in real time, and the time interval cannot be met in time so that the next process can be performed. In addition, it is necessary to temporarily store all the extracted data in order to reproduce it many times.

このため、メモリ容量はたとえば5秒間の音響信号を10
msec間隔で分析し、32個ずつの周波数,位相,振幅の主
スペクトルデータを抽出した場合、(5000/10)×32×
3=48000ワードのメモリとなる。次にアドレス制御回
路302は合成に用いる速度とは無関係な速度で、所定の
時刻の周波数,位相,振幅より成る32個の主スペクトル
データを各時刻毎に周波数メモリ305,位相メモリ304,振
幅メモリ303にそれぞれ転送し、分析における速度と同
じ速度で音響の再合成が可能となるように、一時各デー
タが各メモリにストアされる。次に周波数メモリ305に
は、角速度を決定するN個(以下32個を例にとる)の周
波数データωq(q=1,2,…32)がストアされ、所定の
時間間隔で更新される。このデータは、所定のクロツク
(たとえば分析装置におけるサンプリングクロツク)周
波数でまた32の時分割で常に周波数累算器306により累
算され、時刻tにおける位相角ωqtが算出される。次に
この出力ωqtは、位相加算器307に入力し、もう一方の
入力である位相メモリ304より出力される位相データθq
(q=1,2,…32)と加算され、(ωqt+θq)が算出さ
れる。次に位相加算器307の出力により一周期の正弦波
をデイジタル的に記憶する正弦波テープル308より正弦
波値SIN(ωqt+θq)なる値を読み出し、この出力に振
幅メモリ303より出力される振幅データαq(q=1,2,…
32)と乗算器309で乗算され、αqsin(ωqt+θq)なる
合成スペクトルが出力される。次に乗算器309より出力
される32個の時分割された合成スペクトル波形データを
スペクトル累算器310で累算し、 なるデイジタル合成波形を算出し、さらにDAC311により
デイジタル−アナログ変換された後、サウンドシステム
400から音響出力される。
For this reason, the memory capacity is, for example, 10 seconds for audio signals of 5 seconds.
When analysis is performed at msec intervals and main spectrum data of 32 units of frequency, phase, and amplitude are extracted, (5000/10) × 32 ×
3 = 48000 words of memory. Next, the address control circuit 302 has a frequency memory 305, a phase memory 304, and an amplitude memory at each time of 32 main spectrum data consisting of frequency, phase and amplitude at a predetermined time, at a speed irrelevant to the speed used for synthesis. Each piece of data is transferred to 303 and each piece of data is temporarily stored in each memory so that the sound can be resynthesized at the same speed as in the analysis. Next, the frequency memory 305 stores N pieces of frequency data ω q (q = 1, 2, ... 32) for determining the angular velocity (hereinafter, 32 pieces are taken as an example), and is updated at predetermined time intervals. . This data is always accumulated by the frequency accumulator 306 at a predetermined clock (eg sampling clock in the analyzer) and in 32 time divisions to calculate the phase angle ω q t at time t. Next, this output ω q t is input to the phase adder 307, and the phase data θ q output from the phase memory 304, which is the other input.
(Q = 1, 2, ... 32) is added to calculate (ω q t + θ q ). Next, the sine wave value SIN (ω q t + θ q ) is read from the sine wave table 308 that digitally stores one cycle of the sine wave by the output of the phase adder 307, and is output from the amplitude memory 303 to this output. Amplitude data α q (q = 1,2, ...
32) is multiplied by the multiplier 309, and a combined spectrum of α q sin (ω q t + θ q ) is output. Next, the 32 time-division synthesized spectrum waveform data output from the multiplier 309 is accumulated by the spectrum accumulator 310, After the digital composite waveform is calculated, and digital-analog conversion is performed by the DAC311, the sound system
Sound is output from 400.

このように、合成装置300は32個の主スペクトル成分を
時分割で合成再生可能で、さらに各スペクトル毎に任意
の周波数,位相,振幅を設定できる装置である。
As described above, the synthesizing device 300 is a device capable of synthesizing and reproducing 32 main spectrum components in a time division manner, and further capable of setting an arbitrary frequency, phase and amplitude for each spectrum.

さらに、抽出データメモリ301とアドレス制御回路302は
合成の速度とは無関係な速度でデータ転送が行なわれる
ため、装置303〜311を複数系列でもつことにより容易に
複音化することも可能である。
Further, since the extracted data memory 301 and the address control circuit 302 perform data transfer at a speed that is independent of the speed of synthesis, it is possible to easily produce a polyphonic sound by having a plurality of devices 303 to 311.

以上は別提案の音響解析装置の構成と動作を説明したも
のであり、これにより従来解析や合成が困難であつた非
調和性の音を含む自然音全体の音の解析を比較的少ない
N個ずつの主スペクトル成分値を用い高速処理が可能と
なるものである。
The above is a description of the configuration and operation of the separately proposed acoustic analysis device. As a result, it is possible to analyze a relatively small number of N natural sounds including anharmonic sounds High-speed processing is possible by using each main spectrum component value.

ここで用いた合成装置は、音色変化が前述のように10ms
ec毎に設定されると、たとえば5秒間分のメモリ容量は
4800ワードのメモリが必要となる。この音色変化を高速
にすると、音色の忠実度は高くなり変化によるノイズも
低減される。しかし、その結果大きなメモリ容量が必要
となる。そこで、音色変化をスムーズにしてノイズを低
減する手法として、主スペクトルデータの変化10msecに
対し、さらに分割してその間を補間するような処理を合
成段階で補間回路のみで行なうことを考えたものであ
る。
The synthesizer used here has a timbre change of 10 ms as described above.
If it is set for each ec, for example, the memory capacity for 5 seconds
It requires 4800 words of memory. When this timbre change is made faster, the timbre fidelity is increased and the noise due to the change is also reduced. However, as a result, a large memory capacity is required. Therefore, as a method of smoothing tone color changes and reducing noise, it is considered to perform processing such as subdividing the main spectrum data change of 10 msec and interpolating between them with only the interpolation circuit at the synthesis stage. is there.

すなわち、本発明では上述の処理を実現したもので、主
スペクトルすなわち周波数,位相,振幅データのそれぞ
れに対し主スペクトルの2時刻の変化するデータの両方
を用い、2時刻の間を補間することにより、忠実度が高
くノイズの低減した合成を行なうものである。
That is, in the present invention, the above-described processing is realized. By using both the main spectrum, that is, the data of the main spectrum, that is, the frequency, the phase, and the amplitude data that change at two times of the main spectrum, and interpolating between the two times. , Which has high fidelity and reduces noise.

第7図は本発明の要部である補間回路を含む合成装置30
0′の実施例の構成説明図である。抽出データメモリ301
にはM組の主スペクトルデータ(1組がN個ずつの周波
数,位相,振幅データより構成される)が記憶されてい
る。たとえば0〜5.12秒までを10msec毎に分析した結
果、つまり512組のデータが記憶されている。アドレス
制御回路302は、時刻Tと時刻T−1における2通りの
主スペクトルデータを、振幅メモリ303,位相メモリ304,
周波数メモリ305に10msecで変化しながら転送を行な
う。各メモリには時刻Tの主スペクトルデータStと時刻
T−1の主スペクトルデータSr-1が記憶されている。次
に各メモリの2時刻のデータは各々振幅補間回路313,位
相補間回路314,周波数補間回路315に与えられる。一
方、10msecの時間間隔をさらにm個に分割し、Δtの時
間間隔を発生し、たとえば64個に分割しΔt=0.156mse
cの時間間隔で補間を行なうための時刻をカウントする
補間カウンタ312の出力Kが前記各補間回路313,314,315
に入力している。ここで、各補間回路は主スペクトルデ
ータS=St-1+(St−St-1)・Kまたは=K・St+(1
−K)・St(Kは補間カウンタ312の出力値)を満足す
るような演算回路で構成されている。
FIG. 7 shows a synthesizing device 30 including an interpolation circuit, which is an essential part of the present invention.
It is a structure explanatory view of the example of 0 '. Extracted data memory 301
Stores M sets of main spectrum data (one set consists of N pieces of frequency, phase, and amplitude data). For example, as a result of analyzing 0 to 5.12 seconds every 10 msec, that is, 512 sets of data are stored. The address control circuit 302 converts the two types of main spectrum data at time T and time T-1 into the amplitude memory 303, the phase memory 304,
Transfers to the frequency memory 305 while changing in 10 msec. The main spectrum data S t at time T and the main spectrum data S r-1 at time T -1 are stored in each memory. Next, the two-time data in each memory are given to the amplitude interpolation circuit 313, the phase interpolation circuit 314, and the frequency interpolation circuit 315, respectively. On the other hand, the time interval of 10 msec is further divided into m pieces, and the time interval of Δt is generated. For example, it is divided into 64 pieces and Δt = 0.156 mse.
The output K of the interpolation counter 312 that counts the time for performing the interpolation at the time interval of c is the above-mentioned interpolation circuits 313, 314, 315.
Are typing in. Here, each interpolation circuit has main spectrum data S = S t-1 + (S t −S t-1 ) · K or = K · S t + (1
-K) · S t (K is an operational circuit that satisfies the output value) of the interpolation counter 312.

この補間された主スペクトルデータSは第6図で説明し
た通り なるデイジタル合成波形を算出し音響出力される。第8
図は補間回路を代表して、振幅補間回路313の詳細説明
図である。なお、位相補間回路314と周波数補間回路315
も全く同様の回路で構成されている。
This interpolated main spectrum data S is as described in FIG. Then, a digital synthesized waveform is calculated and output as sound. 8th
The figure is a detailed explanatory diagram of the amplitude interpolation circuit 313 as a representative of the interpolation circuit. The phase interpolation circuit 314 and the frequency interpolation circuit 315
Is also composed of exactly the same circuit.

アドレス制御回路302は、時刻Tの主スペクトルデータ
を転送するためのアドレスを発生する(T)アドレスカ
ウンタ302Bと、時刻T−1の主スペクトルデータを転送
するためのアドレスを発生する(T−1)アドレスカウ
ンタ302Cと、両アドレスを選択するセレクタ302Aによつ
て構成され、抽出データメモリ301をアドレス制御し、1
0msec毎に2時刻の主スペクトルデータを読み出す。次
に2時刻の主スペクトルデータのうち、振幅データは時
刻Tの振幅データAを記憶する振幅メモリ303Aと時刻T
−1の振幅データBを記憶する振幅メモリ303Bより成る
振幅メモリ303に転送される。
The address control circuit 302 generates an address for transferring the main spectrum data at time T (T) and an address counter 302B for generating the main spectrum data at time T-1 (T-1). ) An address counter 302C and a selector 302A that selects both addresses are used to control the address of the extraction data memory 301,
The main spectrum data at two times is read every 0 msec. Next, of the main spectrum data at two times, the amplitude data is the amplitude memory 303A that stores the amplitude data A at time T and the time T.
It is transferred to the amplitude memory 303 including the amplitude memory 303B for storing the amplitude data B of -1.

次に振幅データAと振幅データB,さらに補間カウンタ31
2より出力される2時刻の間をさらに細かくカウントし
た値Cが振幅補間回路313にそれぞれ入力する。振幅補
間回路313は、減算器313AにてA−Bが算出され、次に
乗算器313Bにて(A−B)×Cが算出され、次に加算器
313CにてB+(A−B)×Cが算出される。この振幅補
間回路313は累算−乗算器(たとえばTRW社製TDC1010J)
によつて演算するようにA×C+B×(1−C)を実行
するようなものでも良い。
Next, the amplitude data A and the amplitude data B, and the interpolation counter 31
The value C obtained by further finely counting the two times output from 2 is input to the amplitude interpolation circuit 313. In the amplitude interpolation circuit 313, the subtractor 313A calculates AB, the multiplier 313B calculates (AB) × C, and then the adder.
B + (A−B) × C is calculated at 313C. The amplitude interpolation circuit 313 is an accumulator-multiplier (for example, TDC1010J manufactured by TRW).
It is also possible to execute A × C + B × (1-C) so as to calculate according to

以上のように、主スペクトル成分値のN個はそのままに
して2値間を補間する上記の補間回路により十分きめ細
かい忠実度のよいノイズの少ない合成出力が得られるも
のである。
As described above, the above-mentioned interpolating circuit for interpolating between two values while leaving the N main spectrum component values as it is, a sufficiently finely tuned composite output with little noise can be obtained.

(6)発明の効果 以上説明図したように、本発明によれば、前記提案によ
る分析装置と、主スペクトル抽出装置と、合成装置とよ
り成る音響解析装置において、前記合成装置が提案例に
おける時刻Tと時刻T−1までの間をさらに分割して補
間する手段を設け、該補間手段により補間された主スペ
クトルデータを前記正弦波発生器に入力することによ
り、周波数,位相,振幅のそれぞれが補間されて音響波
形を合成するものである。これにより、提案例の粗い主
スペクトル成分値の場合に比べて細かい音響波形を合成
することができ、情報量を余り増加することなく補間し
た分割数倍の細かい音響波形が得られ忠実度を高めかつ
ノイズを減少することが可能となる。
(6) Effects of the Invention As described above, according to the present invention, in the acoustic analysis device including the analysis device according to the proposal, the main spectrum extraction device, and the synthesis device, the time of the synthesis device is the time in the proposed example. By providing means for further dividing and interpolating between T and time T-1 and inputting the main spectrum data interpolated by the interpolating means to the sine wave generator, each of the frequency, phase, and amplitude can be obtained. It is interpolated to synthesize an acoustic waveform. As a result, it is possible to synthesize a fine acoustic waveform compared to the case of the coarse main spectrum component value of the proposed example, and to obtain a fine acoustic waveform of the number of divisions that is interpolated without increasing the amount of information so much and improve the fidelity. And it is possible to reduce noise.

【図面の簡単な説明】[Brief description of drawings]

第1図は本発明の基本となる提案例の概略説明図、第2
図は波形に関する用語の説明図、第3図,第4図,第6
図はそれぞれ第1図の要部の詳細説明図、第5図は第4
図に関連する動作波形図、第7図,第8図は本考案の要
部の詳細説明図を示し、図中100は分析装置、101はAD
C、102は波形メモリ、103はFFT演算回路、200は主スペ
クトル抽出装置、201はスペクトルエンベロープメモ
リ、202はn次ピークレベル検出回路、203は振幅抽出回
路、204は周波数抽出回路、205は位相抽出回路、300,30
0′は合成装置、301は抽出データメモリ、302はアドレ
ス制御回路、303は振幅メモリ、304は位相メモリ、305
は周波数メモリ、306は周波数累算器、307は位相加算
器、308は正弦波テーブル、309は乗算器、310はスペク
トル累算器、311はDAC、312は補間カウンタ、313は振幅
補間回路、314は位相補間回路、315は周波数補間回路、
400はサウンドシステムを示す。
FIG. 1 is a schematic explanatory view of a proposed example which is the basis of the present invention, and FIG.
The figure is an explanation of terms related to waveforms, Fig. 3, Fig. 4, Fig. 6
FIG. 4 is a detailed explanatory view of the main part of FIG. 1, and FIG.
FIG. 7 and FIG. 8 are operation waveform diagrams related to the figure, and are detailed explanatory views of the main part of the present invention, in which 100 is an analyzer and 101 is an AD.
C and 102 are waveform memory, 103 is FFT operation circuit, 200 is main spectrum extraction device, 201 is spectrum envelope memory, 202 is nth peak level detection circuit, 203 is amplitude extraction circuit, 204 is frequency extraction circuit, and 205 is phase. Extraction circuit, 300,30
0'is a synthesizer, 301 is an extraction data memory, 302 is an address control circuit, 303 is an amplitude memory, 304 is a phase memory, 305
Is a frequency memory, 306 is a frequency accumulator, 307 is a phase adder, 308 is a sine wave table, 309 is a multiplier, 310 is a spectrum accumulator, 311 is a DAC, 312 is an interpolation counter, 313 is an amplitude interpolation circuit, 314 is a phase interpolation circuit, 315 is a frequency interpolation circuit,
400 indicates a sound system.

フロントページの続き (56)参考文献 特開 昭56−51795(JP,A) 特開 昭57−182798(JP,A) 特開 昭58−147798(JP,A) 電子技術Vol.22−13(1980−12) P.22−25 日本音響学会講演論文集(S54.10)3 −2−3P.557−558 日本音響学会講演論文集(S55.5)3 −4−9P.639−636 日本音響学会講演論文集(S55.10)1 −1−1P.359−360 日野「スペクトル解析」(1977−10− 1)朝倉書店P.42−43 マーケル,グレイ著、鈴木訳「音声の線 形予測」(S55−3−25)P.198−232Continuation of front page (56) References JP-A-56-51795 (JP, A) JP-A-57-182798 (JP, A) JP-A-58-147798 (JP, A) Electronic Technology Vol. 22-13 (1980-12) P. 22-25 Proceedings of ASJ (S54.10) 3-2-3P. 557-558 Proceedings of the Acoustical Society of Japan (S55.5) 3-4-9P. 639-636 Proceedings of the Acoustical Society of Japan (S55.10) 1-1-1-1. 359-360 Hino “Spectral Analysis” (1977-10-1) Asakura Shoten P. 42-43 Markel and Gray, Suzuki Translation, "Linear Prediction of Speech" (S55-3-25) p. 198-232

Claims (1)

【特許請求の範囲】[Claims] 【請求項1】時間的に変化する楽音信号のディジタルサ
ンプル値を所定の時間窓で切取り、これをずらして1ブ
ロック宛の周波数スペクトルをフーリエ変換手法により
所定の時間間隔で時系列的に周波数スペクトルを算出す
る分析装置と、 該周波数スペクトルの数より十分少ないN個以内の正弦
波成分を表わし周波数,位相,振幅データより成る主ス
ペクトル成分値を、該周波数スペクトルのエンベロープ
よりピークレベルの順に抽出し、前記時系列に従い抽出
する主スペクトル抽出装置と、 周波数,位相,振幅を個別に制御可能なN個分のディジ
タル正弦波発生器を有し、前記主スペクトル抽出装置よ
り抽出された主スペクトル成分値の各値に基づき前記時
系列に従い時間的に変化する楽音信号を合成する合成装
置と、 からなり調和又は非調和な楽音を分析,合成する音響解
析装置において、 前記合成装置は、 ある時刻においてN個の周波数,位相,振幅より成る1
組の主スペクトル成分値を、時系列的に各値を複数組記
憶する主データメモリと、 時刻Tにおける1組の主スペクトル成分値を記憶し、前
記合成装置の合成速度に合せて該主スペクトル成分値の
各値を読み出す第1メモリと、 単位時刻前のT−1における1組の主スペクトル成分値
を記憶し、前記合成装置の合成速度に合せて該主スペク
トル成分値の各値を読み出す第2メモリと、 前記主データメモリをアドレスし時刻Tと時刻T−1に
おけるN個の主スペクトル成分値を合成に用いる速度よ
り高速度で順次読みだし第1および第2メモリに転送す
るアドレス制御手段と、 時刻T−1から時刻Tまでの分割した時刻を出力する補
間カウンタと、 該出力を受け前記第1および第2メモリより合成速度に
合せて読み出された主スペクトル成分値の各値を時刻T
−1から時刻Tまで直線的に補間するための補間手段
と、 を具え、該補間手段により補間された主スペクトル成分
値データを前記正弦波発生器に入力することによりノイ
ズの少ない調和又は非調和な時間的に変化する楽音を合
成することを特徴とする音響解析装置の合成装置。
1. A digital sample value of a musical tone signal which changes with time is cut out in a predetermined time window, and the frequency spectrum addressed to one block is shifted in a time series at a predetermined time interval by a Fourier transform method. And a main spectrum component value consisting of frequency, phase, and amplitude data representing N or less sine wave components, which is sufficiently smaller than the number of the frequency spectra, and extracted from the envelope of the frequency spectrum in the order of peak level. A main spectrum extracting device for extracting in accordance with the time series, and N digital sine wave generators capable of individually controlling frequency, phase, and amplitude, and main spectrum component values extracted by the main spectrum extracting device. A synthesizer for synthesizing a musical tone signal that changes with time in accordance with the time series based on each value of Analysis harmony musical tones, the acoustic analyzer for synthesis, the synthesis device comprises at a certain time N frequency, phase, than the amplitude 1
A main data memory that stores a plurality of sets of main spectrum component values in time series, and a set of main spectrum component values at time T, and stores the main spectra in accordance with the synthesis speed of the synthesizer. A first memory for reading each value of the component value, and a set of main spectral component values at T-1 before the unit time is stored, and each value of the main spectral component value is read in accordance with the combining speed of the combining device. Address control for sequentially addressing the second memory and the main data memory, and sequentially reading the N main spectrum component values at time T and time T-1 at a speed higher than the speed used for combining and transferring to the first and second memories Means, an interpolation counter for outputting divided times from time T-1 to time T, and a main spectrum component read out from the first and second memories in accordance with the combining speed Time T each value of
An interpolating means for linearly interpolating from -1 to time T, and by inputting the main spectrum component value data interpolated by the interpolating means to the sine wave generator, harmonic or anharmonic with less noise. A synthesizer for an acoustic analysis device, which synthesizes various musical sounds that change over time.
JP58199816A 1983-10-25 1983-10-25 Sound analyzer synthesizer Expired - Lifetime JPH0715640B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP58199816A JPH0715640B2 (en) 1983-10-25 1983-10-25 Sound analyzer synthesizer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP58199816A JPH0715640B2 (en) 1983-10-25 1983-10-25 Sound analyzer synthesizer

Publications (2)

Publication Number Publication Date
JPS6091227A JPS6091227A (en) 1985-05-22
JPH0715640B2 true JPH0715640B2 (en) 1995-02-22

Family

ID=16414114

Family Applications (1)

Application Number Title Priority Date Filing Date
JP58199816A Expired - Lifetime JPH0715640B2 (en) 1983-10-25 1983-10-25 Sound analyzer synthesizer

Country Status (1)

Country Link
JP (1) JPH0715640B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005044867A1 (en) 2003-11-11 2005-05-19 Huntsman Advanced Materials (Switzerland) Gmbh Initiator systems for polymerisable compositions

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002049001A1 (en) * 2000-12-14 2002-06-20 Sony Corporation Information extracting device
JP5001664B2 (en) * 2007-01-25 2012-08-15 株式会社小野測器 Signal converter
JP5057535B1 (en) * 2011-08-31 2012-10-24 国立大学法人電気通信大学 Mixing apparatus, mixing signal processing apparatus, mixing program, and mixing method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5853351B2 (en) * 1979-10-03 1983-11-29 日本電信電話株式会社 speech synthesizer
JPS5919360B2 (en) * 1981-05-06 1984-05-04 日本電信電話株式会社 speech synthesizer
JPS58147798A (en) * 1982-02-27 1983-09-02 松下電工株式会社 Voice synthesizer

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
マーケル,グレイ著、鈴木訳「音声の線形予測」(S55−3−25)P.198−232
日本音響学会講演論文集(S54.10)3−2−3P.557−558
日本音響学会講演論文集(S55.10)1−1−1P.359−360
日本音響学会講演論文集(S55.5)3−4−9P.639−636
日野「スペクトル解析」(1977−10−1)朝倉書店P.42−43
電子技術Vol.22−13(1980−12)P.22−25

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005044867A1 (en) 2003-11-11 2005-05-19 Huntsman Advanced Materials (Switzerland) Gmbh Initiator systems for polymerisable compositions

Also Published As

Publication number Publication date
JPS6091227A (en) 1985-05-22

Similar Documents

Publication Publication Date Title
US11094307B2 (en) Electronic musical instrument and method of causing electronic musical instrument to perform processing
US20060111908A1 (en) Data synthesis apparatus and program
JP2006215204A (en) Voice synthesizer and program
Bonada et al. Sample-based singing voice synthesizer by spectral concatenation
JPH079591B2 (en) Instrument sound analyzer
JP3482685B2 (en) Sound generator for electronic musical instruments
JPH0715640B2 (en) Sound analyzer synthesizer
JPH11259066A (en) Musical acoustic signal separation method, device therefor and program recording medium therefor
US6208969B1 (en) Electronic data processing apparatus and method for sound synthesis using transfer functions of sound samples
GB2294799A (en) Sound generating apparatus having small capacity wave form memories
JP3404850B2 (en) Sound source device
JP2722482B2 (en) Tone generator
JP2705063B2 (en) Music signal generator
JPH05119782A (en) Sound source device
JP2727089B2 (en) Electronic musical instrument sound generator
JP4012410B2 (en) Musical sound generation apparatus and musical sound generation method
Izmirli A hierarchical constant Q transform for partial tracking in musical signals
JPS6088326A (en) Sound analyzer
JP2784399B2 (en) Tone generator
JPS6088327A (en) Main spectrum extractor for sound analyzer
JPH1031496A (en) Musical sound generating device
JP3304709B2 (en) Music signal analysis method
JP3304701B2 (en) Music synthesizer
JPS6097398A (en) Sound analyzer
JP2990897B2 (en) Sound source device