%创建fs表示音频信号的采样频率
%创建t表示音频信号的采样时间
>> t=0:(1/fs):20
%创建y表示每个采样时刻的音频信号幅值
>> y=sin(1.8*2*pi*t)+sin(2.1*2*pi*t)
%绘制y对t的图,可以看到“拍频现象”
>> plot(t,y)
%傅里叶变换返回信号的频率信息
>> yfft=fft(y)
%使用numel函数返回数组中的元素数目
>> n=numel(y)
%创建f变量表示与yfft中的值对应的频率
>> f=0:(fs/n):(fs*(n-1)/n)
fft的输出值为复数,绘制它们的模对f的图。可以看到两个峰值,这两个峰值与之前创建的两个正弦波的频率相对应,由于这些峰值非常接近,因此信号会出现拍频现象。
>> plot(f,abs(yfft))
网友评论