編集履歴

回答編集履歴

日本語が下手すぎる点を訂正

2018/02/16 10:54

投稿

スコア18394

test CHANGED Viewed

@@ -6,13 +6,13 @@
 - 最大振幅からいきなり振幅0にしても同じ
-- 同時撥音数が変化した際にampを突然変化させると波に連続性がなくなるため音割れする
-同時撥音数というのは和音を再生する際には結構変化するので、本式のシンセサイザー（あるいはミキサー）には波形を調整するコンプレッサーとかリミッターと呼ばれるものが備わっていると思います。これはご質問のコードにあるように単に範囲外をclampするという方法では実現できません。そういう方式では波形がひずむため結局音割れします。
-波形を直接生成するような場合、こうした音の特性に注意する必要があると思います。自然界の音は[エンベロープ](https://www.g200kg.com/jp/docs/dic/envelope.html)に従って振幅が徐々に変化するのが普通です。MIDIのシンセサイザーにはピンキリがあると思いますが、どんな単純なものでも最小限のエンベロープ制御を行っており波形が不連続にならないように（音割れしないように）なっていると思います。
+- 同時発音数が変化した際にampを突然変化させると波に連続性がなくなるため音割れする
+同時発音数というのは和音を再生する際には結構変化するので、本式のシンセサイザー（あるいはミキサー）には波形を調整するコンプレッサーとかリミッターと呼ばれるものが備わっています。これはご質問のコードにあるように単に範囲外をclampするという方法では実現できません。そうやってしまうと波形がひずみ結局音割れします。
+波形を直接生成するような場合、こうした音の特性に注意する必要があります。自然界の音は[エンベロープ](https://www.g200kg.com/jp/docs/dic/envelope.html)に従って振幅が徐々に変化するのが普通ですので、ソフトウェアで波形を生成する場合でも最小限のエンベロープ制御を行い波形が不連続にならないようにするとよいでしょう。
@@ -44,7 +44,7 @@
-なお、前述したコンプレッサーとかリミッターをソフトウェア的に実現するのは少々面倒だと思います。振幅がオーバーフローしそうになったときだけ滑らかに波形を調整しないといけないのですが、**滑らかに調整**というのがなかなかめんどくさいのです。そこで、とりあえずは同時撥音数を適当な値に固定してしまい(例えば8とか)、演奏する際にその限界を超えて発音するような事態を避けるのが簡単と思います。同時撥音数というのは同時にNOTE_ONとなっている数の最大数よりも大きくなる点に注意してください。NOTE_OFFしてもリリースタイムの期間の間振幅を徐々に小さくしながら音が鳴り続けるため「音が完全に消えるまでには若干時間がかかる」からです。
+なお、前述したコンプレッサーとかリミッターをソフトウェア的に実現するのは少々面倒だと思います。振幅がオーバーフローしそうになったときだけ滑らかに波形を調整しないといけないのですが、**滑らかに調整**というのがなかなかめんどくさいのです。そこで、とりあえずは同時発音数を適当な値に固定してしまい(例えば8とか)、演奏する際にその限界を超えて発音するような事態を避けるのが簡単でしょう。同時発音数というのは同時にNOTE_ONとなっている数の最大数よりも大きくなる点に注意してください。NOTE_OFFしてもリリースタイムの期間の間振幅を徐々に小さくしながら音が鳴り続けるため「音が完全に消えるまでには若干時間がかかる」からです。
@@ -54,7 +54,7 @@
-下記コードでは0.4秒ごとにC4～B#4までの音を順に自動演奏していますが、その限りでは雑音はでないようでした。しかし・・・エンベロープを滑らかにしたつもりかつ同時撥音数10固定なんですが、自動演奏で和音にしたとたん演奏開始時点で雑音が・・・
+下記コードでは0.4秒ごとにC4～B#4までの音を順に自動演奏していますが、その限りでは雑音はでないようでした。しかし・・・エンベロープを滑らかにしたつもりかつ同時発音数10固定なんですが、自動演奏で和音にしたとたん演奏開始時点で雑音が・・・

誤記訂正・コードコメント訂正

2018/02/16 10:54

投稿

KSwordOfHaste

スコア18394

test CHANGED Viewed

@@ -44,7 +44,7 @@
-なお、前述したコンプレッサーとかリミッターをソフトウェア的に実現するのは少々面倒だと思います。振幅がオーバーフローしそうになったときだけ滑らかに波形を調整しないといけないのですが、**滑らかに調整**というのがなかなかめんどくさいのです。そこで、とりあえずは同時撥音数を適当な値に固定してしまい(例えば8とか)、演奏する際にその限界を超えて発音するような事態を避けるのが簡単と思います。同時撥音数というのは同時にNOTE_ONとなっている数の最大数よりも大きくなる点に注意してください。NOTE_OFFしてもリリースタイムの期間の間神父が徐々に小さくなるため「音が完全に消えるまでには若干時間がかかる」からです。
+なお、前述したコンプレッサーとかリミッターをソフトウェア的に実現するのは少々面倒だと思います。振幅がオーバーフローしそうになったときだけ滑らかに波形を調整しないといけないのですが、**滑らかに調整**というのがなかなかめんどくさいのです。そこで、とりあえずは同時撥音数を適当な値に固定してしまい(例えば8とか)、演奏する際にその限界を超えて発音するような事態を避けるのが簡単と思います。同時撥音数というのは同時にNOTE_ONとなっている数の最大数よりも大きくなる点に注意してください。NOTE_OFFしてもリリースタイムの期間の間振幅を徐々に小さくしながら音が鳴り続けるため「音が完全に消えるまでには若干時間がかかる」からです。
@@ -138,7 +138,7 @@
         # 以下のように、3度の和音にしたら同時発音数=10でも音割れした。(Win10 64bit)
-        # なぜそうなるか分かりませんでした... ＞＜
+        # =>PyCharmのデバッグ実行では音割れしたが通常実行にすると音割れしなくなった
         # ---------------------------------------------

追記２

2018/02/16 10:49

投稿

KSwordOfHaste

スコア18394

test CHANGED Viewed

@@ -62,6 +62,12 @@
+追記２：PyCharm上で実行していたのですが、ふと思いついて「デバッグ実行」ではなく普通の「実行」にしたら和音にしても雑音でなくなりました!
+想像ですが、PyCharmでデバッグ実行していたためインタープリタの実行スピードが音声再生においついてなかったのだと思います・・・
 ```Python
 import math

コード追記

2018/02/16 09:44

投稿

KSwordOfHaste

スコア18394

test CHANGED Viewed

@@ -45,3 +45,355 @@
 なお、前述したコンプレッサーとかリミッターをソフトウェア的に実現するのは少々面倒だと思います。振幅がオーバーフローしそうになったときだけ滑らかに波形を調整しないといけないのですが、**滑らかに調整**というのがなかなかめんどくさいのです。そこで、とりあえずは同時撥音数を適当な値に固定してしまい(例えば8とか)、演奏する際にその限界を超えて発音するような事態を避けるのが簡単と思います。同時撥音数というのは同時にNOTE_ONとなっている数の最大数よりも大きくなる点に注意してください。NOTE_OFFしてもリリースタイムの期間の間神父が徐々に小さくなるため「音が完全に消えるまでには若干時間がかかる」からです。
+---
+追記：実際に単純なエンベロープにしてやってみた例を挙げてみます。少々長いのでちょっと気が引けますが。自分の環境ではMIDI-INがすぐに使えないので、Key_inputの部分は自動演奏に置き換えています。質問者さんの環境ではこの関数を書き替えずに実行すれば演奏の通りに音が鳴らせると思います。（すみませんが、自動演奏に変えたため元のコードからpygame.midiの部分のコードは削除してあります）
+下記コードでは0.4秒ごとにC4～B#4までの音を順に自動演奏していますが、その限りでは雑音はでないようでした。しかし・・・エンベロープを滑らかにしたつもりかつ同時撥音数10固定なんですが、自動演奏で和音にしたとたん演奏開始時点で雑音が・・・
+単なるバグなのかも知れませんし、Pythonのスピードとか自分の配慮が不足している点が何かあるのかも知れません・・・
+```Python
+import math
+import pyaudio
+import struct
+import time
+# PyAudio制御関連=============================
+SAMPLE_RATE = 44100
+if __name__ == '__main__':
+    p = pyaudio.PyAudio()
+    # ストリームを開く
+    stream = p.open(format=pyaudio.paInt16, channels=1, rate=SAMPLE_RATE, output=1, stream_callback=None)
+def End():   # 終了処理(必ず入れる)
+    stream.close()
+    p.terminate()
+# Key管理================================
+key_list = []  # 押されてるキーを保持するリスト
+start_time = time.time()  # MIDI-INなしで自動演奏するための小細工
+def Key_input():
+    # MIDI-INなしで自動演奏するための小細工
+    # MIDI-INが使えるならオリジナルのコードそのままでよいはず・・・
+    global key_list
+    t = (time.time() - start_time)
+    if 3 < t < 20:
+        key = 60 + int((t - 3) / 0.4) % 12
+        key_list = [[key, 110]]
+        # 以下のように、3度の和音にしたら同時発音数=10でも音割れした。(Win10 64bit)
+        # なぜそうなるか分かりませんでした... ＞＜
+        # ---------------------------------------------
+        # key_list = [[key, 120], [key + 4, 100]]
+    else:
+        key_list = []
+    return
+# 周波数とか作るところ=======================
+class Oscillator:
+    def __init__(self, key, max_level, sustain_level, attack_time, decay_time, release_time):
+        self.key = key
+        freq = 440.0 * (2 ** ((key - 69) / 12))
+        self.dt = math.pi * 2 * freq / SAMPLE_RATE
+        self.max = max_level
+        self.sustain = sustain_level
+        self.attack_t = int(attack_time * SAMPLE_RATE)
+        self.decay_t = int(decay_time * SAMPLE_RATE)
+        self.release_t = int(release_time * SAMPLE_RATE)
+        self.released = None
+        self.time = 0
+        self.amp = 0
+    def calc_amp(self):
+        if self.released is None:
+            # まだNOTE_ON状態
+            if self.time < self.attack_t:
+                return Oscillator.interpolate(self.time, self.attack_t, 0.0, self.max)
+            else:
+                return Oscillator.interpolate(self.time - self.attack_t, self.decay_t, self.max, self.sustain)
+        else:
+            # NOTE_OFF状態
+            return Oscillator.interpolate(self.time - self.released, self.release_t, self.sustain, 0.0)
+    # 振幅の直線補間用
+    @staticmethod
+    def interpolate(cur_t, max_t, start_level, end_level):
+        t = cur_t / max_t if cur_t < max_t else 1
+        return start_level * (1 - t) + end_level * t
+    def wave(self):
+        self.amp = self.calc_amp()
+        w = self.amp * math.sin(self.time * self.dt)
+        self.time += 1
+        return w
+    def note_off(self):
+        if self.released is None:
+            self.sustain = self.amp  # attack/decay中にNOTE_OFFされた場合への配慮
+            self.released = self.time
+            return True
+        else:
+            return False
+    def is_playing(self):
+        return self.released is None or self.time < self.released + self.release_t
+class Synthesizer:
+    def __init__(self, polyphony=8):
+        self.oscillators = {}
+        self.polyphony = polyphony
+        buffer_size = 1048 // 2  # 44100 1sec 4410 0.1sec 441 0.01sec
+        latency = buffer_size / SAMPLE_RATE
+        self.wave = [0] * int(latency * SAMPLE_RATE)
+        self.pack_format = 'h' * len(self.wave)
+        print("polyphony=", self.polyphony, "latency=", latency)
+    def generate_wave(self, keys):
+        """実際の波形を生成"""
+        self.update_oscillators(keys)
+        wave = self.wave
+        # print(keys)
+        amp = 32767 / self.polyphony  # 同時発音数を固定とする
+        for i in range(len(self.wave)):
+            s = int(amp * sum(map(Oscillator.wave, self.oscillators.values())))
+            wave[i] = s  # max(-32767, min(s, 32767))  # clampしない
+        wave = struct.pack(self.pack_format, *wave)
+        return wave
+    def update_oscillators(self, keys):
+        """NOTE_ONのキー集合によりオシレーター群を更新する"""
+        # もはや発音しないオシレーターを削除
+        osc_exists = len(self.oscillators) > 0
+        self.oscillators = {osc.key: osc for osc in filter(Oscillator.is_playing, self.oscillators.values())}
+        if osc_exists and len(self.oscillators) == 0:
+            print("all oscillators removed")
+        # 新たにNOTE_ONとなったオシレーターを追加
+        new_keys = filter(lambda key: key not in self.oscillators, keys)
+        for key in new_keys:
+            max_level = keys[key] / 127
+            sustain_level = max_level * 0.9
+            attack_time = 0.3
+            decay_time = 0.1
+            release_time = 0.1
+            self.oscillators[key] = Oscillator(key, max_level, sustain_level, attack_time, decay_time, release_time)
+            print("note on", key)
+        # NOTE_OFFとなったオシレーターの状態を変更(リリース制御のため)
+        for osc_key in self.oscillators:
+            if osc_key not in keys:
+                if self.oscillators[osc_key].note_off():
+                    print("note off", osc_key)
+# 再生================================
+def play(stream, data):  # 再生用関数、ストリームと波形データを引数に
+    # チャンク単位でストリームに出力し音声を再生
+    chunk = 1048
+    sp = 0  # 再生位置ポインタ
+    buffer = data[sp:sp+chunk]
+    while buffer:
+        stream.write(buffer)
+        sp = sp + chunk
+        buffer = data[sp:sp+chunk]
+# ==========================================
+if __name__ == '__main__':
+    synthesiser = Synthesizer(10)
+    while True:
+        Key_input()  # キーの取得
+        keys = {k[0]: k[1] for k in key_list}
+        wave = synthesiser.generate_wave(keys)
+        play(stream, wave)  # 波形の再生
+    End()  # このままでは決して呼び出されませんが...
+```

文章不備訂正

2018/02/16 09:33

投稿

KSwordOfHaste

スコア18394

test CHANGED Viewed

@@ -12,11 +12,11 @@
-こうした音の特性に注意する必要があるのではないでしょうか。自然界の音は[エンベロープ](https://www.g200kg.com/jp/docs/dic/envelope.html)に従って振幅が徐々に変化するのが普通です。MIDIのシンセサイザーにはピンキリがあると思いますが、どんな単純なものでも最小限のエンベロープの配慮により波形が不連続になることによる音割れを避けるようになっていると思います。
+波形を直接生成するような場合、こうした音の特性に注意する必要があると思います。自然界の音は[エンベロープ](https://www.g200kg.com/jp/docs/dic/envelope.html)に従って振幅が徐々に変化するのが普通です。MIDIのシンセサイザーにはピンキリがあると思いますが、どんな単純なものでも最小限のエンベロープ制御を行っており波形が不連続にならないように（音割れしないように）なっていると思います。
-基本体なエンベロープのしかたとしては以下のような感じに考えるとよいと思います。
+なお、基本体なエンベロープ制御は以下のような感じに考えるとよいと思います。