DTMF产生传输及检测软件仿真课程设计

1、袃燕山大学课程设计指导书膂 DTMF通信系统设计与仿真薁燕山大学信息科学与工程学院膇通信电子教研室羃 1、任务、性质与意义薂用数学模型的方法对信息传输的过程进行仿真，是分析、设计及改进通信系统性能的重要手段。本次课程设计的任务是综合运用通信原理、数字信号处理等课程所学知识，实现一个双音多频(DTMF)信号传输过程的软件仿真。罿本次课程设计的性质是通信原理、数字信号处理等课的综合实践环节。它的意义在于，通过本次课程设计可以从理论上加深对有关通信、数字信号处理原理的综合理解，了解DTMF信号原理和用数字信号处理的方法实现 DTMF信号的产生、传输与检测的过程。还可以从实践能力上掌握一门新的计算机编

2、程工具MATLAB语言。羅 2、基本原理肃 2.1 双音多频(DTMF)信号羃 双音多频 (DTMF)信号就是用两个不同频率的音频信号的组合来传输信息的一种通信信号形式，可以用数学式表示为螁s(t ) A(cos 2 f H t cos2 f L t) （ 1）羈式中 f H和 f L 分别是高频组和低频组的两个信号频率。对这两个频率的要求是：1、位于人而能够听到的音频频率范围内；2、频率间隔近似为乐音音阶的频率间隔（按十二平均律，相邻各个音的频率呈等比关系，其比值约为 1.059463 ，例如 C3130.8Hz，C4=261.7Hz，D4293.7 ）。这种信号在电话线路上传输，而且声音悦

3、耳。在电话系统中，高低频组的频率已经标准化，如图1 所示。肀 DTMF最早应用是在电话通信系统中作为控制信令使用的。在电话通信系统中，最早的控制信令是使用直流脉冲信号。 脉冲拨号方式是由接在旋转拨盘位置上的开关或电子开关控制二线电路的通与断。每拨出一位数字，都有唯一一串电流脉冲相对应，图2 中给出了数字3的电流脉冲序列。每个脉冲周期通常为100ms，其中有40%的占空时间，在人工控制条件下，两个相继数字的时间间隔可以从0.5S 到数秒之间变化不等。由此可见，利用拨号脉冲拨号，发送一个10 位长的号码大约需要7 秒的时间。腿 图 2 脉冲拨号示意图螇 当采用DTMF拨号方式时，号码的每一个数字用

4、一对音频表示，话机中有8 个单音频，分为两组，通过拨号盘选拨号码时，各位数字由触键开关输入，其对应的某个频率对（高频和低频）同时传输，图3 所示为各频率对应位置。对DTMF信号所规定的指标为：传送速率为每秒10 个数字，即每个数字100ms。每个数字传送过程中，信号必须存在至少45ms，且不得多于55ms，100ms 里其余时间是静音。另外ITU( 国际电信联盟 ) 还规定， DTMF信号在每个频率点上允许有1.5%的偏移，任何超出给定3.5%的信号认为无效，拒绝承认接收。还规定，在最坏检测条件下，信噪比不得低于15dB。显然。 DTMF的拨号方式比双脉冲拨号方式可快至10 倍，明显缩短了拨号

5、时间。而且DTMF拨号方式发出的信号抗干扰能力大大高于脉冲拨号方式。因此，目前普遍使用DTMF拨号方式。节采用 DTMF信号，用话音频率发送数字，可以避免占用额外的信道，又比脉冲拨号方式节约时间。在通信、测量、控制、自动服务等领域有着广泛的应用。在基于电话的各种信息服务系统中，广泛使用了DTMF信号来传送按键操作信息。利用这种按键信息，人们可以直接通过电话查询所需要的信息以及进行各种远程控制。在实际系统中， 首次拨号由程控交换机识别，完成主叫与被叫之间的接续；二次乃至多次拨号的识别以及操作由用户系统自己完成。 DTMF信号还可以用来在话音信道上传送各种类型的控制指令，例如利用电话控制家用电器的

6、启停，传送远方的状态监测信息。目前DTMF信号已经不仅仅限于在电话系统中应用，在测量、控制、遥测遥控等各个领域都有应用。蒁 2.2DTMF信号的硬件产生与检测袁 通常是用硬件产生和检测DTMF信号。已经有多种专用的产生和检测DTMF信号的集成电路器件。具有代表性的 DTMF发送、接收器的型号为MT8860,MT8862,MT8870,MT8872等。薆 DTMF信号发送器电路原理如图3 所示，它主要包括：(1) 晶体振荡器外接晶体( 通常采用3.58MHz) 与片内电路构成振荡器、经分频产生参考信号。(2) 键控可变速率时钟产生电路它是一种可变分频比的分频器，通常由 n 级移位寄存器与键控反馈

7、逻辑单元组成。(3) 正弦波产生电路它出正弦波编码器与D A变换器构成。通常可变速时钟信号先经过5 位移位寄存器，产生组5 位移位代码，再由可编程逻辑阵列(PLA) 将其转换成二进制代码，送到D A 变换器转换成台阶形正弦波。显然台阶的宽度等于时钟频率的倒数，这样形成的正弦波频率必然对应于时钟的速率和按键的号码。(4) 混合电路将键盘所对应产生的行、列正弦波信号( 即低高群f L、f H单音 ) 相加，混合成双音信号而输出。(5) 附加功能单元，如含有单音抑制、输出控制( 禁止 ) 、双键向按无输出等控制电路。袂 DTMF接收器主要包括DTMF分组滤波器和DTMF译码器，其基本原理如图4 所示

8、。 DTMF接收信号先经高、低群带通滤波进行 f/fH区分，然后过零检测、比较，得到相应于DTMF的两路 f /fH信号输出。该两路信号经译码、LL锁存、缓冲，恢复成对应于16 种 DTMF信号音对的4 比特二进制码 (L 1 一 L4) 。荿蚆 2.3 双音多频信号软件产生与检测芃 近年来 DTMF也应用在交互式控制中，诸如语言菜单、语言邮件、电话银行和ATM终端等。将DTMF信令的产生与检测集成到任一含有数字信号处理器(DSP)的系统中，是一项较有价值的工程应用。肁 为了产生DTMF信号， DSP用软件产生两个正弦波叠加在一起后发送，软件具体实现方式可以有三种方法：莈 （ 1）用软件模拟D

9、SP的方式产生两个数字正弦波振荡器并把输出合成起来，建立所希望的双音频。螆 （ 2）通过产生的某频率的方波后再附加滤波器滤出二次以上的谐波，产生DTMF信号。蚄 （ 3）查表输出产生DTMF信号。蕿 软件检测DTMF信号主要有以下方法：肇 （ 1）利用离散傅立叶变换(DFT) 及其快速算法(FFT) ，求取 DTMF信号频谱的峰值点，进而利用峰值点的频率判断发送的数字。这种方法的缺点是计算量大，不利于实时处理。但是如果使用MATLAB软件仿真则可以用很简单的几条语句实现。袆（ 2）对于 DTMF信号的波形， 波形过零点数量决定了低频信号的频率，而极值点数量决定了高频信号的频率。因此通过软件分别

10、计算出高频频率和低频频率，从而实现DTMF的解码。袁 这里选择采样频率为8K，已知每一个信号的频率至少持续40ms 的时间，我们对DTMF采样 160 个点，即20ms，这样总会有一个20ms的采样点全部落在40ms的发送范围内，对采样的点进行6 次求和取平均，去掉噪声成分，然后分别计算过零点的数目a 和极值点的数目b，对应频率可以计算为：芀 既然可以求出 DTMF信号的一组频率值， 那么就可以得出对应电话号码值。 这种算法比较适合于单片机硬件实现，但是其抗噪声性能较差。袆（ 3）计算接收到的DTMF信号在 8 个既定频率的能量， 进而确定是否有DTMF信号到达以及收到的是哪一对信号，此外通过

11、计算接收信号的总能量来防止误判。采用Geortzel算法可进一步提高计算效率，Goertzel算法实质是一个两极点的IIR 滤波器， 对应每一个频点有一个匹配滤波器，在抽样时刻检测。这种方法理论上属于最佳接收方法，而且算法有利于用DSP实现。在本次课程设计中要求使用这种方法。羆 DTMF解码即是在输入信号中搜索出有效的行频和列频。芁 2.4 数字正弦波振荡器原理蚈 数字正弦波振荡器的功能是利用数字信号处理的方法产生正弦波信号cos( t) 的抽样序列cos( nT ) 。基本原理是：羈设数字正弦波振荡器要产生的波形为肆 y( n)cos( nT )u(nT )1 (e j nT e j nT

12、)u(nT ) （ 1）2蚂则其 Z 变换为莀 上式的分子与分母同时乘以z 2蚇 y(n) 可以看作是一个单位冲击函数激励了一个线性系统的冲击响应。这个线性系统的传递函数就是Y(z) 。由 Y(z) 可以写出该线性系统的差分方程：肆式中0Ts 2 f / f s 为数字角频率。肃 令输入 x(n) 为单位冲击函数，即袈则可以得到下面的递推方程蒆 当 n<0 时， y(n)=0膆 当 n=0 时， y(0)=0膀 当 n=1 时，y(1)2 cos(0) ( 0) cos(0)x(1) cos(0)y薀 当 n=2 时， y(2)2 cos( 0 ) y(1)y( 0)芅当 n>2

13、时， y(n) 2 cos(0 ) y(n 1)y(n 2) (2)芅 实现时，首先将每个频率常数存在一个表中，用来初始化给定键的振荡器；再按照(2) 进行迭代计算，就得到正弦序列的输出；将两个这样的输出叠加即可得到双音频的输出。每个数字输出的样本数取决于数字速率与采样速率。双音频中的每一个音调之后都是同样长度的无音周期，从而可以检测到按键的释放。薁 2.5 数字匹配滤波器原理（ Geortzel 算法）肈计算接收到的DTMF信号在8 个既定频率的能量，进而确定是否有DTMF信号到达以及收到的是哪一对信号，此外通过计算接收信号的总能量来防止误判。为满足检测8 个DTMF频点（偏差1.5%），需

14、计算256 点FFT，因只对 8 个频点感兴趣，故可以直接计算8 个频点附近的DFTN1芈 X (k)x(n) e j 2 kn / N , k k1, k2 ,k3 , k4 , k5 , k6 ,k 7 , k8n0莅 计算数字信号的频谱可以采用DFT及其快速算法FFT，而在实现DTMF解码时，采用Goertzel算法要比FFT更快。通过FFT 可以计算得到信号所有谱线，了解信号整个频域信息，而对于DTMF信号只用关心其8 个行频 /列频及其二次谐波信息即可（二次谐波的信息用于将DTMF信号与声音信号区别开）。此时 Goertzel算法能更加快速地在输入信号中提取频谱信息。Goertzel

15、算法实质是一个两极点的IIR 滤波器，其算法原理框图如图5。羂图 5Geortzel算法信号流图螀 对应每一个频点有一个匹配滤波器，第k 个频率匹配滤波器的传递函数为肇其中蒅 H k1( z)112 ， H k 2 ( z)1 WNk z 1 ，2zcos kz1f kkfkk2 kN222， WN ejcos k j sin k莃kffsNN膈对应的前、后向差分方程为：螆 后一个方程无需全部计算，只需要在n=N 时刻，即最后的输出时刻计算。薅 如果考虑同步相位误差，也可以在n=N时刻附近计算一段时间内的输出。螄前向差分方程也可以改写为便于计算的递归形式：罿 并令初始值为vk ( n)0,n0

16、 。在 n=N时刻既可以得到所需要的X(k) 。衿 由于在 DTMF检测中，输入的信号是实数序列，并不需要检测出8 个行频 / 列频的相位，只需要计算出其幅度平方即可。因此计算|X(k)|2 如下：蚅 这里 vk ( N ), vk ( N1) 分别表示状态变量在N 和 N-1 时刻的值。羀 在输入信号中检测DTMF信号，并将其转换为实际的数字，这一解码过程本质是连续的过程，需要在输入的数据信号流中连续地搜索 DTMF信号频谱的存在。 整个检测过程分两步： 首先采用 Goertzel 算法在输入信号中提取频谱信息；接着作检测结果的有效性检查。蚁 2.6 仿真系统框图莁 DTMF螂 高聿DTMF

17、信 号 的斯噪信 号 检产生声莆电话信道测螅判断输出薇螁 图 6 仿真系统框图膅软件流程图（ 1）（ 2） 袅输入键号（ 3）（ 4） 膃 生成 DTMF信号（时域、频域显示）（ 5）（ 6） 艿 产生高斯噪声（时域、频域显示）（ 7）（ 8） 膈 信道传输（高斯白噪声信道，可选作带限滤波器信道）（ 9）（ 10）羄接收端随机延迟（ 30 个点以内）（ 11）（ 12）芀匹配滤波接收（ 13）（ 14）羁 DTMF信号到达检测（ 15）（ 16）羇 DTMF信号识别（时域、频域显示）（ 17）（ 18）肄 DTMF键号识别蚁 3. 工作内容与思考题葿 3.1 工作内容蚆 1、设计计算正弦波数字

18、振荡器、数字匹配滤波器的参数；膄 2、用 MATLAB语言设计一个利用 DTMF信号传输学号代码的仿真系统。观察下列信号的波形：学号代码的非归零基带信号、 DTMF信号、信道加入噪声、接收信号、匹配滤波器输出、译码判决结果；肂 3、利用仿真的系统进行抗干扰性能等方面的试验分析，试验分析的具体内容详见思考题。膁 3.2 思考题蝿 1 分析 DTMF信号的频谱，为什么 DTMF信号能够在电话信道传输，要求的带宽为多少？芄 2 设信噪比分别为 20dB，10dB，0dB，-6dB ，-10dB ， -20dB ，分析在理想（信道无失真）情况下，仿真计算误码率与信噪比的关系，并与理论分析结果比较。蒃

19、3 仿真分析匹配滤波器中心频率误差对DTMF信号检测有什么影响？薈 4 仿真分析 DTMF信号的时间长度对信号的检测有什么影响？蒈 5 如果在匹配滤波器前插入 FIR 带通滤波器，通频带在 600-1800Hz 之间。用窗函数法设计 FIR带通滤波器，并分析插入带通滤波器对检测性能有什么影响？芄 4 设计要求袃 1、自学复习通信原理、数字信号处理的有关内容，掌握算法原理；莀 2、掌握 MATLAB语言的使用；芆 3、完整实现仿真系统，发送自己的学号，接收自己的学号；莄 4、记录观察到的各种波形；羀 5、撰写课程设计报告。螈 报告统一使用B5 纸，统一报告封面， 报告篇幅以3000 字为宜。 报

20、告中要有设计目的；设计基本内容及步骤；每一步中应用的原理方法；试验仿真结果及相关波形；对试验结果进行的性能分析。报告正文用手写，图、程序可以打印作为附件。肅 5编程工具 MATLAB简介（亦可以采用 C语言、 C+、VC、 VB等编程语言实现）蒄 5.1 基本运算与函数莁 在 MATLAB下进行基本数学运算，只需将运算式直接打入提示号（>>）之后，并按入Enter 键即可。例如：蒀 >>(5*2+1.3-0.8)*10/25肈 ans=薃 4.2000螂 MATLAB会将运算结果直接存入一变量ans，代表 MATLAB运算后的答案 （ Answer），并显示其数值于屏幕

21、上。 （为简便起见， 在下述各例中， 我们不再印出MATLAB的提示号。）我们也可将上述运算式的结果设定给另一个变量x：羈 x=(5*2+1.3-0.8)*102/25袇 x=蚃 42芃 此时 MATLAB会直接显示x 的值。若不想让MATLAB每次都显示运算结果，只需在运算式最後加上分号（；）即可。蚀 由上例可知，MATLAB认识所有一般常用到的加（+）、减（ - ）、乘（ * ）、除（ / ）的数学运算符号，以及幂次运算（ ）。MATLAB将所有变量均存成double 的形式，所以不需经过变量声明。MATLAB同时也会自动进行存储器的使用和回收，而不必像C 语言，必须由使用者一一指定。蚆变

22、量也可用来存放向量或矩阵，并进行各种运算，如下例的行向量（Rowvector ）运算：螃 x=1352;蚄 y=2*x+1膈 y=虿 37115袃 5.2 控制命令螁 5.2.1 for 语句，其基本形式为：袀for 变量 =矩阵；蒈运算式；袃end膂 其中变量的值会被依次设定为矩阵的每一行，来执行介于for和 end 之间的运算式。因此，若无意外情况，运算式执行的次数会等于矩阵的行数。薁 举例来说，下列命令会产生一个长度为6 的调和数列：膇 x=zeros(1,6);%x 是一个 16 的零矩阵羃 fori=1:6,薂 x(i)=1/i;罿 end羅 5.2.2while 语句肃基本形式为：

23、羃while条件式；螁运算式；羈end膃也就是说，只要条件示成立，运算式就会一再被执行。例如先前产生调和数列的例子，我们可用圈改写如下：肀 x=zeros(1,6);%x 是一个 16 的零矩阵腿 i=1;螇 whilei<=6,节 x(i)=1/i;蒁 i=i+1;袁 end薆 5.2.3 条件语句 if,.,end薆基本形式为：袂 if 条件式；荿运算式；蕿 end蚆例：芃 ifrand(1,1)>0.5,肁 disp('Givenrandomnumberisgreaterthan0.5.');莈 end螆 5.2.4 MATLAB文件蚄若要一次执行大量的MAT

24、LAB命令，可将这些命令存放于一个扩展名为m的文件。包含MATLAB命令的文件都以m为扩展名，因此通称m文件，只要直接键入文件名，就可执行其所包含的命令。严格来说，m文件可再细分为命令集（Scripts）及函数（Functions）。蕿命令集的效果和将命令逐一输入完全一样，因此若在命令集可以直接使用工作空间的变量，而且在命令集中设定的变量，也都在工作空间中看得到。肇 函数则需要用到输入参数和输出参数来传递变量，这就像是C 语言的函数。在调用函数时，你只能经由输入参数来控制函数的输入，经由输出参数来得到函数的输出，但所有的暂时变量都会随着函数的结束而消失。袆 5.3 MATLAB常用的数学函数袁

25、 abs(x) ：标量的绝对值或向量的长度芀 angle(z) ：复数 z 的相角 (Phaseangle)袆 sqrt(x) ：开平方羆 real(z) ：复数 z 的实部芁 imag(z) ：复数 z 的虚部蚈 conj(z) ：复数 z 的共轭复数羈 round(x) ：四舍五入至最近整数肆 fix(x) ：无论正负，舍去小数至最近整数蚂 sign(x)：符号函数。当x<0 时， sign(x)=-1；当 x=0 时， sign(x)=0;当 x>0 时， sign(x)=1。莀 rem(x,y) ：求 x 除以 y 的余数蚇 gcd(x,y) ：整数 x 和 y 的最大公因

26、数肆 lcm(x,y) ：整数 x 和 y 的最小公倍数肃 exp(x) ：自然指数袈 pow2(x) ： 2 的指数蒆 log(x)：以 e 为底的对数，即自然对数或膆 log2(x)：以 2 为底的对数膀 log10(x) ：以 10 为底的对数薀 sin(x) ：正弦函数芅 cos(x) ：余弦函数芅 tan(x) ：正切函数薁 asin(x) ：反正弦函数肈 acos(x) ：反馀弦函数芈 atan(x) ：反正切函数莅 min(x): 向量 x 的元素的最小值羂 max(x): 向量 x 的元素的最大值螀 sort(x):对向量 x 的元素进行排序（Sorting）肇 length(

27、x): 向量 x 的元素个数蒅 sum(x): 向量 x 的元素总和莃 cumsum(x): 向量 x 的累计元素总和膈 dot(x,y): 向量 x 和 y 的内积螆 cross(x,y): 向量 x 和 y 的外积薅 5.4MATLAB绘图螄 MATLAB不但擅长于矩阵相关的数值运算，还适合于各种科学绘图。这里简介几种基本平面绘图命令。罿 plot是绘制一维曲线的基本函数，但在使用此函数之前，需要先定义曲线上每一点的x 及 y 座标。下例可画出一条正弦曲线：衿 closeall;x=linspace(0,2*pi,100);%100个点的 x 座标蚅 y=sin(x);% 对应的 y 座标

28、羀 plot(x,y);蚁 MATLAB基本绘图函数薇 plot:x轴和 y 轴均为线性刻度（Linearscale）蚅 loglog:x轴和 y 轴均为对数刻度（Logarithmicscale）莁 semilogx:x 轴为对数刻度， y 轴为线性刻度聿 semilogy:x 轴为线性刻度， y 轴为对数刻度莆 若要画出多条曲线，只需将座标对依次放入plot函数即可：螅 plot(x,sin(x),x,cos(x);螂若要改变颜色，在座标对後面加上相关字串即可：螁 plot(x,sin(x),'c',x,cos(x),'g');膅 若要同时改变颜色及图线型态

29、（Linestyle），也是在座标对後面加上相关字串即可：袅 plot(x,sin(x),'co',x,cos(x),'g*');膃 plot 绘图函数的叁数艿字符颜色字符图线型态膈y黄色.点羄k黑色o圆芀w白色xx羁b蓝色+羇肄蚁葿g绿色*r红色-实线c亮青色:点线m锰紫色-.点虚线蚆 用 subplot来同时画出数个小图形于同一个视窗之中：例如：膄 subplot(2,2,1);plot(x,sin(x);肂 subplot(2,2,2);plot(x,cos(x);膁 5.5MATLAB帮助蝿 MATLAB功能强大，命令众多，很难记全，但是MATLAB提供

30、了友好的帮助界面。因此利用MATLAB提供的帮助可以方便使用，边用边学。 MATLAB帮助界面有两种形式：命令行帮助和网页浏览帮助。这里仅简单介绍命令行帮助。在 MATLAB命令（ command）窗口中键入 helpXXX， MATLAB就会给出相应命令的功能、语法解释。芄例如，在MATLAB命令（ command）窗口中键入helpplot蒃 MATLAB就会给出如下提示：薈 PLOTLinearplot.蒈PLOT(X,Y)plotsvectorYversusvectorX.IfXorYisamatrix,thenthevectorisplottedversustherowsorcolu

31、mnsofthematrix,whicheverlineup.IfXisascalarandYisavector,length(Y)disconnectedpointsareplotted.芄 PLOT(Y)plotsthecolumnsofYversustheirindex.袃 IfYiscomplex,PLOT(Y)isequivalenttoPLOT(real(Y),imag(Y).莀InallotherusesofPLOT,theimaginarypartisignored.Variouslinetypes,plotsymbolsandcolorsmaybeobtainedwithPL

32、OT(X,Y,S)whereSisacharacterstringmadefromoneelementfromanyorallthefollowing3columns:芆 bblue.point-solid莄 ggreenocircle:dotted羀 rredxx-mark-.dashdot螈 ccyan+plus-dashed肅 mmagenta*star蒄 yyellowssquare莁 kblackddiamond蒀 vtriangle(down)肈triangle(up)薃<triangle(left)螂>triangle(right)羈 ppentagram袇 hhex

33、agram蚃Forexample,PLOT(X,Y,'c+:')plotsacyandottedlinewithaplusateachdatapoint;PLOT(X,Y,'bd')plotsbluediamondateachdatapointbutdoesnotdrawanyline.芃PLOT(X1,Y1,S1,X2,Y2,S2,X3,Y3,S3,.)combinestheplotsdefinedbythe(X,Y,S)triples,wheretheX'sandY'sarevectorsormatricesandtheS'sares

34、trings.蚀Forexample,PLOT(X,Y,'y-',X,Y,'go')plotsthedatatwice,withasolidyellowlineinterpolatinggreencirclesatthedatapoints.ThePLOTcommand,ifnocolorisspecified,makesautomaticuseofthecolorsspecifiedbytheaxesColorOrderproperty.ThedefaultColorOrderislistedinthetableaboveforcolorsystemswher

35、ethedefaultisblueforoneline,andformultiplelines,tocyclethroughthefirstsixcolorsinthetable.Formonochromesystems,PLOTcyclesovertheaxesLineStyleOrderproperty.蚆PLOTreturnsacolumnvectorofhandlestoLINEobjects,onehandleperline.TheX,Ypairs,orX,Y,Striples,canbefollowedbyparameter/valuepairstospecifyadditiona

36、lpropertiesofthelines.螃 6. 程序与波形示例蚄部分程序代码如下：膈 %Generalparameter虿 clearall;袃 t_symbol=0.1;螁 f_sample=8000;薆 t_sample=1/8000;膄 frequency_LL=697770852941;虿 frequency_HH=1209133614771633;袈 symbols='1''2''3''a''4''5''6''b''7''8&#

37、39;'9''c''#''0''*''d'芇 ID_string=input('pleaseinputyourIdenticalString','s')羂 SNR_dB=input('pleaseinputSNRindB');蝿 SNR=power(10,SNR_dB/10);芈 signal_amplitude=1;螅 power_of_signal=1;蚁 sigma=sqrt(1/SNR);蝿 LL=length(ID_string);虿 tot

38、al_ss=;蒇 fornn=1:LL螄 row,col=identify_ID(ID_string(nn);衿 ss1,ss2=generate_DTMF_signal(row,col);袆 ss=ss1+ss2;袅 total_ss=total_ssss;蒃 end羈 figure芇 noise=sigma*randn(1,length(total_ss);蚇 zz=total_ss+noise;节 subplot(2,1,1);莂 plot(zz);蚈 subplot(2,1,2);肅 plot(abs(fft(zz);芅 zzz=zz;蒂 %detectsignalstartingpo

39、sition聿 NN1=t_symbol/t_sample;螇 NN2=NN1/2;肄 NN=length(zzz)/NN1;蒂 vector=1,NN1,-500,500;蒀 %符号检测芄 fornn=0:NN-1% 对 N个符号检测袃 figure薂 fork=1:4% 对低频组 4 个匹配滤波器运算薇 ss(1)=1;羆 ss(2)=0;蚁 coef=2*pi*frequency_LL(k)/f_sample;蚂 b0=sin(coef);羇 a1=2*cos(coef);蒄 forn=3:NN1蚄 ss(n)=a1*ss(n-1)-ss(n-2)+b0*zzz(n+nn*NN1);螂 end莈 %匹配滤波器输出信号的模值膆 XX(k)=0;蒃 fortt=1:5,袂 XX(k)=XX(k)+ss(NN2-tt)*ss(NN2-tt);蝿 end薄 XX(k)膂 subplot(4,2,k);羁 plot(ss,'r');袆 AXIS(vector);芆 end羁 fork=1:4% 对高频组 4 个匹配滤波器运算