山内_backup の履歴(No.30)

B4 前期授業スケジュール

	月曜日	火曜日	水曜日	木曜日	金曜日
1-2					研究会
3-4					研究会
5-6			集まる日
7-8		研究会(サブ)
9-10		研究会(サブ)
11-12

章立て

やることリスト

やること
本論　表紙、章立て更新
参考文献

提案システムの流れ　図で書く
提案システムの有効性　初歩的な例を図表で書く

ポスター
本論から文章と図表
図表の中身作成

＊＊＊＊リアルタイムに処理
https://qiita.com/ha161553/items/80f3056544a3d4ae3352

https://qiita.com/kinopee0120/items/7bbcda4afa7e501d6272

リアルタイムhttps://skimie.com/articles/6a3bfa82712f59cb6b5a6c10d7

import sys
import scipy.io.wavfile
import math
import numpy as np
import matplotlib.pyplot as plt

import math
from scipy import signal
from scipy import fftpack

class GetMonoInterval:
    #音名
    vecInterval=[]
    #周波数
    vecFreq=[]
    #周波数：音名
    mapInterval=[]
    #ファイル名
    filename=""
    #生データ
    data=[]
    #サンプリング周期
    rate=0
    #時間軸
    time=[]
    #FFT処理後データ
    fftData=[]
    #周波数リスト
    fftFreq=[]
    
    #コンストラクタ
    def __init__(self,file):
        #音名マップ
        #音名の配列作成
        self.vecInterval = ["A0","A#0","B0"]
        for i in range(1,8):
            self.vecInterval.append("C" + str(i))
            self.vecInterval.append("C#" + str(i))
            self.vecInterval.append("D" + str(i))
            self.vecInterval.append("D#" + str(i))
            self.vecInterval.append("E" + str(i))
            self.vecInterval.append("F" + str(i))
            self.vecInterval.append("F#" + str(i))
            self.vecInterval.append("G" + str(i))
            self.vecInterval.append("G#" + str(i))
            self.vecInterval.append("A" + str(i))
            self.vecInterval.append("A#" + str(i))
            self.vecInterval.append("B" + str(i))
        self.vecInterval.append("C8")
        #音程の周波数の配列作成
        self.vecFreq.append(27.5)
        semtiones = pow(2,1/12)
        for i in range(0,87):
            self.vecFreq.append(self.vecFreq[i]*semtiones)
        #二つの配列で「周波数:音名」dictを作成
        self.mapInterval = dict(zip(self.vecFreq,self.vecInterval))
        
        #ファイル読み込み
        self.filename = file
        
        #音源読み込み
        self.rate, self.data = scipy.io.wavfile.read(self.filename)
        #横軸（時間）を作成
        self.time = np.arange(0,self.data.shape[0]/self.rate,1/self.rate)    
        
        #FFT
        self.fftData = np.abs(np.fft.rfft(self.data))
        self.fftFreq = np.fft.fftfreq(self.fftData.shape[0],1.0/self.rate)
        
    #そのままの波形を表示
    def displayFigUneditWav(self):

#データプロット
        fig = plt.figure()
        ax = fig.add_subplot(111)
        #軸ラベル
        ax.set_xlabel("Time[s]")
        ax.set_ylabel("Amplitude")
        #プロット
        ax.plot(self.time,self.data,label=self.filename)

plt.legend()
        #グリッド
        plt.grid(axis="both")

#表示
        plt.show()
            
    #FFTを処理した波形を表示
    def displayFigFFTWav(self,peak):

#ピーク値を取得
        fft_max_data = signal.argrelmax(self.fftData,order=peak)

#データプロット
        fig = plt.figure()
        ax = fig.add_subplot(111)
        #軸ラベル
        ax.set_xlabel("Frequency[Hz]")
        ax.set_ylabel("Gain")
        #プロット
        ax.loglog(self.fftFreq,self.fftData,".-",label=self.filename)
        ax.plot(self.fftFreq[fft_max_data],self.fftData[fft_max_data],"ro",label="peak")
        ax.legend()
        #グリッド
        ax.grid(axis="both")
        #表示
        plt.show()

#最もゲインの大きいピークの周波数を返す
        maxPeakFreq = self.getPeakFreq(self.fftFreq[fft_max_data],self.fftData[fft_max_data])
        #その周波数に最も近い音名を返す
        intervalName = self.searchInterval(maxPeakFreq)
        return intervalName,self.fftFreq[fft_max_data]
        
    
   #searchInterval(周波数)
    def searchInterval(self,freq):
        """
        @freq  :音程を知りたい周波数
        @return:音程を返す
        """
        #平均率で最も近い周波数を検索
        nearfreq = self.getNearestValue(self.vecFreq,freq)
        #その音名
        nameInterval = self.mapInterval[nearfreq]
        return nameInterval
        
    #最も近い値を返す
    def getNearestValue(self,list, num):
        """
        概要: リストからある値に最も近い値を返却する関数
        @param list: データ配列
        @param num: 対象値
        @return 対象値に最も近い値
        """
        # リスト要素と対象値の差分を計算し最小値のインデックスを取得
        idx = np.abs(np.asarray(list) - num).argmin()
        return list[idx]
        
    #最もゲインの大きいピークの周波数を返す(周波数配列、ゲイン配列)
    def getPeakFreq(self,peakFreq,peakGain):
        mapMaxPeak = dict(zip(peakGain,peakFreq))
        maxPeakValue = np.amax(peakGain)
        return mapMaxPeak[maxPeakValue]

maou = GetMonoInterval("ou-A.wav")
gare = GetMonoInterval("gare-A-mono.wav")
koe  = GetMonoInterval("koe-A.wav")

# maou.displayFigUneditWav()
# koe.displayFigUneditWav()

maouInterval, maouFreq = maou.displayFigFFTWav(1000)
# maouInterval

# gareInterval, gareFreq = gare.displayFigFFTWav(900)
# gareInterval

gareInterval, gareFreq = gare.displayFigFFTWav(900)
gareInterval

maouFreq[1:5]
gareFreq[0:4]

https://skimie.com/articles/6a3bfa82712f59cb6b5a6c10d7

研究会(山内)

卒論・中間発表(山内)

勉強用(山内)

B4 前期授業スケジュール

章立て

やることリスト

参考文献(山内)

専門ゼミ(山内)