大学毕业课程设计-数字信号处理课程设计.docVIP

数字信号处理 课程设计报告 信息科学与工程学院 通信0804 张驰 0909082328 目录 课程设计目的和任务 设计原理与方案 设计实现的过程 结果显示 附程序清单 个人心得 一、课程设计原理方案 1．所学知识 2．提高综合运用所学知识独立分析和解决问题的能力； （一）．信号分析 1、单位抽样序列 2）矩形序列 ???? 三角波序列 反三角波序列 5）Gaussian序列 ??? 6） 正弦序列 ?? ? 进行采样，可得到测试序列 。令A=50,采样周期T=1ms,即fs=1000Hz,f0=62.5,a=100。 2. 对上述信号完成下列信号分析 1）对三角波序列和反三角波序列，作N=8点的FFT，观察比较它们的幅频特性，说明它们有什么异同？绘出两序列及其它们的幅频特性曲线。 在和的尾部补零，作N=16点的FFT，观察它们的幅频特性发生了什么变化？分析说明原因。 ??2)、观察高斯序列，固定信号中的参数p=8，令q分别等于2，4，8，观察它们的时域和幅频特性，了解当q取不同值时，对信号序列的时域幅频特性的影响；固定q=8，令p分别等于8，13，14，观察参数p变化对信号序列的时域及幅频特性的影响，观察p等于多少时，会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。 3）对于正弦序列 4)、观察衰减正弦序列的时域和幅频特性，绘出幅频特性曲线，改变采样频率fs,使fs=300Hz， 观察此时的频谱的形状和谱峰出现位置？说明产生现象的原因。 3．设有一连续时间信号s(t),其由20Hz、220Hz和750Hz的正弦信号叠加而成，分析确定采样频率及数据分析长度，计算并绘出信号的频谱，指出各个频率份量。 （二）、数字滤波器设计 矩形序列：采样点N=31，取前5点为0，中间21点为1，后5点为0的矩形序列。MATLAB实现为： 三角波序列,反三角波序列及Gaussian序列 MATLAB实现为： n1=0:3;y1=n1+1;n2=4:7;y2=8-n2; n1=0:3;y1=4-n1;n2=4:7;y2=n2-3; n=0:15;p=5;q=8;x=exp(-1*(n-p).^2/q); 正弦序列：给出的抽样频率fs=64Hz;采样点N=16， ：其MATLAB实现如下： 衰减正弦序列：抽样频率A=50fs=1000hz，f0=62.5,a=100。对其取采样N=50。其MATLAB实现如下： A=50;a=100;f0=62.5;T=1/1000;n=1:50; x=A*exp(-1*a*n*T).*sin(2*pi*f0*n*T); 对三角波序列和反三角波序列，作N=8点的FFT对以上序列进行FFT变换，观察它的频谱图。 MATLAB实现在以上的基础上加上：即可。 求幅度特性： 求相位特性： 数字滤波器的设计与实现：在本次设计中采用了基于双线性变换IIR巴特沃思滤波器，为了方便随意输入或更改滤波器的技术指标，在设计中特意将指标设计成在滤波器界面可以随意输入。设计过程为首先将数字指标预畸变成模拟指标，根据模拟指标可以求得过渡模拟滤波器，再将模拟滤波器直接转换成数字滤波器。 其MATLAB实现过程主要如下： 低通滤波器：fp=input(输入通带截止频率fp:);fs=input(输入阻带截止频率fs:);Fs=input(输入抽样频率Fs:); rp=input(输入通带最大衰减系数rp:);rs=input(输入阻带最小衰减系数rs:); wp=2*Fs*tan(pi*fp/Fs);ws=2*Fs*tan(pi*fs/Fs); [n,wc]=buttord(wp,ws,rp,rs,s)；[z,p,G]=buttap(n); [bp,ap]=butter(n,wc,s)；[bs,as]=lp2lp(bp,ap,wp); [bz,az]=bilinear(bs,as,Fs);[h,w]=freqz(bz,az,Fs); 高通滤波器MATLAB实现过程同上类似，只是高通在： 上不同，低通与带通时在‘ftype’上为缺省，高通与带阻是在‘ftype’上分别为’high’,’stop’; 带通数字滤波器数字指标不一样： fpl=input(输入通带上截止频率fpl:);fpu=input(输入通带下截止频率fpu:);fsl=input(输入阻带上截止频率fsl:);fsu=input(输入阻带下截止频率fsu:); Fs=input(输入抽样频率Fs:); rp=input(输入通带最大衰减系数rp:);rs=input(输