信号与系统（2）---录制和播放采样率关系对于声音信号的影响

本博文是“信号与系统”系列博文的第二篇，主要依托于“信号与系统”课程设计，后续还会继续推出第三篇，希望大家持续关注！！

---

前言

上篇博文介绍了采样率对于声音信号采集的影响，并且从声音信号播放和时频域波形分析两个角度阐释了采样率对于声音信号采集具体会产生哪些影响。本博文继续以“采样率”为抓手，研究采样率的另外一个问题---“录制和播放采样率关系对于声音信号的影响”。

---

下面开始我们的正文部分，由于在上篇的博文中已经介绍过采样的相关知识以及“奈奎斯特采样定理”，因此本篇博文直接从代码入手研究“录制和播放采样率关系对于声音信号影响”。

一、特定采样率声音信号录制

研究“录制和播放采样率关系对声音信号影响”前提条件是以特定采样率录制声音信号。

recorder = audiorecorder;

Fs2 = 44100;              %设置特定的采样率为44100Hz（生活中常见的采样频率） 
nBits = 16;
NumChannels = 1;

recorder = audiorecorder(Fs2,nBits,NumChannels);  % 调用电脑麦克风，开始声音录制


disp('strat speaking:');
recordblocking(recorder,7);   % 录制7s的声音信号
disp('stop speaking');

二、特定采样率声音信号播放

前面代码实现以特定采样率 Fs = 44100Hz 采集一段时长为7秒的语音信号，接下来实现以特定的采样率播放该录制的声音文件。项目研究的目的在于探讨录制和播放采样率对于声音信号的影响，因此在播放采样率设置时设置“一半”“两倍”“相同”三种情况进行研究。

myrecording = getaudiodata(recorder);         %getaudiodata--读取声音信号信息
audiowrite('data_record.wav',myrecording,44100)
info = audioinfo('data_record.wav')
audiowrite('data_record1.wav',myrecording,22050)
info1 = audioinfo('data_record1.wav')
audiowrite('data_record2.wav',myrecording,88200)
info2 = audioinfo('data_record2.wav')

% sound(myrecording,44100) % 播放时采样率设置为与录制时相同
% sound(myrecording,22050) % 播放时采样率设置为录制时的一半
% sound(myrecording,88200) % 播放时采样率设置为录制时的二倍

3.结果分析

针对于结果分析板块，与上篇博文相同，从两个不同的角度进行结果的分析，首先便是从声音信号播放的角度直观感受录制和播放采样率关系对于声音信号的影响，另外便是借助时频域分析更加清晰深刻的分析其对声音信号的影响。

（1）声音信号播放角度

在matlab中，借助于sound命令即可以特定的频率播放声音信号，本项目中分别设置“相同”“一半”“二倍”三种情况分析采样率关系对于声音信号的影响。
声音信号录制的采样率为 Fs = 44100Hz，如果播放时采样率为44100Hz，即录制和播放时采样率相同，此时，播放录音文件，听到的为正常声音信号，清晰洪亮。
但此时如果将播放时的采样率设置为Fs = 22050Hz，即录制和播放采样率关系为“一半”，此时声音信号播放速度变缓，声音变得更加低沉；此时如果将播放时的采样频率设置为Fs = 88200Hz，即录制和播放采样率的关系为二倍，此时声音信号播放速度加快，声音音调变高，略显刺耳。

（2）时频域分析角度

上面是从声音信号播放的角度借助于人耳直观感受录制和播放采样率关系对于声音信号的影响，接下来借助于时频域波形进行可视化分析，以另外的方式阐释其对于声音信号的影响。 时频域可视化代码：

% 数据可视化代码
[y,f] = audioread('data_record.wav');
t = 0:seconds(1/f):seconds(info.Duration);
t = t(1:end-1);
subplot(1,3,1)
plot(t,y)
title('采样频率相同')

[y1,f1] = audioread('data_record1.wav');
t1 = 0:seconds(1/f1):seconds(info1.Duration);
t1 = t1(1:end-1);
subplot(1,3,2)
plot(t1,y1)
title('采样频率变为“一半”')

[y2,f2] = audioread('data_record2.wav');
t2 = 0:seconds(1/f2):seconds(info2.Duration);
t2 = t2(1:end-1);
subplot(1,3,3)
plot(t2,y2)
title('采样频率变为“二倍”')

录制和播放采样率关系对声音信号的影响
通过时频域波形图可以明显看出录制和播放采样率不同对于音频文件长度有影响，换句话说，录制和播放采样率关系决定音频文件的播放速度。
通过对时频域波形的观察，明显可以看到当播放时采样率设置为录制时的一半时，声音文件的音频长度变为源文件二倍，当播放时采样率设置为录制时的二倍时，声音文件的音频长度变为源文件一半，但是如果观察三个音频文件的时域波形，发现波形完全一致，通过上述可以得到，录制和播放采样率关系对于音频文件内容没有影响，改变的为音频文件的长度，换句话说改变音频文件播播放的速度。

---

4.现象原理说明

为方便理解项目现象的原理，在这里对采样率的定义做一下补充，所谓“采样率”，即 每秒从连续信号中提取并组成离散信号的采样个数 。
声音文件录制时，以特定的采样率录制音频文件，只要满足“奈奎斯特采样定理”就可以采集到完整，丰富的声音信号。
而在声音信号播放时，按照自己的理解，相当于在已有音频文件的基础上进行二次采样，由于声音文件的内容是固定的，即音频文件样本数目是固定的，当我们以之前一半的采样率进行二次采样时，相当于每秒钟采集到的样本数目是减半的，因此，如果采集同样大小的音频样本，所需要的时间变为之前二倍，但是由于音频内容未发生改变，所以仅仅改变的是音频文件的播放速度，实现0.5倍速播放或者2倍速播放。

总结

本博文依托“信号与系统”课程设计，承接上篇文章，继续针对采样率展开研究探讨，通过matlab编码以及时频域分析阐释“录制和播放采样率对于声音信号的影响”，论证倍速播放的实现原理，并借此加深对于采样率的理解。
本博文依托课程设计，如有偏驳，敬请指正，共同学习！！

项目已上传至GitHub，GitHub地址：
https://github.com/booue/signals-processing