2021-06-20から1日間の記事一覧

スペクトログラム

import os import numpy as np from scipy import signal from scipy.io import wavfile import wave import matplotlib.pyplot as plt def WaveAnalyzer(filename): path = os.getcwd() fp = path+'/'+filename print(fp) wavefile = wave.open(fp,"rb") ra…

カイザー窓の振幅スペクトル その2

clear; close all; clc; beta=2.3052; N=2^5; w_kaiser=kaiser(N,beta); amp=2.0*fft(w_kaiser,2048)/N; freq = linspace(-0.5,0.5,length(amp)); mag = 20*log10(abs(fftshift(amp/max(abs(amp))))); plot(freq,mag); title('Kaiser window: N=32'); figure…