介绍:

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