利用MATLAB仿真多普勒效应

1、无 利用利用 MATLABMATLAB 仿真仿真多普勒效应多普勒效应 某某某某某某 摘摘 要：要：分析多普勒效应特性，建立数学模型，利用 MATLAB 软件对其进行仿真试验，进行定量分析， 根据仿真试验结果绘制出听者接收到的信号的频率变化曲线以及用信号处理工具箱函数 spectrogram绘制的时间频率图，并生成相应的声音。 关键词关键词：Doppler effect MATLAB/ Simulink 0 0 研究背景研究背景 多普勒效应是由生在德国的奥地利物理学家多普勒(Johann Doppler 1802 一 1853)发现的。1845 年，荷兰气象学家巴依斯巴洛(h.d.Buys Ba

2、llot)测得了声音的多普勒效应。一辆汽车在我们身旁急驰而过， 车上喇叭的音调有一个从高到低的突然变化； 站在铁路旁边听火车的汽笛声也能够发现， 火车迅速迎面而来时音调较静止时为高， 而火车迅速离去时则音调较静止时为低。这是日常生活中的一个多普勒效应的例子。在天文、通信等领域还有众多的例子。 当波源或观察者相对于媒质运动时， 或者说波源和观察者有相对运动时， 观察者接受到的震动频率与波源震动频率不同的现象，称为多普勒效应。 对于多普勒效应的讨论， 一般仅限于声源和听者在同一直线上运动的情况。 当声源和听者不在同一直线上运动时， 接收频率变化比较复杂， 听者接收到的信号波形方程也难以用解析式表示

3、。 MATLA 具有强大的数值计算和仿真功能以及图形技术。本文试图从 MATLAB 编程的角度出发，应用 MATLAB 的 Simulink 仿真试验方法，建立仿真的试验环境，对声源和听者不在同一直线上运动的情况下产生的多普勒效应特性进行分析， 产生极好的模拟， 实现多普勒效应的验证，绘制出听者接收到的信号的频率变化曲线以及用信号处理工具箱函数 spectrogram绘制的时间频率图，并生成相应的声音。 1 1 基本原理基本原理 以下公式描述了多普勒效应现象的各个物理量之间的定量关系： fsvf/cos10 （1-1） 其中， f0 是声源发出的声音的频率; v是听者与声源的相对运动速度; 为

4、速度矢量与声源和听者的连线夹角; vs 为声音在空气中传播的速度，f 是听者听到的声音频率。 2 2 理论理论模型模型 多普勒效应-设声源距离听者的水平距离为0 x，以 v的速度沿水平方向向听者直线驶 无 来，其轨迹与听者的最小距离为垂直距离 y ， 声源的频率为 f0，试求听者接受到的信号的 频率变化曲线以及用信号处理工具箱函数 sp ectrogram绘制的时间频率图，并生成相应的 声音。 3 3 仿真试验仿真试验 下面是一段用 MATLAB 的 M 文件依据公式（1-1）原理编写的程序（1-1） 。它会产生一个名为 dopp.wav 的表现多普勒效应的声音文件， 描述火车向一个距离铁路

5、30 米 （垂直距离） 、距火车 150（水平距离）的听者开来时他听到的声音。将 u1 送入图 3-1 所示的 Simulink 仿真系统，再一次听到该声音，并且看到用频谱仪表现的声音频率随时间的变化的情况。表3-1 和表 3-2 给出了仿真系统中两个模块的主要参数设置。 表表 3 3- -1 From Multimedia File1 From Multimedia File（波形文件波形文件）的主要参数的主要参数 模块名称 From Multimedia File 位置 DSP System toolbox / Signal Processing Sources 参 数 名 称 参 数 值

6、 File name(文件名) Dopp.wav Samples per output frame（输出每帧中的样值数） 1024 表表 3 3- -2 2 To Audio DeviceTo Audio Device（扬声器）的主要参数（扬声器）的主要参数 模块名称 To Audio Device 位置 DSP System toolbox / Signal Processing Sinks 参 数 名 称 参 数 值 Device Default Queue duration(seconds)(队列延迟时间) 1 Automatically determine Buffer size（自动

7、确定缓冲区长度） 选中 图图 3 3- -1 1 显示声音信号（多普勒效应）的仿真系统模型显示声音信号（多普勒效应）的仿真系统模型 无 程序程序 1 1- -1 1 % 多普勒效应程序 程序程序难点难点分析分析 1 在计算机中所能存储的数字序列，也即模拟信号必须通过取样和量化后，变成相应的数字信号，才能被计算机存储和处理。对音频信号来说，实现模拟音频信号与数字音频信号之间的转化模块就是声卡， MATLAB 可以方便地对声卡进行诸如采样频率等输入/输出参数的配置。 2 播放火车静止时的汽笛声。使用命令“sound”将数字序列以设定的采样速率输出到声卡，通过声卡转化为模拟音频信号。 “sound”

8、用法是： sound（y,Fs） ； 其中：y 为取值范围必须在-1，,1区间的 n 行 1 列的数字序列（单声道输出） 。Fs为设定的采用速率，一般声卡支持 5000Hz 到 441000Hz。 x0=150; %m v=50; %车速 m/s y0=30; %m c=330; %声音在空气中的速度 m/s f0=2000; %声音的频率 Hz fs=8000; %采样频率 t=0:1/fs:6; %规定 t 的范围和步进 x_t=-x0+v.*t; %火车与观察者的水平距离 y=y0; %火车与观察者的垂直距离 r=sqrt(x_t.2+y.2); %火车与观察者间距离 costheta=

9、x_t./r; %声速矢量与火车和听者间连线夹角的余弦 f=f0./(1+v.*costheta./c); %听者接收到的频率 fmin=min(f); fmax=max(f); voltage=(f-fmin)./(fmax-fmin)*2-1; %归一化调频电压在 -1 +1 之间 1 signal=0.5*vco(voltage,fmin fmax,fs); %压控震荡器，由输入电压控制输出信号频率 u=sin(2*pi*f0*t); %声源发出的信号 sound(u,fs); %播放声源发出的信号 2 pause(5); %暂停 3 sound(signal,fs); %播放听者接收到

10、的信号 4 wavwrite(signal,fs,dopp.wav); %音频数据存盘 5 figure(1); plot(t,f); xlabel(Time/s); ylabel(Freq/Hz); %做出听者接收到的信号的频率变化曲线 figure(2); spectrogram(signal,kaiser(256,5),220,512,fs,yaxis); %做出时间频率图 6 axis(0 6,1000,spectro3000) 无 3 此处暂停必不可少，且时间要足够，以便再打开声音系统，这个量与计算机硬件有关。 4 播放静止的听者听到运动的火车的汽笛声。 5 MATLAB 可 以 记

11、 录 的 音 频 信 号 直 接 保 存 为 “ *.wav ” 格 式 。 利 用 命 令“wavwrite(y,Fs,Filename);”就可以将向量 y 存储为取样频率为 Fs 的 wav 音频 文件。 6 另一个函数发生器是压控振荡器 （VCO） ， 它由输入向量来决定某个震荡频率的信号。此处，绘制生成信号的光谱图。 频域测量仪器频域测量仪器（频谱仪（频谱仪 Spectrum Scope） 计算机仿真中的频谱仪应用的是数字信号处理中的快速傅里叶变换（FFT）技术，它能完成数据流从时域到频域的变换。 将时域的数据流取出一段，即 FFT size(快速傅里叶变换的长度)确定为 N，通常要

12、求 N是 2 的幂次方。同时需要设置相应长度的 Buffer Size(缓存器)。N 的大小（即时窗的长短）决定了频谱仪的分辨率，N 越大，频率的分辨率就越高，但相应的计算时间也愈长。 希望所研究的谱线内容出现在频谱仪显示窗的中间部分， 能看到在频率轴上谱线的低端和高端的情况，以便于观察和分析。要做到这一点，将输入信号的采样频率取为期望的频率显示窗最大值的两倍即可。 注意，频谱仪的采样频率与被测信号的采样频率要一致。 4 4 仿真结果仿真结果及分析及分析 （1 1）仿真结果）仿真结果 图 4-1 所示：M 文件程序运行结束时绘出的听者接收到信号的频率变化曲线。 图 4-2 所示：用信号处理工具

13、箱函数 spectrogram 绘制的时间频率图。 图 4-3 所示：某个时刻仿真系统中的频谱仪上显示的功率谱。 图图 4 4- -1 1 听者接收到的信号的频率变化曲线听者接收到的信号的频率变化曲线 无 图图 4 4- -2 2 时间频率图（时间频率图（SweptSwept- -Frequency WaveformsFrequency Waveforms） T T=0s=0s 无 T T=3s=3s T T= =6 6s s 图图 4 4- -3 3 频谱仪显示频谱仪显示 无 结果结果分析分析 傅里叶变换是把各种形式的信号用正弦信号表示， 因此非正弦信号进行傅里叶变换， 会得到与原信号频率不

14、同的成分都是原信号频率的整数倍。 这些高频信号是用来修饰频率与原信号相同的正弦信号，使之趋近于原信号的。所以说，频谱上频率最低的一个峰（往往是幅度上最高的） ，就是原信号频率。 由上可知，听者接收到的信号的频率先高于原来的汽笛声，后低于原来的汽笛声。 5 5 结束语结束语 在 MATLAB/Simulink 下，从数学模型到计算机仿真模型的转换非常容易。本文利用Simulink 结合 M 文件编程的方法， 使其同时具备图形界面的直观性和字符界面的强大功能。事实上，所有 Simulink 的模块以及系统构建、仿真参数等均可通过编程语句实现。 针对声源和听者不在同一直线上运动情况， 利用 MATL