介绍:
(1)C 和弦录音存储在一个名为 Cchord.mat 的文件中。此文件包含两个变量:
y:录音信号
fs:采样频率
此任务使用 numel 函数返回数组中的元素数目。
任务:加载文件 Cchord.mat。创建一个名为 n 的变量,表示 y 中的元素数目。然后使用 n 创建等间距向量 t,该向量以 0 开头,以 n-1 结尾,元素之间的间距为 1。
(2)t
现在有正确的点数,但它需要表示音频信号的采样时间。您可以使用采样频率 fs
将向量转换为时间(以秒为单位)。
任务:将 t
除以 fs
。将输出赋给同一个变量 t
。然后绘制 y
对 t
的图。
(3)在绘图中,请注意 y
是周期性的,但它不是简单的正弦波。它由具有不同频率的多个正弦波组成。傅里叶变换将返回信号的频谱信息。主频的位置将显示和弦中包含的音符。
可以使用 fft
函数来计算向量的离散傅里叶变换。
fft(y)
fft
的输出值为复数。您可以使用 abs
函数来获得幅值。
任务: 创建一个名为 yfft
的变量,表示 y
的离散傅里叶变换的绝对值。
(4)在任务 1 和 2 中,您已计算了信号 y 的时间向量 t。同样,您需要为您的 FFT 向量 yfft 计算频率向量 f。
任务: 创建等间距向量 f,该向量以 0 开头,以 n-1 结尾,元素之间的间距为 1。
(5)向量 f
现在包含 n
个点。要将这些点转换为频率,您可以将整个向量乘以采样频率 (fs
),然后除以点数 (n
)。f
将包含从 0
到 fs
的频率。主频位于 f
的开头位置。您可以使用 xlim
函数放大所关注的区域。
xlim([xmin xmax])
任务: 将 f
乘以 fs/n
。将输出赋给同一个变量 f
。 使用 x 的限值 0
和 1000
绘制 yfft
对 f
的图