matlab 项目 - 音频频率

C 和弦录音存储在一个名为 Cchord.mat 的文件中。此文件包含两个变量：

y：录音信号
fs：采样频率

此任务使用 numel 函数返回数组中的元素数目。

任务

加载文件 Cchord.mat。

创建一个名为 n 的变量，表示 y 中的元素数目。然后使用 n 创建等间距向量 t，该向量以 0 开头，以 n-1 结尾，元素之间的间距为 1。

Task 1

load Cchord.mat;
n = numel(y);
t = 0:n-1;

t 现在有正确的点数，但它需要表示音频信号的采样时间。您可以使用采样频率 fs 将向量转换为时间（以秒为单位）。

任务

将 t 除以 fs。将输出赋给同一个变量 t。然后绘制 y 对 t 的图。

Task 2

t = t / fs;
plot(t,y)

在绘图中，请注意 y 是周期性的，但它不是简单的正弦波。它由具有不同频率的多个正弦波组成。

傅里叶变换将返回信号的频谱信息。主频的位置将显示和弦中包含的音符。

您可以使用 fft 函数来计算向量的离散傅里叶变换。

fft(y)

fft 的输出值为复数。您可以使用 abs 函数来获得幅值。

任务

创建一个名为 yfft 的变量，表示 y 的离散傅里叶变换的绝对值。

Task 3

yfft = abs(fft(y))

在任务 1 和 2 中，您已计算了信号 y 的时间向量 t。同样，您需要为您的 FFT 向量 yfft 计算频率向量 f。

任务

创建等间距向量 f，该向量以 0 开头，以 n-1 结尾，元素之间的间距为 1。

Task 4

f = 0:n-1

向量 f 现在包含 n 个点。要将这些点转换为频率，您可以将整个向量乘以采样频率 (fs)，然后除以点数 (n)。

f 将包含从 0 到 fs 的频率。主频位于 f 的开头位置。您可以使用 xlim 函数放大所关注的区域。

xlim([xmin xmax])

任务

将 f 乘以 fs/n。将输出赋给同一个变量 f。

使用 x 的限值 0 和 1000 绘制 yfft 对 f 的图。

Task 5

f = f*fs/n;
plot(f,yfft);
xlim([0 1000])