产生一个连续信号包含低频中频高频分量对其进行采样分别设计三种高通_低通_带通滤波器对信号进行滤波处理
通信原理课程设计——DSB调制解调系统设计与仿真通信原理
通信原理课程设计设计题目：DSB调制解调系统设计与仿真通信原理班级：学生姓名：学生学号：指导老师：目录引言 (3)1、课程设计目的 (3)2、课程设计要求 (3)一、DSB调制解调模型的建立 (4)1、DSB信号的模型 (4)2、DSB信号调制过程分析 (4)3、高斯白噪声信道特性分析 (6)4、DSB解调过程分析 (9)5、DSB调制解调系统抗噪声性能分析 (10)二、仿真过程 (13)三、心得体会 (15)四、参考文献 (15)引言本课程设计用于实现DSB信号的调制解调过程。
信号的调制与解调在通信系统中具有重要的作用。
调制过程是一个频谱搬移的过程，它是将低频信号的频谱搬移到载频位置。
解调是调制的逆过程，即是将已调制的信号还原成原始基带信号的过程。
信号的接收端就是通过解调来还原已调制信号从而读取发送端发送的信息。
因此信号的解调对系统的传输有效性和传输可靠性有着很大的影响。
调制与解调方式往往决定了一个通信系统的性能。
双边带DSB信号的解调采用相干解调法，这种方式被广泛应用在载波通信和短波无线电话通信中。
1、课程设计目的本课程设计是实现DSB的调制解调。
在此次课程设计中，我们将通过多方搜集资料与分析，来理解DSB调制解调的具体过程和它在MATLAB中的实现方法。
预期通过这个阶段的研习，更清晰地认识DSB的调制解调原理，同时加深对MATLAB这款通信仿真软件操作的熟练度，并在使用中去感受MATLAB的应用方式与特色。
利用自主的设计过程来锻炼自己独立思考，分析和解决问题的能力，为我们今后的自主学习研究提供具有实用性的经验。
2、课程设计要求（1）熟悉MATLAB中M文件的使用方法，掌握DSB信号的调制解调原理，以此为基础用M文件编程实现DSB信号的调制解调。
（2）绘制出SSB信号调制解调前后在时域和频域中的波形，观察两者在解调前后的变化，通过对分析结果来加强对DSB信号调制解调原理的理解。
（3）对信号分别叠加大小不同的噪声后再进行解调，绘制出解调前后信号的时域和频域波形，比较未叠加噪声时和分别叠加大小噪声时解调信号的波形有何区别，由所得结果来分析噪声对信号解调造成的影响。
《机械工程测试技术基础》期末试题及答案
第一章 信号及其描述（一）填空题1、 测试的基本任务是获取有用的信息，而信息总是蕴涵在某些物理量之中，并依靠它们来传输的。
这些物理量就是 信号 ，其中目前应用最广泛的是电信号。
2、 信号的时域描述，以 时间 为独立变量；而信号的频域描述，以 频率 为独立变量。
3、 周期信号的频谱具有三个特点： 离散性 ， 谐波性 ， 收敛性 。
4、 非周期信号包括 准周期 信号和 瞬变周期 信号。
5、 描述随机信号的时域特征参数有 均值 、 均方值 、 方差 。
6、 对信号的双边谱而言，实频谱（幅频谱）总是 关于Y 轴 （偶） 对称，虚频谱（相频谱）总是 关于原点（奇） 对称。
（二）判断对错题（用√或×表示）1、 各态历经随机过程一定是平稳随机过程。
（ √ ）2、 信号的时域描述与频域描述包含相同的信息量。
（ √ ）3、 非周期信号的频谱一定是连续的。
（ × ）4、 非周期信号幅频谱与周期信号幅值谱的量纲一样。
（ × ）5、 随机信号的频域描述为功率谱。
（ √ ）（三）简答和计算题1、 求正弦信号t x t x ωsin )(0=的绝对均值μ|x|和均方根值x rms 。
2、 求正弦信号)sin()(0ϕω+=t x t x 的均值x μ，均方值2x ψ，和概率密度函数p(x)。
3、 求指数函数)0,0()(≥>=-t a Ae t x at 的频谱。
4、 求被截断的余弦函数⎩⎨⎧≥<=T t T t t t x ||0||cos )(0ω的傅立叶变换。
5、 求指数衰减振荡信号)0,0(sin )(0≥>=-t a t e t x at ω的频谱。
第二章 测试装置的基本特性（一）填空题1、 某一阶系统的频率响应函数为121)(+=ωωj j H ，输入信号2sin )(t t x =，则输出信号)(t y 的频率为=ω ，幅值=y ，相位=φ 。
2、 试求传递函数分别为5.05.35.1+s 和2224.141n n ns s ωωω++的两个环节串联后组成的系统的总灵敏度。
高通滤波器
1 低通滤波器法
举例：观察有高斯噪声Lenna图像的傅立叶谱和不 同半径下的谱图像的信号能量。
E _ T 1.5387 1015 E _ 5 1.3886 1015 E _ 5 E _ T 0.9025 E _10 1.41911015 E _10 E _ T 0.9223 E _ 20 1.4346 1015 E _ 20 E _ T 0.9323 E _ 50 1.4483 1015 E _ 50 E _ T 0.9412
2 高通滤波器法
-4
4
0.8 0.6 0.4 0.2 -2 -2 0 v u 2 2
-4
4
0.8 0.6 0.4 0.2
-4
-2
0
2 t
4
2 高通滤波器法
3）巴特沃思高通滤波器（BHPF）
n阶巴特沃思（Butterworth）高通滤波器 D0 1 2 2 u v n 11 ， 阶巴特沃思高通滤波器 1 H u, v 2 D0 1 2 2 u v H u, v 1
1 低通滤波器法
有高斯噪声的Lenna图像
D0=5
1 低通滤波器法
D0=10
D0=20
1 低通滤波器法
D0=50
有高斯噪声的原Lenna图像
1 低通滤波器法
问题：
（1）模糊
对于半径为5，包含了全部90%的能量。但严重的 模糊表明了图片的大部分边缘信息包含在滤波器滤 去的10%能量之中。随着滤波器半径增加，模糊的 程度就减少。 模糊产生的原理：根据卷积定理
2 t
4
3 带通和带阻滤波器法
4）巴特沃斯带通滤波器
电子信息工程笔试题
1、基尔霍夫定理的内容是什么？（仕兰微电子）2、平板电容公式(C=εS/4πkd)。
（未知）3、最基本的如三极管曲线特性。
（未知）4、描述反馈电路的概念，列举他们的应用。
（仕兰微电子）5、负反馈种类（电压并联反馈，电流串联反馈，电压串联反馈和电流并联反馈）；负反馈的优点（降低放大器的增益灵敏度，改变输入电阻和输出电阻，改善放大器的线性和非线性失真，有效地扩展放大器的通频带，自动调节作用）（未知）6、放大电路的频率补偿的目的是什么，有哪些方法？（仕兰微电子）7、频率响应，如：怎么才算是稳定的，如何改变频响曲线的几个方法。
（未知）8、给出一个查分运放，如何相位补偿，并画补偿后的波特图。
（凹凸）9、基本放大电路种类（电压放大器，电流放大器，互导放大器和互阻放大器），优缺点，特别是广泛采用差分结构的原因）。
（未知）10、给出一差分电路，告诉其输出电压Y+和Y-,求共模分量和差模分量。
（未知）11、画差放的两个输入管。
（凹凸）12、画出由运放构成加法、减法、微分、积分运算的电路原理图。
并画出一个晶体管级的运放电路。
（仕兰微电子）13、用运算放大器组成一个10倍的放大器。
（未知）14、给出一个简单电路，让你分析输出电压的特性（就是个积分电路），并求输出端某点的rise/fall时间。
(Infineon笔试试题)15、电阻R和电容C串联，输入电压为R和C之间的电压，输出电压分别为C上电压和R上电压，要求制这两种电路输入电压的频谱，判断这两种电路何为高通滤波器，何为低通滤波器。
当RC<<T时，给出输入电压波形图，绘制两种电路的输出波形图。
（未知）16、有源滤波器和无源滤波器的原理及区别?（新太硬件）17、有一时域信号S=V0sin(2pif0t)+V1cos(2pif1t)+V2sin(2pif3t+90),当其通过低通、带通、高通滤波器后的信号表示方式。
（未知）18、选择电阻时要考虑什么？（东信笔试题）19、在CMOS电路中，要有一个单管作为开关管精确传递模拟低电平，这个单管你会用P管还是N管，为什么？（仕兰微电子）20、给出多个mos管组成的电路求5个点的电压。
通信原理简答题答案2（个人整理）
通信原理简答题答案2（个⼈整理）第⼀章绪论1-2何谓数字信号？何谓模拟信号？两者的根本区别是什么？答：数字信号：电信号的参量值仅可能取有限个值。
模拟信号：电信号的参量取值连续。
两者的根本区别是携带信号的参量是连续取值还是离散取值。
1-3何谓数字通信？数字通信偶哪些优缺点？答：利⽤数字信号来传输信息的通信系统为数字通信系统。
优点：抗⼲扰能⼒强，⽆噪声积累传输差错可控；便于现代数字信号处理技术对数字信息进⾏处理、变换、储存；易于集成，使通信设备微型化，重量轻；易于加密处理，且保密性好。
缺点：⼀般需要较⼤的传输带宽；系统设备较复杂。
1-4 数字通信系统的⼀般模型中各组成部分的主要功能是什么？答：信源编码：提⾼信息传输的有效性（通过数字压缩技术降低码速率），完成A/D转换。
信道编码/译码：增强数字信号的抗⼲扰能⼒。
加密与解密：认为扰乱数字序列，加上密码。
数字调制与解调：把数字基带信号的频谱搬移到⾼频处，形成适合在信道中传输的带通信号。
同步：使收发两端的信号在时间上保持步调⼀致。
1-5 按调制⽅式，通信系统如何分类？答：基带传输系统和带通传输系统。
1-6 按传输信号的特征，通信系统如何分类？答：模拟通信系统和数字通信系统。
1-7 按传输信号的复⽤⽅式，通信系统如何分类？答：FDM,TDM,CDM。
1-8 单⼯、半双⼯及全双⼯通信⽅式是按什么标准分类的？解释他们的⼯作⽅式。
答：按照消息传递的⽅向与时间关系分类。
单⼯通信：消息只能单向传输。
半双⼯：通信双⽅都能收发消息，但不能同时进⾏收和发的⼯作⽅式。
全双⼯通信：通信双⽅可以同时收发消息。
1-9 按数字信号码元的排列顺序可分为哪两种通信⽅式？他们的适⽤场合及特点？答：分为并⾏传输和串⾏传输⽅式。
并⾏传输⼀般⽤于设备之间的近距离通信，如计算机和打印机之间的数据传输。
串⾏传输使⽤与远距离数据的传输。
1-10 通信系统的主要性能指标是什么？—答：有效性和可靠性。
1-11 衡量数字通信系统有效性和可靠性的性能指标有哪些？答：有效性：传输速率，频带利⽤率。
数字信号处理讲义第7章滤波器的设计方法
第７章滤波器的设计方法教学目的1．掌握由连续时间滤波器设计离散时间IIR滤波器的方法，包括冲激响应不变法，双线性变换法等；2．了解常用的窗函数，掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法；3．掌握FIR滤波器的逼近原理与设计方法。
教学重点与难点重点：本章是本课程的重中之重，滤波器的设计是核心内容之一。
1．连续时间滤波器设计离散时间IIR滤波器的方法，包括冲激响应不变法，双线性变换法等；2．常用的窗函数，掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法；3．掌握FIR滤波器的逼近原理与设计方法。
难点：1.冲激响应不变法，双线性变换法2.用窗函数法设计FIR滤波器FIR滤波器的逼近原理与设计方法基本概念7.0.1 选频滤波器的分类数字滤波器是数字信号处理的重要基础。
在对信号的过滤、检测与参数的估计等处理中, 数字滤波器是使用最广泛的线性系统。
数字滤波器是对数字信号实现滤波的线性时不变系统。
它将输入的数字序列通过特定运算转变为输出的数字序列。
因此，数字滤波器本质上是一台完成特定运算的数字计算机。
我们已经知道，一个输入序列x(n)，通过一个单位脉冲响应为h(n)的线性时不变系统后，其输出响应y(n)为∑∞-)(y))()()(n(nn=m*=xmhnhx将上式两边经过傅里叶变换，可得式中，Y (e j ω)、X (e j ω)分别为输出序列和输入序列的频谱函数， H (ejω)是系统的频率响应函数。
可以看出，输入序列的频谱X (e j ω)经过滤波后，变为X (e j ω)H (e j ω)。
如果|H (e j ω)|的值在某些频率上是比较小的，则输入信号中的这些频率分量在输出信号中将被抑制掉。
因此，只要按照输入信号频谱的特点和处理信号的目的，适当选择H (ej ω)，使得滤波后的X (e j ω)H (e j ω)符合人们的要求，这就是数字滤波器的滤波原理。
和模拟滤波器一样，线性数字滤波器按照频率响应的通带特性可划分为低通、高通、带通和带阻几种形式。
基于MATLAB的信号分析与处理
山东建筑大学课程设计说明书题目：基于MATLAB的信号分析与处理课程：数字信号处理课程设计院（部）：信息与电气工程学院专业：通信工程班级：通信学生姓名：学号：指导教师：完成日期：目录目录 (1)摘要 (2)正文 (3)1设计目的和要求 (3)2设计原理 (3)3设计内容 (4)3.1源程序代码 (4)3.2程序执行的结果．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． (7)3.3调试分析过程描述 (12)3.4结果分析 ................................... 错误！未定义书签。
总结与致谢 (14)参考文献 (15)摘要随着科学技术的飞速发展，人们对信号的要求越来越高。
然而，学好《数字信号处理》这门课程是我们处理信号的基础。
MATLAB是一个处理信号的软件，我们必须熟悉它的使用。
本次课程设计利用MATLAB软件首先产生成低频、中频、高频三种频率信号，然后将三种信号合成为连续信号，对连续周期信号抽样、频谱分析，并设计低通、带通、高通三种滤波器对信号滤波，观察滤出的信号与原信号的关系，并分析了误差的产生，通对数字信号处理课程的理论知识的综合运用。
从实践上初步实现对数字信号的处理。
关键词：MATLAB；连续信号；采样定理；滤波器；频谱分析；正文1设计目的和要求（1）、产生一个连续信号，该信号中包含有低频、中频、高频分量，对其进行采样，用MATLAB绘制它们的时域波形和频域波形，对其进行频谱分析；（2）、根据信号频谱分析的结果，分别设计合适的低通、带通、高通滤波器，用MATLAB绘制其幅频及相频特性图；（3）、用所设计的滤波器对信号进行滤波处理，对滤波后的信号进行FFT 频谱分析，用MATLAB绘制处理过程中的各种波形及频谱图，比较滤波前后的时域波形及频谱，对所得结果和滤波器性能进行分析，阐明原因，得出结论；（4）学会使用MATLAB对信号进行分析和处理；2设计原理理论上信号的采样要符合奈奎斯特采样定律，就是采样频率要高一点，一般为被采信号最高频率的2倍，只有这样，才能保证频域不混叠，也就是采样出来数字信号中包含了被采信号的所有信息，而且没有引入干扰。
《现代信号处理》课程设计任务书(09级)
中南大学本科生课程设计任务书课程名称现代信号处理指导教师赵亚湘学院信息科学与工程学院专业班级通信工程0901-0905班0()()sin()()anTa x n x nT Ae nT u nT -==Ω中 南 大 学课程设计任务书一、课程设计目的：1．全面复习课程所学理论知识，巩固所学知识重点和难点，将理论与实践很好地结合起来。
2．提高综合运用所学知识独立分析和解决问题的能力； 3．熟练使用一种高级语言进行编程实现。
二、课程设计内容1. 给定模拟信号：)()sin()(0t u t Ae t at a x Ω=-，式中128.444=A，α=，s rad /2500π=Ω。
对()a t x 进行采样，可得采样序列1） 选择采样频率s f =1 kHz ，观测时间50=p T ms ，观测所得序列()x n 及其幅频特性|()|jw X e 2） 改变采样频率s f =300Hz ，观测此时|()|jw X e 的变化 3） 令采样频率s f =200Hz ，观测此时|()|jw X e 的变化要求分析说明原理，绘出相应的序列及其它们对应的幅频特性曲线，指出|()|jw X e 的变化，说明为什么？2. 已知Gaussian 序列固定序列()x n 中的参数p=8，令q 分别等于2，4，8，观察它们的时域和幅频特性，了解当q 取不同值时，对信号序列的时域及幅频特性的影响；固定q=8，令p 分别等于8，13，14，观察参数p 变化对信号序列的时域及幅频特性的影响，观察p 等于多少时，会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。
并理论分析说明产生这些变化的原因 3. 一个连续信号含两个频率分量，经采样得2(),015()0,n p q en x n --⎧⎪≤≤=⎨⎪⎩其它x(n)=sin2π*0.125n+cos2π*(0.125+Δf)n n=0,1……,N -1已知N=16，Δf 分别为1/16和1/64，观察其幅频特性；当N=128时，Δf 不变，其结果有何不同，为什么？分析说明原因，并打印出相应的幅频特性曲线4. 产生一个淹没在噪声中的信号()x t ，例如由50Hz 和120Hz 的正弦信号以及一个零均值的随机噪声叠加而成。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、题目：产生一个连续信号，包含低频，中频，高频分量，对其进行采样，进行频谱分析，分别设计三种高通，低通，带通滤波器对信号进行滤波处理，观察滤波后信号的频谱。
二、设计内容：1.设计要求:编程生成连续的时域信号,对其进行采样,分析幅频曲线。
分析频谱后,根据自己生成的信号,设计对应的低通滤波器,带通滤波器,高通滤波器,将产生的信号通过滤波器,得到对应频率的信号,然后分析对应频率的频谱,和预期的进行对比,然后调试,设计出较为经济且又能满足要求的滤波器2.设计原理:滤波器的设计有多种方法,可以设计IIR滤波器或FIR滤波器，其中IIR滤波器可以根据巴特沃斯滤波器或契比雪夫滤波器的原理来设计；FIR滤波器可以根据窗函数法、频率采样法和契比雪夫逼近法等等。
实验中有涉及到切比雪夫滤波器的设计。
这里我选用了巴特沃斯滤波器的设计原理和和窗函数法的设计原理。
3.滤波器设计产生一个连续信号，包含低频，中频，高频分量，对其进行采样，进行频谱分析，分别设计低通，带通，高通滤波器对信号进行滤波处理，观察滤波前后信号的频谱。
三、实验程序及结果分析：1.输入信号产生程序如下：clear,clcfigure(1)f1 = 100;f2 = 1000;f3 = 2000;fs = 6000;t = 0:1/fs:0.1;ft = sin(2*pi*f1*t)+3*sin(2*pi*f2*t)+cos(2*pi*f3*t)+0.3*randn(size(t)); subplot(211)plot(t,ft)title('输入信号的波形')xlabel('时间t');ylabel('幅度');F = fft(ft);OMEGA = 0:fs/length(F):fs-fs/length(F);subplot(212)plot(OMEGA,abs(F))title('输入信号的频谱');xlabel('频率');ylabel('幅度');n = 0:length(t)-1;w1 = 2*pi*f1/fs;w2 = 2*pi*f2/fs;w3 = 2*pi*f3/fs;fn = sin(w1*n)+3*sin(w2*n)+cos(w3*n)+0.3*randn(size(n)); pausefigure(2)plot(n,fn)title('离散后的输入信号')xlabel('时间n');ylabel('幅度');分析：这次实验我用的输入信号是含有三种频率成分的信号，频率和幅度分别是100Hz，1000Hz，2000Hz和1V，3V，1V，采样频率为6000Hz。
从得到的波形上看，该波形比较不平滑，这是因为在产生该信号时对该信号加0.3*rand(size(t）)的高斯噪声。
对其频谱分析,同样可以看出其频谱有很多毛刺，这是由于高斯噪声所引起的.2.低通滤波器及其单位脉冲响应频谱和滤出的信号的波形及其频谱，程序如下：figure(3)wp1 = 2.1*w1;ws1 = 1.0*w2;wn1 = (wp1+ws1)/2;N1 = ceil(8*pi/wn1);n1=0:N1-1;wc1=(ws1+wp1)/2;a=(N1-1)/2;m=n1-a+eps;hd1=sin(wc1*m)./(pi*m);window1 = (bartlett(N1))';h1=hd1.*window1;w=0:0.01:pi;H1=freqz(h1,1,w);dbH1=20*log10(abs(H1)/max(abs(H1)));subplot(311)plot(dbH1);title('FIR低通滤波器的单位冲击响应')f1 = filter(h1,1,ft);subplot(312)plot(t,f1);title('低频波的波形');Hw1 = fft(f1);subplot(313)plot(OMEGA,abs(Hw1))title('低频波的频谱')分析：低通滤波器是采用的是FIR滤波器的中的窗函数法设计的，同代截止频率取的是2.1ω1，阻带截止频率是0.65ω2，通带衰减αp=3db，αs=25。
观察低频信号的波形，可以看出，滤出的波还是比较好的，当然还可以将阻带衰减频率变得更小一些，将衰减系数变得更大一些。
观察滤波的效果还是很好地，对于频谱，滤出的频谱还是很干净的。
而波形上的毛刺，那是由于高斯噪声带来的影响，这个是无法完全消除的。
3.带通滤波器及其单位脉冲响应频谱和滤出的信号的波形及其频谱，其程序如下：figure(4)wp21 = 0.7*w2; wp22 = 1.6*w2;ws21 = 0.55*w2;ws22 = 1.45*w2;wn2 = wp21-ws21;N2 = ceil(8*pi/wn2);n2 = 0:N2-1;wc21 = (wp21+ws21)/2;wc22 = (wp22+ws22)/2;a = (N2-1)/2;m2 = n2-a+eps;hd2=sin(wc22*m2)./(pi*m2)-sin(wc21*m2)./(pi*m2);window2 = (hamming(N2))';h2 = hd2.*window2;H2 = freqz(h2,1,w);dbH2=20*log10(abs(H2)/max(abs(H2)));subplot(311)plot(dbH2)title(' FIR带通滤波器的单位取样响应')f2 = filter(h2,1,ft);subplot(312);plot(t,f2);title('中频波的波形');subplot(313);plot(OMEGA,abs(fft(f2)));title('中频波频谱');分析：带通滤波器的设计也是采用FIR滤波器的中的窗函数法设计的，这里的上通带截止频率为0.7ω2,，下通带截止频率为1.6ω2，上阻带截止频率选的是0.55ω2，下阻带截止频率为1.45ω2。
这里的通带宽带为0.9ω2，阻带宽度为0.15ω2，原则上选用其他的窗函数也是可以的，但是事实上，选用三角窗或是矩形窗，设计的滤波器，其效果不好，滤出来的波形不是很理想，因此用的是汉明窗，滤出的效果如图三的第二个图。
可以看出其波形是个很好地正弦波，其频谱如第三个图，频率成分比较的单一，说明滤波器的设计还是很成功的。
4.高通滤波器及其单位脉冲响应频谱和滤出的信号的波形及其频谱，其程序如下：figure(5)wp3 = 0.85*w3;ws3 = 0.7*w3;wn3 = wp3-ws3;N3 = ceil(12*pi/wn3);n3 = 0:N3-1;wc3 = (wp3+ws3)/2;a = (N3-1)/2;m3 = n3-a+eps;window3 = (blackman(N3))';hd3 = sin(pi*m3)./(pi*m3)-sin(wc3*m3)./(pi*m3);h3 = hd3.*window3;H3=freqz(h3,1,w);dbH3=20*log10(abs(H3)/max(abs(H3))); subplot(311)plot(w,dbH3)f3 = filter(h3,1,ft);title('FIR高通滤波器的单位采样响应') subplot(312)plot(t,f3)title('高频波的波形')subplot(313)plot(OMEGA,abs(fft(f3)))title('高频波的频谱')通带截止频率为0.85ω3，阻带截止频率选的是0.7ω3。
阻带宽度为0.15ω3，这里选用的是布莱克曼窗。
这里选用的布莱克曼窗是为了有较大的阻带，观察图形，其波形相对来说还是比较好的，其频谱如第三个图，频率成分比较的单一但是还是有一些纹波，这也是由于高斯噪声产生的，但总体的波形还是比较的接近正弦波。
四、实验总结：通过这次试验加深了我对确定性信号和系统的分析方法、相关算法、系统实现等相关知识的理解以及把握，借助于数字滤波器的设计及实现，掌握数字系统的分析以及设计方法。
通过本次课程设计使我们把该课程的理论与工程应用的紧密结合,进一步理解了信号处理的内涵和实质。
在整个过程中还是碰到不少问题的，从刚开始的没有思路，经过交流以及查阅相关知识和课本内容，再结合以前做过的信号与系统的实验，渐渐找到了方法。
总的来说，在整个实验过程中对信号与系统这个课程又有了进一步的认识，把以前学的一些知识都融会贯通的用了起来，受益不少。
对以前学懂的知识有了更深一步的巩固，对一些不大明白的地方也又有了新的认识，把理论与实践结合起来了。
附录：完整源代码%1.输入信号产生程序如下：clear,clcfigure(1)f1 = 100;f2 = 1000;f3 = 2000;fs = 6000;t = 0:1/fs:0.1;ft = sin(2*pi*f1*t)+3*sin(2*pi*f2*t)+cos(2*pi*f3*t)+0.3*randn(size(t));subplot(211)plot(t,ft)title('输入信号的波形')xlabel('时间t');ylabel('幅度');F = fft(ft);OMEGA = 0:fs/length(F):fs-fs/length(F);subplot(212)plot(OMEGA,abs(F))title('输入信号的频谱');xlabel('频率');ylabel('幅度');n = 0:length(t)-1;w1 = 2*pi*f1/fs;w2 = 2*pi*f2/fs;w3 = 2*pi*f3/fs;fn = sin(w1*n)+3*sin(w2*n)+cos(w3*n)+0.3*randn(size(n));pausefigure(2)plot(n,fn)title('离散后的输入信号')xlabel('时间n');ylabel('幅度');pause%2.低通滤波器及其单位脉冲响应频谱和滤出的信号的波形及其频谱，程序如下：figure(3)wp1 = 2.1*w1;ws1 = 1.0*w2;wn1 = (wp1+ws1)/2;N1 = ceil(8*pi/wn1);n1=0:N1-1;wc1=(ws1+wp1)/2;a=(N1-1)/2;m=n1-a+eps;hd1=sin(wc1*m)./(pi*m);window1 = (bartlett(N1))';h1=hd1.*window1;w=0:0.01:pi;H1=freqz(h1,1,w);dbH1=20*log10(abs(H1)/max(abs(H1)));subplot(311)plot(dbH1);title('FIR低通滤波器的单位冲击响应')f1 = filter(h1,1,ft);subplot(312)plot(t,f1);title('低频波的波形');Hw1 = fft(f1);subplot(313)plot(OMEGA,abs(Hw1))title('低频波的频谱')pause%3.带通滤波器及其单位脉冲响应频谱和滤出的信号的波形及其频谱，其程序如下：figure(4)wp21 = 0.7*w2; wp22 = 1.6*w2;ws21 = 0.55*w2;ws22 = 1.45*w2;wn2 = wp21-ws21;N2 = ceil(8*pi/wn2);n2 = 0:N2-1;wc21 = (wp21+ws21)/2;wc22 = (wp22+ws22)/2;a = (N2-1)/2;m2 = n2-a+eps;hd2=sin(wc22*m2)./(pi*m2)-sin(wc21*m2)./(pi*m2);window2 = (hamming(N2))';h2 = hd2.*window2;H2 = freqz(h2,1,w);dbH2=20*log10(abs(H2)/max(abs(H2)));subplot(311)plot(dbH2)title(' FIR带通滤波器的单位取样响应')f2 = filter(h2,1,ft);subplot(312);plot(t,f2);title('中频波的波形');subplot(313);plot(OMEGA,abs(fft(f2)));title('中频波频谱');pause%4.高通滤波器及其单位脉冲响应频谱和滤出的信号的波形及其频谱，其程序如下：figure(5)wp3 = 0.85*w3;ws3 = 0.7*w3;wn3 = wp3-ws3;N3 = ceil(12*pi/wn3);n3 = 0:N3-1;wc3 = (wp3+ws3)/2;a = (N3-1)/2;m3 = n3-a+eps;window3 = (blackman(N3))';hd3 = sin(pi*m3)./(pi*m3)-sin(wc3*m3)./(pi*m3); h3 = hd3.*window3;H3=freqz(h3,1,w);dbH3=20*log10(abs(H3)/max(abs(H3)));subplot(311)plot(w,dbH3)f3 = filter(h3,1,ft);title('FIR高通滤波器的单位采样响应')subplot(312)plot(t,f3)title('高频波的波形')subplot(313)plot(OMEGA,abs(fft(f3)))title('高频波的频谱')。