双音多频通信设计Matlab仿真.docVIP

双音多频（DTMF）通信设计的MATLAB仿真 摘要：讨论以MATLAB作为仿真工具产生DTMF信号，并用FFT算法、DFT算法、卷积法及迭代法来对DTMF信号进行解码。 关键词：FFT；DFT；频谱分析；卷积；滤波；差分方程；MATLAB 0 引 言 双音多频(DTMF:Double Tone Multi-Frequency)是按键电话通信，也广泛应用于电子邮件和银行系统中。用户可从电话发送DTMF信号来选择菜单进行操作。DTMF信号容易用软件产生和解码。MATLAB是一个高度集成的软件系统，通过交互式的命令（语句）可以十分简便地实现许多复杂的数值计算。 本文采用MATLAB作为仿真工具产生DTMF信号，并用FFT算法、DFT算法、卷积法及迭代法来对DTMF信号进行解码，由此得出：时域和频域是研究信号的两个窗口，其中信号处理大都在时域中进行，而信号分析往往在频域中比较方便直观。且数字信号处理技术中的DFT、FFT、卷积、滤波、差分方程这几个概念之间有内在联系。 1 DTMF信号的产生 DTMF是数字音频信号，在DTMF通信系统中共有8个频率，分为4个高频音和4个低频音，用一个高频音和一个低频音的组合表示一个信号，这样共有16种组合，分别代表16种信号，如表1所示： fH(Hz) fL(Hz) 1209 1336 1477 1633 697 1 2 3 A 770 4 5 6 B 852 7 8 9 C 941 * 0 # D 表1 DTMF信号组合表 例如，当按下数字键“1”时，则产生低频697Hz和高频1209Hz这两个正弦信号的迭加。由于语音信号的最高频率为4KHz,根据奈奎斯特取样定理，取样频率fs应大于或等于原信号最高频率fc的两倍，即 fs≥2fc (1) 才能保证取样后的信号不失真，所以电话音频信号在数字信号处理时，取样频率fs为 2×4k=8kHz,这里，每个数字信号持续时间为100ms,后面加上100ms的间隔时间（用0表示）。上述DTMF信号产生方法如下： 建立拨号数字的表矩阵，用查表法（查表1）求用户所按数字键对应的高、低频音。 为简化起见，仅允许选择“0-9”这十个键，在开始时还可拨空信号。 产生相应的DTMF信号及间隔时间。 由于fs=8kHz,各信号持续时间为100ms,因此在程序中每个信号取800点，间隔时间也取800点，结果存入数组中。 （3）画图并监听产生的DTMF信号。 程序如下： %[程序1]， DTMF信号的产生 clear %清除内存 TAB=[941 1336;697 1209;697 1336;697 1477;770 1209;770 1336;770 1477;852 1209;852 1336;852 1477];%拨号数字表矩阵 n=input(n=); %DTMF信号的个数 l=input(n0=);%空信号点数 for i=1:n k=input(0~9);%输入的数字键 fL=TAB(k+1,1); %k对应的低频音 fH=TAB(k+1,2); %k对应的高频音 n1=800;fs=8000;%产生k对应的DTMF信号，取样频率8kHz,每个信号共1600点，其中， j=0:1:n1-1; 前800点为信号持续时间，后800点为间隔时间，结果存入数组out中。 x=sin(2*pi*fL*j/fs)+sin(2*pi*fH*j/fs); out(1600*(i-1)+1+l:1600*i-800+l)=x; out(1600*i-799+l:1600*i+l)=0; end out=out./2; %out数组中每个数据除以2 subplot(211); plot(out); %绘图并监听DTMF信号。 sound(out,fs) wavwrite(out,fs,D2.wav);%out另存入声音文件 程序结果： 例如：若要产生数字键“1”所对应的DTMF信号，且开始时有200点空信号，则程序运行如下： 询问： 输入： n= 1 n0= 200 1 得图1： 图1 DTMF信号图 图1中，前200点为空信号，第201-1000点为键“1”对应的DTMF信号，第1001-1800点为间隔时间。 另外，在程序中为使