信号与系统大作业报告.docVIP

信号与系统大作业报告 课程名称：信号与系统 设计题目：低通、带通和高通三种滤波器的设计 姓名：xxx 学号：xxxxxxxx 专业班级：通信工程xxx班 实验内容 （1）内容描述：输入信号及采样频率设定：输入由三个频率分量组合的载波，比如输入模拟信号为 其中三个频率分量为f1=900Hz, f2=1800Hz, f3=3600Hz, 采样频率为fs=10800Hz （2）基本要求：设计三个滤波器，实现单独提取单个频率分量的功能，即实现低通、带通和高通。 二、实验目的 加深对matlab相关语法以及程序编写的认识与理解。 更好地熟悉和认识低通、带通、高通三种滤波器的原理与功能。 初步建立自我开发matlab程序的意识和概念，增强进一步学习matlab的兴趣。 三、实验原理 我把三个滤波器分成了三个独立的程序分别进行设计，所采用的方法不同。不过三个程序运行结果所呈现出来的图的规格是一样的，都是在一个绘图窗口中显示2*2=4个图，它们分别是原信号时域图、滤波器频域图、滤波后的时域信号图和滤波前后信号频谱图。 低通滤波器 对于低通滤波器，我采用了契比雪夫cheby2滤波器。首先用subplot和plot这两个函数画了原信号时域图；然后设置低通滤波器的截止频率，用cheby2和freqz这两个主要的函数设计IIR滤波器，并求得滤波器的幅频响应，画出滤波器频域图；接着用filter函数对原信号进行滤波，画出滤波后的时域信号图；最后测量滤波前后信号的长度，用fft函数对这两个信号进行傅立叶变换，再用虚线和实现分别画出原信号频谱图和滤波后信号的频谱图。 带通滤波器 对于带通滤波器，我也采用了契比雪夫cheby2滤波器。首先用subplot和plot这两个函数画了原信号时域图；然后设置带通滤波器的截止频率，值得强调的是带通滤波器的截止频率与低通滤波器的截止频率不同，前者是个区间，后者是个点。用cheby2和freqz这两个主要的函数设计IIR滤波器，并求得滤波器的幅频响应，画出滤波器频域图；接着用filter函数对原信号进行滤波，画出滤波后的时域信号图；最后测量滤波前后信号的长度，用fft函数对这两个信号进行傅立叶变换，再用虚线和实现分别画出原信号频谱图和滤波后信号的频谱图。 高通滤波器 对于高通滤波器，我采用了巴特沃斯butter滤波器。首先用subplot和plot这两个函数画了原信号时域图；然后设置高通滤波器的通带和阻带截止频率，用butter和freqz这两个主要的函数设计滤波器，并求得滤波器的幅频响应，画出滤波器频域图；接着用filter函数对原信号进行滤波，画出滤波后的时域信号图；最后测量滤波前后信号的长度，用fft函数对这两个信号进行傅立叶变换，再用虚线和实现分别画出原信号频谱图和滤波后信号的频谱图。 四、主要函数介绍 契比雪夫cheby2滤波器的调用格式：[b a]=cheby2(n,Rs,Wn,ftype) b:滤波器传递函数分子系数向量 a:滤波器传递函数分母系数向量 n:滤波器的阶数 Rs:滤波器阻带衰减 Wn:滤波器截止频率 ftype:滤波器类型，即低通(low)、带通(bandpass)、高通(high)、带阻(stop) 巴特沃斯butter滤波器的调用格式：[b a]=butter(n,Wn,ftype) b:滤波器传递函数分子系数向量 a:滤波器传递函数分母系数向量 n:滤波器的阶数 Wn:滤波器截止频率 ftype:滤波器类型 频率响应函数的调用格式：[h w]=freqz(b,a,n) b:滤波器传递函数分子系数向量 a:滤波器传递函数分母系数向量 n:复频率响应的采样点数 h:滤波器的复频率响应 w:频率点，该w是归一化以后的频率再乘以π 滤波函数的调用格式：xf=filter(b,a,x) b:滤波器传递函数分子系数向量 a:滤波器传递函数分母系数向量 x:被滤波的信号 五、运行结果 低通滤波器 从图上可以看出该低通滤波器已成功把s1=cos(2*pi*900*t)从原信号 s=cos(2*pi*900*t)+cos(2*pi*1800*t)+cos(2*pi*3600*t)提取出来了。 带通滤波器 从图上可以看出该低通滤波器已成功把s2=cos(2*pi*1800*t)从原信号 s=cos(2*pi*900*t)+cos(2*pi*1800*t)+cos(2*pi*3600*t)提取出来了。 高通滤波器 从图上可以看出该低通滤波器已成功把s3=cos(2*pi*3600*t)从原信号 s=cos(2*pi*900*t)+cos(2*pi*1800*t)+cos(2*pi*3600*t)提取出来了。 六、遇到的困难以及相应的解决方法 由于我对于matlab软件应用还不是很熟悉，此次滤波器设计困扰