大学电子信息课程设计.docVIP

二○一○～二○一一学年第 一 学期 信息科学与工程学院 课程设计报告书 课程名称： 数字信号处理课程设计 班 级： 电子信息工程2008级 6 班 学 号： 200804135167 姓 名： 石深刚 指导教师： 刘老师 2011 年 1 月 05日 信号与系统课程设计 一、课程设计目的： 通过本课程设计使学生进一步巩固数字信号处理的基本概念、理论、分析方法和实现方法；使学生掌握的基本理论和分析方法方面的知识得到进一步扩展；使学生能有效地将理论和实际紧密结合；增强学生软件编程实现能力和解决实际问题的能力。 要求学生能够熟练地用Matlab等编写数字信号处理的应用程序；熟悉用FFT对连续信号和离散信号进行谱分析；编程实现IIR数字滤波器和FIR数字滤波器；了解各种窗函数对滤波器特性的影响等，进一步明确数字信号处理的工程应用。 (1) 已知某测试系统观测到的信号，其中 ， ， n(t) 为高斯白噪声，试设计滤波系统，输入为 ，输出为和。 要求：确定采样频率、采样长度，对进行采样； 设计数字滤波器； 通过计算机仿真对滤波器的性能进行分析因此，截取信号的持续时间最好为0.1s的整数倍，取为0.2s。 信号x(t)的最高频率为Ｆh=200Hz，因此采样频率或采样间隔需满足 Fs≥2Ｆh=400Hz,T≤1/2Fh=2.5ms 取采样频率和采样间隔分别为：fs=500Hz, T=2ms N=0.2s/2ms=100 令t=nT，ｎ=０,1,2,…,99,代入原信号，得 x(n) = 20sin(70π×n×2.5×0.001)+ 40sin(200π×n×2.5×0.001) =20sin(0.175π×n)+40sin(0.5π×n) MATLAB源程序： fs=500; T=0.2; N=100; n=0:1:N-1; t=n/fs; s1=20.*sin(70*pi*t); s2=40.*sin(200*pi*t); n=0.1*randn(1,N); x=s1+s2+n; figure(1) plot(x); w=0:0.001:pi; Wp=90*pi; Ws=180*pi; Rp=1; Rs=15; [N,Wn]=buttord(Wp,Ws,Rp,Rs,s); [B,A]=butter(N,Wn,s); [bz,az]=bilinear(B,A,fs); y=filter(bz,az,x); figure(2) freqz(bz,az); figure(3) plot(y); Wp1=[40*pi 100*pi]; Ws1=[60*pi 80*pi]; Rp1=3; Rs1=30; w=0:1:50; [N1,Wn1]=buttord(Wp1,Ws1,Rp1,Rs1,s); [B1,A1]=butter(N1,Wn1,stop,s); [bz1,az1]=bilinear(B1,A1,fs); y1=filter(bz1,az1,x); figure(4) freqz(bz1,az1); figure(5) plot(y1); (2) 针对含有5Hz、15Hz和30Hz的混合正弦波信号，试设计滤波系统，5Hz和30Hz的正弦分量，阻带的最小衰减不小于50dB。要求：确定采样频率、采样长度，对进行采样；窗设计数字滤波器； 通过计算机仿真对滤波器的性能进行分析取采样频率和采样间隔分别为：fs=120Hz,T=10ms N=0.4 /T=40 MATLAB源程序： fs=120; T=0.4 N=40; n=0:1:N-1; t=n/fs; x=sin(5*2*pi*t)+sin(15*2*pi*t)+sin(30*2*pi*t); figure(1) plot(x); w=0:0.01:pi; Wc=[10*2/fs;22.5*2/fs]; h=fir1(N-1,Wc,hamming(N)); figure(2); freqz(h,1,w); figure(3); y=filter(h,1,x); plot(t,y); 实验小结： 经过一周的课程设计，我学到了很多东西。对于以前不理解的知识，通过试验的学习得到了理解，学会的知识也得到了进一步深化。 这学期开设的数字信号处理课程是信号与系统课程的延续，带着对信号与系统学习的兴趣，我满怀信心的开始了对数字信号处理这门课程的学习。 因为对信号与系统这门课程学习的还算透彻，，所以以为数字信号处理这门课程也应该不在话下，但事实上并非如此。信号与系统