通信技术方向课程设计题目：双音多频信号检测
物联网工程学院电子信息工程专业
学号0703070106
学生姓名时雅茹
二〇一〇年六月
一、原理介绍
双音多频（Dual Tone Multi Frequency, DTMF ）信号是音频电话中的拨号信号，由美国AT&T 贝尔公司实验室研制，并用于电话网络中。这种信号制式具有很高的拨号速度，且容易自动监测识别，很快就代替了原有的用脉冲计数方式的拨号制式。这种双音多频信号制式不仅用在电话网络中，还可以用于传输十进制数据的其它通信系统中，用于电子邮件和银行系统中。这些系统中用户可以用电话发送DTMF 信号选择语音菜单进行操作。
DTMF 信号系统是一个典型的小型信号处理系统，它要用数字方法产生模拟信号并进行传输，其中还用到了D/A 变换器；在接收端用A/D 变换器将其转换成数字信号，并进行数字信号处理与识别。为了系统的检测速度并降低成本，还开发一种特殊的DFT 算法，称为戈泽尔(Goertzel)算法，这种算法既可以用硬件（专用芯片）实现，也可以用软件实现。下面首先介绍双音多频信号的产生方法和检测方法，包括戈泽尔算法，最后进行模拟实验。
二、内容及结论
1、双音多频（DTMF ）信号的组成
在电话中，数字0~9的中每一个都用两个不同的单音频传输，所用的8个频率分成高频带和低频带两组，低频带有四个频率：679Hz,770Hz,852Hz 和941Hz ；高频带也有四个频率：1209Hz,1336Hz,1477Hz 和1633Hz.。每一个数字均由高、低频带中各一个频率构成，例如1用697Hz 和1209Hz 两个频率，信号用
)2sin()2sin(21t f t f ππ+表示，其中Hz f 6791=，Hz f 12092=。这样8个频率形成
16种不同的双频信号。具体号码以及符号对应的频率如表1所示。表中最后一列在电话中暂时未用。
表1 双频拨号的频率分配
列 行 1209Hz
1336Hz
1477Hz 633Hz 697Hz 1 2 3 A 770Hz 4 5 6 B
852Hz 7 8 9 C 942Hz
*
#
D
DTMF 信号在电话中有两种作用，一个是用拨号信号去控制交换机接通被叫的用户电话机，另一个作用是控制电话机的各种动作，如播放留言、语音信箱等。
2、电话中的双音多频（DTMF ）信号的产生与检测 （1）双音多频信号的产生
假设时间连续的 DTMF 信号用)2sin()2sin()(21t f t f t x ππ+=表示，式中2
1f f 和是按照表1选择的两个频率，1f 代表低频带中的一个频率，2f 代表高频带中的一个频率。显然采用数字方法产生DTMF 信号，方便而且体积小。下面介绍采用数字方法产生DTMF 信号。规定用8KHz 对DTMF 信号进行采样，采样后得到时域离散信号为
)8000/2sin()8000/2sin()(21n f n f n x ππ+=
形成上面序列的方法有两种，即计算法和查表法。用计算法求正弦波的序列值容易，但实际中要占用一些计算时间，影响运行速度。查表法是预先将正弦波的各序列值计算出来，寄存在存储器中，运行时只要按顺序和一定的速度取出便可。这种方法要占用一定的存储空间，但是速度快。
因为采样频率是8000Hz ，因此要求每125ms 输出一个样本，得到的序列再送到D/A 变换器和平滑滤波器，输出便是连续时间的DTMF 信号。DTMF 信号通过电话线路送到交换机。
（2）基于Goertzel 算法的双音多频信号检测
在接收端，要对收到的双音多频信号进行检测，检测两个正弦波的频率是多少，以判断所对应的十进制数字或者符号。显然这里仍然要用数字方法进行检测，因此要将收到的时间连续 DTMF 信号经过A/D 变换，变成数字信号进行检测。检测的方法有两种，一种是用一组滤波器提取所关心的频率，根据有输出信号的2个滤波器判断相应的数字或符号。另一种是用DFT （FFT ）对双音多频信号进行频谱分析，由信号的幅度谱，判断信号的两个频率，最后确定相应的数字或符号。当检测的音频数目较少时，用滤波器组实现更合适。FFT 是DFT 的快速算法，但当DFT 的变换区间较小时，FFT 快速算法的效果并不明显，而且还要占用很多内存，因此不如直接用DFT 合适。下面介绍Goertzel 算法，这种算法的实质是直接计算DFT 的一种线性滤波方法。这里略去Goertzel 算法的介绍，可以直接调用MATLAB 信号处理工具箱中戈泽尔算法的函数Goertzel ，计算N 点DFT 的几个感兴趣的频点的值。 3、MATLAB 工具箱函数goertzel Goerztel 函数的调用格式为：
Xgk=goertzel(xn,K)
xn 是被变换的时域序列，用于DTMF 信号检测时，xn 就是DTMF 信号的205个采样值。
K 是要求计算的DFT[xn]的频点序号向量，用N 表示xn 的长度，则要求1≤K ≤N 。由表2可知，如果只计算DTMF 信号8个基频时，
K=[18，20，22，24，31，34，38，42]，
如果同时计算8个基频及其二次谐波时，
K=[18，20，22，24，31，34，35，38，39，42，43，47，61，67，74，
82]。
Xgk 是变换结果向量，其中存放的是由K 指定的频率点的DFT[x(n)]的值。设X(k)= DFT[x(n)]，则()(()), 1,2,,length()Xgk i X K i i K =
=。
4、检测DTMF 信号的DFT 参数选择
用DFT 检测模拟DTMF 信号所含有的两个音频频率，是一个用DFT 对模拟信号进行频谱分析的问题。根据第三章用DFT 对模拟信号进行谱分析的理论，确定三个参数：（1）采样频率s F ，（2）DFT 的变换点数N ，（3）需要对信号的观察时间的长度p T 。这三个参数不能随意选取，要根据对信号频谱分析的要求进行确定。这里对信号频谱分析也有三个要求： （1）频率分辨率，（2）谱分析的频谱范围，（3）检测频率的准确性。 （1）、频谱分析的分辨率。
观察要检测的8个频率，相邻间隔最小的是第一和第二个频率，间隔是73Hz ，要求DFT 最少能够分辨相隔73Hz 的两个频率，即要求Hz F 73min =。DFT 的分辨率和对信号的观察时间p T 有关，ms F T p 7.1373/1/1min === 。考虑到可靠性，留有富裕量，要求按键的时间大于40ms 。 （2）、频谱分析的频率范围
要检测的信号频率范围是697～1633Hz ，但考虑到存在语音干扰，除了检测这8个频率外，还要检测它们的二次倍频的幅度大小，波形正常且干扰小的正弦波的二次倍频是很小的，如果发现二次谐波很大，则不能确定这是DTMF 信号。这样频谱分析的频率范围为697～3266Hz 。按照采样定理，最高频率不能超过折叠频率，即Hz F s 36225.0≥，由此要求最小的采样频率应为7.24KHz 。因为数字电话总系统已经规定s F ＝8KHz ，因此对频谱分析范围的要求是一定满足的。按照ms T p 7.13m in =，s F ＝8KHz ，算出对信号最少的采样点数为
110min min ≈?=s p F T N 。
（3）、检测频率的准确性
这是一个用DFT 检测正弦波频率是否准确的问题。序列的N 点DFT 是对序列频谱函数在0～π2区间的N 点等间隔采样，如果是一个周期序列，截取周期序列的整数倍周期，进行DFT ，其采样点刚好在周期信号的频率上，DFT 的幅度最大处就是信号的准确频率。分析这些DTMF 信号，不可能经过采样得到周期序列，因此存在检测频率的准确性问题。
DFT 的频率采样点频率为N k k /2πω=（k =0,1,2,---,N -1），相应的模拟域采样点频率为N k F f s k /=（k =0,1,2,---,N -1），希望选择一个合适的N ，使用该公式算出的k f 能接近要检测的频率，或者用8个频率中的任一个频率
'k f 代入公式
'/k s f F k N =中时，得到的k 值最接近整数值，这样虽然用幅度最大点检测的
频率有误差，但可以准确判断所对应的DTMF 频率，即可以准确判断所对应的数字或符号。经过分析研究认为N ＝205是最好的。按照s F ＝8KHz ，N ＝205，算出8个频率及其二次谐波对应k 值，和k 取整数时的频率误差见表2。 表2 8个基频 Hz 最近的整数k 值 DFT 的 k 值 绝对误差
二次谐波 Hz
对应的 k 值
最近的 整数k 值
绝对误差
697 17.861 18 0.139 1394 35.024 35
0.024 770 19.531 20 0.269 1540 38.692 39 0.308 852 21.833 22 0.167 1704 42.813 43 0.187 941 24.113 24 0.113 1882 47.285 47 0.285 1209 30.981 31 0.019 2418 60.752 61 0.248 1336 34.235 34 0.235 2672 67.134 67 0.134 1477 37.848 38 0.152 2954 74.219 74 0.219 1633 41.846
42
0.154
3266
82.058 82
0.058
通过以上分析，确定s F ＝8KHz ，N ＝205，ms T p 40≥。
三、仿真程序运行结果
DTMF信号的产生与识别仿真实验在MATLAB环境下进行，编写仿真程序，运行程序，送入8位电话号码，程序自动产生每一位号码数字相应的DTMF信号，并送出双频声音，再用DFT进行谱分析，显示每一位号码数字的DTMF信号的DFT 幅度谱，安照幅度谱的最大值确定对应的频率，再安照频率确定每一位对应的号码数字，最后输出8位电话号码。
1、程序分四段：第一段（2—7行）设置参数，并读入8位电话号码；第二段（9—20行）根据键入的8位电话号码产生时域离散DTMF信号，并连续发出8位号码对应的双音频声音；第三段（22—25行）对时域离散DTMF信号进行频率检测，画出幅度谱；第四段（26—33行）根据幅度谱的两个峰值，分别查找并确定输入8位电话号码。
程序清单如下：
% DTMF双频拨号信号的生成和检测程序
%clear all;clc;
tm=[1,2,3,65;4,5,6,66;7,8,9,67;42,0,35,68]; % DTMF信号代表的16个数
N=205;K=[18,20,22,24,31,34,38,42];
f1=[697,770,852,941]; % 行频率向量
f2=[1209,1336,1477,1633]; % 列频率向量
TN=input('键入8位电话号码= '); % 输入8位数字
TNr=0; %接收端电话号码初值为零
for l=1:8;
d=fix(TN/10^(8-l));
TN=TN-d*10^(8-l);
for p=1:4;
for q=1:4;
if tm(p,q)==abs(d); break,end % 检测码相符的列号q end
if tm(p,q)==abs(d); break,end % 检测码相符的行号p
end
n=0:1023; % 为了发声，加长序列
x = sin(2*pi*n*f1(p)/8000) + sin(2*pi*n*f2(q)/8000);% 构成双频信号
sound(x,8000); % 发出声音
pause(0.1)
% 接收检测端的程序
X=goertzel(x(1:205),K+1); % 用Goertzel算法计算八点DFT 样本
val = abs(X); % 列出八点DFT向量
subplot(4,2,l);
stem(K,val,'.');grid;xlabel('k');ylabel('|X(k)|') % 画出DFT(k)幅度
axis([10 50 0 120])
limit = 80; %
for s=5:8;
if val(s) > limit, break, end % 查找列号
end
for r=1:4;
if val(r) > limit, break, end % 查找行号
end
TNr=TNr+tm(r,s-4)*10^(8-l);
end
disp('接收端检测到的号码为：') % 显示接收到的字符
disp(TNr)
2、运行结果：
（1）、运行程序，根据提示键入8位电话号码12345678，然后回车。
（2）、回车后可以听见8位电话号码对应的DTMF信号的声音，并输出相应的8幅频谱图如下图所示。
（3）、左上角的第一个图在k=18和k=31两点出现峰值，所以对应第一位号码数字1。最后显示检测到的电话号码12345678。
四、总结
这学期我们学习了通信原理课程这门课程并开设了实验，同时在短学期中进行了为期一周的课程设计，目的是让我们更深刻的理解通信原理这门课，并且把所学的理论与实践相联系。这次的课程设计共有个八题目，我完成的是第六个——双音多频信号检测。
我拿到题目后，上网查阅老师给的参考资料，将英文翻译成中文，了解了双音多频信号检测的原理。参考资料所给内容比较简单且Matlab的程序有缺省和错误，所以我又学习通信原理中相关知识并查阅一些资料，学习Matlab的使用方法，最终写出程序，进行了调试、仿真，并完成报告。
这次课程设计，我学会了很多东西。熟悉了仿真软件Matlab的使用方法，学会了Matlab程序的编写，完成仿真。我还意识到英文水平的重要性，我们学习电信专业，会需要了解国外先进的技术及最新，也需要阅读经典书籍和各种文献、资料，因此需要扎实的英语阅读和翻译功底，尤其是相关的专业词汇。此外，这次课程设计还加深了对我通信原理理论知识的理解，了解通信原理知识和其它课程的联系，例如交换技术、数字信号处理等，及其在实际生活中的应用，增强了自己查阅资料的动手能力。
课程设计报告 课程名称信号与系统课程设计指导教师 设计起止日期 学院信息与通信工程 专业电子信息工程 学生 班级/学号 成绩 指导老师签字
目录 1、课程设计目的 (1) 2、课程设计要求 (1) 3、课程设计任务 (1) 4、课程设计容 (1) 5、总结 (11) 参考文献 (12) 附录 (12)
1、课程设计目的 “信号与系统”是一门重要的专业基础课，MATLAB作为信号处理强有力的计算和分析工具是电子信息工程技术人员常用的重要工具之一。本课程设计基于MATLAB完成信号与系统综合设计实验，以提高学生的综合应用知识能力为目标，是“信号与系统”课程在实践教学环节上的必要补充。通过课设综合设计实验，激发学生理论课程学习兴趣，提高分析问题和解决问题的能力。 2、课程设计要求 （1）运用MATLAB编程得到简单信号、简单信号运算、复杂信号的频域响应图； （2）通过对线性时不变系统的输入、输出信号的时域和频域的分析，了解线性时不变系统的特性，同时加深对信号频谱的理解。 3、课程设计任务 （1）根据设计题目的要求，熟悉相关容的理论基础，理清程序设计的措施和步骤； （2）根据设计题目的要求，提出各目标的实施思路、方法和步骤； （3）根据相关步骤完成MATLAB程序设计，所编程序应能完整实现设计题目的要求； （4）调试程序，分析相关理论； （5）编写设计报告。 4、课程设计容 （一）基本部分 （1）信号的时频分析 任意给定单频周期信号的振幅、频率和初相，要求准确计算出其幅度谱，并准确画出时域和频域波形，正确显示时间和频率。 设计思路： 首先给出横坐标，即时间，根据设定的信号的振幅、频率和初相，写出时域波形的表达式；然后对时域波形信号进行傅里叶变化，得到频域波形；最后使用plot函数绘制各个响应图。 源程序： clc; clear; close all; Fs =128; % 采样频率 T = 1/Fs; % 采样周期 N = 600; % 采样点数 t = (0:N-1)*T; % 时间，单位：S x=2*cos(5*2*pi*t);
第1章绪论 双音多频DTMF（Dual Tone Multi-Frequency）信令，就是用两个频率——行频和列频来表示机键盘上的一个数字。双音多频信号是音频中的拨号信号，由美国AT&T贝尔公司实验室研制，并用于网络中。这种信号制式具有很高的拨号速度，且容易自动监测识别，很快就代替了原有的用脉冲计数方式的拨号制式。逐渐在全世界围使用在按键式机上，这种双音多频信号制式不仅用在网络中，还可以用于传输十进制数据的其它通信系统中，用于电子和银行系统中。这些系统中用户可以用发送DTMF信号选择语音菜单进行操作。作为实现快速可靠传输的一种技术，它具有很强的抗干扰能力和较高的传输速度，因此，可广泛用于通信系统中。但绝大部分是用作的音频拨号。另外，它也可以在数据通信系统中广泛地用来实现各种数据流和语音等信息的远程传输。近年来DTMF也应用在交互式控制中，诸如语言菜单、语言、银行和ATM 终端等。通过软件产生与检测DTMF 信令，是一项较有价值的工程应用。这是一种技术，就是机上的一个按键按下去时，机向交换机同时发送两个频率的信号，告诉交换机按的是哪个按键，以前采用脉冲方式，速度慢，一共有8个频率的音频信号，分为2组，每组4个，两两组合共可以代表16个按键，分别代表0-9 、#、*等按键。
第2章双音多频（DTMF）信号的设计 2.1设计目的及意义 双音多频信号（DTMF）是系统中机与交换机之间的一种用户信令，通常用于发送被叫。双音多频信号是贝尔实验室发明的，其目的是为了自动完成长途呼叫。 拨号有两种，脉冲和音频，所谓音频也称双音多频（DTMF）信号的拨号方式，双音多频既是拨号时每按一个键，有两个音频频率叠加成一个双音频信号，十二个按键由七个音频频率区分。在使用双音多频信号之前，系统中使用一连串的断续脉冲来传送被叫，称为脉冲拨号。脉冲拨号需要电信局中的操作员手工完成长途接续。双音多频的拨号键盘是4×4的矩阵，每一行代表一个低频，每一列代表一个高频。每按一个键就发送一个高频和低频的正弦信号组合，比如'1'相当于697和1209赫兹(Hz)。交换机可以解码这些频率组合并确定所对应的按键。本次课设的目的就是通过学习和掌握现代交换原理的基础上，设计一个双音多频检测模块并对电路进行仿真，综合应用所学知识，进行一次比较全面的训练，为今后的学习和工作积累经验。 此外，该题目还涵盖了《通信原理》、《电路分析》、《交换原理》等主要课程的知识点，学生通过该题目的设计过程，可以初步掌握DTMF编/解码技术原理和相关电路设计、开发原理，得到系统的训练，提高解决实际问题的能力。 2.2双音多频（DTMF）信号的组成 双音多频信号作为实现快速可靠传输的一种技术，它具有很强的抗干扰能力和较高的传输速度，因此，广泛应用于通信系统中。近年来，双音多频信号也应用在交互式控制中，如在语言控制、语言、银行和ATM 终端等的应用。 音频拨号当人们按下某一个按键时，会产生一组特定的双音信号，称为双音多频信号，交换机会对该信号进行处理，根据两个单音频率来识别所按下的。将拨号盘上的数字0～9 和两个标有“*”和“#”的特殊按钮进行频率分配，如图所示。包括两个频率低频段包括的频率是697Hz，770Hz，852Hz 和941Hz，称为行频。高频段含的频率 1209Hz，1336Hz，1477Hz，1633Hz，称为列频，它们可频率分配构成16 种频率组合，每一种组合由一对正弦频率信号唯一确定。其中第四列1633Hz 对应的按键目前并没有使用，留待将来扩展使用。
抽样定理的应用 摘要 抽样定理表示为若频带宽度有限的，要从抽样信号中无失真地恢复原信号，抽样频率应大于2倍信号最高频率。抽样频率小于2倍频谱最高频率时，信号的频谱有混叠。抽样频率大于2倍频谱最高频率时，信号的频谱无混叠。 语音信号处理是研究用数字信号处理技术和语音学知识对语音 信号进行处理的新兴学科，是目前发展最为迅速的学科之一，通过语音传递信息是人类最重要，最有效，最常用和最方便的交换信息手段，所以对其的研究更显得尤为重要。 Matlab语言是一种数据分析和处理功能十分强大的计算机应用 软件，它可以将声音文件变换成离散的数据文件，然后用起强大的矩阵运算能力处理数据。这为我们的本次设计提供了强大并良好的环境！ 本设计要求通过利用matlab对模拟信号和语音信号进行抽样，通过傅里叶变换转换到频域，观察波形并进行分析。 关键词：抽样Matlab
目录 一、设计目的： (2) 二、设计原理： (2) 1、抽样定理 (2) 2、MATLAB简介 (2) 3、语音信号 (3) 4、Stem函数绘图 (3) 三、设计内容： (4) 1、已知g1(t)=cos(6πt),g2(t)=cos(14πt),g3(t)=cos(26πt),以抽样频率 fsam=10Hz对上述三个信号进行抽样。在同一张图上画出g1(t)，g2(t)，g3(t)及其抽样点,对所得结果进行讨论。 (4) 2、选取三段不同的语音信号，并选取适合的同一抽样频率对其进 行抽样，画出抽样前后的图形，并进行比较，播放抽样前后的语音。 (6) 3、选取合适的点数，对抽样后的三段语音信号分别做DFT，画图 并比较。 (10) 四、总结 (12) 五、参考文献 (13)
信号与系统课程设计报告 实验题目：信号的运算与处理 内容简介： 设计一个信号，对其进行信号运算和处理，利用Matlab仿真。 课设方式： 利用电子技术、电路理论和信号与系统的知识学习验证信号的运算和处理，如延时、相加、微分、抽样等。自已设计信号及运算方式，并利用Matlab仿真。 分析计算结果。 课程设计要求： 独立完成； 完成信号设计（任意信号均可）及其某种运算（任意运算均可，也可多做几种，或做组合运算）的验证； 学会利用Matlab仿真；提交课程设计报告。 例如： 设计一个信号为f(t)=3sin2t 对其做微分运算得到f/(t) , 用MATLAB 编程实现计算过程，画出f(t)和f/(t)
本次课程设计本人选的信号运算是： 设计一个信号为y1=y(x)=sin2x,对其作微分运算得到dy1,用MATLAB对其实现运算过程，后画出y1,dy1,y1+dy1的图像 实验步骤（操作过程） 1、 首先打开MATLAB软件，在其命令窗口直接输入以下程序，对y(x)进 行微分运算。得到dy1 clear >> syms x y1; >> y1=sin(2*x); >> dy1=diff(y1,'x') dy1 =2*cos(2*x) 运算过程如下图所示： 2、 接着便是对其进行验证，点击fire，新建一个文件，输入以下程序（绘制出y1=sin2x, dy1=2cos2x, 以及y1+ dy1=sin2x+2cos2x。的波形）
3、保存文件，后缀名为.m,随后按F5执行输出输出图形。实验结果如下图所示 、
结果分析 如图所示绿色波形为y1=sin2x，蓝色为dy1=2cos2x,红色波形为y1+dy1。仿真结果与运算结果一致。 实验心得体会（调试过程） 总的来说，这次课程设计难度并不是太高，而我选取的正玄信号也是较为简单常用的一种函数，对其进行微分运算之后，得到了余弦函数，其仿真结果波形也如上所示，与预期一致。在设计过程中，还是出现了几个小问题的，一个是变量的定义，之前没有定义x，直接取范围结果出错了，还有一个是注意各种函数的调用以及运算格式，还是希望能在之后再接再厉，掌握好matlab软件！（附上调试过程图片） 左边为文件、历史窗口，底下是命令窗口，最右下角为实验仿真波形，中间为运算程序，绘图画图程序。
实验一、数字信号处理在双音多频拨号系统中的应用 一、实验目的 1．了解双音多频信号的产生、检测、包括对双音多频信号进行DFT 时的参数选择等。 2．初步了解数字信号处理在是集中的使用方法和重要性。 3．掌握matlab 的开发环境。 二、实验原理 双音多频（Dual Tone Multi Frequency, DTMF ）信号是音频电话中的拨号信号，由美国AT&T 贝尔公司实验室研制，并用于电话网络中。这种信号制式具有很高的拨号速度，且容易自动监测识别，很快就代替了原有的用脉冲计数方式的拨号制式。这种双音多频信号制式不仅用在电话网络中，还可以用于传输十进制数据的其它通信系统中，用于电子邮件和银行系统中。这些系统中用户可以用电话发送DTMF 信号选择语音菜单进行操作。 DTMF 信号系统是一个典型的小型信号处理系统，它要用数字方法产生模拟信号并进行传输，其中还用到了D/A 变换器；在接收端用A/D 变换器将其转换成数字信号，并进行数字信号处理与识别。为了系统的检测速度并降低成本，还开发一种特殊的DFT 算法，称为戈泽尔(Goertzel)算法，这种算法既可以用硬件（专用芯片）实现，也可以用软件实现。下面首先介绍双音多频信号的产生方法和检测方法，包括戈泽尔算法，最后进行模拟实验。下面先介绍电话中的DTMF 信号的组成。 在电话中，数字0-9的中每一个都用两个不同的单音频传输，所用的8个频率分成高频带和低频带两组，低频带有四个频率：679Hz,770Hz,852Hz 和941Hz ；高频带也有四个频率：1209Hz,1336Hz,1477Hz 和1633Hz.。每一个数字均由高、低频带中各一个频率构成，例如1用697Hz 和1209Hz 两个频率，信号用)2sin()2sin(21t f t f ππ+表示，其中Hz f 6791=，Hz f 12092=。这样8个频率形成16种不同的双频信号。具体号码以及符号对应的频率如表4.1所示。表中最后一列在电话中暂时未用。 表4.1 双频拨号的频率分配
实验：入侵检测系统（Snort）的安装与配置 一、实验目的 学会WINDOWS下SNORT的安装与配置 二、实验环境 WinXP虚拟机 三、实验步骤与结果 一．在“我的电脑”中C盘中建立文件夹“zhangxiaohong” 二．安装WinPcap，运行WinPcap_4_1_2.zip，默认安装。 三．安装mysql，运行mysql-5.0.22-win32.zip，选择自定义安装选择安装路径 C:\zhangxiaohong\mysql 下，安装时注意：端口设置为3306（以后要用到），密码本实验设置成123 四．安装apache 1.运行apache_ 2.2.4-win32-x86-no_ssl.zip，安装到c:\zhangxiaohong\Apache 2.安装Apache，配置成功一个普通网站服务器 3.出现Apache HTTP Server 2.0.55的安装向导界面，点“Next”继续 4.确认同意软件安装使用许可条例，选择“I accept the terms in the license agreement”，点“Next”继续 5.将Apache安装到Windows上的使用须知，请阅读完毕后，按“Next”继续 6.选择安装类型，Typical为默认安装，Custom为用户自定义安装，我们这里选 择Custom，有更多可选项。按“Next”继续 7.出现选择安装选项界面，如图所示，左键点选“Apache HTTP Server 2.0.55”，
选择“This feature, and all subfeatures, will be installed on local hard drive.” 8.即“此部分，及下属子部分内容，全部安装在本地硬盘上”。点选 “Change...”，手动指定安装目录。 9.我这里选择安装在“C:\zhangxiaohong\Apache”，各位自行选取了，一般建议 不要安装在操作系统所在盘，免得操作系统坏了之后，还原操作把Apache配置文件 也清除了。选“OK”继续。 10.返回刚才的界面，选“Next”继续。 11.好了现在我们来测试一下按默认配置运行的网站界面，在IE地址栏打 “.0.1”，点“转到”，就可以看到如下页面，表示Apache服务器已安装成功。 12. 五．安装和配置PHP53、安装winpcap 1.解压php-5. 2.5-Win32到c:\zhangxiaohong\php 2.添加gd图形库支持 复制c:\zhangxiaohong\php\php5ts.dll和c: \zhangxiaohong\php\libmysql.dll文件到 C:\Windows\system32 复制c: \zhangxiaohong\php\php.ini-dist到C:\Windows文件夹并重命名为php.ini， 修改php.ini，分别去掉“extension=php_gd2.dll”和“extension=php_mysql.dll”前的分号， 3.并指定extension_dir="c:\zhangxiaohong\php\ext"， 4.同时复制c:\zhangxiaohong\php\ext下的php_gd2.dll与php_mysql.dll到C:\Windows\system32 在C:\zhangxiaohong\apache\conf\httpd.conf中添加 LoadModule php5_module c:/zhangxiaohong/php/php5apache2.dll AddType application/x-httpd-php .php AddType application/x-httpd-php-source .phps AddType application/x-httpd-php .html AddType application/x-httpd-php .htm 5.重启Apache服务 在C:\zhangxiaohong\apache\htdocs目录下新建webinf.php（文件内容为：）并使用访问测试是否能够显示当前Apache服务器的信息，如果能够显示表明Apache和php工作基本正常 六．安装snort 1.运行Snort_2_9_0_5_Installer.exe 安装在C:\zhangxiaohong\Snort下即可， 运行C:\zhangxiaohong\Snort\bin\snort.exe或者在DOS中找到该位置， 如果安装Snort成功会出现一个可爱的小猪 2.并按照以下修改C:\zhangxiaohong\Snort\etc\snort.conf文件
语音信号的数字滤波 ——利用双线性变换法实现IIR数字滤波器的设计一．课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法，掌握利用数字滤波器对信号进行滤波的方法，掌握数字滤波器的计算机仿真方法，并能够对设计结果加以分析。 二．设计方案论证 1.IIR数字滤波器设计方法 IIR数字滤波器是一种离散时间系统，其系统函数为 假设M≤N，当M＞N时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数和，它 是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。 2.用双线性变换法设计IIR数字滤波器 脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到Ｚ平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用z=e sT转换 平面的-π/T～π到Z平面上。也就是说，第一步先将整个S平面压缩映射到S 1 /T一条横带里；第二步再通过标准变换关系z=e s1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图1所示。 图1双线性变换的映射关系 为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的-π/T到π/T段上，可以通过以下的正切变换实现
一、课程设计题目 去除干扰蜂鸣音 1.目的：掌握信号时频域分析方法，正确理解采样定理，准确理解滤波器的概念。 2.内容：提供一个包含某人说话语音片段的声音文件，但该语音信号被一个包含有几个谐波分量的蜂鸣信号干扰了。 用Matlab 的wavread 命令读取该声音文件。注意，该命令可以同时得到声音文件的采样率和采样位宽，请查阅Matlab 的帮助文件。 （1） 用快速傅立叶变换（FFT ）计算并画出声音信号的频谱，列写出蜂鸣信号的谐波频率。 （2） 思考如何将这些蜂鸣音去除？将去除了蜂鸣音的语音片段播放出来，仔细聆听并写下语音片 段中人物所说的话。注意：由于只能播放实信号，因此记得提取信号的实部。 Matlab 命令：wavread, wavplay, fft, fftshift, fir1, filter, plot, figure. 二、设计思路 用waveread()函数读取音频和其采样率和采样位宽，对读取的音频信号使用fft()函数进行快速傅立叶变换并绘出得到的频谱。观察频谱分析噪声（蜂鸣信号）的谐波频率分布，选择合适的滤波模式将噪声信号的谐波滤去，便可以得到去除噪声后的人声。 设计滤波器的频域特性便成了除去噪声并留下原声的关键，我们注意到所学的采样定理以及一维sinc 函数（辛格函数）x x x Sinc ππ) sin()(=，然而汉宁窗可以看作是3个矩形时间窗的频谱之和，或者 说是 3个)(x Sinc 型函数之和，而括号中的两项相对于第一个谱窗向左、右各移动了π/T ，从而使旁瓣互相抵消，消去高频干扰和漏能。它适用于非周期性的连续信号。利用它的选择特性使用fir1()建立一个汉宁窗函数，并用filter()函数进行滤波，去除噪声部分。最后用play()函数播放音频检查效果。 三、设计过程 1.音频的读取和分析 先将原始音频文件读入， [audio0, Fs, nbits] = wavread('C:\Users\Administrator\Desktop\signal\buzz.wav')；%按路径读取音频存入audio0变量，并用Fs 变量记录采样率，nbits 变量记录采样位宽。 其中， Fs=11025 #采样率为11025Hz nbits=32 #采样带为32 p0=audioplayer(audio0,Fs);%将audio0载入音频播放器 play(p0)；%并进行播放 subplot(2,1,1);%分屏绘图 plot(audio0);%绘制原始音频时域图,如下图所示 title('时域');%标注题目
实验三双音多频信号的合成与检测 一实验目的 1．理解电话拨号音的合成与检测的基本原理； 2．深入理解信号频谱分析理论中相关参数的作用和意义； 3．了解频谱分析在实际工程中的应用实例。 二实验基础 双音多频(dual-tone multifrequency, DTMF信号的产生及检测在现代通信系统中有着广泛的应用，家用电话、移动电话以及公共程控交换机(PBX都采用DTMF 信号发送和接收电话拨号号码。本实验要求利用信号的时域分析和频域分析的基本理论实现DTMF 的合成和检测。 1. DTMF信号合成 DTMF 信号由低频组和高频组两组频率信号构成。按键电话上每个按键都由对应的两个频率组成，如表4.1。当按下某个键时，所得到的按键信号是由相应两个频率的正弦信号叠加而成。设x(n为DTMF 信号，产生方式为： x (n =sin (ωH n +sin (ωH n 式中：ωH = f s DTMF 信号的标准是：在传送过程中每个按键字占用100ms ，其中信号必须持续至少40ms ，且不得多于55ms ，100ms 里的其余时间为静音（无信号）。 表4.1按键频率对应表
2. DTMF信号检测 ,ωL = f s f s =8KHz 。 DTMF 信号的检测是将信号的两个频率提取出来，从而确定接收到的DTMF 对应的按键。利用DFT 对DTMF 信号进行N 点的频谱分析，N 的选取决定了频率分辨率以及捕捉N 个样值所需要的时间。根据谱峰出现的频率点位置m 就可以确定DTMF 信号的频率f k： f k =kf s /N 这样计算出的DTMF 信号频率可能与实际的DTMF 信号频率有一定的差别，但可以通过加大N 的选取来减小这种频率差异。然而从另外一方面来考虑，虽然加大N 值会减小检测频率误差，但这势必会带来捕捉N 个样值所需要的时间增加，从而会对检测的效果造成一定影响。 由DTMF 信号频率所具有的特性不难发现要选取一定的N 值使得计算出的频率和真实的DTMF 信号的频率相一致几乎不可能，而实际中也并不需要计算出来的频率值与其真实频率相一致，只需偏差保持在±1.5%即可认为是DTMF 信号的真实频率。国际上通用N=205点或N=106点。当N=205点时，各个频率所对应的DFT 结果X[k]中的序号k 如表4.2。N=106时对应表4.3。
数字信号处理课程设计 姓名：刘倩 学号：201014407 专业：信息与计算科学 实验一：常见离散信号产生和实现 一、实验目的： 1、加深对常用离散信号的理解； 2、掌握matlab 中一些基本函数的建立方法。 二、实验原理： 1.单位抽样序列 在MATLAB 中可以利用zeros()函数实现。 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： 2．单位阶越序列 在MATLAB 中可以利用ones()函数实现。 3．正弦序列 在MATLAB 中 4．复指数序列 在MATLAB 中 5．指数序列 在MATLAB 中
实验内容：由周期为10的正弦函数生成周期为20的余弦函数。 实验代码： n=0:30; y=sin(0.2*pi*n+pi/2); y1=sin(0.1*pi*n+pi/2); subplot(121) stem(n,y); xlabel ('时间序列n');ylabel('振幅');title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); subplot(122) stem(n,y1); xlabel ('时间序列n');ylabel('振幅'); title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); 实验结果： 实验二：离散系统的时域分析 实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。实验原理：离散系统 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号， 记系统单位冲激响应 则系统响应为如下的卷积计算式：
当N k d k ,...2,1,0==时，h[n]是有限长度的（n ：[0，M]），称系统为FIR 系统；反之，称系统为IIR 系统。 在MATLAB 中，可以用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应。 实验内容：用MATLAB 计算全解 当n>=0时，求用系数差分方程y[n]+y[n-1]-6y[n-2]=x[n]描述的一个离散时间系统对阶跃输入x[n]=8μ[n]的全解。 实验代码： n=0:7; >> [y,sf]=filter(1,[1 1 -6],8*ones(1,8),[-7 6]); >> y1(n+1)=-1.8*(-3).^n+4.8*(2).^n-2; >> subplot(121) >> stem(n,y); >> title('由fliter 函数计算结果'); >> subplot(122) >> stem(n,y1); >> title('准确结果'); 实验结果： 结果分析：有图可得由fliter 函数得出的结果与计算出的准确结果完全一致。 实验三FFT 算法的应用
南通大学电子信息学院信号与系统课程设计报告 班级： 姓名： 学号： 指导老师： 设计时间： 2014—2015学年第一学期
一、连续信号的时域分析 二、 1. 信号的产生 (1)阶跃函数 function [t,y1]=jieyue(t1,t2,t0) dt=0.01; ttt=t1:dt:t0-dt; tt=t0:dt:t2; t=t1:dt:t2; n=length(ttt); nn=length(tt); u=zeros(1,n); uu=ones(1,nn); y1=[u,uu]; return 冲激函数 function [t,y2]=chongji(t1,t2,t0) dt=0.01; t=t1:dt:t2; n=length(t); y2(1:n)=0; y2(1,(t0-t1)/dt+1)=1/dt; （2）调用上述函数产生信号）2-t ε（，）（4-t δ，-t e ）（t ε，-6s ≤t ≤6s,并画出波形。 Command Window subplot(3,1,1); [t1,y1]=jieyue(-6,6,2);
stairs(t1,y1); axis([-6 6 0 1.5]); subplot(3,1,2); [t2,y2]=chongji(-6,6,4);plot(t2,y2); subplot(3,1,3); [t3,y3]=jieyue(-6,6,0); y3=exp(-(t3)).*y3;plot(t3,y3); 波形如下图所示： （3）根据f(t)画出f(2t)和f(1-0.5t)的波形 t=-3:0.01:3; y=tripuls(t,4,0.6); subplot(3,1,1); plot(t,y);
通信技术方向课程设计题目：双音多频信号检测 物联网工程学院电子信息工程专业 学号0703070106 学生姓名时雅茹 二〇一〇年六月
一、原理介绍 双音多频（Dual Tone Multi Frequency, DTMF ）信号是音频电话中的拨号信号，由美国AT&T 贝尔公司实验室研制，并用于电话网络中。这种信号制式具有很高的拨号速度，且容易自动监测识别，很快就代替了原有的用脉冲计数方式的拨号制式。这种双音多频信号制式不仅用在电话网络中，还可以用于传输十进制数据的其它通信系统中，用于电子邮件和银行系统中。这些系统中用户可以用电话发送DTMF 信号选择语音菜单进行操作。 DTMF 信号系统是一个典型的小型信号处理系统，它要用数字方法产生模拟信号并进行传输，其中还用到了D/A 变换器；在接收端用A/D 变换器将其转换成数字信号，并进行数字信号处理与识别。为了系统的检测速度并降低成本，还开发一种特殊的DFT 算法，称为戈泽尔(Goertzel)算法，这种算法既可以用硬件（专用芯片）实现，也可以用软件实现。下面首先介绍双音多频信号的产生方法和检测方法，包括戈泽尔算法，最后进行模拟实验。 二、内容及结论 1、双音多频（DTMF ）信号的组成 在电话中，数字0~9的中每一个都用两个不同的单音频传输，所用的8个频率分成高频带和低频带两组，低频带有四个频率：679Hz,770Hz,852Hz 和941Hz ；高频带也有四个频率：1209Hz,1336Hz,1477Hz 和1633Hz.。每一个数字均由高、低频带中各一个频率构成，例如1用697Hz 和1209Hz 两个频率，信号用 )2sin()2sin(21t f t f ππ+表示，其中Hz f 6791=，Hz f 12092=。这样8个频率形成 16种不同的双频信号。具体号码以及符号对应的频率如表1所示。表中最后一列在电话中暂时未用。 表1 双频拨号的频率分配 列 行 1209Hz 1336Hz 1477Hz 633Hz 697Hz 1 2 3 A 770Hz 4 5 6 B 852Hz 7 8 9 C 942Hz * # D DTMF 信号在电话中有两种作用，一个是用拨号信号去控制交换机接通被叫的用户电话机，另一个作用是控制电话机的各种动作，如播放留言、语音信箱等。
遵义师范学院计算机与信息科学学院 实验报告 （2013—2014学年第1 学期） 课程名称：网络安全实验 班级： 学号： 姓名： 任课教师： 计算机与信息科学学院
实验报告
闸门，在不影响网络性能的情况下能对网络进行监测，从而提供对内部攻击、外部攻击和误操作的实时保护。这些都通过它执行以下任务来实现： ·监视、分析用户及系统活动 · 系统构造和弱点的审计 · 识别反映已知进攻的活动模式并向相关人士报警 · 异常行为模式的统计分析 · 评估重要系统和数据文件的完整性 ·操作系统的审计跟踪管理，并识别用户违反安全策略的行为。 对一个成功的入侵检测系统来讲，它不但可使系统管理员时刻了解网络系统（包括程序、文件和硬件设备等）的任何变更，还能给网络安全策略的制订提供指南。更为重要的一点是，它应该管理、配置简单，从而使非专业人员非常容易地获得网络安全。而且，入侵检测的规模还应根据网络威胁、系统构造和安全需求的改变而改变。入侵检测系统在发现入侵后，会及时作出响应，包括切断网络连接、记录事件和报警等。 入侵检测系统所采用的技术可分为特征检测与异常检测两种: 特征检测(Signature-based detection) 又称Misuse detection ，这一检测假设入侵者活动可以用一种模式来表示，系统的目标是检测主体活动是否符合这些模式。它可以将已有的入侵方法检查出来，但对新的入侵方法无能为力。其难点在于如何设计模式既能够表达“入侵”现象又不会将正常的活动包含进来。 异常检测(Anomaly detection)的假设是入侵者活动异常于正常主体的活动。根据这一理念建立主体正常活动的“活动简档”，将当前主体的活动状况与“活动简档”相比较，当违反其统计规律时，认为该活动可能是“入侵”行为。异常检测的难题在于如何建立“活动简档”以及如何设计统计算法，从而不把正常的操作作为“入侵”或忽略真正的“入侵”行为 Snort入侵检测系统： Snort简介：在1998年，Martin Roesch先生用C语言开发了开放源代码(Open Source)的入侵检测系统Snort.直至今天，Snort已发展成为一个多平台 (Multi-Platform),实时(Real-Time)流量分析，网络IP数据包(Pocket)记录等特性的强大的网络入侵检测/防御系统(Network Intrusion Detection/Prevention System),即NIDS/NIPS.Snort符合通用公共许可(GPL——GUN General Pubic License),在网上可以通过免费下载获得Snort,并且只需要几分钟就可以安装并开始使用它。snort基于libpcap。 Snort有三种工作模式：嗅探器、数据包记录器、网络入侵检测系统。嗅探器模式仅仅是从网络上读取数据包并作为连续不断的流显示在终端上。数据包记录器模式把数据包记录到硬盘上。网络入侵检测模式是最复杂的，而且是可配置的。我们可以让snort 分析网络数据流以匹配用户定义的一些规则，并根据检测结果采取一定的动作。 Snort原理：Snort能够对网络上的数据包进行抓包分析，但区别于其它嗅探器的是，它能根据所定义的规则进行响应及处理。Snort 通过对获取的数据包，进行各规则
数字信号处理课程设计报告《应用Matlab对信号进行频谱分析及滤波》 专业： 班级： 姓名： 指导老师： 二0 0五年一月一日
目录 设计过程步骤（） 2.1 语音信号的采集（） 2.2 语音信号的频谱分析（） 2.3 设计数字滤波器和画出其频谱响应（） 2.4 用滤波器对信号进行滤波（） 2.5滤波器分析后的语音信号的波形及频谱（） ●心得和经验（）
设计过程步骤 2.1 语音信号的采集 我们利用Windows下的录音机，录制了一段开枪发出的声音，时间在1 s内。接着在C盘保存为WAV格式，然后在Matlab软件平台下．利用函数wavread对语音信号进行采样，并记录下了采样频率和采样点数，在这里我们还通过函数sound引入听到采样后自己所录的一段声音。通过wavread函数和sound的使用，我们完成了本次课程设计的第一步。其程序如下： [x,fs,bite]=wavread('c:\alsndmgr.wav',[1000 20000]); sound(x,fs,bite); 2.2 语音信号的频谱分析 首先我们画出语音信号的时域波形；然后对语音信号进行频谱分析，在Matlab中，我们利用函数fft对信号进行快速傅里叶变换，得到信号的频谱特性性。到此，我们完成了课程实际的第二部。 其程序如下： n=1024; subplot(2,1,1); y=plot(x(50:n/4)); grid on ; title('时域信号') X=fft(x,256); subplot(2,1,2); plot(abs(fft(X))); grid on ; title('频域信号'); 运行程序得到的图形：
一、 1.正弦信号 A = input('input A=') ;% 给正弦信号的幅度A赋值 w = input('input w=') ; % 给正弦信号的频率w赋值 theta =input('input theta='); % 给正弦信号的初始相位theta 赋值disp(['这个信号是周期信号']) T=2*pi/w t = 0 : 0.01 : 3*T ; % 定义时间点 ft = A * sin( w * t + theta ) ; % th计算函数值 plot( t ,ft ) ; % 画图 title( '正弦信号' ) ; % 为图像加标题注释 grid on ; % 在图上画方格
2.复指数信号 j00 = sqrt( - 1 ) ; % 定义复数j a = input('input a='); % 复指数信号赋值w = input('input w='); K = input('input K='); if a==0 disp('这是一个周期信号') T=2*pi/w else if a>0 disp('这不是一个周期信号') else disp('这不是一个周期信号') end end t = -1.5*abs(a) : 0.01 : 1.5*abs(a) % 定义时间点 ft = K*exp( ( a + j00 * w ) * t ) ; subplot( 2 , 2 , 1 ) ; plot( t , real( ft ) ) ; title( '实部' ) ; %画图subplot( 2 , 2 , 2 ) ; plot( t , imag( ft ) ) ; title( '虚部' ) ; subplot( 2 , 2 , 3 ) ; plot( t , abs( ft ) ) ; title( '模' ) ; subplot( 2 , 2 , 4 ) ; plot( t , angle( ft ) ) ; title( '相角' ) ;
数字信号处理 课 程 设 计 院系：电子信息与电气工程学院 专业：电子信息工程专业 班级：电信班 姓名： 学号： 组员：
摘要 滤波器设计在数字信号处理中占有极其重要的地位，FIR数字滤波器和IIR 滤波器是滤波器设计的重要组成部分。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR和IIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。 关键词数字滤波器 MATLAB 窗函数法巴特沃斯
目录 摘要 (1) 1 引言 (1) 1.1课程设计目的 (1) 1.2 课程设计内容及要求 (1) 1.3课程设计设备及平台 (1) 1.3.1 数字滤波器的简介及发展 (1) 1.3.2 MATLAB软件简介 (2) 2 课程设计原理及流程 (4) 3.课程设计原理过程 (4) 3.1 语音信号的采集 (4) 3.2 语音信号的时频分析 (5) 3.3合成后语音加噪声处理 (7) 3.3.1 噪声信号的时频分析 (7) 3.3.2 混合信号的时频分析 (8) 3.4滤波器设计及消噪处理 (10) 3.4.1 设计IIR和FIR数字滤波器 (10) 3.4.2 合成后语音信号的消噪处理 (13) 3.4.3 比较滤波前后语音信号的波形及频谱 (13) 3.4.4回放语音信号 (15) 3.5结果分析 (15) 4 结束语 (15) 5 参考文献 (16)
信号与系统课程设计 课程名称：信号与系统 题目名称：滤波器的设计与实现 学院：电气与电子工程学院 专业班级：电气工程及其自动化 学号：3 学生：宗喜 指导教师：黄劲 2015年12 月20 日
目录 一、设计要求 (2) 二、设计原理 (2) 三、设计思路 (3) 四、设计容 (3) A、一阶有源滤波电路 (3) B、二阶有源滤波电路 (5) 1、二阶低通滤波电路 (5) 2、二阶高通滤波电路 (6) 3、二阶带通滤波电路 (8) C、用仿真软件设计滤波器 (10) 1、给定性能参数设计滤波器 (10) a、二阶低通滤波器 (10) b、二阶高通滤波器 (11) c、二阶带通滤波器 (12) 2、不同阶数滤波器性能比较 (12) D、滤波器的Matlab设计仿真 (13) 1、二阶低通滤波器 (13) 2、二阶高通滤波器 (14) 五、参考文献 (16)
一、设计要求 自已设计电路系统，构成低通滤波器、高通滤波器和带通滤波器。利用Matlab或其他仿真软件进行仿真。 有源滤波器由是有源元件和无源元件(一般是R和C)共同组成的电滤波器。和无源滤波器相比，它的设计和调整过程较简便，此外还能提供增益。因此，本课程设计中选择了二阶有源滤波器作为主要研究对象。 1、自行设计电路图，确定前置放大电路，有源滤波电路，功率放大电路的方案， 并使用绘图软件（Electronics Worrkbench）画出设计电路，包括低通、高通和带通。 2、所设计的滤波器不仅有滤波功能，而且能起放大作用，负载能力要强。 3、根据给定要求和电路原理图计算和选取单元电路的元件参数。 4、用Matlab或其他仿真软件（FilterLab）对滤波器进行仿真，记录仿真结果。 二、设计原理 1、电容器C具有通高频阻低频的性能。 2、由源滤波器由放大电路部分和滤波电路部分组成。 3、仿真软件可以将滤波器的性能直观的表现出来。 4、各种滤波器的幅频特性：
第27卷 第3期 吉首大学学报(自然科学版) V ol.27 N o.3 2006年5月 Journal of Jishou University (Natural Science Edition ) M ay.2006 文章编号:1007-2985(2006)03-0043-05 双音多频信号检测在DSP 中的实现 Ξ 李义府,彭卫韶 (中南大学信息科学与工程学院,湖南长沙 410083) 摘 要:双音多频DT MF 信号是音频电话中的拨号信号,将DT MF 信号的检测集成到含有数字信号处理器(DSP )的系统中,是一项较有价值的研究课题.笔者设计出TI 公司浮点DSP 芯片C6711中的实现方案,通过20个并行的哥兹柔信号滤波器能成功地使双音多频信号的检测变得迅速和简单,采用谐音检测可以显著地提高检测系统的准确性. 关键词:哥兹柔滤波器;双音多频;浮点DSP ;Bellcore 标准中图分类号:T N914.3;T N911.7 文献标识码:A 按键式电话拨号广泛采用双音多频信号,近年来双音多频信号(DT MF )逐渐应用于工程信号发生与检 测系统中,并与DSP 、FPG A 相互促进,具有广泛的应用前景[1-2] .一个有效的音频信号由一个行频信号和一个列频信号叠加而成.例如,要表示“4”这个音频信号,可由一个770H z 的行频信号和一个1209H z 的列频信号叠加而成.电话音频拨号使用的正弦音频叠加信号如表1所示: 表1 电话机键盘的频率矩阵 行频组ΠH z 列频组ΠH z 1209136614771633697123A 770456B 8527 89C 941 3 # D 1 算法与滤波器 1.1算法 图1 离散通频带的分布 由于在实现DT MF 解码时,采用哥兹柔算法(G o 2 ertzel Alg orithm Theory )要比FFT 更快,因为通过FFT 可以计算得到信号所有频谱线,但处理DT MF 信号只考虑其中的8个频率及其二次谐波信息,运用G A T 能更加快速的从输入信号中提取频谱信息,所以使用C6177浮点DSP 进行信号检测不失为一种可行的技 术解决方案[3-4] .通过对信号作离散傅立叶变换得到其离散通频带.离散通频带的数目用字母N 表示,这些通频带段在频域中均匀分布如图1所示. Ξ 收稿日期:2006-02-26 基金项目:国家自然科学基金资助项目(69974043);湖南省自然科学基金资助项目(99JJ Y 20062) 作者简介:李义府(1946-),男,湖南省长沙市人,中南大学信息科学与工程学院教授,主要从事电子线路和故障诊断应用研究.