停滞しがちな VCA の製作と並行してエンベロープジェネレータの試作を開始しました。今回のものはエンベロープは計算で生成します。マイクロコントローラは前回と同じく CAN コントローラ内蔵の STM32C092KCT を使います。
マイクロコントローラでエンベロープジェネレータを実装することにはいくつか期待する利点があります。
- CAN バスからゲート信号を拾うことができるのでゲート信号のためのパッチングをしなくても動く
- ベロシティー対応にするが、アナログのゲート信号にベロシティー情報を乗せると精度が問題になるが CAN バスからゲート信号を拾えば精度が問題になることはない
- 計算さえ間に合えばどんな非線形曲線も生成できる。ADSR にこだわる必要もなくもっとステップ数を増やしたりもできる
- パラメータをセーブ・ロードできる。これは Analog3 の重要な特徴
まずは Analog3 共通の CAN 周りの基本機能を実装しました。この部分はゆくゆく Github の analog3 プロジェクトに共用のコードとして置く予定です。次にエンベロープ生成アルゴリズムの実装です。このマイクロコントローラは別に信号処理用というわけではないので、FPU などはついておらず、演算速度を稼ぐためには固定小数点演算で信号生成をしたいところです。固定小数点演算のやり方は忘れてしまったので過去の記事を読んで復習します。
曲線を求める計算は next = r * prev を繰り返すだけで良くて単純でした。問題は、ノブの位置から減衰率を計算する方法です。まずは愚直にノブの位置と上昇・減衰の時定数 τ の関係を指数関数にして、さらに時定数から減衰率 r を求めることにします。まず時定数の定義から減衰曲線は以下の通り
減衰率 r は、エンベロープの更新の時間差 Δt の間にどれぐらいの減衰するか決める値なので、
時定数をノブの位置 k の指数関数とすると
これを減衰率を求める式に当てはめると
ややこしい計算式ですが、実装が固まってくるまでこれを素直に計算することにします。時定数を求める式に現れる定数 α と β はエンベロープジェネレータが快適に使える範囲に決めないといけません。これは試行錯誤で決めるよりほかにありません。ノブの値は実際には ADC の結果を 16 ビット整数に読み込むので範囲は 0 から 65535 までです。gnuplot の力を借りて、α = 0.00015、β = 7.5 とすると、大まかに良い感じの時定数範囲が得られました。
gnuplot> plot [k=0:65535] exp(k * 0.00015 - 7.5)
ここから減衰率は
gnuplot> plot [r=0:65535] exp(-1.e-3/(exp(r * 0.000115 - 7.5)))
時定数が長くなるにつれ減衰率は 1.0 に近づくことが見て取れます。最初の版では、この計算式を使って減衰率を求めることにしました。
（続く）
