JP3097167B2 - 楽音合成装置 - Google Patents

楽音合成装置

Info

Publication number
JP3097167B2
JP3097167B2 JP03104805A JP10480591A JP3097167B2 JP 3097167 B2 JP3097167 B2 JP 3097167B2 JP 03104805 A JP03104805 A JP 03104805A JP 10480591 A JP10480591 A JP 10480591A JP 3097167 B2 JP3097167 B2 JP 3097167B2
Authority
JP
Japan
Prior art keywords
signal
resonance
pressure
frequency
model
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 - Fee Related
Application number
JP03104805A
Other languages
English (en)
Other versions
JPH04311995A (ja
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.)
Yamaha Corp
Original Assignee
Yamaha Corp
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 Yamaha Corp filed Critical Yamaha Corp
Priority to JP03104805A priority Critical patent/JP3097167B2/ja
Priority to US07/864,632 priority patent/US5272275A/en
Publication of JPH04311995A publication Critical patent/JPH04311995A/ja
Application granted granted Critical
Publication of JP3097167B2 publication Critical patent/JP3097167B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H5/00Instruments in which the tones are generated by means of electronic generators
    • G10H5/007Real-time simulation of G10B, G10C, G10D-type instruments using recursive or non-linear techniques, e.g. waveguide networks, recursive algorithms
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H1/00Details of electrophonic musical instruments
    • G10H1/02Means for controlling the tone frequencies, e.g. attack or decay; Means for producing special musical effects, e.g. vibratos or glissandos
    • G10H1/06Circuits for establishing the harmonic content of tones, or other arrangements for changing the tone colour
    • G10H1/16Circuits for establishing the harmonic content of tones, or other arrangements for changing the tone colour by non-linear elements
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H2250/00Aspects of algorithms or signal processing methods without intrinsic musical character, yet specifically adapted for or used in electrophonic musical processing
    • G10H2250/131Mathematical functions for musical analysis, processing, synthesis or composition
    • G10H2250/155Graham function, i.e. mathematical description of the fluid dynamics of air flowing through a gap, where there is a given pressure differential on either side of the gap, e.g. to model air velocity in wind instruments for physical modeling sound synthesis
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H2250/00Aspects of algorithms or signal processing methods without intrinsic musical character, yet specifically adapted for or used in electrophonic musical processing
    • G10H2250/315Sound category-dependent sound synthesis processes [Gensound] for musical use; Sound category-specific synthesis-controlling parameters or control means therefor
    • G10H2250/461Gensound wind instruments, i.e. generating or synthesising the sound of a wind instrument, controlling specific features of said sound
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H2250/00Aspects of algorithms or signal processing methods without intrinsic musical character, yet specifically adapted for or used in electrophonic musical processing
    • G10H2250/471General musical sound synthesis principles, i.e. sound category-independent synthesis methods
    • G10H2250/511Physical modelling or real-time simulation of the acoustomechanical behaviour of acoustic musical instruments using, e.g. waveguides or looped delay lines
    • G10H2250/515Excitation circuits or excitation algorithms therefor
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H2250/00Aspects of algorithms or signal processing methods without intrinsic musical character, yet specifically adapted for or used in electrophonic musical processing
    • G10H2250/471General musical sound synthesis principles, i.e. sound category-independent synthesis methods
    • G10H2250/511Physical modelling or real-time simulation of the acoustomechanical behaviour of acoustic musical instruments using, e.g. waveguides or looped delay lines
    • G10H2250/535Waveguide or transmission line-based models
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S84/00Music
    • Y10S84/09Filtering
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S84/00Music
    • Y10S84/10Feedback

Landscapes

  • Physics & Mathematics (AREA)
  • Nonlinear Science (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Electrophonic Musical Instruments (AREA)

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、金管楽器のシミュレー
トに用いて好適な楽音合成装置に関する。
【0002】
【従来の技術】自然楽器の発音メカニズムをシミュレー
トすることにより得られたモデルを動作させ、これによ
り、自然楽器の楽音を合成する方法が知られている。特
に、クラリネット等の管楽器の最も基本的なモデルとし
ては、リードの弾性特性をシミュレートした非線形増幅
素子と、共鳴管をシミュレートした双方向伝送回路とを
接続した閉ループ構造のモデルが知られている。このモ
デルでは、非線形増幅素子から信号が出力されると、こ
の信号は後退波信号が加算された後、進行波信号として
双方向伝送回路に入力される。
【0003】次に、この進行波信号は双方向伝送回路の
終端部で反射され、双方向伝送回路を逆方向に伝搬され
る。しかる後に反射波信号は進行波信号に加算され、非
線形増幅素子(励振回路)に帰還される。このように、
非線形増幅素子と双方向伝送回路とからなる閉ループ回
路によって、管楽器における空気圧力波の伝播がシミュ
レートされる。
【0004】また、実際の管楽器には、音高操作用の
孔、いわゆるトーンホールが設けられているものもある
が、このトーンホールをも含めて管楽器をシミュレート
したモデルも知られている。このモデルでは、トーンホ
ールに対応し、各双方向伝送回路間に信号散乱ジャンク
ション(以下、ジャンクションと略す)と呼ばれる信号処
理回路が介挿される。そして、各ジャンクションによ
り、隣接する双方向伝送回路からの各入力信号に対し係
数乗算等の演算処理が行われ、演算結果が隣接する双方
向伝送回路に供給される。この演算処理における乗算係
数等は当該トーンホールの開閉状態に対応して切り換え
られる。
【0005】この場合、非線形増幅素子に帰還される信
号は、各ジャンクションにおいて折り返された成分の総
和となる。しかも、上述したように、各ジャンクション
における演算用の乗算係数は当該トーンホールの開閉状
態に対応して切り換えるので、結局、非線形増幅素子か
ら双方向伝送回路側を見た場合の伝送量周波数特性はト
ーンホールの開閉状態に対応して切り換えられる。この
伝送量周波数特性は、非線形増幅素子の出力信号が開放
状態のトーンホールに対応したジャンクションにおいて
折り返されて非線形増幅素子に帰還されるまでの遅延時
間に対応した周波数(1次)、およびそのほぼ整数倍の各
周波数(高次)に共振周波数を有する多峰性の特性とな
る。なお、この種の技術は、例えば特開昭63−401
99号公報に開示されている。
【0006】
【発明が解決しようとする課題】ところで、上述の技術
においては、主に木管楽器のシミュレーションを行うこ
とを主眼としていたため、金管楽器(リップリード楽
器)のシミュレーションを行うことが困難であった。す
なわち、金管楽器は、楽音の種々のパラメータが演奏者
の唇の状態に大きく支配されるとともに、楽音の立上が
り時に独特の不安定さがあり、さらに、楽音に高調波成
分が多く含まれる等の特徴を有するため、このような特
徴を楽音に反映させなければ不自然な楽音になりがちで
あった。本発明は上述した事情に鑑みてなされたもので
あり、金管楽器の楽音を忠実にシミュレートすることが
できる楽音合成装置を提供することを目的としている。
【0007】
【課題を解決するための手段】上記課題解決するため、
請求項1に記載の発明にあっては、演奏者の息圧に関す
る息圧情報を発生する息圧情報発生手段と、前記演奏者
の唇の押圧力に関する押圧力情報を発生する押圧力情報
発生手段と、前記押圧力が大となると楽音信号の周波数
を高くする周波数制御手段と、前記楽音信号の周波数に
乱数成分を付加する乱数成分付加手段と、前記息圧が大
となると前記乱数成分のレベルを大とする乱数レベル制
御手段とを具備することを特徴としている。また、請求
項2に記載の発明にあっては、演奏者の息圧に関する息
圧情報を発生する息圧情報発生手段と、前記演奏者の唇
の押圧力に関する押圧力情報を発生する押圧力情報発生
手段と、前記息圧および押圧力によって共振特性が変化
する第1の共振手段と、前記共振手段から出力された信
号波を遅延させて前記第1の共振手段に帰還させる遅延
手段とを具備することを特徴としている。また、請求項
3に記載の発明にあっては、さらに楽音周波数近傍の共
振周波数を有する第2の共振手段と、低周波の共振周波
数を有する第3の共振手段とを具備することを特徴とし
ている。
【0008】
【作用】請求項1に記載の構成にあっては、息圧情報発
手段によって発生された息圧情報に基づき演奏者の息
圧が大となると乱数レベル制御手段によって息圧信号に
付加される乱数成分が大となる。従って、実際の金管楽
器において息圧を大とした場合に発生する楽音のピッチ
の揺らぎがシミュレートされる。また、請求項2に記載
の発明にあっては、息圧情報発生手段および押圧力情報
発生手段によって発生された息圧情報および押圧力情報
に基づき第1の共振手段の共振特性が変化する。従っ
て、該息圧および押圧力に応じて楽音信号に変化が付与
される。また、請求項3に記載の発明にあっては、さら
に楽音周波数近傍の共振周波数を有する第2の共振手段
と、低周波の共振周波数を有する第3の共振手段とを具
備するから、第2の共振手段と遅延手段とによって楽音
周波数が決定されるとともに、第3の共振手段によって
楽音の発生当初における不安定な振動状態がシミュレー
トされる。
【0009】
【実施例】以下、図面を参照してこの発明の一実施例に
ついて説明する。A.実施例の前提理論 A.1 金管楽器の物理モデル 本実施例の電子金管楽器は、金管楽器の発音機構をシミ
ュレートするものである。そこで、実際の金管楽器の物
理モデルを図2(イ)、(ロ)を参照し説明する。な
お、図2(イ)、(ロ)に示す物理モデルは、マッキン
タイヤー氏らによる木管楽器のモデル(M.E. Mcintyre,
R.T. Schumacher, J.Woodhouse “On theoscillations
of musical instruments", J.Acoust.Soc.Am.74(5),No
vember 1983 0001-4966/83/111325-21$00.88,1983 Aco
ustic Society of America)に、金管楽器(リップリー
ド楽器)における唇の動作を適用したものである。
【0010】図2(イ)において20は金管楽器の共鳴
管、21は共鳴管21に挿入されたマウスピースであ
る。3は演奏者の唇であり、マウスピース21に押し当
てられている。演奏者がマウスピース21に息を吹き込
むと、唇3の直下部分の圧力が変化し、唇3の非線形特
性により流量fが発生する。この流量fによる圧力変化
は、この瞬間における圧力後退波R1が加算され圧力進
行波F1となって共鳴管20の終端部20eに向かって
伝搬される。
【0011】次に、圧力進行波F1は共鳴管20の各部
で反射されつつ伝搬され、時間の経過とともに進行波F
2、F3に変化する。次に、圧力進行波F3は終端部20
eで反射され反射波R2として唇3に向って伝搬する。
この反射波R2も共鳴管20の各部で反射されつつ伝搬
され、時間の経過とともに反射波R3、R4に変化し、唇
3の直下に帰還される。この反射波R4と、この瞬間に
おける進行波F4との加算結果が、唇3の直下部分にお
ける圧力qとなる。また、口腔内部の空気圧(吹奏圧)
pとマウスピース1内部の空気圧(反射波Rに基づく圧
力)qとの差(以下、圧力差phという)が大きい程、
大きな流入速度が得られる。また、図2(ロ)に同図
(イ)におけるA−A'断面図を示す。図において網目
を付した部分の面積がスリット(開口面積)Sである。
【0012】ここで、グレアム(Graham)の法則によれ
ば、p≧qの場合において、単位面積を単位時間に流れ
る流量(空気速度v)は、次式(A1)で表わされる。v=
√{2(p−q)/ρ}≡L(ph)……(A1)ただし、
ρは空気密度である。また、体積流量fは空気速度vに
スリットSを乗じたものに等しい。ここで、式(A1)の関
係をグラフにしたものを図14に示す。
【0013】A.2 唇3の物理モデル 図2(イ)において、反射波R4が唇3に達すると、こ
れが唇3の振動にフィードバックされ、唇3の定常的発
振に寄与する。そのメカニズムについては、以下説明す
る二つのモデルが代表的である。外側モデル マーチン氏の論文(Daniel W Martin: Lip Vibrations i
n a Cornet Mouthpiece, J.Acoust Soc.Am.,Vol 13,194
2)およびフレッチャー氏の論文(N.H.Fletcher: Exci
tation Mechanisms in Woodwind and Brass Instrument
s, Acustica,Vol 43,1979)等によれば、唇3の振動
は、吹奏圧力によって唇3がマウスピース21の内部で
押し開かれることによって振動が起こると説明されてい
る。これらの説によるならば、唇3の物理モデルは図3
に示すようになる。このモデルにおいては、共鳴管20
の長手方向に移動可能にピストン23が設けられ、この
ピストン23がスプリング24を介して共鳴管20に取
付けられている。そして、このピストン23が唇3に相
当し、図示の位置より左右に移動すると、元の位置の戻
る方向にスプリング24によって付勢される。また、ピ
ストン23と共鳴管20との間隙25が最も小となる位
置の断面積がスリットS(図2(ロ)参照)に相当す
る。
【0014】このモデルによれば、口腔28から印加さ
れる吹奏圧pとマウスピース内部の空気圧qとの差(p
−q)が大となれば大となるほど間隙25が広がる。ま
た、空気圧qがある程度以上大きい場合には間隙25が
完全になくなる一方、吹奏圧pがある程度以上大きい場
合にはピストン23の移動が間隙25の大きさの変化に
ほとんど影響を与えないことも判る。このように、スリ
ットSは圧力差(p−q)の関数S(p−q)として表
すことができ、この関数をスリット関数という。そし
て、外側モデルにおけるスリット関数Sは、図5(イ)
に示すような特性を有する。
【0015】上側モデル 一方、吉川氏の論文(吉川茂 管楽器の発音問題におけ
る問題点,音楽音響研究会資料 1986年3月)によれ
ば、吹奏圧によって生ずる吹奏流圧によって上下唇間に
ベルヌイ圧が生じ、その結果上下唇は閉じる方向に作用
することが発振原理になると説明されている。この説に
よるならば、唇3の物理モデルは図4に示すように構成
される。図においてピストン26はスプリング27を介
して吊り下げられており、共鳴管20内のベルヌイ圧に
よって上下方向に移動する。ここで、ベルヌイ圧とは、
流速によって生ずる圧力であって流速に対して垂直な方
向に発生し、その値は流速が大となるほど大となる。
【0016】なお、図4に示す物理モデルの構築にあっ
ては、以下のような仮定を設けている。(a)ベルヌイ圧
を受ける唇の面積は唇開度にはあまり影響されないこと
とする。(b)唇は上下にある場合も上側だけにある場合も
大差がないこととする(上下の唇がほぼ対称的に動くも
のと仮定する)。上述した仮定(a)、(b)は、上記吉川氏
の論文によるものであるが、実験等により妥当性が証明
されている。次に、ベルヌイ圧PBを数式で示すと下式
(A2)のようになる。PB=−(ρ/2)|2(p−q)
/ρ| ……(A2) 式(A2)によれば、圧力差(p−q)とベルヌイ圧PB
の関係は図5(ロ)に示すようになる。そして、ベルヌ
イ圧PBとスリットSとの関係は、図5(イ)と同様に
なるから、圧力差(p−q)とスリット関数Sとの関係
は図5(ハ)に示すようになる。図5(ハ)によれば、
圧力差(p−q)が大となるほど空気の流速が大とな
り、スリット関数Sが小となることが判る。すなわち、
上側モデルによるならば、外側モデルとは全く逆の結論
が導き出される。
【0017】実際の金管楽器にあっては、上側モデルに
よる作用と外側モデルによる作用とによって発振が行わ
れ、その比率は前者が大であると考えられる。その詳細
を図12を参照してさらに詳細に説明する。図12にお
いて、上下唇は外側モデルに基づいて振動する外側リー
ド部3aと上側モデルに基づいて振動する上側リード部
3bとから成っている。ここで、外側リード部3aは比
較的大きな質量を有するため、大きな時定数を有する。
また、上側リード部3bは外側リード部3aの上に小さ
く付いており、外側リード部3aが押し開かれるように
上下に駆動されると、上側リード部3bもこれに伴って
動く。ここで、外側リード部3aは大きな時定数を有す
ることにより上側リード3bに比較してきわめてゆっく
りと動くものと考えられる。すなわち、振動周波数が3
0ヘルツ前後であると考えられ、楽音のピッチ自体に直
接影響を与えるものではない。従って、上側リード部3
bがピッチ周波数の振動を行い、この振動状態における
唇3の初期スリットSが外側リード部3aによって与え
られるものと考えられる。これにより、演奏状態におけ
る唇3の実際のスリットSは、外側リード部3aの変位
と上側リード部3bの変位との和によって与えられるも
のと考えられる。
【0018】A.3 物理モデル音源の全体構成 上述したような物理モデルに基づいて管楽器の楽音を発
生する物理モデル音源としては、図6に示す構成が一般
的である。図において100は非線形部であり、第2図
における唇3等をシミュレートし、唇3を介して発生す
る圧力進行波を示す信号を、ジャンクション101を介
して線形部102に供給する。線形部102は共鳴管2
0をシミュレートするものであり、第2図において説明
した圧力進行波および圧力反射波を示す信号が伝搬され
る。また、ジャンクション101は、唇3における圧力
進行波と圧力反射波の相互干渉をシミュレートするもの
である。そして、線形部102を伝搬する圧力進行波あ
るいは圧力反射波が楽音信号として出力される。以下、
金管楽器および唇の物理モデルに基づいた非線形部10
0および線形部101の基本構成を順次説明する。
【0019】A.4 非線形部の基本構成 外側モデルに基づく非線形部 まず、外側モデルに基づく非線形部100の基本構成を
図7を参照して説明する。図において103は加算器で
あり、吹奏圧を示す圧力信号pとマウスピース内部の圧
力を示す圧力信号qとが供給されると、両者の差(q−
p)を演算し演算結果を圧力差信号phとして出力す
る。次に、104はAVR(自動電圧調整器)、105
はDCF(直流フィルタ)であり、それぞれ回路の寄生
発振を防止するために設けられたものである。すなわ
ち、AVR104は圧力差信号phの振幅の過剰な変動
を抑制する一方、DCF105は圧力差信号phの高周
波成分を除去する。圧力差信号phはこれらAVR10
4およびDCF105を順次介してDCF106および
グレアム関数テーブル108に供給される。
【0020】グレアム関数テーブル108は、圧力差信
号phに対して式(A1)の演算(グレアム関数)を施し、
その結果を速度信号vとして乗算器109に供給する。
一方、DCF106においては圧力差信号phの直流成
分が抽出され、この直流成分がスリット関数テーブル1
07に供給される。スリット関数テーブル107は、図
5(イ)に示す特性に従って、圧力差信号phに対応す
るスリットSを出力する。次に、速度信号vおよびスリ
ットSは乗算器109において乗算され、乗算結果が流
量信号fとして乗算器110に供給される。
【0021】乗算器110においては、流量信号fと定
数zとが乗算される。ここで、定数zは、第2図(イ)
の物理モデルにおいてマウスピース21および共鳴管2
0の空気流に対する抵抗、すなわち空気流量と空気圧と
の比例定数である。したがって、乗算器110からは、
空気圧を示す信号が出力される。そして、この信号は、
ジャンクション101内部の加算器101aにおいて反
射波信号qiと加算される。なお、反射波信号qiは図2
(イ)のモデルにおいて唇3に到達する圧力反射波
1、R4を示す信号である。この結果、第2図(イ)の
モデルにおける圧力進行波F1、F4に対応する進行波信
号qoがジャンクション101から出力される。
【0022】進行波信号qoは管体線形部102に供給
され、その結果、所定時間経過後に管体線形部102か
ら圧力反射波を示す反射波信号qiが出力される。反射
波信号qiは加算器101aに供給されるとともに乗算
器101cを介して二倍され、加算器101bに供給さ
れる。加算器101bにおいては、二倍された反射波信
号qiと乗算器110から出力された空気圧を示す信号
とが加算され、加算結果がマウスピース内部の圧力を示
す圧力信号qとして出力される。そして、この時点にお
ける圧力信号pから圧力信号qが減算され、減算結果が
新たな圧力差信号phとして出力される。このように、
圧力差信号phがAVR104を介して連続的に出力さ
れ、管楽器における圧力波の伝搬がシミュレートされ
る。
【0023】上側モデルに基づく非線形部 次に、非線形部100を上側モデルで構成した場合のブ
ロック図を図8に示す。図8の構成は大略図7と同様で
あるが、図7におけるDCF106およびスリット関数
テーブル107の代りに、全波整流回路111および共
振回路112が設けられている。ここで、全波整流回路
111は、DCF105を介して供給された圧力差信号
hの絶対値を演算し、演算結果の符号を「−」(マイ
ナス)に変換して出力する。これにより、式(A2)におけ
る絶対値符号の内部の演算が実現される。また、共振回
路112は唇3における周波数特性をシミュレートする
ための回路であり、詳細は後述する。
【0024】混合モデルに基づく非線形部 (i)混合モデルの全体構成 次に、上記外側モデルと上側モデルとを混合して非線形
部100を構成した場合のブロック図を図9に示す。図
9においては、上側モデルを実現するために全波整流回
路111、共振回路112およびスリット関数テーブル
107が設けられており、外側モデルを実現するため
に、共振回路114およびスリット関数テーブル115
が設けられている。
【0025】スリット関数テーブル107から出力され
たスリットSは、乗算器113を介して定数rが乗算さ
れ、乗算結果が加算器117に供給される。また、スリ
ット関数テーブル115から出力されたスリットSは、
乗算器116を介して定数(1ーr)が乗算され、乗算
結果が加算器117に供給される。そして、加算器11
7においては、両乗算器から供給された信号が加算さ
れ、加算結果が乗算器109に供給される。
【0026】これにより、上側モデルによって求められ
るスリットと外側モデルによって求められるスリットと
が定数rによって重み付けされることが判る。なお、上
記以外の構成は、図7および図8に示す構成と同様であ
る。なお、図9に示すブロック図は、例えば図10に示
すように簡略化することも可能である。
【0027】(ii)混合モデルにおける各共振回路の構成 次に、上述した共振回路112,114の特性について
説明する。まず、共振回路112は上側モデルに基づく
唇3の共振状態をシミュレートするための回路であり、
バンドパスフィルタ的な共振特性を有している。すなわ
ち、共振回路112は、その共振周波数Fu付近にピッ
チ周波数Fpが設定され、その振幅上昇比Quはきわめ
て大である。従って、通過帯域以外の信号はほとんど除
去されるものと考えることができる。ここで、共振回路
112を実際の電子楽器に使用する場合を考えると、音
作りのために振幅上昇比Quを適宜変更する場合がある
と考えられる。従って、共振回路112は、振幅上昇比
Quを変更した場合においてもピッチ周波数Fp付近の
ゲインの変化が小となるように構成すると好適である。
また、この回路に付与するパラメータとしては、上述し
た共振周波数Fpおよび振幅上昇比Quで充分である。
【0028】一方、共振回路114は外側リード部3a
(図12参照)をシミュレートするための回路である
が、上述したように、外側リード部3aは上側リード部
3bに比較してきわめてゆっくりと動作するため、共振
回路114の時定数は共振回路112に比較してきわめ
て低い値に設定する必要がある。なお、本発明者の解析
によれば、時定数を30msec程度(この場合、共振
周波数は30Hz程度になる)に設定すると好適であ
る。ここで、共振回路112,114は、例えば図11
に示すような2次のDCフィルタによって構成すること
が可能である。図11において211〜215は加算器
であり、216〜218は乗算器である。乗算器216
には乗算係数として「1/Q」(但しQは振幅上昇比で
ある)が供給され、乗算器217,218には乗算係数
として共振周波数Fが供給される。また、219,22
0は遅延回路であり、それぞれ加算器213,214と
の組合せによって積分回路を構成している。
【0029】ここで、図11に示すフィルタ回路は信号
入力端を適宜変更することにより共振回路112,11
4の何れにも使用可能に構成されている。すなわち、図
示の共振回路は、入力端221から信号が入力されると
共振回路114として動作する一方、入力端222から
信号が入力されると共振回路112として動作する。以
下、共振回路112の伝達関数をHu(ω)、共振回路
114の伝達関数をHo(ω)と表記する。
【0030】(iii)混合モデルの動作 次に、図9による混合モデルの動作を図13を参照して
説明する。まず、図13(イ)に示すように、息圧信号
pが大となると、共振回路114の出力信号は徐々に大
となり、多少のリンギングを伴いながら一定値に収束す
る。これは、図12において、口腔内の圧力が大となる
と、外側リード部3aが徐々に押し広げられることによ
りスリットSが徐々に大となることに相当する。一方、
共振回路112は高い共振周波数を有することにより、
管体線形部102との相互作用によって発振現象を呈す
る。従って、その出力信号波形は、同図(ハ)に示すよ
うに、ほぼ一定振幅の正弦波状になる。両共振回路の出
力信号は乗算器113,116を介して重み付けされた
後、加算器117を介して出力される。この加算器11
7の出力信号は唇3のスリットSに相当する信号であ
り、図13(ニ)に示すような信号波形を有する。
【0031】このように、図9に示す混合モデルによれ
ば、共振回路112および114を組合せたことによ
り、息圧信号pが付与された当初において加算器117
の出力信号波形が乱れる。ここで、実際の金管楽器にお
いても、息を吹込む際(アタック時)に不安定な挙動を
示すことが知られている。従って、図9の混合モデルに
よれば、この金管楽器特有のアタック時の不安定さがシ
ミュレートされることが判る。
【0032】ローカルフィードバック排除型モデル
(i)ローカルフィードバック系に係る問題 上記〜において説明した各モデルは、流量信号fが
空気流に対する抵抗zと乗算され、乗算結果が直ちに加
算器101bを介して加算器103にフィードバックさ
れる。従って、このフィードバックループ(ローカルフ
ィードバックループ)においてリミットサイクル(高い
周波数での異常発振)が発生するという危険が有り、こ
れを防ぐために回路中にDCF105、AVR104等
を挿入する必要が有った。しかし、DCF105やAV
R104によってもリミットサイクルを充分に抑制する
ことは困難な場合があり、むしろ音色の幅を狭くすると
いう問題を生ずることもある。従って、非線形部からロ
ーカルフィードバックループを除去することができれ
ば、リミットサイクルを抑制しつつDCF105やAV
R104を取除くことができる。
【0033】(ii)ローカルフィードバックループを除去
するための方策 まず、ローカルフィードバックループを除去するために
は、ローカルフィードバックループを有する系の伝達関
数を解析し、解析された伝達関数と等価であるとともに
しローカルフィードバックの無い系を採用すればよい。
その原理を説明するため、図15にローカルフィードバ
ックループを有する系(図8〜図10)を簡略化したモ
デルを示す。このモデルにおいては、簡単のため、乗算
器109に供給されるスリット信号Sは一定であること
としている。ここで、スリット信号Sを一定とした場合
において、乗算器110の出力信号xに対する乗算器1
01cの出力信号yの応答関係を解析的に解くことによ
り、図16に示すように、図15と等価であってローカ
ルフィードバックループの無い系が構成される。また、
図17に、図15のグレアム関数テーブル108の特性
(曲線A)と、図16のグレアム関数テーブル108a
(曲線B)のグラフを示す。
【0034】ここで、図15、16は等価であるが、そ
の等価性はある一定のスリット信号Sについてのみ成立
するものである。従って、スリット信号Sが種々の値を
取った場合における信号x,yの応答特性を求め、これ
に基づいてテーブル108aにおけるグレアム関数を適
宜修正する必要がある。まず、ローカルフィードバック
ループが存在する場合において、スリット信号Sが変化
した場合のグレアム関数の変化を図18に示す。図18
によれば、スリット関数が大となるに従って、グレアム
関数は一定のゲイン特性に収束するような飽和特性を有
することが判る。
【0035】次に、ローカルフィードバックが無い場合
において、スリット信号Sが変化した場合のグレアム関
数の変化を図19に示す。図19によれば、スリット信
号Sに対してグレアム関数が単調な比例特性を有してい
ることが判る。ここで、金管楽器、特にトランペット等
の楽音は高次倍音を多く含む鋭い音に特徴が有り、これ
をシミュレートするためには図19に示すグレアム関数
の変化を実現することが不可欠である。しかしながら、
上述した方策によって単にローカルフィードバックルー
プを除去すると、オープンループ系(管体線形部102
までも含めたループ系)に対して大きな利得を与える必
要が生じ、これによってオープンループ系において異常
発振を生ずる場合がある。本発明者の解析によれば、こ
の異常発振は、圧力差信号phが「0」近傍付近にある
場合にグレアム関数の勾配がきわめて大となること(図
19参照)に起因するものであることが判明した。従っ
て、ローカルフィードバックループを除去する場合に
は、何らかの方策によって圧力差信号phが「0」近傍
にある場合のグレアム関数の勾配を抑制とすると一層好
適である。
【0036】(ii)グレアム関数の勾配を抑制するための
方策 グレアム関数の勾配を抑制するための一つの方策として
は、図20に示すような回路構成を採用することが考え
られる。図20において124は関数テーブルであり、
テーブル108a(図16参照)に設定されたグレアム
関数を圧力差信号phで除算して成る関数が設定されて
いる。また、126は乗算器であり、関数テーブル12
4の出力信号とスリット関数Sとを乗算して出力する。
125はリミット回路であり、乗算器126の出力信号
が所定値以上になると、この信号を該所定値に抑制して
出力する回路であり、図21(イ)に示すような入出力
特性を有している。次に、128は乗算器であり、リミ
ット回路125の出力信号と、乗算器127を介して供
給された圧力差信号phとを乗算して出力する。
【0037】上記構成において、リミット回路125が
リミット動作を行わない場合には、図20に示す回路は
図16に示す回路と等価である。これは、図20におけ
る関数テーブル124の出力信号は図16におけるグレ
アム関数テーブル108aの出力信号を圧力差信号ph
で除算した値に等しいが、この信号は乗算器128を介
して圧力差信号phと乗算されるからである。一方、リ
ミット回路125がリミット動作を行った場合には、グ
レアム関数が図22に示す直線Aによって抑制された場
合に等しい入出力関係を有する。これにより、グレアム
関数の勾配を抑制しつつ、図19に示された変化に近似
したグレアム関数の変化を実現することができる。ここ
で、図21(イ)のリミット関数を用いた場合のグレア
ム関数の変化を図22に示す。なお、リミット回路12
5の入出力特性は、例えば図21(ロ)に示すようなな
めらかな特性であってもよい。この場合、グレアム関数
の勾配は漸近的に抑制されるから、一層自然な楽音を発
生させることができるものと考えられる。
【0038】A.5 管体線形部の構成 次に、管体線形部102の各種モデルの構成を説明す
る。K−Lラティスモデルに基づく管体線形部 管体線形部を実現するモデルの一例として、図23に示
すようなK−Lラティスモデルが知られている。ここ
で、金管楽器の共鳴管20(図2参照)は一般的には円
錐形であり、K−Lラティスモデルにおいては、この円
錐形をm段(mは「2」以上の自然数)の多段円筒状に
近似することによって共鳴管20の圧力波の伝搬をシミ
ュレートしている。
【0039】図において81−0〜81−mは遅延回路
であり、多段円筒状の管体における各段の振動の伝搬遅
延をシミュレートするものである。次に、83は反転回
路であり、管体の終端部における圧力波の反射をシミュ
レートする。この反射の際の損失はローパスフィルタ8
0によってシミュレートされる。次に、82−0〜82
−mは管体各部における圧力波の反射をシミュレートす
るジャンクションである。上記各構成要素は全体として
ループ状に接続され、図2図(イ)のモデルにおける本
体部20をシミュレートする。なお、このようなK−L
ラティスモデルに基づく管体線形部は、木管楽器をシミ
ュレートするために使用される周知の管体線形部と同様
であり、例えば、特願平1−101308、特願平1−
259735、特願平1−258229、特願平2−2
80196あるいは特願平2−2800197等におい
て、本出願人によって種々の回路が開示されている。
【0040】円筒近似型モデルに基づく管体線形部 また、円錐管たる共鳴管20は、二つの円筒管を並列に
接続することによっても近似することも可能である。ま
ず、図24に示すように、共鳴管たる円錐ホーン30と
点音源29とが存在する場合、円錐ホーン30の入力イ
ンピーダンスZ1は、下式(A3)で示される。
【数1】 式(A3)においてkは点音源の周波数、r1は点音源29
から円錐管30までの距離、hは円錐管30の長さ、r
2はr1およびhの和、ρは音波伝搬媒体たる空気の密
度、cは音速である。また、Z2、Z3は各々近似された
円筒の特性インピーダンスである。
【0041】式(A3)から円錐管30は二つの円筒管を並
列に接続したものと等価であることが判るが、このモデ
ルを用いて構成された管体線形部のモデルを図25に示
す。図においてフィルタ130および遅延回路131等
によって一方の円筒管がシミュレートされる。すなわ
ち、該円筒管における圧力波の伝搬遅延が遅延回路13
1によってシミュレートされ、該円筒管内部における損
失がフィルタ130によってシミュレートされる。
【0042】また、フィルタ130および遅延回路13
1は、図示のように、乗算器134,135および加算
器136を介して加算器140に接続されている。ま
た、同様に、フィルタ132および遅延回路133等に
よって他方の円筒管がシミュレートされる。そして、加
算器140によって、両円筒管の並列に接続されること
がシミュレートされる。
【0043】ここで、エネルギー保存の法則に鑑みれ
ば、乗算器142,135,138の乗算係数b1,b
2,b3は「−1」になるものと考えられ、乗算器14
1,134,137の乗算係数f1,f2,f3は
「0.5」になるものと考えられる。しかし、これらの
値以外にも系として正常に動作する組合せは多数存在す
るから、これらの乗算係数として種々の値を設定するこ
とによって音色の幅を持たせることも可能である。
【0044】管体線形部の設計方法 次に、上述したモデルに基づいて、実際に管体線形部を
設計する場合の好適な手順を説明する。(i)基本的な設
計方法 まず、実際の金管楽器にあっては、任意のピッチを実現
するために管長を調節しつつ唇の緊張のさせ具合によっ
てモードを選ぶ。ここで、実際の金管楽器において管長
が可変である部分は直管になっているから、管長を調節
することは、物理モデルにおいては、例えばシフトレジ
スタを可変遅延回路として使用しすることによってシミ
ュレートすることができる。一方、管長が可変である部
分の前後の部分、すなわちマウスピース部やベル部にあ
っては、楽器の種類毎に様々な形状を有している。この
部分をシミュレートしてヴォイシングするためには、例
えば図41に示すように、マウスピース部やベル部の形
状をディスプレイに表示し、表示された形状に近似した
多段円筒管によってシミュレートするとよい。
【0045】上述した点に鑑みて、K−Lラティスモデ
ルと円筒近似型モデルとを組合せて成る管体線形部のモ
デルを図31に示す。図において遅延回路40、42の
遅延時間は上記各円筒の長さに応じて設定する。また、
積分回路44、遅延回路45,46は可変管長部分をシ
ミュレートしており、管長に応じた遅延時間DL1が与え
られる。また、遅延回路40,42によって、マウスピ
ース部を近似的に構成する各円筒における遅延時間がシ
ミュレートされる。さらに、ベル部においては、円錐型
で近似したモデルが使用されており、フレアの具合から
遅延回路53,58における遅延時間DLL,DLSが設定
される。また、ベル部においては、管形状を変更した場
合に、K−Lラティスモデルの各格子の係数、遅延時間
が影響されるのみならず、放射特性も影響を受ける。す
なわち、ベルを大きくすることにより低い周波数まで放
射されるようになるから、ローパスフィルタ54,59
のカットオフ周波数も下げる必要がある。なお、金管楽
器の物理モデルは、例えば図32に示すように、円筒近
似型モデルを使用せずK−Lラティスモデルのみ使用し
ても構成し得ることは言うまでもない。
【0046】(ii)PLLを用いたオートチューニング方
ところで、上述した設計方法によれば、共振周波数Fu
に関してはキーコードから決定できるが、遅延時間DL1
に対応する管長については実際にウエーブガイドを用い
て管の近似を行った後でなければ困難である。また、音
作りの段階においては、管体線形部あるいは非線形部に
係る種々のパラメータを変更する必要が生じるが、パラ
メータを変更した場合に音程を合せることはきわめて煩
雑であった。そこで、本発明者は、実際に管を含めた系
を発振させ、その発振周波数が所望のピッチになるよう
に管長すなわち遅延時間DL1を調整し得る方法を開発し
た。以下、その手順を説明する。
【0047】まず、所望のピッチ周波数Fpに対して、
共振周波数FuはFu=Ku・Fp ……(A4) なる式で与えられる(ただしKuは定数であり、通常
「1」より僅かに大きな値になる)。なお、リップコン
トローラ情報が無いと決らないパラメータも多いが、そ
のような値に対しては、予めデフォルト値を決めてお
く。また、ランダムに変動するパラメータに対しては、
その変動幅に応じた乱数を付与しておく。次に、図31
に示されたモデルの総遅延時間DLTを求め、これに基づ
いて遅延時間DL1を求める関係式は以下の通りとなる。
LT=DL1+DL2+DL3+DLS+DLL ∴DL1=DLT−(DL2+DL3+DLS+DLL) ……
(A5)
【0048】次に、最低音(モードNO.=2,ピストンは
全開)からチューニングをとっていくとすると、最低音
での遅延時間DLTの近似値(DLT,max)は次のように
して求まる。まず、最低音(Fp,min)は第2モード
であるから、(DLT,max)=2/(Fp,min×T)
……(A6) となる。この値を式(A5)に代入すると、最低音における
遅延時間たる(DL1,max)が求められる。求められた
(DL1,max)に対して、(Fp,min)を目標周波数と
してPLLをかけると、このピッチ周波数での遅延時間
L1の最適値が求められる。次の目標ピッチが上記最低
音に対して半音高い音であるとすると、ピッチは約1.
06倍になる。すなわち、総管長DLTは約1/1.06
倍されることになる。従って、先に求められた最低音で
の遅延時間DL1に基づいて該最低音における遅延時間D
LTを式(A5)から求め、これを1/1.06倍して再度式
(A5)から遅延時間DL1を求め、求めた遅延時間DL1に対
して再度PLLをかけるとよい。以下同様に、順次高い
ピッチに対応する遅延時間DL1を求めてゆくことによ
り、各ピッチに対して精度の高い初期値が得られるか
ら、PLLをかけることによって最適な遅延時間DL1
設定することが可能である。
【0049】このように、順次ピッチを上げながらチュ
ーニングを進めてゆくと、当該モードでの再高音に達す
る。すなわち、ピストンの切換ではなくモードの切換に
よってピッチを上げる必要が生じる。モードを切換える
と、管長は元の長さに戻ることとなるが、この場合の遅
延時間DL1は先に求めた最低音における遅延時間DL1
ほぼ同値になるものと考えられる。従って、最低音にお
ける遅延時間DL1をモードが切換えられた後のDL1の初
期値に設定し、この初期値に対してPLLをかけること
によって同様に最適な遅延時間DL1を求めることができ
る。
【0050】A.6 エンベロープ発生方法の改善 本実施例のようないわゆる物理モデル音源においては、
楽音のエンベロープの発生要因は、外部から与えられる
各種パラメータに依るところが大きい。すなわち、外側
リード部における各種パラメータもエンベロープの発生
に寄与するが、エンベロープ特性の大半は制御ソフトウ
エアによって与えられるエンベロープ信号あるいは演奏
者によって与えられる息圧p等によって決定される。と
ころで、本発明者の観察によれば、実際の金管楽器の演
奏時においては、息圧pに不感帯やヒステリシス特性が
存在することが判明した。すなわち、音の立上がり時に
おいては息圧pが所定値pthになるまで音が発生しない
一方、音の立下がり時においては息圧pが所定値ptl
なるまで音が消えないことが判明した。なお、所定値p
thは所定値ptlよりも常に大である。
【0051】これら所定値pth,ptlは種々のパラメー
タや音域によって変化するが、例えば図42(ロ)に示
すように、キーコードKCに対応して両者の概略値を予
め求めておき、この概略値に基づいて楽音を制御する
と、実際の金管楽器の楽音を一層忠実に再現することが
できる。すなわち、図42(イ)に示すように、エンベ
ロープ信号が所定値pthを越えるまでは楽音を発生しな
い一方、エンベロープ信号が所定値ptl以下になると楽
音の発生を停止するようにすればよい。
【0052】B.一実施例の構成 次に、本実施例の構成を図面を参照し説明する。全体
構成 まず、本実施例の全体構成を図1に示す。図において1
は金管楽器を模したリップコントローラであり、演奏者
によって入力された演奏情報をMIDI信号に変換して
出力する。6はMIDIインターフェースであり、リッ
プコントローラ1から入力されたMIDI信号をバス1
3を介して各部に供給する。 次に、2は中央処理装置
(CPU)であり、読出し専用メモリ(ROM)4に設
定された処理プログラムに基づいて、他の構成要素を制
御するように構成されている。また、ROM4には、制
御プログラムの他に、処理において使用される各種のデ
ータおよびテーブル等も格納されている。
【0053】次に、8はスイッチ回路であり、音色等の
制御のために種々のスイッチ、操作子等が設けられてい
る。次に、7は表示装置であり、CPU2の制御に基づ
いて各種の情報が表示可能に構成されている。次に、9
はDSP(デジタル・シグナル・プロセッサ)であり、
CPU2からバス13を介して供給された演奏情報に基
づいて楽音信号を発生する。発生された楽音信号はDA
C(デジタル・オーディオ・コンバータ)を介してアナ
ログ信号に変換され、サウンドシステム11を介して発
音される。
【0054】次に、5はRAM(読出し/書込み記憶装
置)であり、CPU2の動作に必要なワークメモリとし
て設けられている。なお、RAM5はデータの消失を防
ぐためにバッテリ・バックアップ型の物を使用すると好
適である。また、12は不揮発性RAMから構成された
RAMカートリッジであり、バス13に着脱可能に設け
られるとともに、RAM5と同様にCPU2によって読
出し/書込み自在になっている。RAMカートリッジ1
2は音色データの格納用に使用すると好適である。
【0055】リップコントローラ1の構成 次に、リップコントローラ1の構成を図26を参照して
説明する。図においてリップコントローラ1は細長筒状
に構成されており、その一端に吹込み口31が設けられ
ている。また、吹込み口31にはアンブシュアセンサ3
2が設けられている。、このアンブシュアセンサ32は
唇3(図示せず)の押圧力を検出し、検出結果をアンブ
シュア信号embとして出力する。
【0056】次に、33は圧力センサであり、吹込み口
31を介して吹込まれた演奏者の息圧を検出して出力す
る。また、34、35および36はそれぞれ第1ピスト
ン、第2ピストンおよび第3ピストンであり、実際の金
管楽器における第1〜第3ピストンと同様の作用を行
う。すなわち、第1ピストン34が押下されると楽音の
音程が全音下げられ、第2ピストン35が押下されると
音程が半音下げられ、第3ピストン36が押下される
と、音程が1音半だけ下げられる。
【0057】次に、37、38および39はさらに大き
く音程を下げるために設けられたピストンであり、ピス
トン37が押下されると音程が半オクターブ下げられ、
ピストン38が押下されると音程が1オクターブ下げら
れ、ピストン39が押下されると、音程が2オクターブ
下げられる。ここで、ピストン37〜39の組合せによ
って設定できるオクターブの変位を下表1に示す。
【0058】
【表1】
【0059】表1において最上欄はピストン3桁の数字
が記載されているが、これはピストン37〜39のオン
/オフ状態を示しており、“0”がオフ状態、“1”が
オン状態に対応している。そして、表1の最下欄には各
状態におけるモード番号が示されている。なお、実際の
金管楽器にあっては、モード番号は唇の緊張度で指定で
きるのであるが、唇の緊張度を測定することは若干困難
である一方、上述のようにピストン37〜39で指定し
ても楽音のリアリティーは損われないものと考えられる
ので、本実施例にあっては上述のようにモード番号を指
定することとした。
【0060】ROM4のソフトウエア構成 次に、ROM4に記憶されている制御プログラムの全体
構成を図29を参照して説明する。図において60はタ
スクスケジューラであり、所定周期毎にタイマ割込がか
けられることによって駆動され、他のプログラムを時分
割で実行する。次に、61はMIDIデバイスドライバ
であり、MIDIインターフェース6(図1参照)を介
して受信したMIDI信号を内部で使用される演奏情報
に変換して他のプログラムに供給する。また、62はL
CE(LCエミュレータプログラム)であり、リップコ
ントローラ1の代りに鍵盤が演奏情報入力装置として使
用された場合に、この鍵盤の出力するMIDI信号を、
リップコントローラ1の出力するMIDI信号と等価な
信号に変換するプログラムである。次に、63はALC
(アルゴリズム制御プログラム)であり、MIDIデバ
イスドライバ61あるいはLCE62から供給された演
奏情報に基づいてキーオン信号KONを出力するととも
に、演奏情報に基づいて楽音制御のための各種パラメー
タを求め、これらパラメータをDSP9に供給する。な
お、上記各プログラムの詳細は動作とともに詳述する。
【0061】DSP9の構成 次に、DSP9の構成を説明する。まず、DSP9は高
速のCPUおよびこのCPUの制御プログラムを格納し
たROMおよびワーク用のRAM等から構成されている
(図示せず)。そして、上記ROMに格納された制御プ
ログラムは図6の物理モデル音源における非線形部10
0、ジャンクション101および管体線形部102をそ
れぞれ実現するプログラムを有している。
【0062】ここで、管体線形部102を実現するプロ
グラムのアルゴリズムの構成は図31に示す管体線形部
のモデルと等価であり、非線形部100およびジャンク
ション101を実現するプログラムのアルゴリズムは図
27のブロック図に示すように構成されている。なお、
図27に示すアルゴリズムは、図9および図20に示す
モデルを組合せて構成されたものであり、ローカルフィ
ードバック排除型モデルとして構成されている。また、
図27のアルゴリズムはソフトウエアによって実現され
るが、これに含まれる機能ブロックはハードウエアにお
ける回路素子と等価な動作を行うものであるから、説明
の便宜上、ハードウエア回路として説明する。
【0063】図において共振回路112,114は、圧
力差信号phに基づいて、それぞれ上側モデルおよび外
側モデルにおける唇3の振動状態をシミュレートする。
両モデルにおける重み付けは、それぞれ乗算器113,
116を介して供給される定数u−gainおよびo−
gainによって決定される。なお、これら図27にお
いて使用される各種定数はALC63(図29参照)に
よって決定され、DSP9に供給される(詳細は後述す
る)。
【0064】さらに、乗算器116の出力信号には、オ
フセット量o−offsetが加算される。これは、実
際に金管楽器を演奏する場合において、圧力差信号ph
が「0」の場合においても、唇3がある程度の開口面積
を有することに鑑みてである。 次に、乗算器113の
出力信号と加算器150の出力信号とは加算器117に
おいて加算され、スリット関数テーブル107に供給さ
れる。これにより、両モデルの重み付けを反映したスリ
ット信号Sがスリット関数テーブル107から出力され
る。また、圧力差信号phは、図20と同様に関数テー
ブル124に供給され、この圧力差信号phに対応する
グレアム関数値をさらに圧力差信号phで除算して成る
値の信号が出力される。この信号は乗算器126を介し
てスリット信号Sと乗算され、リミット回路125に供
給される。リミット回路125は図21(ロ)に示すな
めらかなリミット関数が設定されている。
【0065】次に、リミット回路125の出力信号は、
乗算器128を介して圧力差信号phと乗算され、乗算
結果が流量信号fとして出力される。ここで、実際の金
管楽器においては、流量fに乱流成分が存在するため、
図27示す構成においては乱流成分をシミュレートとす
る構成要素151〜158が設けられている。まず、1
58は乱数発生器であり、乱数信号を出力する。次に、
151はHPF(ハイパスフィルタ)であり、乱数信号
の直流成分を除去する。152はLPF(ローパスフィ
ルタ)であり、主に音作りのために設けられ、乱数信号
の周波数スペクトルを設定する。
【0066】次に、157は乗算器であり、乱数のレベ
ルを設定するための係数r−gainを該乱数信号に乗
算する。次に、153は関数テーブルであり、流量信号
fが供給されると、|f|0.5なる値を有する信号を出
力する。そして、この信号は乗算器155を介して乱数
信号に乗算され、この乗算結果が乱流成分として加算器
154を介して流量信号fに加算される。ここで、関数
テーブル153を設けた理由は、本発明者の解析によれ
ば、乱流成分の振幅が流量fの平方根に比例する関係に
あるためである。次に、乱流成分の付加された流量信号
fは、乗算器110を介して定数zが乗算され、さらに
加算器101aを介して管体線形部102を実現するプ
ログラムに供給される。
【0067】C.第1の実施例の動作 次に、第1の実施例の動作を説明する。キーオン時の
動作 一般の電子楽器に使用されるキーボード型の操作子にあ
っては、何れかのキーを押鍵するとキーオン信号が出力
される。一方、本実施例に使用されている管楽器型のリ
ップコントローラ1においては、特にキーオン信号を出
力しない。従って、CPU2(すなわち図29における
ALC63)はキーコード関連情報(モード情報、ピス
トン情報等)を受信した時にキーオン信号を受信したも
のとみなして以下の動作を行う。まず、上記キーコード
関連情報に基づいてキーコードが求められる。ピストン
とモードキーとによってキーコードを算出するには、ま
ずモードキーによるオクターブ情報を得、これに対して
下表2に示すようにピストンキーによるシフトを追加す
るとよい。
【0068】
【表2】
【0069】次に、CPU2は求めたキーコードに対応
するピッチ周波数Fpを求め(ピッチ周波数Fpはキー
コードから一義的に求まる)、さらに式(A4)によって共
振周波数Fuを求める。さらに、当該キーコードに対応
する遅延時間DL1を管体線形部に供給する。ここで、な
んらかのキーコードが入力されることにより遅延時間D
L1’が一度設定され、さらに別のキーコードが入力され
ることによって別の遅延時間DL1”を供給する場合に
は、ポルタメントを付与すると一層好適である。すなわ
ち、遅延時間DL1をDL1’からDL1”に徐々に変化させ
ることにより、実際の金管楽器に近似したピッチ変動を
付与することができるからである。ここで、例えばトラ
ンペット等をシミュレートする場合にはポルタメント時
間を短くし、トロンボーン等をシミュレートする場合に
はポルタメント時間を長くすることが好適なことは言う
までもない。
【0070】また、キーオン時にはLPF54,59
(図31参照)の特性、あるいは出力ミックス比(詳細
は後述する)等もキーコードに応じて設定する。ここ
で、キースケールは、1オクターブあたり1ブレークポ
イントのようなものでも良いが、実際の金管楽器の特性
を考えると図34に示すように各ブレークポイントが段
差を有しているから、ブレークポイントの段差はモード
の切換わりに一致させることが好適である。また、モー
ドの切換わりにおいて楽音が急変することをヴォイシン
グで吸収するため、キースケールされるパラメータは全
て同様に変化される。
【0071】他のリップコントローラ情報を受信した
ときの動作 次に、キーコード関連情報以外のリップコントローラ情
報を受信した場合の動作を説明する。ここで、キーコー
ド関連情報以外のリップコントローラ情報としては、息
圧信号pとアンブシュア信号embとがある。まず、息
圧信号pを受信すると、CPU2は息圧信号pにキース
ケールされた定数を乗算してDSP9に出力する。一
方、アンブシュア信号embを受信した場合には、図3
5に示すアルゴリズムを介して各種定数を求め、これら
定数をDSP9に出力する。その動作を以下、詳細に説
明する。
【0072】まず、アンブシュア信号embは乗算器7
0〜76を介して比例定数が乗算され、乗算結果が加算
器77〜84を介して各定数u−gain,o−gai
n,o−offset,Qo,Fo,Qu,Fuとして
出力される。これら各定数は、上述したように、DSP
9に供給される。なお、各乗算器70〜76に供給され
る比例定数は実験等により求められる。
【0073】図35において「K.S.」とあるのはキース
ケーリングを表す定数である。これらの演算はCPU2
において実行する代りにDSP側において実行してもよ
いが、非線形部のアルゴリズムが各定数を受信する段階
までには、各定数には補完演算が施され、なめらかな変
化を有する状態にしておく必要がある。すなわち、図示
の例にあっては、各定数が変化する際に加算器77〜8
3を介してキースケーリングが施されるため、各定数は
なめらかな変化が付与されつつ出力される。
【0074】さらに、図示の構成にあっては、乱数発生
器60によって乱数信号が発生され、この乱数信号がL
PF(ローパスフィルタ)61、乗算器62、乗算器6
3〜69を順次介して加算器77〜83に供給される。
これは、金管楽器の有する不安定なアンブシュアを実現
するためである。また、乗算器62にあっては比例定数
として圧力信号pが供給されるが、これは一般に口内圧
が大きいほどアンブシュアが不安定になり、ピッチ等の
安定度が低下することに鑑みてである。
【0075】また、ピッチ情報はリップコントローラ1
を使用する場合には「1」に固定する。一方、リップコ
ントローラ1の代りに鍵盤等を用いた場合には任意の値
が設定される。この場合には、上述した遅延時間DL
L1の値に上記任意の値を乗算して出力する。実際の金
管楽器にあっては、トロンボーン等を除いてピッチベン
ドをかけることができないが、このようにすることによ
ってシンセサイザー的ななめらかなベンドやビブラート
が得られる。逆に、実際の金管楽器のような表情あるビ
ブラートをかける場合には、アンブシュアの変調を行う
方が好適であると考えられる。なお、図35において、
LPF61はバンドパスフィルタであってもよいことは
言うまでもない。
【0076】DSP9における動作 CPU2によって各定数u−gain,o−gain,
o−offset,Qo,Fo,Qu,Fu等が求めら
れると、これらの定数に応じてDSP9における非線形
部100の特性が決定される。そして、息圧信号pがD
SP9に供給されると、非線形部100、管体線形部1
02およびジャンクション101を介して圧力波信号が
伝搬され、金管楽器における発音のメカニズムがシミュ
レートされる。そして、この圧力波信号は、DAC10
(図1参照)を介してアナログ信号に変換され、サウン
ドシステム11を介して発音される。
【0077】D.変形例 本発明は上述した実施例に限定されるものではなく、例
えば以下のように種々の変形が可能であることは言うま
でもない。D.1 LCエミュレータ(LCE)を用い
た変形例 まず、上記実施例におけるリップコントローラ1は管楽
器を模した物であったが、その代りに一般的なキーボー
ドを用いてもよい。この場合には、キーボードから出力
される信号を変形することによってリップコントローラ
1から出力される信号を模擬(エミュレート)する必要
がある。ここで、キーボードから出力されるは主にキー
オン、キーオフ、キーコード等であって息圧信号pある
いはアンブシュア信号emb等は出力されない。また、
キーボードにはモジュレータ、ベンドホイールあるいは
ペダル等各種の操作子が付属している場合があるが、基
本的にはそのような操作子の補助なしにエミュレートで
きれば好適である。
【0078】エンベロープの発生方法 図36にLCE62とALC63との関係を示す。AL
C63の入力はリップコントローラからの出力か、ある
いはLCE62の出力か何れか一方が選択される。そし
て、LCE62はALC63の出力するピストン、モー
ドキー、息圧信号pおよびアンブシュアembの4種類
の情報をエミュレートして出力する。図37はLCE6
2の発生するエンベロープ信号と、息圧信号pおよびア
ンブシュアembとの関係を示す図である。図において
階段状の実線はLCE62の発生するエンベロープ信号
の波形であり、破線はハードウエアによる補間をうけて
ALC63に入力されるエンベロープ信号波形である
(なお、補間器はハードウエア側でアルゴリズムの前段
に設けられている)。
【0079】キーオンがあると、LCE62には周期的
な起動がかかる。なお、その周期は通常2〜8msec
程度である。息圧信号pの立上がりには上述したように
不感帯が存在するため、そのレベルpthまではハードウ
エアの補間定数を「0」にして立上げ、その後通常の補
間レートに設定し、エンベロープのステートをADSR
(アタック、ディケイ、サスティン、リリース)に順次
遷移してゆく。息圧信号pは音量を制御するため、その
エンベロープのアタック部のレートやレベルはベロシテ
ィや操作子情報から制御される。同様に、サスティン部
についてもペダルの操作やアフタタッチによって制御を
受ける。
【0080】アンブシュアembのエンベロープが多少
正(唇を閉じる方向)の値から開始するのは、アルゴリ
ズムの性質上、多少o−offsetなどを小さめの値
から立上げた方がスムーズなアタックが得られるからで
ある。これについても、タッチ、操作子などからアタッ
ク部のレート、レベル等を制御できた方が好適である。
次に、キーオフ信号を受信すると、できるだけ速やかん
びリリースに移行させるために周期をリセットする。そ
して、このまま音が消えるまで周期的に起動をかける
が、次のキーオンまでそのままでも良いし、適当な(レ
ベルが小となった)ところで止めてもよい。図37に示
す例にあってはリリース時にはアタック時とは逆方向に
アンブシュアembをオフセットしているが、これは必
須ではなく、演奏表現の一例として示したものである。
【0081】息圧信号p、アンブシュアembおよび
ピッチ情報の発生方法 ところで、実際の息圧信号p、アンブシュアembがキ
ーオンベロシティによって決定されるエンベロープのみ
によって動くのでは、楽音が単調になりがちである。従
って、この場合、楽音をアフタタッチ情報、モジュレー
ションあるいはペダル情報等によって修飾すると一層好
適である。その一例を図38に示す。図においてエンベ
ロープジェネレータ60の出力信号は、低周波発振器6
1およびエンベロープジェネレータ62を介して出力さ
れる操作子出力(操作子は図示せず)によってオフセッ
トが加えられる。なお、低周波発振器61の出力信号あ
るいはエンベロープジェネレータ62の各パラメータ等
も操作子によって可変できるようにすれば一層好適であ
る。そして、息圧信号p、アンブシュアembおよびピ
ッチ情報に関しても基本的に同様な方法によってエンベ
ロープを発生することができる。
【0082】各種コントローラは色々な組合せで使える
ように、図39のように任意のコントローラ入力が指示
できるようにする。図において71は加算マトリクス回
路であり、交点の「O」印は加算ポイントを示す。すな
わち、各種の操作子72〜74の操作情報が楽音に影響
を与えるように構成されており、楽音に影響を与える操
作子が適宜選択可能になっている。また、図示の回路に
あってはTEG(タッチエンベロープジェネレータ)7
0も操作子72〜74と同様に扱われる。次に、TEG
70の構成を図40に示す。図において80,82は加
算器、81は乗算器、83は遅延回路であり、乗算器8
1には補間定数αが供給される。補間定数αは、キーオ
ン後数100msec程度は小さくし、その後アフタタ
ッチに追従させるため大きくすると好適である。なお、
遅延回路83の遅延時間は数msec〜数10msec
程度にするとよい。
【0083】D.2 他の変形例 図31においては、遅延回路58を、インダクタンス要
素を付与するフィルタに変更してもよい。あるいは、図
31および図32において、遅延回路をオールパスフィ
ルタに変更してもよい。このように構成することによ
り、各モードの調和性がずれるため、特定のモードで発
振させたときにその発振波形の倍音成分が管体のモード
からずれることになり、結果として倍音の少ないくすん
だ音色になる。実際の金管楽器においても、ホルン系統
の金管楽器では同じ現象により倍音成分の減少が起こり
特有の音色を発生するため、このような金管楽器の音色
を忠実に再現することが可能になる。
【0084】
【発明の効果】以上説明したように、この発明の楽音合
成装置によれば、金管楽器の楽音を忠実にシミュレート
することが可能である。
【図面の簡単な説明】
【図1】 一実施例の構成を示すブロック図である。
【図2】 金管楽器の動作の原理を示す図である。
【図3】 金管楽器の動作の原理を示す図である。
【図4】 金管楽器の動作の原理を示す図である。
【図5】 金管楽器における各物理量の特性を示す図で
ある。
【図6】 金管楽器の物理モデルの全体構成を示すブロ
ック図である。
【図7】 非線形部100のブロック図である。
【図8】 非線形部100のブロック図である。
【図9】 非線形部100のブロック図である。
【図10】 非線形部100のブロック図である。
【図11】 共振回路112,114のブロック図であ
る。
【図12】 唇3の振動状態を示す図である。
【図13】 図9における各部の波形図である。
【図14】 グレアム関数の概形を示す図である。
【図15】 非線形部100のブロック図である。
【図16】 非線形部100のブロック図である。
【図17】 グレアム関数の概形を示す図である。
【図18】 グレアム関数の概形を示す図である。
【図19】 グレアム関数の概形を示す図である。
【図20】 非線形部100のブロック図である。
【図21】 リミット回路21の特性を示す図である。
【図22】 変形されたグレアム関数の概形を示す図で
ある。
【図23】 管体線形部102のブロック図である。
【図24】 円錐管30と点音源29との位置関係を示
す図である。
【図25】 円錐管の物理モデルのブロック図である。
【図26】 リップコントローラ1の側面図である。
【図27】 非線形部100のブロック図である。
【図28】 共振回路112,114の特性図である。
【図29】 ROM4に設定された制御プログラムのブ
ロック図である。
【図30】 DSP9の全体構成を示すブロック図であ
る。
【図31】 管体線形部102のブロック図である。
【図32】 管体線形部102のブロック図である。
【図33】 管体線形部102のチューニング方法を示
す図である。
【図34】 金管楽器のブレークポイントを示す図であ
る。
【図35】 ALC63のアルゴリズムの構成を示すブ
ロック図である。
【図36】 ROM4に設定された制御プログラムのブ
ロック図である。
【図37】 LCE62の動作を示す図である。
【図38】 一実施例の変形例の概要を示す図である。
【図39】 一実施例の変形例の概要を示す図である。
【図40】 一実施例の変形例の概要を示す図である。
【図41】 円筒管近似したマウスピース部およびベル
部を示す図である。
【図42】 金管楽器のヒステリシス特性を示す図であ
る。
【符号の説明】
中央処理装置(押圧力情報発生手段、息圧情報発生
手段、周波数制御手段、第1の共振手段、遅延手段、乱
数成分付加手段、乱数レベル制御手段、第2の共振手
段、第3の共振手段)32アンブシュアセンサ(押
圧力情報発生手段)33圧力センサ(息圧情報発生
手段)100非線形部(周波数制御手段、第1の共
振手段)102管体線形部(遅延手段)158
乱数発生器(乱数成分付加手段)153関数テーブ
ル(乱数レベル制御手段)112共振回路(第2の
共振手段)114共振回路(第3の共振手段)
フロントページの続き (58)調査した分野(Int.Cl.7,DB名) G10H 1/053 G10H 1/00 G10H 7/00 513 G10H 7/08

Claims (3)

    (57)【特許請求の範囲】
  1. 【請求項1】 演奏者の息圧に関する息圧情報を発生
    る息圧情報発生手段と、 前記演奏者の唇の押圧力に関する押圧力情報を発生する
    押圧力情報発生手段と、 前記押圧力が大となると楽音信号の周波数を高くする周
    波数制御手段と、 前記楽音信号の周波数に乱数成分を付加する乱数成分付
    加手段と、 前記息圧が大となると前記乱数成分のレベルを大とする
    乱数レベル制御手段とを具備することを特徴とする楽音
    合成装置。
  2. 【請求項2】 演奏者の息圧に関する息圧情報を発生
    る息圧情報発生手段と、 前記演奏者の唇の押圧力に関する押圧力情報を発生する
    押圧力情報発生手段と、 前記息圧および押圧力によって共振特性が変化する第1
    の共振手段と、 前記共振手段から出力された信号波を遅延させて前記第
    1の共振手段に帰還させる遅延手段とを具備することを
    特徴とする楽音合成装置。
  3. 【請求項3】 前記第1の共振手段は、楽音周波数近傍
    の共振周波数を有する第2の共振手段と、低周波の共振
    周波数を有する第3の共振手段とを具備することを特徴
    とする請求項2に記載の楽音合成装置。
JP03104805A 1991-04-10 1991-04-10 楽音合成装置 Expired - Fee Related JP3097167B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP03104805A JP3097167B2 (ja) 1991-04-10 1991-04-10 楽音合成装置
US07/864,632 US5272275A (en) 1991-04-10 1992-04-07 Brass instrument type tone synthesizer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP03104805A JP3097167B2 (ja) 1991-04-10 1991-04-10 楽音合成装置

Publications (2)

Publication Number Publication Date
JPH04311995A JPH04311995A (ja) 1992-11-04
JP3097167B2 true JP3097167B2 (ja) 2000-10-10

Family

ID=14390644

Family Applications (1)

Application Number Title Priority Date Filing Date
JP03104805A Expired - Fee Related JP3097167B2 (ja) 1991-04-10 1991-04-10 楽音合成装置

Country Status (2)

Country Link
US (1) US5272275A (ja)
JP (1) JP3097167B2 (ja)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3160981B2 (ja) * 1991-12-13 2001-04-25 ヤマハ株式会社 電子楽器用音源の制御装置
JP2707911B2 (ja) * 1992-03-03 1998-02-04 ヤマハ株式会社 楽音合成装置
US5841054A (en) * 1996-04-06 1998-11-24 Yamaha Corporation Musical tone synthesizing apparatus having competibility of tone color parameters for different systems
JP3045375B2 (ja) * 1996-05-18 2000-05-29 ヤマハ株式会社 楽音合成装置
WO1999051314A1 (en) * 1998-04-03 1999-10-14 Glenn White A toy
FR2792125B1 (fr) * 1999-04-08 2001-06-08 France Telecom Procede de simulation de la propagation non lineaire d'une onde acoustique, notamment dans un resonateur
AU750620B2 (en) * 2000-08-28 2002-07-25 Glen Jacobs Didgeridoo containing an inserted sound emitting device
FR2846768B1 (fr) * 2002-10-31 2005-07-08 Centre Nat Rech Scient Procede de simulation et de synthese numerique d'un phenomene oscillant
ITMC20030032A1 (it) * 2003-03-28 2004-09-29 Viscount Internat Spa Metodo e dispositivo elettronico per riprodurre il suono delle canne ad anima dell'organo liturgico, sfruttando la tecnica della modellazione fisica degli strumenti acustici
US6995307B2 (en) * 2003-06-30 2006-02-07 S&D Consulting International, Ltd. Self-playing musical device
JP6493689B2 (ja) * 2016-09-21 2019-04-03 カシオ計算機株式会社 電子管楽器、楽音生成装置、楽音生成方法、及びプログラム

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4736663A (en) * 1984-10-19 1988-04-12 California Institute Of Technology Electronic system for synthesizing and combining voices of musical instruments
US4984276A (en) * 1986-05-02 1991-01-08 The Board Of Trustees Of The Leland Stanford Junior University Digital signal processing using waveguide networks
US5371317A (en) * 1989-04-20 1994-12-06 Yamaha Corporation Musical tone synthesizing apparatus with sound hole simulation
JP2679299B2 (ja) * 1989-10-03 1997-11-19 ヤマハ株式会社 楽音合成装置
JPH0782342B2 (ja) * 1989-10-04 1995-09-06 ヤマハ株式会社 楽音合成装置
JPH0774956B2 (ja) * 1989-12-06 1995-08-09 ヤマハ株式会社 楽音波形信号形成装置
JPH0778679B2 (ja) * 1989-12-18 1995-08-23 ヤマハ株式会社 楽音波形信号形成装置
JPH0774957B2 (ja) * 1989-12-19 1995-08-09 ヤマハ株式会社 楽音波形信号形成装置

Also Published As

Publication number Publication date
US5272275A (en) 1993-12-21
JPH04311995A (ja) 1992-11-04

Similar Documents

Publication Publication Date Title
JP3097167B2 (ja) 楽音合成装置
US5508473A (en) Music synthesizer and method for simulating period synchronous noise associated with air flows in wind instruments
US5519167A (en) Musical tone synthesizing apparatus
Gibiat et al. Period doubling occurences in wind instruments musical performance
JP2727841B2 (ja) 楽音合成装置
JPH0778679B2 (ja) 楽音波形信号形成装置
US5521328A (en) Electronic musical instrument for simulating wind instrument musical tones
US5371317A (en) Musical tone synthesizing apparatus with sound hole simulation
JP2722947B2 (ja) 楽音波形信号形成装置
JP3149708B2 (ja) 楽音合成装置
JPH0452697A (ja) 楽音合成装置
JP3459948B2 (ja) 楽音合成装置
Cook Implementation of single reed instruments with arbitrary bore shapes using digital waveguide filters
JP3223683B2 (ja) 楽音波形信号合成装置
JP2568759B2 (ja) 楽音合成装置
JP2765479B2 (ja) 楽音波形信号形成装置
JP2824555B2 (ja) 音圧推定装置
Scavone Time-domain synthesis of conical bore instrument sounds
JPH03210600A (ja) 楽音制御装置
JP2770676B2 (ja) 電子楽器
Réveillac Synthesizers and Subtractive Synthesis 1: Theory and Overview
JP2743781B2 (ja) 楽音信号発生装置
JP3484737B2 (ja) 楽音特性制御装置
JP3413890B2 (ja) 楽音合成装置
JP2725497B2 (ja) 楽音合成装置

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20000711

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313532

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20070811

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080811

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090811

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees