音の周波数特性を描くMatlabプログラムを改善したので、メモ的に置いておきます。今回の目玉は「1/3オクターブごとの平均値を描くようにした」です。500-2000Hz帯の平均値が0dBになるようにスケーリングしています。 function [f, MX] = magnitude(x, Fs) %[f, MX] = magnitude(x, Fs) % Calculate and correctly scale the output of the FFT function to obtain % a meaningful power versus frequency plot. This function plots the % result if no output variables are specified. % % Input: % x - signal % fs - sampl