基于matlab的电子琴信号分析

                     基于matlab的电子琴信号分析

交流电频率是50Hz，当经过电力变压器附近时，我们总会听到“呜呜”的声音，这个信号就是几乎标准的50Hz纯音，夏天蚊子从耳边飞过时候的“嗡嗡”声，也几乎是标准的高频纯音信号。电子琴的频率范围从27.5Hz到4186Hz，远远覆盖了从变压器到蚊子振动翅膀的频率，但是从直观上感觉，这个信号好像和电子琴的任何一个琴键发出的信号都有很大差别。这是为什么？是我们的听觉欺骗了我们，还是电子琴的频率不是单一的标准频率？

任何一本乐理知识上都指出了中央C的音频是261.6Hz，现在我们来分析一下这个信号到底是不是纯音信号。

 

用matlab生成一个标准的261.6Hz的频率来实际感受一下,程序如下：

clear

fs=8000;

N=8000; 

n=0:N-1;

t=n/fs; 

x=sin(2*pi*261.6*t);

sound(x,fs);

运行程序后发现这个信号和实际的中央C信号的听觉效果完全不一样，这是什么原因造成的呢？

用录音设备录制中央C的声音，采样频率设定为8000.保存到E盘根目录，我这边命名为E:/DO.wav。现在用matlab进行信号的fft分析得出其频谱。程序如下：

 

clear

fs_do=8000;%采样频率

x_do=audioread('E:/do.wav'); %信号读入

x_do=x_do';

N_do=8000; %采样点数

n_do=0:N_do-1;

t_do=n_do/fs_do; %以上两行生成时间序列

y_do=fft(x_do,N_do);%fft

mag_do=abs(y_do); %模值

ang_do=angle(y_do);%相角

f_do=n_do*fs_do/N_do; %频率序列

plot(f_do,mag_do);

title('电子琴DO FFT ')

axis([-100 2000 0 200]);

sound(x_do,8000);%声音试听

运行后波形如图（1）：

              

 把采集到的声音信号用matlab进行