北大心理学系.利用Matlab生成听觉刺激.pdf

第五章 利用Matlab 生成听觉刺激 【内容提要】 1. 声音的特性； 2. 用Matlab 产生纯音、噪音、乐音和复合音（和音） 3. 构造声音的移动； 4. 声音的淡入和淡出； 5. 声音的其它处理（滤波和声场分析） 一、声音的特性 如前一章提到的，Psychtoolbox 是一组Matlab 函数（M 文件）和用C 语言写的可执行的mex 文件构成，这使得我们非常方便去定制心理学实验所需要的刺激事件。利用 Matlab 和 Psychtoolbox 工具包，我们也能方便地产生各种类型的声音（包括噪音），复杂和虚拟的声 场，并可以对声音的频谱特性进行修饰，以符合心理学实验的需要。 在介绍具体的方法之前，先对声音的基本特性要有清晰的认识。对于一个以正弦波或与余弦 波为基础的纯音，它包含四个要件：声音的时长、频率、采样率和比特率。赫兹是频率单位， 记为Hz，指每秒钟周期性变化的次数。分贝是用来表示声音强度的单位，记为dB。声音是 一种波，我们以极短的时间间隔把波形变成一系列的数字，也就是模拟信号到数字信号的转 换，简称A/D 转换，就是对声音进行采样，每秒钟采样的次数称为采样频率。声音的采样率 与声卡的缓存容量有关。常用的采样频率有 11.025kHz、22.05kHz 和44.1kHz。在采样时， 采样点之间的时间间隔越小，即每秒钟采样的次数越多，采样频率越大，采样就越细腻逼真， 所以理论上讲应该是采样频率越高音质越好。但由于人耳听觉分辨率毕竟有限，最大分辨率 大约是20kHz 左右。比特率则表示记录音频数据每秒钟所需要的平均比特值，比特是计算机 中最小的数据单位，指一个0 或者 1 的数，通常我们使用 Kbps （1024bps）作为单位。CD 中的数字音乐比特率为1411.2Kbps，也就是记录1 秒钟的CD 音乐，需要1411.2×1024 比 特的数据，对于大部分人来说，192Kbps 的MP3 音质已经非常不错了，但是如果你的耳朵比 较灵敏，自然要求也就相应提高。目前大部分MP3 播放机所能够支持的256Kbps 的比特率， 基本上已经能够满足要求了，如果能够支持更高当然更好。影响声音的大小的物理要素是振 幅，电脑上的声音必须也要能精确表示乐曲的轻响，所以一定要对声波的振幅有一个精确的 描述比特”就是这样一个单位，x 比特就是指把波形的振幅划为2 的x 次方个等级，根据模 拟信号的轻响把它划分到某个等级中去，就可以用数字来表示了。比特率越高，越能细致地 反映声音的轻响变化。为了体现正常的声音信息，16bit 为基本的需求，较好的cd 使用的 是20bit 甚至24bit。 图 1 4 比特的声音采样示意图 声音总可以被分解为不同频率不同强度正弦波的叠加。这种变换（或分解）的过程，称为傅 里叶变换。因此，一般的声音总是包含一定的频率范围。人耳可以听到的声音的频率范围在 20 到2 万赫兹 （Hz）之间。高于这个范围的波动称为超声波，而低于这一范围的称为次声 波。 音阶是按照特定规律，将音高顺次排列的一组音列。音阶的概念可以追溯至中古时代的素歌 的调式，用以突出乐曲的某种特质或音质，后来经过不断的演变，使音阶的概念逐步确立。 一般人较常认识的音阶有大调音阶和小调音阶。都是根据特定的排列方式所制造出来的，不 过，大、小二调亦可以通过把次序的先后改变，变成了新的音阶变种，而这类变种亦富有古 代调式的元素在内。各种较常见的音阶包括半音音阶、八声音阶、七声音阶、六声音阶、五 声音阶、爵士音阶等。 表征声音强度的一个单位为声压。特定介质下的声压是指是指声波通过某种媒质时，由振动 所产生的压强改变量，一般会考虑在不同时间或空间下，声压的均方根（RMS）为其平均值。 例如空气中声压均方根为1Pa （94dbSPL）的声音，表示其实际的压强会在（1 atm Pa） 及（1 atm Pa）之间变化，即在101323.6Pa 及101326.4 Pa 之间变化。若以压强的 观点来看，上述声压造成的压强变化很小，但若频率在声频内，此此音却是震耳欲聋，可能 会造成听力损害的程度。由于人耳可以传感的声音振幅范围较广，声压一般会表示为对数尺 度，以分贝表示的声压级（SPL）来表示。声压级（SPL）可以用L 表示．定义如下： 其中，p 为声压的均方根值， 为参考声压，一般用的参考声压是以ANSI S1.1