JPS63311399A - 音合成方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 Ａ、産業上の利用分野 本発明は人間の音声又は楽器の音の合成方法に関する。

Ｂ０発明の概要 本発明は人間の音声又は楽器の音の合成方法において、 音響管の半径を演算し、その演算結果における音響伝搬
の現象が、送電線等の電気回路における系統過渡現象に
近似していることに着目し、音響管の半径を演算し、そ
の演算出力をサージインピーダンス変化に対応させてサ
ージインピーダンスを可変させたことにより、 自然な音質の音の合成を得るようにしたものである。

また、管楽器等の音響管の長さの変化による音の合成は
、インピーダンス結合部に遅延回路を設けて、その長さ
変化に対応した遅延定数を変化させて合成できるように
した音の合成方法である。

Ｃ１従来の技術 音声合成やミュ、−シックシンセサイザー（電子楽器）
等の所謂音を人７エ的に合成して出力する電子装置は、
最近になって１ないし数チップの音声４認識や音声合成
のＬＳＩが音声情報処理と半導体の大規模集積回路技術
により低価格で実現されるようになり、その使用目的、
制約条件により種々の方式が提案されている。この音声
合成には、人間の発生した生の音声を録音しておき、こ
れを適当に結合して文章に編集する録音編集方式と、人
間の声を直接的には利用せず、人間の音声のパラメータ
だけを抽出し、音声合成過程で、そのパラメータを制御
して人工的に音声信号を作り出す方法がある。

このパラメータ方式で良質な合成音が得られることで広
く利用されているパーコール（ＰＡＲＣｏｒｔ）方式が
ある。

音声を電子計算機で扱う場合、音声波形をある周期毎に
サンプリングして各サンプリング点での音声信号の値を
アナログ／ディジタル変換し、その値を０と１の符号で
表示して行われるが、アナログ信号に忠実な記録をする
には、ビット数を増やす必要があるが音声合成信号は大
変多くのメモリーを必要とする。

そこで、この情報量を極力少なくするために各種の高能
率な符号化法が研究開発されている。

その方法の１つとして、１つの音声信号の情報に対し、
最低限１ビツトとした方式で、デルタ変調方式がある。

この方式は、１ビツトの使い方として、次にくる音声信
号値が現在の値より高いか低いかを判定して、高ければ
符号“１”、低ければ符号“０”を与え音声信号の符号
化を行うもので、実際のシステム構成としては一定の振
幅ステップ飛（デルタ）を定めておき、誤差が蓄積され
ないように今までの符号化によって得られる音声の値と
、入力してくる音声信号との残差信号に対して、符号化
を行う。

このような構成は予測コード化といわれ、線形予測法（
何個か前のサンプル値から予測する）およびパーコール
方式（線形予測法の予測係数の代わりにパーコール係数
にといわれる偏自己相関関数を用いる）がある。

Ｄ６発明が解決しようとする問題点 前述のように予測コード化を用いたものは、音と音との
継ぎ目に相当する調音結合が難しいという問題がある。

即ち第１０図は横軸に音声発生の時間ｔをとり、縦軸に
パーコール係数ｋをとったもので、例えば母音から子音
を経て母音に至る発声において、母音の定常から過渡を
経て子音に至りまた母音の過渡を経て母音の定常音に至
る過程で母音と母音の継ぎ目の音が跡切れ、人間が聞い
たときに自然な感じを与えない。

また楽器音合成の場合は、音階の継ぎ目が重要であるが
合成手法が実際の楽器の音発生の原理と異なるため、や
はり自然な感じが無く、特に残響音において顕著にあら
れれる。これら両者において自然な音に近付けるために
は、これを構成するメモリや、演算器等の電子部品を多
く必要とし装置が高価になる等の問題がある。

Ｅ１問題点を解決するための手段 曹 以上の点に鑑み、本発明は人間の音の発生又は楽器の楽
音は人間の口腔又は音響管の長さや断面積等の形状変化
によって作り出される。そこで、これら音響管の音波の
伝達を表す進行波現象を音響管等価回路で解析し、音響
管の半径を演算し、その後円周率と半径の２乗を乗算し
て断面積を求め、その断面積がサージインピーダンスに
反比例することに着目し、サージインピーダンスを変化
させることで音響管の半径を演算して求めた断面積を模
擬的に変化させ、サージインピーダンスを連続的変化す
ることで調音結合をスムーズに行うことができるように
して人間の発声と同様な音の合成を容易となし音声の自
然性の向上を図ったものである。

また楽器音の長さの変化は、遅延回路の段数で模擬し、
断面積変化と相俟ってより自然な楽器音を簡単に実現で
きるようにした。

Ｆ１作用 人間は口腔を動かすことにより、音を発声し、管楽器は
音響管の長さや形状を変化させることによって楽音を作
る。本発明は音響管（人間の声帯から口唇までの声道も
１つの音響管とみなすことができる）の半径ｒを演算し
、その後、πｒ３から断面積を求め、その断面積を等価
回路のサージインピーダンスに１対１に対応させている
ので、このサージインピーダンスを変化させれば音響管
の断面積を変化させたと同じとなる。このサージインピ
ーダンスの変更は、電気技術的に極めて簡単にできるの
で、人間の音の発生と全く同様な音の合成ができ、特に
従来の問題点とされた音と音の継ぎ目゛にあたる調音結
合もサージインピーダンスを連続的に変化することで良
好に行われ、自然に近い音の発声ができる。

また、音響管の長さを変えることは、音波の進行波を遅
らせることであるから電気回路的には遅延回路（メモリ
）の段数を変えることに相当し遅延回路の定数を調整す
ることにより極めて簡単に模擬できる。従って半径を演
算して求めた断面積変化と相俟ってより自然な楽器音も
簡単に実現できる。

Ｇ、実施例 まず、本発明の実施例を述べる前に、本発明の背景につ
いて説明する。

音声を口から外に放射されるには、音源が必要で、この
音源は声帯によって作り出される。一方声帯は一２枚の
ヒダを開閉することによって呼気を断続的に止める働き
があり、その断続によってパフと呼ばれる空気流が発生
し、声帯を緊張させるとこのヒダに張力が加わりヒダの
開閉の周波数が高くなり、周波数の高いパフ音が発生す
る。そして呼気流を大きくすると大きな音となる。

この音源波が声道のような円筒状の音響管を通過すると
、開放端から音波は共振現象によりある成分が強調され
、ある成分が減弱し複雑な母音の波形が作り出される。

音源が同じ波形をもっていても、口唇から放射されるま
でに通過する声道の形によって影響を受ける。即ち、声
道の形状が一定であれば音源のピッチや強度を変えても
スペクトル包絡はあまり変化しない。声道は母音によっ
て極めて複雑な形状を示すが、声道があまり変化しない
部分と大きく変化する部分に分けて考えることができる
。例えば第１図のように長さと断面積がＡｔ、Ａｔとそ
れぞれ異なるような２っ゛の音響管が接続したものと仮
定することができる。−第１図は音響管モデル図、第２
図はその等価回路図で、断面積がＡｔ、Ａｔとそれぞれ
異なる２つの音響管を接続した場合である。′ この音響管の接続する面に着目すると、音波の流れは断
面積の異なる場合、その異なる面で音波の一部が反射す
るという現象を生ずる。この現象は、電気回路でインピ
ーダンスの異なる線路にインパルス電流を流したときの
過渡現象と同じであるら 音声の発生は、前述したように声帯による音源の断続に
よって行われるがこれは電気的には、−インパルスが断
続的に印加されることと等価となる。

音は気体、液体、固体のいずれでも伝わる一種の振動で
あるが、電気回路的には、抵抗の無い無損失のＬＣ分布
回路に対応させることができる。

そしてこ・の等価回路の電気的インピーダンス（Ｖ／Ｉ
）は、ＪＬ／Ｃとなるので、音波の場合に置き換えると
音波の速度、空気密度ρと音速Ｃを掛けたρＣとなり、
音場におけるインピーダンス即ち音響インピーダンスは
気体の質量と音速だけに依存する。

断面積の異なる音響管が連設されていると、その境界面
で反射が起こる。これは電気的なサージインピーダンス
に模擬することができる。即ち、第１図のような音響管
の断面積の異なるブロックの接続された等価回路は第２
図に置き換えられる。

ここで、空気密度をρ、音速をＣとすれば、各音響管の
音響アドミッタンスｙｌ、　Ｙｔは次のように与えられ
る。

但しＺｔ、Ｚｔは音響インピーダンスである。

次に隣接ブロックよりの伝搬電流源を１．、夏ｔとし、
これにより決定される電流分布ａｌｏ　λ、及び１１＋
１！および接合点の電圧をｅとするとまたａ＋＝ｉｔ＋
Ｉｔ、　　ａｔ＝Ｉ＊＋Ｉｔｉｔ＝ａ＋　　　Ｉｔ、　
　　ｉ＊＝ａｔ　　　Ｉとなり、次のステップのための
隣接ブロックの伝搬電流源１１′、Ｉｔ′は、 １１’＝ｌｔ＋ａ＋＋　　Ｉ＊’＝ｉｌ＋ａｔとなる。

上式で１ｎ＝ａ＋　−■ｌおよび１ｔ＝ａｔ　　Ｉｔを
代入してＩｔ′＝２ａｔ−Ｉｔ、Ｉｔ’＝２ａｔ−Ｉｔ
としてもよい。

上記の方式において断面積Ａの時間に対する補間状況を
第５図に示す。この第５図は最も演算の簡単な直線補間
を示している。

第３図は音響管の電気回路等価モデル図で、その（ア）
図は声帯から口唇までの声道を１つの音響管とみなした
音響管モデル図、（イ）図はその電気回路モデル図、（
つ）図は進行波等価モデル図を示している。

第３図を説明するに先立ち、人間の母音はどのようにし
て作られるかを説明する。

第４図は音声発生時の声道の断面積変化を模擬したもの
で、その（ア）図は、「ア」の発声の場合で喉の奥が狭
く口唇が開いた状態で肺から押し出される呼気で声帯が
呼気を断続的に開閉して声道（音響管）の中で反射を繰
り返して出てくる音波が「ア」の音声波形となって出て
くる。「イ」は（イ）図のように喉の方が広く口唇の先
が狭いと「イ」の音声波形が出力される。

このように口の恰好で周波数が決まり、口の恰好を模擬
すれば「ア」なり「イ」が発声される。

口の恰好は音響管の断面積で模擬でき、また音響管の断
面積の変化は、サージアドミッタンスの変化で模擬でき
、サージアドミッタンスの変化は、電軍回路上極めて容
易に可変できる。第３図（ア）は断面積Ａ＋、Ａ１・・
Ａｎと異なる断面積をもった音響管を接続して声道を模
擬したものである。同図（イ）はその音響インピーダン
スを電気回路のＬＣ回路に置き換えたもので、各音響管
を１個のＬＣ線路とし、全体を集中線路のｎ−１の電気
回路としたものである。また第３図（つ）は進行波等価
モデル図で、各音響管の音響インピーダンスｚ、、Ｚｌ
・・・Ｚｎは、音響管の断面積に反比例（音響アドミッ
タンスは比例）し、音波の速度に比例するので となる。なお、同図でＺ。は音源インピーダンス。

ＺＬは放射インピーダンスを示し、またブロック間の矢
印は、進行波と反射波を表している。

今「ア」という音声を発声させる場合は、第４図の口唇
の先の断面積に相当する断面積ＡＩのところで「ア」の
口の恰好を与えて、インパルスＰを断続的に印加するこ
とで、「ア」の音が得られ、また「ア」から「イ」の音
を発声させる場合は、同図（イ）に示すように断面積を
Ａ、′に狭め「イ」の口の恰好を与えることで「イ」が
得られる。

インパルスＰが連続して断続的に与えられ、断面積全体
を「イ」の口の恰好に変化させる場合、声道は第３図に
示すｎ個の音響管によって模擬しているので、これらの
各断面積を「ア」から動かして口の恰好を「アーイ」と
連続的に変えることになる。この音響管の断面積を変え
るということは、サージインピーダンスを徐々に変える
ことによって行われる。

上記のように断面積をＡ１からＡ　１　’に連続的に変
えた直線補間の例は第５図に示しである。この第５図に
おける直線補間の場合でも合成としてはかなり良くなる
けれども、人間の口腔の動きは断面積の直線的な動きに
対応しているのではなく、「アゴ」や「舌」等の上下関
係に対応しているため、合成音としては第５図に示す直
線補間では多少不自然になることがある。

そこで、本発明では音響管の半径ｒを補間対象として演
算し、この演算結果から断面積ＡをＡ＝πｒ１にて計算
し、この計算結果を前記（１）式、（２）式に代入する
ことによって、（■）式。

（２）式の面積比、すなわち係数は次式となる。

なお％ｒｌ＋ｒｌは隣接する音響管の半径である。

上述のように、音響管の半径ｒを演算し、その結果から
断面積を得るようにしたので、定常状態の「ア」、「イ
」の音が得られることは勿論であるが、その中間の音、
即ち、音と音との間の音が自然に近い音質かつ品質で得
ることができる。第６図は上記本発明による音響管の半
径ｒを補間対象として演算した結果の説明図で、図中一
点鎖線で示す曲線は半径ｒを演算したときの断面積の変
化特性曲線である。この曲線から音と音との間の音が自
然に近い音質で得られるようになる。

上記のものは隣接する２つの音響管の場合であるが、隣
接する多分岐管部の演算においては次の一般式を用いて
行う。

ｒ１′ Σ　　ｒ雪、 （但し、ｒｋはに＝１．２・・・ｎ、ｎ分岐管部のそれ
ぞれの音響管の半径、ｒｌはｉ番目分岐管の半径である
） 次に音波の伝搬速度を考えると、これは長さＱでＬＣを
持った電線路にインパルスを印加した時の過渡現象に似
ている。

即ち第７図に示すようにＬＣを有する線路を等価的に表
すと第８図のようになる。ここで両端部からみたサージ
インピーダンスＺ。ｌ＋　Ｚｏｔは、Ｚｏｔ＝ＪＬ／Ｃ
、Ｚｏｔ＝ＪＬ／Ｃとなる。

ここで相手から到達してきた進行波を等価的な電流源と
考えると、 となり電流は中間にｎ個の遅延回路ブロックＺがあれば
、ｎ時間後に出力される。即ち左側の回路で発生したも
のがτ時間後右側に到達したということになる。

となる。但し、ディジタル計算においては、電圧または
電流を細分割するのでＶ、、Ｖ、は計測時刻ｔにおける
電圧、τは経過時間を示している。

第８図では、Ｌ、Ｃ回路にインパルスを印加すれば、１
時間後に出力管側に出る。そしてτ時間前到達されたも
のは相手にも到達しているということを等価的に表して
いる。線路の長さａを１にするということは、遅延ブロ
ックｎを正規化してｌにすることで計算し易くなる。ｅ
を３ｃｍに刻む場合は遅延ブロックのｎを３ブロツクに
すればよい。

第３図（ア）を人間の声道は男性で約１７ｃｍなので、
１ｃｍ刻みで１７本の音響管で模擬すれば、Ａ、から入
った波形は、半周期の電流を１０に分割しそのΔＬを１
０μｓｅｃとすれば、１７０μｓｅｃかかってＡｎ側か
ら出てくる。楽器のトロンポンを考えると、トロンポン
は音響管の長短によって楽音を変える。本発明によれば
、トロンポンの「ア」の音からトロンポンの「二」の音
のパラメータを２つ持てば良いトロンポンの「ア」の音
はトロンポンの「ド」の音からトロンポンの上の音とい
う２つのパラメータがあればよい。その中間音は、遅延
回路の遅延ブロックを変えることによって自由に調音す
ることができる。

即ち第８図における遅延回路Ｚｌ→２．．２．→Ｚ１を
可変することで同じ面積をもつ音響管の長今を変化する
ことに対応させることができる。

次に、第３図（つ）の進行波等価モデルの演算処理を第
９図のフローチャートに基づいて説明する。

音響管ＡＩ（管の半径ｒ＋）にインパルスが入力されろ
と、コンピュータよりなる演算処理装置は、ステップＳ
、にてメモリよりＡＩのａＯ＾、ｉｏ＾。

■い、Ｅを取り出す。取り出した値をもとに、ステップ
Ｓ、では、 ａ　ＯＡ’＝　ｆ　−（Ｂ　、　　Ｉ　ＩＩＡ）ｆ　ｏ
Ａ””　ａ　ＧＡ’　　　［ＧＡの演算を行う。この演
算値ＩＬＯＡＺ　　ＩＯＡ’およびステップＳ３でメモ
リより導入された管Ａｔ（半径ｒｔ）の値ａｌ”＋　ａ
ｌ＾、　　ｌ　ｌＢ＋　　ｌ　ｌＡ＋　　Ｉ　ｌＢ＋　
　Ｉ　ＩＡを用いてステップＳ、では、 ａ　＋ａ’＝　Ｓ　ｔａ−（Ｉ　＋ｅ　＋　Ｉ　ＩＡ）
ａ＋Ａ’　＝ＳＩＡ（１ｌＢ＋　Ｉ　ＩＡ）Ｌ　１．’
＝　ａ　ｔ８’　　１　ｔａｊ　ｌＡ’−ａ　ＩＡ’　
　Ｉ　＋８ Ｉ　＋ａ’−ｉ　ＧＡ’　　ａ　ＧＡ’の演算を行う。

ステップＳ、ではＳ４にて求められた１ＩＢＺ　ａｌａ
′を用いて １０Ａ’　＝　ｉ　ｌＢ＋　ＦＬ　＋１！１を演算する
。また一方、Ｓ、にて求められた値ｉ＋Ａ’、ａ＋Ａ′
と、ステップＳ８においてメモリより導入された管Ａｓ
（半径ｒ３）の値ａ　＊ｓ、λ宏＾。

ｆ　１ｅ、　ｆい、Ｉｔｍ、１いとを用いてステップＳ
７にて次の演算が行われる。

ａ　２ｅ’＝　Ｓ　鵞ａ　（１＊＊＊　Ｉ　ｔＡ）λ＊
Ａ’＝Ｓい（Ｉ　＊＊＊　Ｉ　ｔＡ）ｉ　＊ａ’＝　ａ
　ｔａ’　　ｆ　ｔｓｉ　＊Ａ’＝　ａ　！Ａ’−１ａ
ｈ ｌ　＊ａ’＝　Ｉ　ＩＡ＋　ａｌＡ’ ステップＳ０ではＳ７にて求められた１ｔｌｋ’＋　ｌ
ｔＡ’を用いて １１Ａ’＝　ｉ　ｔａ’＋　ａ　！Ａ’の演算が行われ
る。以下同様にして模擬された音響管の半径ｒを演算し
、この演算結果からπｒ２を計算して断面積が求められ
、ステップ５ｎ−ｔでは、 ａ　ｎＢ’：　ｆ　（Ｉ　ｎｓ） ｌ　ｎＢ’＝　ａ　ｎＢ’　　１　ｎＢＩ　ｎｓ”　ｌ
　ｎ−ＩＡ　＋　ａ　ｎ−ｔＡの演算を行う。その結果
を用いてステップＳＬ１では、 Ｉ　ｎ−ＬＡ’＝　１　、、’＋　ａ　ｎＢ’の演算を
行う。すなわち、音響管の半径ｒ、〜ｒｎを演算し、そ
の演算結果からπｒ、！〜πｒｎ’を計算し、得られた
断面積に対応した等価回路の最終段（ｎ段）における演
算結果の出力がＤ／Ａ変換されて図示省略されたスピー
カに出力され、スピーカより音声として出力される。

したがって演算処理装置は音響管Ａｌ−Ａｌ％に対応し
た演算を行うものであるから、この演算処理装置は音響
管のＡ１〜Ａ、個々の等価回路を流れる各部の電流値お
よび関数ｆ、Ｓ１！１，５ＩＡ（ｉ＝１゜２・・・ｎ−
１）をテーブルとして有しているメモリと、当該等価回
路の各部の電流値を演算する第１の演算手段と、当該等
価回路とは相隣接する等価回路の電流値を用いて当該等
価回路の電流値を演算する第２の演算手段とを備えてい
る。

Ｈ、発明の効果 以上のように本発明は、音響管の半径ｒを演算し、その
演算結果からπｒ３を計算し、得られた断面積が等価回
路のサージインピーダンスにｌ対Ｉに対応していること
に着目し、音響管の半径を演算した結果により求められ
た断面積の変化をサージインピーダンスの変化で模擬す
るようにしたものであるから、パラメータとしては、音
響管の半径を演算するだけを持てばよいので、極めて簡
単にでき、従来のように音声合成のための多くのメモリ
を必要としない。また、調音結合も、サージインピーダ
ンスの連続的な変化によってスムーズに行われ、自然音
に近い発声が得られる。

また、楽音も遅延回路のメモリーの段数を変えることに
より容易に得ることができ、半径を演算した演算値の変
化と相俟ってより自然な楽器音を簡単に実現出来る。

【図面の簡単な説明】

第１図〜第９図は本発明を説明するための図で、第１図
は音響管モデル図、第２図は音響管等価回路図、第３図
は音響管の電気回路等価モデル図、第４図は声道の変化
説明図、第５図は音響管断面積の時間に対する補間説明
図、第６図は本発明による音響管の半径ｒの時間に対す
る補間説明図、第７図は音声伝搬を電気的に模擬した電
気回路図、第８図は第７図の等価回路図、第９図は本発
明の音合成をコンピュータ処理するプログラムの一例を
示すフローチャート図、第１０図は従来のパーコール合
成による調音結合説明図を示す。 Ａ１．Ａｔ・・・Ａｎ・・・音響管、ｒ＋、ｒｔ・・・
ｒｎ・・・音響管の半径、Ｙｌ、Ｙｔ・・・音響管の音
響アドミッタンス、ＣＩ、ＣＩ・・・Ｃ１・・・電気回
路モデルの静電容５１、Ｌ　Ｉｎ　Ｌ　ｔ・・・Ｌｎ・
・・同上のりアクタンス、Ｚｌ＋Ｚｅ・・Ｚ、１・・・
サージインピーダンス。 第４図 菖１ 第５図　　　　第６図

Claims (3)

【特許請求の範囲】
1. （１）音響管を用いて模擬する音合成において、隣合う
   ２つの音響管を伝搬電流源と音響断面積に反比例したサ
   ージインピーダンスの並列回路として扱い、音響管の形
   状変化に伴い音響管の半径を演算し、その後円周率と前
   記演算により求められた半径とを乗算して断面積を求め
   る際に、次式の係数 ｒ＾２＿ｉ／（ｒ＾２＿ｉ＋ｒ＾２＿ｉ＿＋＿１）、ｒ
   ＾２＿ｉ＿＋＿１／（ｒ＾２＿ｉ＋ｒ＾２＿＋＿１）（
   但し、ｒ＿１、ｒ＿ｉ＿＋＿１は隣合う２つの音響管の
   半径、ｉは１、２、３・・・ｎである）を用いることを
   特徴とする音合成方法。
2. （２）音響管を用いて模擬する音合成において、隣合う
   ２つの音響管を伝搬電流源と音響断面積に反比例したサ
   ージインピーダンスの並列回路として扱い、音響管の形
   状変化に伴い音響管の半径を演算し、その後円周率と前
   記演算により求められた半径とを乗算して断面積を求め
   る際に、次式の係数 ｒ＾２＿ｉ／（ｒ＾２＿ｉ＋ｒ＾２＿ｉ＿＋＿１）、ｒ
   ＾２＿ｉ＿＋＿１／（ｒ＾２＿ｉ＋ｒ＾２＿＋＿１）（
   但し、ｒ＿ｉ、ｒ＿ｉ＿＋＿１は隣合う２つの音響管の
   半径、ｉは１、２、３・・・ｎである）を用いると共に
   隣へ新しく伝搬すべき電流源をサージインピーダンスに
   流れ込む電流と、当該音響管に流れ込む電流の和、若し
   くはサージインピーダンスに流れ込む電流の２倍の値と
   当該音響管の電流源との差を用いることを特徴とする音
   合成方法。
3. （３）音響管を用いて模擬する音合成において、隣合う
   ２つの音響管を伝搬電流源と音響断面積に反比例したサ
   ージインピーダンスの並列回路として扱い、音響管の形
   状変化に伴い音響管の半径を演算し、その後円周率と前
   記演算により求められた半径とを乗算して断面積を求め
   る際に、次式の係数 ｒ＾２＿ｉ／（ｒ＾２＿ｉ＋ｒ＾２＿ｉ＿＋＿１）、ｒ
   ＾２＿ｉ＿＋＿１／（ｒ＾２＿ｉ＋ｒ＾２＿ｉ＿＋＿１
   ）（但し、ｒ＿ｉ、ｒ＿ｉ＿＋＿１は隣合う２つの音響
   管の半径、ｉは１、２、３・・・ｎである）を用い、前
   記伝搬電流源の伝搬に遅延回路を持ち、且つこの遅延回
   路を可変とすることで同面積を持つ音響管の長さを変化
   することに対応させたことを特徴とした音合成方法。