信号系统课程设计.docx

信号与系统课程设计 PAGE \* MERGEFORMAT 18 目 录一、实验要求…………………………………………………………………2二、实验原理…………………………………………………………………21、用FFT对信号进行谱分析……………………………………………2 2、模拟低通滤波器的设计………………………………………………2 3、用双线性法设计数字滤波器…………………………………………5 4、用窗函数法设计数字滤波器…………………………………………65、输入信号的取样，滤波器的参数设计………………………………9三、实验结果…………………………………………………………………9 1、输入信号波形……………………………………………………… 10 2、输入信号频谱……………………………………………………… 10 3、低通滤波器幅频响应……………………………………………… 11 4、输出的低频信号波形……………………………………………… 11 5、输出信号的低频分量频谱图……………………………………… 12 6、带通滤波器幅频响应……………………………………………… 12 7、输出信号中频分量波形…………………………………………… 13 8、输出信号的中频分量频谱图……………………………………… 13 9、高通滤波器幅频响应……………………………………………… 14 10、输出信号高频分量波形……………………………………………14 11、输出信号高频分量频谱图…………………………………………15四、实验总结…………………………………………………………………15五、参考文献…………………………………………………………………16六、附录………………………………………………………………………16 一、实验要求产生一个连续信号，包含低频，中频，高频分量，对其进行采样，进行频谱分析，分别设计低通，带通，高通滤波器对信号进行滤波处理，观察滤波前后信号的频谱。二、实验原理 1、用FFT对信号进行谱分析DFT是信号分析与处理中的一种重要变换。因直接计算DFT的计算量与变换区间长度N的平方成正比，当N较大的时，计算量太大，所以有了快速傅立叶FFT的出现。各种语言的FFT子程序都可在有关信号处理的程序库中找到。本实验用到matlab中的函数主要有fft和ifft 。工程实际中，经常遇到连续信号x(t),其频谱函数X(jΩ)也是连续函数，设时域连续信号x(t)持续时间为Tp，最高频率为fc。对x(t)以采样间隔fs≥2fc(fs=1/T)采样得x(t)=x(nT)。设共采样N点，采样间隔为F，F称为频率分辨率。参数fs,Tp,N和F满足如下关系式： F=fs/N=1/NT由于NT=Tp，所以F=1/Tp，说明要提高频率分辨率即减小F就要增加Tp，即增加信号的有效长度。实验给出的信号要求有噪声干扰，可用randn(1,N)函数产生噪声，其中N为采样点数。 2、模拟低通滤波器的设计设计模拟低通滤波器是根据一组技术指标来设计模拟系统函数，使其逼近理想低通滤波器的幅频特性。模拟低通滤波器的幅频特性如图7-1所示。图7-1 模拟低通滤波器的幅频特性模拟低通滤波器的技术指标有、、和。其中，和分别称为通带截止频率和阻带截止频率，是通带（）内的最大衰减系数，是阻带（）内的最小衰减系数。和可分别表示为： 如果处的幅度已归一化到1，即，则和可进一步表示为： (7.1) (7.2)模拟滤波器的技术指标给定后，需要根据这组指标设计模拟系统函数，使其逼近理想滤波器特性，一般是根据幅度平方函数来逼近的。典型的模拟滤波器有巴特沃斯（Butterworth）滤波器、切比雪夫(Chebyshev)滤波器、椭圆（Ellipse）滤波器等。每种典型滤波器都有其不同的特点，下面分别介绍。（1）巴特沃斯原型低通滤波器的设计巴特沃斯低通滤波器的幅度平方函数定义为： (7.3)式中，为正整数，表示滤波器的阶次；为3dB截止频率。Matlab中提供了用于设计巴特沃斯原型低通滤波器的函数，其调用格式如下：[N,Wn]=buttord(Wp,Ws,Rp,Rs,s)：其中，参数Wp和Ws分别是通带边界频率和阻带边界频率，Wp和Ws的单位是rad/s。Rp和Rs分别为通带最大衰减和阻带最小衰减（dB）。返回的参数N和Wn分别为滤波器的阶数和3dB截止频率。对于带通和带阻滤波器，Wp和Ws都是二维向量，向量的第一个元素对应低