介绍:

(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 的图











