《matlab及数字信号处理课程设计报告》由会员分享，可在线阅读，更多相关《matlab及数字信号处理课程设计报告（25页珍藏版）》请在金锄头文库上搜索。
1、数字信号处理专业课程设计报告书 实验报告题目四 :Using the bilinear transformation and a lowpass analog Butterworth prototype, design a highpass digital filter operating at a rate of 20kHz and having passband extending to 5kHz with maximum passband attenuation of 0.5dB,and stopband ending at 4kHz with a minimum stopband at
2、tenuation of 10dB.备注：题目3,4要求：实验报告中要求写出对应的滤波器H(z)，并在H(z)表达式中将共轭极点对组成二阶基本节，以极点在Z平面上分布顺序写出H(z)形式并将各二阶基本节系数以顺序列表。画出幅度频谱图的|H()|及其以(dB)为单位的幅度谱图。二：实验目的1）熟练掌握低通滤波器的设计方法。2）学会利用低通滤波器设计高通滤波器。3）掌握用双线形变换法设计数字高通滤波器的方法。4）熟悉MATLAB提供的各种滤波器设计函数。5）掌握各种关于滤波器的幅度频谱设计函数。三：实验原理本题利用双线性变换法和巴特沃斯低通滤波器来设计数字高通滤波器：双线形变换法是利用s=2*（1
3、-z-1）/T*（1+z-1）将s域转换到z域，从而得到系统函数H（Z）。根据所要设计要求将高通数字滤波器指标转化为低通模拟滤波器技术指标，主要利用双线性变换式=2/ Ttan(W/2)。滤波器设计中主要用到的函数：Buttord函数用来选择巴特沃斯滤波器最小阶数，调用方式如下：n,wn=buttord(Wp,Ws,rp,rs,s) ：返回符合要求性能指标的数字滤波器最小阶数n和巴特沃斯滤波器截止频率wn； n,wn=buttord(Wp,Ws,rp,rs)：（同上）此处Wp,Ws都是归一化频率。z,p,k=buttap(n)：返回设计的巴特沃斯滤波器的零点（z），极点（p）和增益（k），n为
4、滤波器阶数。b,a=zp2tf(z,p,k)：零极点增益滤波器参数转换为传输函数形式，b，a分别为传输函数的分子分母。bt,at=lp2hp(b,a,wn)：模拟低通滤波器参数转化为模拟高通滤波器参数。bd,ad=bilinear(bt,at,Fs); %模拟高通滤波器参数转化为数字高通滤波器参数。Fs为采样频率。sos,g=tf2sos(bd,ad,order)：将传递函数模型转化为二次分式模型sos。order指定sos中行的顺序： Up：首行中所包含的极点离原点最近，离单位圆最远； Down：首行中所包含的极点离原点最远，离单位圆最近；与画图有关的函数：（1）Zplane(b,a)：绘制
5、系统零极点图；（2）求解数值滤波器频率响应函数：1）freqz(b,a,n)：无输出参数，直接在当前命令窗口绘制频率响应的幅频响应（dB形式的）和相频响应。2）H,W=freqz(b,a,n)：返回数字滤波器的n点复频率响应。四：实验步骤简述 （1）先将数字高通滤波器技术指标转化为模拟低通滤波器技术指标。（2）确定滤波器最小阶数.（3）确定零极点增益。（4）由零极点增益确定模拟低通滤波器传输函数。（5）将模拟低通滤波器传输函数转换成模拟高通滤波器传输函数。（6）利用双线性变法将模拟高通滤波器转换成数字滤波器。（7）利用传递函数模型求其二次分式模型sos。（8）绘制相关幅度频谱图。五：程序框图
6、六：源程序：clc; close; clear;Fs=20000;T=1/Fs;%数字高通技术指标转化为模拟低通技术指标 fs=4000; fp=5000;rp=0.5; rs=10;ws=fp*2*pi*T; wp=fs*2*pi*T;Wp=2/T*tan(wp/2); Ws=2/T*tan(ws/2);%*滤波器设计*N,wn=buttord(Wp,Ws,rp,rs,s)%根据技术指标选择数字滤波器最小阶数，巴特沃斯截止频率z,p,k=buttap(N);%返回设计的巴特沃斯滤波器的零点极点和增益b,a=zp2tf(z,p,k);%零极点增益滤波器参数转换为传输函数形式bt,at=lp2h
7、p(b,a,wn);%模拟低通转化为模拟高通bd,ad=bilinear(bt,at,Fs); %模拟高通转化为数字高通%模拟低通和数字高通二阶基本节系数顺序表sos,g=tf2sos(b,a,up)sos,g=tf2sos(bd,ad,up)%*相关设计图*%模拟低通和数字高通零极点图figure(1); subplot(121);zplane(b,a); gtext(模拟低通零极点图);set(gcf,defaultAxesLinestyleOrder,-|:);set(gcf,defaultAxesColorOrder,1 0 0;0 0 1);%按先后顺序对图形线条进行颜色设置subp
8、lot(122);zplane(bd,ad);gtext(数字高通零极点图);pause; figure(2);%滤波器幅频特性图freqz(bd,ad,100);axis(0 1 -40 0);pause; figure(3);H,W=freqz(bd,ad,100);%求出数字滤波器幅频特性plot(W*Fs/2/pi,abs(H);axis(0 8000 0 1.5);xlabel( 频率 f/HZ ); ylabel(归一化幅度 );title(滤波器); grid on;pause;close all;七：程序结果及图表滤波器阶数N=7数字高通滤波器二阶基本节系数按“up（首行中所包
9、含的极点离原点最近，离单位圆最远）”排列:sos = 1.0000 -0.9900 0 1.0000 -0.0783 0 1.0000 -2.0181 1.0182 1.0000 -0.1647 0.0582 1.0000 -2.0044 1.0045 1.0000 -0.1927 0.23771.0000 -1.9874 0.9876 1.0000 -0.2552 0.6396g =0.0279分子系数1.0000 -0.9900 0 1.0000 -2.0181 1.0182 1.0000 -2.0044 1.00451.0000 -1.9874 0.9876分母系数:1.0000 -0.
10、0783 0 1.0000 -0.1647 0.0582 1.0000 -0.1927 0.2377 1.0000 -0.2552 0.6396图（1）：图（2）：以dB为单位的幅度谱图及相位频率图： 图（3）：数字高通滤波器幅度频谱图：八：实验总结 在本次试验中，我熟练的掌握了用模拟低通巴特沃斯滤波器设计数字高通滤波器的方法。本次设计重在设计思维的建立及其所用函数的熟练掌握。在进行具体设计时，主要抓住“给定条件”，“目的要求”展开设计构架。根据“目的要求”：用模拟低通巴特沃斯滤波器设计数字高通滤波器；可知主要设计出传输函数H(Z)。再结合“给定条件”：数字高通滤波器技术指标；设计出系统框图。
11、最终应用相关函数编辑出具体代码，在进行调试。在本次试验中我遇到一个小问题：二阶基本节系数矩阵sos为什么是“4行3列”，而不是“1行3列”，最终我找到了答案：在本题中二次分式形式的传输函数H(Z)是由4个二阶基本节串联而成，因此是“4行3列”。本次实验我掌握了利用双线性变换法设计滤波器的方式方法，相对成功的设计出要求的高通数字滤波器，达到了实验要求及实验目的。题目七：该题目的目的是说明一个PN扩频信号在抑制正弦干扰中的有效性。现考虑下图所示的二进制通信系统，对信号发生器的输出乘上一个二进制（1）PN序列。同一个二进制PN序列用来与解调器输入相乘，因此消除了这个PN序列在期望信号上的影响。信道将
12、传送信号受到一宽带加性噪声序列(n)和一正弦干扰序列i（n）=Asin0n, 00的污损。可以假定AM,这里M是来自调制器的每比特的样本数。用和不用PN序列执行这个仿真系统，并在条件AM下，对不同的M值（如M50，100，500，1000）测量差错率。画出每种情况的差错率曲线，作比较并说明结论。说明这个PN序列对于正弦干扰信号的效果。由此说明为什么PN扩频系统在正弦干扰信号存在下优于常规的二进制通信系统。二：实验目的（1）利用matlab实现一个简单的通信系统。（2）了解扩频序列及其实现方法。 （3）在一个简单的通信系统，信号通过一个有噪声的系统，对加上PN序列和不加PN序列后的差错率比较，说
13、明PN序列能有效的抑制正弦干扰信号，即PN扩频系统在正弦干扰信号存在下优于常规的二进制通信系统。三：实验原理 原始信号S(n)与PN扩频序列相乘（异或实现）后，叠加正弦信号噪声，再加一般噪声后与PN序列相乘恢复出经过系统后的信号，选择不同的调制解调器的每比特的样本数，使恢复信号与原输入信号相比较后得到差错率。根据差错率可以证明PN扩频信号是否有利于通信系统中正弦噪声的抑制。四：实验步骤简述（1）理解PN序列概念，掌握产生二进制PN序列的方法。（2）编写输入信号及正弦干扰信号和噪声信号子函数。（3）编写主函数，按照通信系统模型，通过对子函数的调用实现有PN序列的通信系统和无PN序列的通信系统：1): 有PN扩频序列时：原始信号S(n)与PN扩频序列相乘（异或实现）后，叠加正弦信号噪声，再加一般噪声完成信号的调制和在信道中的传输；无PN扩频序列时：原始信号S(n) 叠加正弦信号噪声，再加一般噪声完成信号的调制和在信道中的传输；2): 设置合适门限，滤除噪声，通过与PN扩频序列相乘（有PN扩频序列时）和不与PN扩频序列相乘（无PN扩频序列时）解调出信号；
