数字信号处理数字滤波器设计及在心电信号滤波中的应用的课程设计
课程设计报告
课程名称数字信号处理课程设计
课题名称数字滤波器设计及在心电信号滤波中的应用
专业通信工程
班级1081
学号201013120103
姓名刘献文
指导教师彭祯谭小兰
2012年9 月10日
湖南工程学院
课程设计任务书
课程名称数字信号处理课程设计
课题名称数字滤波器设计及在心电信号滤波中的应用
专业通信工程
班级1081
学号201013120103
姓名刘献文
指导教师彭祯谭小兰
审批
任务书下达日期2012 年9月 1 日
任务完成日期2012 年9月10日
《数字信号处理》课程设计指导
一、课程设计的性质与目的
《数字信号处理》课程是通信专业的一门重要专业基础课，是信息的数字化处理、存储和应用的基础。通过该课程的课程设计实践，使学生对信号与信息的采集、处理、传输、显示、存储、分析和应用等有一个系统的掌握和理解；巩固和运用在《数字信号处理》课程中所学的理论知识和实验技能，掌握数字信号处理的基础理论和处理方法，提高分析和解决信号与信息处理相关问题的能力，为以后的工作和学习打下基础。
数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应（IIR）滤波器和有限冲激响应（FIR）滤波器。
二、课程设计题目
课程设计题目分为2个方向，每班选做方向一的同学约占2/3,选做方向二的约占1/3。
方向一：数字滤波器设计及在语音信号分析中的应用。
方向二：数字滤波器设计及在心电信号滤波中的应用。
步骤：
1、心电信号采集
心电信号作为心脏电活动在人体体表的表现，信号一般比较微弱，幅度在10μV~5mV，频率为0.05~100Hz。在心电信号的采集、放大、检测及记录过程中，有来自外界的各种干扰。记录一段时间内的人体心电信号波形，要求长度不小于10秒，并对记录的信号进行数字化，保存为数据文件；这里，请同学们使用美国的MIT/BIH心电原始数据，由实验老师给出一定长度的的心电原始数据，数据保存在文件“a01.txt～a10.txt”中，在MATLAB中通过如下语句读取：%从当前路径下的a01.txt文件读取心电原始数据到变量a01中，a01为二维数据，第一列%为心电信号时间，第二列为心电信号幅度。
2、心电信号分析
使用MATLAB绘出数字化后的心电信号的时域波形和频谱图。根据频谱图求出其带宽，并说明心电信号的基本特征。
3、含噪心电信号合成
在MATLAB软件平台下，给原始的心电信号叠加上噪声或干扰，干扰类型分为如下几种：（1）白噪声；（2）工频干扰（50Hz）；（3）谐波干扰（二次、三次谐波为主，分别为100Hz、150Hz）；（4）其它干扰，可设置为低频、高频、带限噪声，或冲激干扰。绘出叠加噪声后的心电信号时域和频谱图，在视觉上与原始心电信号图形对比，绘出其时域波形差，分析频域基本特征变化。
4、数字滤波器设计及滤波,完成以下题目中的一个
给定滤波器的规一化性能指标（参考指标,实际中依据每个同学所叠加噪声情况而定）例如：通带截止频率wp=0.25*pi, 阻通带截止频率ws=0.3*pi; 通带最大衰减Rp=1 dB; 阻带最小衰减Rs=15 dB，每个题目至少设计出5个用不同方法的不同类型滤波器。
题目（1）：采用窗函数法与等波纹法分别设计各型FIR滤波器（低通、高通、带通、带阻中的至少3种类型）来对叠加干扰前后的心电信号进行滤波处理，绘出滤波器的频域响应，绘出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；在相同的性能指标下比较各方法的滤波效果，并从理论上进行分析（或解释）。
题目（2）：采用双线性变换法与脉冲响应不变法，分别利用不同的原型低通滤波器（Butterworth型与切比雪夫I型）来设计各型IIR滤波器（低通、高通、带通、带阻中的至少3种类型），绘出滤波器的频域响应；并用这些数字滤波器对含噪心电信号分别进行滤波处理，比较不同方法下设计出来的数字滤波器的滤波效果，并从理论上进行分析（或解释）。
5、心电信号波形观察、频谱观察
对滤波后的心电信号观察其时域、频域特征变化。绘出滤波后、滤波前、加噪后三个心电信号的差值波形，观察相互间的差异性；同时，分析频谱变化。
学生也可选用持续时间更长的心电原始数据，加上干扰后按上述要求设计滤波器。
三、课程设计要求
1、在一周内学生须上机16小时以上，程序调试完后，须由指导老师在机器上检查运行结果，经教师认可后的源程序可通过打印机输出，并请教师在程序清单上签字。
2、课程设计报告内容和格式：设计题目，设计的详细步骤，设计过程中的结果、图形等，设计总结。
3、每组每人必须独立完成，成绩的考核按设计结果、答辩成绩及课程设计报告来综合评定。成绩分为优、良、中、及格、不及格五级分评定。
4、指导教师：彭祯，谭小兰。
四、设计进度安排
通信工程1081:
1周周一下午，E-412, 任务讲解与布置，学生分组选题，查找相关资料，准备课程设计，学生上机，按任务要求进行课程设计；分组选题；
1周周二下午，E-412,学生上机，按任务要求进行课程设计；分组选题；
1周周三下午，E-412,学生上机，按任务要求进行课程设计；分组选题；
1周周四下午，E-412,学生任务完成，答辩并提交课程设计报告。
附：
课程设计报告装订顺序：封面、任务书、目录、正文、评分、附件（A4大小的图纸及程序清单）。
正文的格式:一级标题用3号黑体,二级标题用四号宋体加粗,正文用小四号宋体;行距为22。
正文的内容:一、课题的主要功能；二、课题的功能模块的划分；三、主要功能的实现；四、程序调试；
五、总结；六、附件（所有程序的原代码，要求对程序写出必要的注释）；七、评分表。
目录
1、课程设计的目的 (1)
2、课程设计的原理 (1)
2.1用窗函数法设计FIR滤波器 (1)
2.2用巴特沃斯法设计IIR低通滤波器 (1)
3、课程设计设计步骤及结果分析 (1)
3.1心电信号采集................... 错误！未定义书签。
3.2源数据的导入 (2)
3.3绘出源心电信号的时域波形图和频谱图并进行分析 (2)
3.4含噪心电信号合成 (3)
3.5时域波形差 (5)
3.6数字滤波器设计及滤波 (6)
3.7心电信号时域波形及频谱的观察和分析 (12)
4、心得体会 (18)
5、附件：（源代码） (18)
6、评分表 (22)
1 课程设计的目的
通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；熟悉设计数字滤波器的方法，掌握利用数字滤波器对信号进行滤波的方法，掌握数字滤波器的计算机仿真，并能够对设计结果加以分析。
2 课程设计的原理
2.1、用窗函数法设计FIR滤波器
根据过渡宽带及阻带衰减要求，选择窗函数的类型并估计窗口长度N，窗函数类型可根据最小阻带衰减As独立选择，因为窗口长度N对最小阻带衰减As没有影响，在确定窗函数类型以后，可根据过渡宽带小于给定指标确定所拟用的窗函数的窗口长度N，设待求滤波器的过渡宽带，它与窗口长度N近似成反比，窗函数类型确定后，其计算公式也确定了，不过这些公式是近似的，得出的窗口长度还要在计算中逐步修正，原则上在保证阻带衰减满足要求的情况下，尽量选择较小的N，在N和窗函数类型确定后，即可调用MATLAB 中的窗函数求出窗函数。
2.2、用巴特沃斯法设计IIR低通滤波器
巴特沃斯低通滤波器的幅度平方函数中，N为滤波器的阶数。在Ω=Ωc附近，随着Ω加大，幅度迅速下降。幅度下降的速度与阶数N有关，N越大，通带愈平坦，过渡带愈窄，过渡带与阻带幅度下降的速度愈快，总的频响特性与理想低通滤波器的误差愈小。
3 课程设计设计步骤及结果分析
3.1、心电信号采集
心电信号作为心脏电活动在人体体表的表现，信号一般比较微弱，幅度在10μV~5mV，频率为0.05~100Hz。在心电信号的采集、放大、检测及记录过程中，有来自外界的各种干扰。记录一段时间内的人体心电信号波形，要求长度不小于10秒，并对记录的信号进行数字化，保存为数据文件；这里，请同学们使用美国的MIT/BIH心电原始数据，由实验老师给出一定长度的的心电原始数据，数据保存在文件“a01.txt～a10.txt”中，在MATLAB中通过如下语句读取：
%从当前路径下的a01.txt文件读取心电原始数据到变量a01中，a01为二维数据，第一列%为心电信号时间，第二列为心电信号幅度。
3.2、源数据的导入
a=load('F:/心电信号数据/a16.txt');%选第十六组数据
3.3、绘出源心电信号的时域波形图和频谱图并进行分析
%时域波形图以及频谱图
a=load('F:/心电信号数据/a16.txt');
t=a(:,1);b=a(:,2);
figure(1);
subplot(211);
plot(t,b);
title('心电信号的时域波形');
xlabel('时间t/s');ylabel('幅值/A');
n=1000;
m=abs(fft(b,n));
fs=100;
f=fs/n*(0:n-1);
subplot(212);
plot(f,m);
title('心电信号的频谱图');
xlabel('频率 f/Hz');ylabel('幅值/db');
0123
45678910
-2
024
6心电信号的时域波形
时间t/s
幅值/A
010203040
5060708090100
50
100
150心电信号的频谱图
频率 f/Hz
幅值/d b
3.4、含噪心电信号合成
1、加入单频正弦干扰 y=0.2*sin(2*40*pi*t); y1=y+b; figure(2); subplot(211); plot(t,y1);
title('加单频噪声时域波形'); xlabel('时间t/s');ylabel('幅值/A'); k=abs(fft(y1,n)); subplot(212); plot(f,k);
title('加单频噪声频谱图'); xlabel('频率 /Hz');ylabel('幅值')
0123
45678910
-2
024
6加单频噪声时域波形
时间t/s
幅值/A
010203040
5060708090100
50
100
150加单频噪声频谱图
频率 /Hz
幅值
选择添加单频正弦信号。通过转置后添加到原始信号b 上，添加后的信号用y1表示。对于变量、函数的操作是矩阵操作，两个信号相加必须长度一致。
2、继续加入白噪声 y2 = awgn(y1,5); figure(3); subplot(211); plot(t,y2);
title('加单频,白噪声时域波形'); xlabel('时间t/s');ylabel('幅值/A'); h=abs(fft(y2,n)); subplot(212); plot(f,h);
title('加单频,白噪声频谱图'); xlabel('频率 /Hz');ylabel('幅值');
0123
45678910
-2
024
6加单频,白噪声时域波形
时间t/s
幅值/A
10
20
30
40
506070
80
90
100
050100150
200加单频,白噪声频谱图
频率 /Hz
幅值
通过调用函数y2 = awgn(y1,5)对已经加了单频噪声的信号进行加噪，加噪后的信号用y2保存。对y2信号的分析，如图所示分别给出了它的时域波形图和频谱图，通过对比，时域波形变化不大，但频谱图有了大幅度的变化。
3.5、时域波形差
s=y2-b; figure(4); subplot(111); plot(t,s);
title('时域波形差');
012345678910
-2
-1.5-1-0.500.511.5
2时域波形差
3.6、数字滤波器设计及滤波
1、设计低通FIR 滤波器 %低通FIR 滤波器 N=50; wc=0.6;
window=blackman(N); hn=fir1(N-1,wc,window); lv=filter(hn,1,y2); freqz(hn,1);
00.10.2
0.30.40.50.60.70.80.91
-4000
-3000-2000-1000
Normalized Frequency (?π rad/sample)
P h a s e (d e g r e e s )
0.1
0.2
0.30.40.50.60.70.80.9
1
-150
-100-500
50
Normalized Frequency (?π rad/sample)
M a g n i t u d e (d B )
2、用低通FIR 滤波器对单频噪声进行滤波 %低通滤波器滤波后图形 subplot(211); plot(t,lv);
title('滤波后时域波形');
xlabel('时间t/s');ylabel('幅值/A'); LV=abs(fft(lv)); subplot(212); plot(f,LV);
title('滤波后频谱图');
xlabel('频率 /Hz');ylabel('幅值')
0123
45678910
-2
024
6滤波后时域波形
时间t/s 幅值/A
10
20
30
40
506070
80
90
100
050100150
200滤波后频谱图
频率 /Hz
幅值
3、设计带阻FIR 滤波器 %带阻FIR 滤波器
wlp=0.2*pi;wls=0.4*pi;wus=1.2*pi;wup=0.6*pi; B=wls-wlp; M=ceil(12*pi/B);
wp=[(wls+wlp)/2/pi,(wus+wup)/2/pi]; hn=fir1(M,wp,'stop',blackman(M+1)); liu=filter(hn,1,y2); freqz(hn,1);
00.10.2
0.30.40.50.60.70.80.91
-4000
-3000-2000-1000
Normalized Frequency (?π rad/sample)
P h a s e (d e g r e e s )
0.1
0.2
0.30.40.50.60.70.80.9
1
-150
-100-500
50
Normalized Frequency (?π rad/sample)
M a g n i t u d e (d B )
4、用带阻FIR 滤波器对单频噪声滤波后图形 %带阻FIR 滤波器滤波后图形 subplot(211); plot(t,liu);
title('滤波后时域波形');
xlabel('时间t/s');ylabel('幅值/A'); LIU=abs(fft(liu)); subplot(212); plot(f,LIU);
title('滤波后频谱图');
xlabel('频率 /Hz');ylabel('幅值');
10
20
30
40
506070
80
90
100
050100150
200滤波后频谱图
频率 /Hz
幅值
0123
45678910
-2
2
4
滤波后时域波形
时间t/s 幅值/A
5、设计低通IIR 滤波器 %低通IIR 滤波器
wp=0.4*pi;ws=0.35*pi;rp=1;rs=10; [N,wc]=buttord(wp,ws,rp,rs,'s'); [B,A]=butter(N,wc,'s'); [Bz,Az]=impinvar(B,A); wen=filter(Bz,Az,y2); freqz(Bz,Az);
00.10.2
0.30.40.50.60.70.80.91
-1500
-1000
-500
Normalized Frequency (?π rad/sample)
P h a s e (d e g r e e s )
0.1
0.2
0.30.40.50.60.70.80.9
1
-150
-100-500
50
Normalized Frequency (?π rad/sample)
M a g n i t u d e (d B )
6、用低通IIR 滤波器对单频噪声滤波后图形 %低通IIR 滤波器滤波后图形 subplot(211); plot(t,wen);
title('滤波后时域波形');
xlabel('时间t/s');ylabel('幅值/A'); WEN=abs(fft(wen)); subplot(212); plot(f,WEN);
title('滤波后频谱图');
xlabel('频率 /Hz');ylabel('幅值');
0123
45678910
-2
024
6滤波后时域波形
时间t/s 幅值/A
010203040
5060708090100
50
100
150滤波后频谱图
频率 /Hz
幅值
3.7、心电信号时域波形及频谱的观察和分析
1、低通FIR 滤波器滤波效果 subplot(331); plot(t,b);
title('心电信号的时域波形'); xlabel('时间t/s');ylabel('幅值/A'); subplot(332); plot(f,m);
title('心电信号的频谱图');
xlabel('频率 f/Hz');ylabel('幅值/db'); subplot(333); plot(t,y1);
title('加单频噪声时域波形'); xlabel('时间t/s');ylabel('幅值/A') subplot(334);
plot(f,k);
title('加单频噪声频谱图');
xlabel('频率 /Hz');ylabel('幅值'); subplot(335);
plot(t,y2);
title('加单频,白噪声时域波形'); xlabel('时间t/s');ylabel('幅值/A'); subplot(336);
plot(f,h);
title('加单频,白噪声频谱图');
xlabel('频率 /Hz');ylabel('幅值'); subplot(337);
plot(t,lv);
title('滤单频波后时域波形');
xlabel('时间t/s');ylabel('幅值/A'); subplot(338);
plot(f,LV);
title('滤单频波后频谱图');
xlabel('频率 /Hz');ylabel('幅值');
510-5
5心电信号的时域波形时间t/s
幅值/A
501000
100
200心电信号的频谱图频率 f/Hz 幅值/d b
0510-5
5
加单频噪声时域波形时间t/s
幅值/A
50100
0100
200
加单频噪声频谱图频率 /Hz
幅值
510-50
5加单频,白噪声时域波形时间t/s
幅值/A
50100
0100
200
加单频,白噪声频谱图
频率 /Hz
幅值
510
-5
5滤单频波后时域波形时间t/s
幅值/A
501000
100
200
滤单频波后频谱图频率 /Hz
幅值
2、带阻FIR 滤波器滤波效果 subplot(331); plot(t,b);
title('心电信号的时域波形'); xlabel('时间t/s');ylabel('幅值/A'); subplot(332); plot(f,m);
title('心电信号的频谱图');
xlabel('频率 f/Hz');ylabel('幅值/db'); subplot(333); plot(t,y1);
title('加单频噪声时域波形');
滤波器设计步骤及实现程序
数字滤波器的设计步骤及程序实现 湖南理工学院信息与通信工程学院 一、IIR 脉冲响应不变法设计步骤 1、已知实际数字指标as s ap p ，，，ωω 2、将数字指标化为原型模拟指标As s Ap p ，，，ΩΩ，可设T=pi, T /ω=Ω 3、求原型模拟滤波器的c N Ω，，其中：??? ???ΩΩ--=)/lg(2)]110/()110lg[(10/10/s p A A s p N N A p cp p 210 /1 10 -Ω= Ω N A s cs s 210 /1 10 -Ω= Ω ][cs cp c ΩΩ∈Ω， 4、根据N 写出归一化原型系统函数)(p G a 5、用c s p Ω=/代入得原型系统函数c s p a a p G s H Ω==/)()( 6、将)(s H a 化为部分分式展开形式∑-=k k a s s A s H )( 7、写出)(z H 的极点T s k k e z =，并写出)(z H 的部分分式展开形式∑--?= 11)(z z A T z H k k 8、将)(z H 化为分子分母形式，验证设计结果。 二、IIR 双线性变换法设计步骤 1、已知实际数字指标as s ap p ，，，ωω 2、将数字指标化为原型模拟指标As s Ap p ，，，ΩΩ，可设T=2, 2 tan 2ω?= ΩT 3、求原型模拟滤波器的c N Ω，，其中：?? ? ???ΩΩ--=)/lg(2)]110/()110lg[(10/10/s p A A s p N N A p cp p 210 /1 10 -Ω= Ω N A s cs s 210 /1 10 -Ω= Ω ][cs cp c ΩΩ∈Ω， 4、根据N 写出归一化原型系统函数)(p G a 5、用c s p Ω=/代入得原型系统函数c s p a a p G s H Ω==/) ()( 6、用11 112--+-?=Z Z T s 代入原型系统函数)(s H a 得1 1 112)()(--+-? ==Z Z T s a s H z H 8、将)(z H 整理成分子分母形式，验证设计结果。
数字信号处理课程设计报告
抽样定理的应用 摘要 抽样定理表示为若频带宽度有限的，要从抽样信号中无失真地恢复原信号，抽样频率应大于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)
二阶有源带通滤波器的设计
设计任务书 一、设计目的 掌握二阶压控电压源有源滤波器的设计与测试方法 二、设计要求和技术指标 带通滤波器：通带增益 up A 2；中心频率：0f =1kHz ；品质因数Q=0.707.要求设计电路具有元件少、增益稳定、幅频响应好等特点。 2、设计内容及步骤 （1）写出电路的传递函数，正确计算电路元件参数，选择器件，根据所选器件画出电路原理图，并用multisim 进行仿真。 （2）安装、调试有源滤波电路。 （3）设计实验方案，完成滤波器的滤波性能测试。 （4）画出完整电路图，写出设计总结报告。 三、实验报告要求 1、写出设计报告，包括设计原理、设计电路、选择电路元器件参数、multisim 仿真结论。 2、组装和调试设计的电路检验该电路是否满足设计指标。若不满足，改变电路参数值，使其满足设计题目要求。 3、测量电路的幅频特性曲线。 4、写出实验总结报告。
前言 随着计算机技术的发展，模拟电子技术已经成为一门应用范围极广，具有较强实践性的技术基础课程。电子电路分析与设计的方法也发生了重大的变革，为了培养学生的动手能力，更好的将理论与实践结合起来，以适应电子技术飞速的发展形势，我们必须通过对本次课程设计的理解，从而进一步提高我们的实际动手能力。 滤波器在日常生活中非常重要，运用非常广泛，在电子工程、通信工程、自动控制、遥测控制、测量仪器、仪表和计算机等技术领域，经常需要用到各种各样的滤波器。随着集成电路的迅速发展，用集成电路可很方便地构成各种滤波器。用集成电路实现的滤波器与其他滤波器相比，其波形质量、幅度和频率稳定性等性能指标，都有了很大的提高。 滤波器在电路实验和设备检测中具有十分广泛的用途。现在我们通过对滤波器器的原理以及结构设计一个带通滤波器。我们通过对电路的分析，参数的确定选择出一种最合适本课题的方案。在达到课题要求的前提下保证最经济、最方便、最优化的设计策略。 RC有源滤波器设计 1.1总方案设计 1.1.1方案框图 图1.1.1 RC有源滤波总框图
二阶低通滤波器课程设计报告(昌航版)
课程设计说明书课程设计名称：模拟电子技术课程设计课程设计题目：二阶低通滤波器的设计学院名称：信息工程学院 专业：电子信息工程班级： 学号：： 评分：教师： 20 12 年 3 月日
模拟电子技术 课程设计任务书 20 10 －20 11 学年 第 2 学期 第 1 周－ 3 周 注：1、此表一组一表二份，课程设计小组组长一份；任课教师授课时自带一份备查。 2、课程设计结束后与“课程设计小结”、“学生成绩单”一并交院教务存档。 题目 二阶低通滤波器的设计 容及要求 （1）分别用压控电压源和无限增益多路反馈两种方法设计电路 （2）截止频率kHz 2f p = (3) 增益2A v = 进度安排 第1周：周一至周三查资料，完成原理图设计及仿真； 第1周：周四至第2周周二，完成系统的制作、调试； 第2周：周三设计结果检查。 学生： 指导时间 指导地点： 楼 室 任务下达 20 年 月 日 任务完成 20 年 月 日 考核方式 1.评阅 □ 2.答辩 □ 3.实际操作□ 4.其它 □ 指导教师 系（部）主任
摘要 低通滤波器是一个通过低频信号而衰减或抑制高频信号的部件。理想滤波器电路的频响在通带具有一定幅值和线性相移，而在阻带幅值应为零。有源滤波器是指由放大电路及RC网络构成的滤波器电路，它实际上是一种具有特定频率相应的放大器。滤波器的阶数越高，幅频特性的速率越快，但RC网络节数越多，元件参数计算越繁琐，电路的调试越困难。根据指标，本次设计选用二阶有源低通滤波器。 该电路主要采用了uA741运放，并且在一阶的基础上增加一节RC网络，加大幅频特性衰减斜率，以达到在给定的频段，让信号无衰减的通过电路，而通带外的其他信号将受到很大的衰减，从而提高滤波效率。 关键词：低通滤波器集成运放uA741 RC网络
IIR数字滤波器的设计流程图讲课讲稿
目录 目录 0 前言 (1) 1.1数字滤波器简介 (1) 1.2使用数字滤波器的原因 (1) 1.3设计的原理和内容 (1) 工程概况 (2) 正文 (2) 3.1 设计的目的和意义 (2) 3.2 目标和总体方案 (2) 3.3 设计方法和内容 (3) 3.4 硬件环境 (3) 3.5软件环境 (3) 3.6IIR数字滤波器设计思路 (3) 3.7 IIR数字滤波器的设计流程图 (3) 3.8 IIR数字滤波器设计思路 (4) 3.9设计IIR数字滤波器的两种方法 (4) 3.10双线性变换法的基本原理 (5) 3.11用双线性变换法设计IIR数字滤波器的步骤 (6) 3.12程序源代码和运行结果 (6) 3.12.1低通滤波器 (6) 3.12.3带通滤波器 (10) 3.12.4带阻滤波器 (13) 3.13结论 (15) 3.13.1存在的问题 (15) 3.13.2解决方案 (16) 致谢 (16)
参考文献 (16) 前言 1.1数字滤波器简介 数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。可以设计系统的频率响应，让它满足一定的要求，从而对通过该系统的信号的某些特定的频率成分进行过滤，这就是滤波器的基本原理。如果系统是一个连续系统，则滤波器称为模拟滤波器。如果系统是一个离散系统，则滤波器称为数字滤波器。 信号通过线性系统后，其输出信号就是输入信号和系统冲激响应的卷积。从频域分析来看，信号通过线性系统后，输出信号的频谱将是输入信号的频谱与系统传递函数的乘积。除非为常数，否则输出信号的频谱将不同于输入信号的频谱，某些频率成分较大的模，因此，中这些频率成分将得到加强，而另外一些频率成分的模很小甚至为零，中这部分频率分量将被削弱或消失。因此，系统的作用相当于对输入信号的频谱进行加权。 1.2使用数字滤波器的原因 数字滤波器具有比模拟滤波器更高的精度，甚至能够实现后者在理论上也无法达到的性能。数字滤波器相比模拟滤波器有更高的信噪比。数字滤波器还具有模拟滤波器不能比拟的可靠性。根据其冲击响应函数的时域特性可将数字滤波器分为IIR（有限长冲击响应）和FIR（无限长冲击响应）。 1.3设计的原理和内容 在windows环境下进行语言信号采集，通过IIR数字滤泼器的设计，数字带滤波器就是用软件来实现上面的滤波过程，可以很好的克服模拟滤波器的缺点，数字带滤波器的参数一旦确定，就不会发生变化。IIR型有较好的通带与阻带特性，所以，在一般的设计中选用IIR 型。IIR型又可以分成Butterworth型滤波器，ChebyshevII型滤波器和椭圆型滤波器等。 IIR数字滤波器的设计一般是利用目前已经很成熟的模拟滤波器的设计方法来进行设计，通常采用模拟滤波器原型有butterworth函数、chebyshev函数、bessel函数、椭圆滤波器函数等。 IIR数字滤波器的设计步骤： （1）按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标； （2）根据模拟滤波器技术指标设计为响应的模拟低通滤波器； （3）很据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器；
数字信号处理课设+语音信号的数字滤波
语音信号的数字滤波 ——利用双线性变换法实现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段上，可以通过以下的正切变换实现
二阶带通滤波器课程设计.
一、制作一个1000Hz 的正弦波产生电路： 图1.1 正弦波产生电路 1.1 RC 桥式振荡电路 RC 桥式振荡电路如图（1.1）所示。这个电路由两部分组成，即放大电路和选频网络。其中，R1、C1和R2、C2为串、并联选频网络，接于运算放大器的输出与同相输入端之间，构成正反馈，以产生正弦自激振荡。R3、W R 及R4组成负反馈网络，调节W R 可改变负反馈的反馈系数，从而调节放大的电压增益，使电压增益满足振荡的幅度条件。RC 串并联网络与负反馈中的R3、W R 刚好组成一个四臂电桥，电桥的对角线顶点接到放大器A1的两个输入端，桥式振荡电路的名称即由此得来。 分析RC 串并联网络的选频特性，根椐正弦波振荡电路的振幅平衡条件，选择合适的放大指标，构成一个完整的振荡电路。 1.2 振荡电路的传递函数 由图（1.1）有 1111 Z R sC =+，2 2222 1Z 1R R C sC =+＝2221R sC R + 其中，1Z 、2Z 分别为图1.1中RC 串、并联网络的阻值。 得到输入与输出的传递函数： F ν(s)= 21 2 1212221121()1 sR C R R C C s R C R C R C s ++++ =12 21122111212 11111()s R C s s R C R C R C R R C C ++++ （1.1） 由式（1.1）得 21212 R R 1 C C =ω 2 1210R R 1 C C = ?ω
取1R ＝2R ＝16k Ω,12C C =＝0.01μF ，则有 1.3 振荡电路分析 就实际的频率而言，可用s j ω=替换，在0ωω=时，经RC 选频网络传输到运放同相端的电压与1o U 同相，这样，放大电路和由Z1和Z2组成的反馈网络刚好形成正反馈系统，可以满足相位平衡条件。 12 2 11221212 ()12v j C R F j j C R j C R C C R R ωωωωω= ++- （1.2） 令2 12101R R C C = ω,且R R R C C C ====2121,，则式（1.2）变为 ) (31 )(00ω ωωωω-+= j j F v （1.3） 由此可得RC 串并联选频网络的幅频响应 2 002)( 31ω ωωω-+= V F （1.4） 相频响应 3 )( arctan 0ω ωωω?--=f （1.5） 由此可知，当 2 12101R R C C = =ωω，或CR f f π21 0= = 时，幅频响应的幅度为最大，即 而相频响应的相位角为零，即 这说明，当2 12101R R C C = =ωω时，输出的电压的幅度最大（当输入电压的幅 度一定，而频率可调时），并且输出电压时输入电压的1/3，同时输出电压与输入
课程设计--二阶低通滤波器电路设计及分析
成都理工大学核技术与自动化工程学院 电子技术课程设计 课题名称：二阶低通滤波器电路设计及分析 指导老师：蒋开明 姓名： 学号： 专业：电气工程及其自动化 日期：2010年6月16日
二阶低通滤波器电路设计及分析 一、课题目的： 1、进一步掌握各种滤波电路的工作原理。 2、了解Multisim10的基本操作，并学会用Multisim10进行仿真设计。 3、学会对比并结合理论分析结果进行仿真软件分析。 4、锻炼实际动手能力，增强对课本知识的理解。 二、软件简介： Multisim是美国国家仪器（NI）有限公司推出的以Windows为基础的仿真工具，适用于板级的模拟/数字电路板的设计工作。它包含了电路原理图的图形输入、电路硬件描述语言输入方式，具有丰富的仿真分析能力。 工程师们可以使用Multisim交互式地搭建电路原理图，并对电路进行仿真。Multisim提炼了SPICE仿真的复杂内容，这样工程师无需懂得深入的SPICE技术就可以很快地进行捕获、仿真和分析新的设计，这也使其更适合电子学教育。通过Multisim和虚拟仪器技术，PCB设计工程师和电子学教育工作者可以完成从理论到原理图捕获与仿真再到原型设计和测试这样一个完整的综合设计流程。 Multisim 10是IIT公司推出Multisim 20006年底又发布最新的版本。Multisim 10提供了全面集成化的设计环境，完成从原理图设计输入、电路仿真分析到电路功能测试
等工作。当改变电路连接或改变元件参数，对电路进行仿真时，可以清楚地观察到各种变化对电路性能的影响。 工具栏菜单栏 元器件栏 仪器仪表栏 电路工作区 状态栏 图一Multisim基本界面 三、原理： 由于一阶低通滤波器的幅频特性下降速率只有-20 dB／10 f，与理想情况相差太大，其滤波效果不佳。为了加快下降速率，使其更接近理想状态，提高滤波效果，我们经常使用二阶RC有源滤波器。采取的改进措施是在一阶的基础上再增加一节RC网络。 电路上半部分是一个同相比例放大电路，由两个电阻R1，Rf和一个理想运算放大器构成。R1与Rf均为16 kΩ。下半部分是一个二阶RC滤波电路，由两个电阻R2，R3及两个电容C1，C2构成。其中R2，R3均为4 kΩ，C1，C2均为0.1μF。电路由一个幅度为1 mV，频率可调的交流电压源提供输入信号，用一个阻值为1 kΩ的电阻作为负载，如图二。
数字滤波器设计步骤
数字信号处理 数字滤波器的设计 学院计算机与电子信息学院 专业电子信息科学与技术班级电子１5－2 班姓名学号 指导教师刘利民
数字滤波器的设计 一、模拟低通滤波器的设计方法 1、B ｕtterw ｏrth 滤波器设计步骤： ⑴。确定阶次N ① 已知Ωc 、Ωs 和As 求Bu ｔt ｅr ｗｏrth DF 阶数N ② 已知Ωc 、Ωs 和Ω=Ωp （3dB p Ω≠-)的衰减A ｐ 求Bu ｔterwort ｈ DF 阶数N ③ 已知Ωｐ、Ωｓ和Ω=Ωp 的衰减A ｐ 和As 求B ｕtte ｒwo ｒth DF 阶数Ｎ /10 /1022(/)101,(/)101p s A A N N p c s c ΩΩ=-ΩΩ=-则：
⑵.用阶次N 确定 ()a H s 根据公式： 1,2,2N ()()a a H s H s -在左半平面的极点即为()a H s 的极点，因而 2,,N 2、切比雪夫低通滤波器设计步骤: ⑴.确定技术指标p Ω p α s Ω s α 归一化： /1p p p λ=ΩΩ= /s s p λ=ΩΩ ⑵．根据技术指标求出滤波器阶数N 及ε: 0.12 10 1δε=- p δα= ⑶.求出归一化系统函数 其中极点由下式求出:
或者由N 和Ｓ直接查表得()a H p 二、数字低通滤波器的设计步骤: 1、 确定数字低通滤波器的技术指标：通带截止频率p ω、通带最大衰减系数 p α、 阻带截止频率ω、阻带最小衰减系数s α。 2、 将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。 巴特沃斯： 切比雪夫：/s s p λ=ΩΩ 0.1210 1δ ε=- p δα=
数字信号处理课设共18页文档
数字信号处理课程设计 姓名：刘倩 学号：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 算法的应用
二阶压控型低通滤波器设计
二阶压控型低通滤波器设计 1. 设计要求 设计一个二阶压控型低通滤波器，要求通带增益为2，截止频率为2KHz ，可以选择0.01uF 电容器，阻值尽量接近实际计算值，电路设计完后，画出频率响应曲线，并采用Multisim 软件进行仿真分析。 2. 设计目的 （1） 进一步掌握滤波器电路的工作原理和参数计算。 （2） 熟练使用Multisim 进行简单的电路设计和仿真。 3. 问题分析与参量计算 3.1 问题的简单分析 二阶压控型低通LPF 电路基本原理图可参照教材P345页（如下） 而题目中已经给出了电容的值，故我们所要做的只是确定电阻阻值以及进行电路合理的相关改善。 实验所选取的运放器是a741，实验是在Multisim 环境仿真完成的。 3.2 计算电路相关参数 （1） 低通滤波器在通带将内电容视为开路，给电路引入负反馈从而满足“虚短”、“虚断”，通带增益 3412up R A R =+ =,则34R R =，取34R R == 10k Ω。 （2） 传递函数：为方便计算，取1212,R R R C C C ====，由“虚短”、“虚断”及叠 加定理，得()() ()()() ()()()677776/1()()[()]0up p p p i U s A U s U s U s sCR U s U s U s U s U s U s sC R R ==+-----= 得到传递函数：62()1()()1(3)()u up i up U s A s A U s A sCR sCR ==+-+ 令s j ω=，取012f RC π=，2f ωπ=，2 001(3)()up u up A A f f j A f f ?=+-- （3） 当f 为截止频率时，200|1(3)()|2up f f j A f f +--=,令0f x f =,则得方程 4210x x --=，解得x ，因为2f kHz =，取0.01C F μ=可解得10.1224R k ≈Ω电阻，由于实际试验中难以的到10.1224k Ω的电阻，故实际试验中用10k Ω的电阻代替之 （4）入10,1p V mv f kHz ==的信号源 最终得到的电路图： 3.3二阶压控电压源低通滤波器（LPF ）的幅频特性 Q=13-Aup =13-2 =1 ,所以Q=1的曲线即为此二阶压控电压源低通滤波器（LPF ）的幅频特性。
FIR数字滤波器设计与软件实现(精)讲解学习
实验二:FIR 数字滤波器设计与软件实现 一、实验指导 1.实验目的 (1掌握用窗函数法设计 FIR 数字滤波器的原理和方法。 (2掌握用等波纹最佳逼近法设计 FIR 数字滤波器的原理和方法。 (3掌握 FIR 滤波器的快速卷积实现原理。 (4学会调用 MA TLAB 函数设计与实现 FIR 滤波器。 2. 实验内容及步骤 (1认真复习第七章中用窗函数法和等波纹最佳逼近法设计 FIR 数字滤波器的原理; (2调用信号产生函数 xtg 产生具有加性噪声的信号 xt ,并自动显示 xt 及其频谱,如图 1所示;
图 1 具有加性噪声的信号 x(t及其频谱如图 (3请设计低通滤波器,从高频噪声中提取 xt 中的单频调幅信号,要求信号幅频失真小于 0.1dB ,将噪声频谱衰减 60dB 。先观察 xt 的频谱,确定滤波器指标参数。 (4根据滤波器指标选择合适的窗函数,计算窗函数的长度 N ,调用 MATLAB 函数 fir1设计一个 FIR 低通滤波器。并编写程序,调用 MATLAB 快速卷积函数 fftfilt 实现对 xt 的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 (5 重复 (3 , 滤波器指标不变, 但改用等波纹最佳逼近法, 调用MA TLAB 函数 remezord 和 remez 设计 FIR 数字滤波器。并比较两种设计方法设计的滤波器阶数。 提示:○ 1MA TLAB 函数 fir1的功能及其调用格式请查阅教材; ○ 2采样频率 Fs=1000Hz,采样周期 T=1/Fs;
○ 3根据图 1(b和实验要求,可选择滤波器指标参数:通带截止频率 fp=120Hz,阻带截 至频率 fs=150Hz, 换算成数字频率, 通带截止频率 p 20.24 p f ωπ =T=π, 通带最大衰为 0.1dB , 阻带截至频率 s 20.3 s f ωπ =T=π,阻带最小衰为 60dB 。 3、实验程序框图如图 2所示,供读者参考。 图 2 实验程序框图 4.信号产生函数 xtg 程序清单(见教材 二、滤波器参数及实验程序清单 1、滤波器参数选取 根据实验指导的提示③选择滤波器指标参数: 通带截止频率 fp=120Hz,阻带截至频率 fs=150Hz。代入采样频率 Fs=1000Hz,换算成 数字频率,通带截止频率 p 20.24 p f
数字信号处理课程规划报告
数字信号处理课程设计报告《应用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 课程设计的目的与作用 (1) 2 设计任务及所用multisim软件环境介绍 (1) 2.1 设计任务 (1) 2.2 Multisim软件环境介绍 (1) 3 电路模型的建立 (2) 4 理论分析及计算 (3) 5 仿真结果分析 (4) 6 设计总结和体会 (4) 7 参考文献 (5)
1 课程设计的目的与作用 目的:根据设计任务完成对二阶带通滤波器的设计，进一步加强对模拟电子技术的理解。了解二阶带通滤波器的工作原理，掌握对二阶带通滤波器频率特性的测试方法。 带通滤波器:其作用是允许某一段频带范围内的信号通过，而将此频带以外的信号阻断。常用于抗干扰设备中，以便接收某一段频带范围内的有效信号，而消除高频段和低频段的干扰和噪声。 2 设计任务及所用multisim软件环境介绍 2.1 设计任务 学会使用Multisim10软件设计二阶带通滤波器的电路，使学生初步了解和掌握二阶带通滤波器的设计、调试过程及其频率特性的测试方法，能进一步巩固课堂上学到的理论知识，了解带通滤波器的工作原理。 2.2 Multisim软件环境介绍 Multisim是美国国家仪器（NI）有限公司推出的以Windows为基础的仿真工具，适用于板级的模拟/数字电路板的设计工作。它包含了电路原理图的图形输入、电路硬件描述语言输入方式，具有丰富的仿真分析能力。 工程师们可以使用Multisim交互式地搭建电路原理图，并对电路进行仿真。Multisim 提炼了SPICE仿真的复杂内容，这样工程师无需懂得深入的SPICE技术就可以很快地进行捕获、仿真和分析新的设计，这也使其更适合电子学教育。通过Multisim和虚拟仪器技术，PCB设计工程师和电子学教育工作者可以完成从理论到原理图捕获与仿真再到原型设计和测试这样一个完整的综合设计流程。
-二阶有源低通滤波器设计-
一题目要求与方案论证 1.1 （设计题题目）二阶有源低通滤波器 1.1.1 题目要求 设计二阶有源低通滤波器。要求截止频率 f 0=1000比通带内电压放大倍数A o=15,品 质因数Q=0.707。分析电路工作原理，设计电路图，列出电路的传递函数，正确选择电路中的参数。 1.1.2 方案论证 （1）：对信号进行分析与处理时, 常常会遇到有用信号叠加上无用噪声的问题, 这些噪声有的是与信号同时产生的, 有的是传输过程中混入的。因此, 从接收的信号中消除或减弱干扰噪声, 就成为信号传输与处理中十分重要的问题。根据有用信号与噪声的不同特性, 消除或减弱噪声, 提取有用信号的过程称为滤波, 实现滤波功能的系统称为滤波器。 滤波器分为无源滤波器与有源滤波器两种： ① 无源滤波器: 由电感L、电容C及电阻R等无源元件组成 ② 有源滤波器: 一般由集成运放与RC网络构成，它具有体积小、性能稳定等优点，同 时，由于集成运放的增益和输入阻抗都很高，输出阻抗很低，故有源滤波器还兼有放大与缓冲作用。利用有源滤波器可以突出有用频率的信号，衰减无用频率的信号，抑制干扰和噪声，以达到提高信噪比或选频的目的，因而有源滤波器被广泛应用于通信、测量及控制技术中的小信号处理。 从功能来上有源滤波器分为： 低通滤波器（LPF）、高通滤波器（HPF、 带通滤波器（BPF、带阻滤波器（BEF、 全通滤波器（APF、。 其中前四种滤波器间互有联系，LPF与HPF间互为对偶关系。当LPF的通带截止频率高于HPF的通带截止频率时，将LPF与HPF相串联，就构成了BPF而LPF与HPF 并联，就构成BEF。在实用电子电路中，还可能同时米用几种不同型式的滤波电路。滤波电路的主要性能指标有通带电压放大倍数AVP通带截止频率fP及阻尼系数Q等。 工作原理： 二阶有源滤波器是一种信号检测及传递系统中常用的基本电路, 也是高阶虑波器的基本组成单元。常用二阶有源低通滤波器的电路型式有压控电压源型、无限增益多路反馈型和双二次型。本次课程设计米用压控电压源型设计课题。
数字信号处理课程设计
数字信号处理 课 程 设 计 院系：电子信息与电气工程学院 专业：电子信息工程专业 班级：电信班 姓名： 学号： 组员：
摘要 滤波器设计在数字信号处理中占有极其重要的地位，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)
二阶压控压源型巴特沃斯低通滤波器设计
利用VCVS型二阶RC有源网络实现巴特沃斯型低通滤波器 的设计 一．二阶压控电压源低通滤波器的构成 下图所示就是压控压源二阶型滤波网络电路： 其传递函数为： 与一般低通滤波传输函数相比： 可得： 截至角频率： 增益因子：
选择性因子： 二阶低通滤波器归一化低通传输函数为: 去归一化低通传输函数为: 令： 得： R2应有实根 得：
二．各参数的设计 由于所需的滤波网络阶次为二阶 因为设计指标里通带截至频率规定： f p =100.1KHz,设运放的电压增益为2，而两 个电容的值最好相同，则令 C C C ==21，带入上式品质因 素公式中，可得： 因为品质因素在数值上等于截止频率时的滤波网络电压增益和通带电压增益只比，则 2 1=Q 则 R R R 2212== （1） 因为 2 121121 R R C C f p π = （2） 则由式（1）（2）可求得 R C 1 10125.16 -??= 由实际电子元器件标称值可以设定： 三．结果的验证 利用Multisim 对设计的电路进行仿真。首先搭建整个电路如下： 2 1R R Q =
其中XFG1是信号发生器，XBP1是波特仪，而XSC1是示波器。我们设计的时候所设定的截止频率是100.1K。所以先选择一个比较低的频率值，看其运放的放大倍数。所以先设定信号源频率为1K，仿真结果如下： 示波器示数： 从图中可以看出在低频段时：通道1的峰值为29.98mv，通道2的峰值为62.029mv，滤波网络的放大倍数可以算得A1=2.069。现在把信号源的频率调到预设截至频率，继续仿真，结果如下：
实验五 IIR数字滤波器设计与滤波(附思考题程序)
实验五 IIR 数字滤波器设计与滤波 1.实验目的 （1）加深对信号采样的理解， （2）掌握滤波器设计的方法； （3）复习低通滤波器的设计。 2．实验原理 目前，设计IIR 数字滤波器的通用方法是先设计相应的低通滤波器，然后再通过双线性变换法和频率变换得到所需要的数字滤波器。模拟滤波器从功能上分有低通、高通、带通及带阻四种，从类型上分有巴特沃兹（Butterworth ）滤波器、切比雪夫（Chebyshev ）I 型滤波器、切比雪夫II 型滤波器、椭圆（Elliptic ）滤波器以及贝塞尔（Bessel ）滤波器等。 典型的模拟低通滤波器的指标如下：,P S ΩΩ分别为通带频率和阻带频率，,P S δδ分别为通带和阻带容限（峰波纹值）。在通带内要求1()1P a H J δ-≤Ω≤,有时指标由通带最大衰减p α和阻带最小衰减s α给出，定义如下：20lg(1)p p αδ=-- 和20lg()s s αδ=- 第二种常用指标是用参数ε和A 表示通带和阻带要求，如图所示： 二者之间的关系为：21/2[(1)1]p εδ-=--和1/s A δ=，根据这几个参数可导出另外两个参数d ，k ，分别称为判别因子和选择性因子。 21d A = - /p s k =ΩΩ
BUTTERWORTH 低通滤波器：幅度平方函数定义为221()1(/)a N c H J Ω=+ΩΩ，N 为滤波器阶数，c Ω为截止频率。当c Ω=Ω 时，有()1/a H J Ω=3DB 带宽。 BUTTERWORTH 低通滤波器系统函数有以下形式： 11111()...() N c a N N N N N k H s s a s a s a k s s --=Ω==++++∏- 由模拟滤波器设计IIR 数字滤波器，必须建立好s 平面和z 平面的映射关系。使模拟系统函数()a H s 变换成数字滤波器的系统函数()H z ，通常采用冲激相应不变法和双线性变换法。冲激相应不变法存在频谱混叠现象，双线性变换法消除了这一线象，在IIR 数字滤波器的设计中得到了更广泛的应用。 s 平面和Z 平面的映射关系为1 121()1s Z s f Z T Z ---==+,将s j =Ω和jw z e =待入数字频率和等效的模拟频率之间的映射关系：tan()2 w Ω=，由于二者不是线性关系，所以称为预畸变。 3.实验内容及其步骤 实验的步骤： （1）给定数字滤波器的幅度相应参数。 （2）用预畸变公式将数字滤波器参数变换为相应的等效模拟滤波器参数。 （3）采用模拟滤波器设计方法设计等效模拟滤波器()a H s （4）采用双线性变换公式把等效模拟滤波器映射为所期望的数字滤波器。 其中第三步中模拟滤波器设计步骤为： 首先，根据滤波器指标求选择因子k 和判别因子d 其次，确定满足技术所需的滤波器阶数N, log log d N k ≥ 再次，设3db 截止频率c Ω
二阶带通滤波器设计
物理学院课程设计任务书 专业：学生姓名：学号：学生班级：题目：二阶带通滤波器的设计 指导教师姓名及职称张晓培 电子线路课程设计 题目: 有源带通滤波器 作者姓名：覃万晴 学号：
学院：机械与船舶海洋工程学院 专业：过程控制自动化 指导教师姓名：张晓培 2016年10月1日 二阶带通滤波器的设计 一、设计要求和意义 1）实验要求：中心频率为1KHZ 2）设计意义：近几年随着冶金、化工、纺织机构等工业使用的各种非线性用电设备而产生的大量的高次谐波，已导致电网上网正常波形发生严重畸变，影响到供电系统的电能质量和用户用电设备的安全经济运行。 3）随着生产技术方式的变化，生产力确实得到较大提高，可同时也受到方方面面的限制。如当人们做出了具体的制度设计需要付诸实践进行试验，试验过程中不可避免地会受到一些偶然随即因素的干扰，为评价新方案的效果，需排除这些随即因素的影响，即需要一个滤波器。经滤波以后对新方案的效果进行检验。 4）有源滤波器一般由集成运放与RC网络构成，它具有体积小、性能稳定等优点，同时，由于集成运放的增益和输入阻抗都很高，输出阻抗很低，故有源滤波器还兼有放大与缓冲作用。 5）利用有源滤波器可以突出有用频率的信号，衰减无用频率的信号，抑制干扰和噪声，以达到提高信噪比或选频的目的，因而有源滤波器被广泛应用于通信、测量及控制技术中的小信号处理。 6）若将低通滤波器和高通滤波器串联，并使低通滤波器的通带截止频率fp2大于高通滤波器的通带截止频率fp1，则频率在fp1 二阶低通滤波器课程设计报告昌航版 课程设计说明书 课程设计名称：模拟电子技术课程设计课程设计题目：二阶低通滤波器的设计学院名称：信息工程学院 专业：电子信息工程班级： 学号：姓名： 评分：教师： 20 12 年 3 月日 模拟电子技术 课程设计任务书 20 10 －20 11 年 第 2 学期 第 1 周－ 3 周 注：1、此表一组一表二份，课程设计小组组长一份；任课教师授课时自带一份备查。 2、课程设计结束后与“课程设计小结”、“学生成绩单”一并交题目 二阶低通滤波器的设计 内容及要求 （1）分别用压控电压源和无限增益多路反馈两种方法设计电路 （2）截止频率kHz 2f p = (3) 增益2A v = 进度安排 第1周：周一至周三查资料，完成原理图设计及仿真； 第1周：周四至第2周周二，完成系统的制作、调试； 第2周：周三设计结果检查。 学生姓名： 指导时间 指导地点： 楼 室 任务下达 20 年 月 日 任务完成 20 年 月 日 考核方式 1.评阅 □ 2.答辩 □ 3.实际操作□ 4.其它 □ 指导教师 系（部）主任 院教务存档。 摘要 低通滤波器是一个经过低频信号而衰减或抑制高频信号的部件。理想滤波器电路的频响在通带内具有一定幅值和线性相移，而在阻带内幅值应为零。有源滤波器是指由放大电路及RC网络构成的滤波器电路，它实际上是一种具有特定频率相应的放大器。滤波器的阶数越高，幅频特性的速率越快，但RC网络节数越多，元件参数计算越繁琐，电路的调试越困难。根据指标，本次设计选用二阶有源低通滤波器。 该电路主要采用了uA741运放，而且在一阶的基础上增加一节RC网络，加大幅频特性衰减斜率，以达到在给定的频段内，让信号无衰减的经过电路，而通带外的其它信号将受到很大的衰减，从而提高滤波效率。 关键词：低通滤波器集成运放uA741 RC网络二阶低通滤波器课程设计报告昌航版