數字信号处理系统设计书—课程设计书论文.docVIP

本科生课程设计论文 题 目：数字信号处理 学生姓名： 学 号： 专 业：通信工程 班 级： 指导教师： 2013年 12 月 27日 内蒙古科技大学课程设计答辩书 姓名 班级 通信2班 学号 指导教师 分数 设计总结：通过本次的课程设计，不仅对课本所学知识有了进一步的理解，而且对MATLAB软件的使用也熟练了许多。在对题目做预期结果分析的同时，巩固了在课堂上所学的诸如“信号Fourier变换”、“滤波器设计”、“功率谱估计”等重要内容。用MATLAB软件对所做课题实现编程的同时，深切的体会到了计算机辅助软件对于信号处理的快捷与方便，MATLAB贴近数学公式化的编程思想另复杂的数学计算简单了许多，尤其对于FFT运算的方法。 尽管开始对于程序的设计不是很明了，但通过查阅课本及相关参考数目，对于MATLAB在信号处理方面的基本语法有所掌握，让我能够编写出一些简单的语句，并且进行调试和纠错。本次设计最大的收获摸过于亲手对课题进行分析并编程实现后，真正的将理论与实际联系起来，正所谓学以至用，数字信号处理的在通信工程系统中的重要意义在设计中得到了体现，以及这门的在实际处理领域的用处。 对于在设计过程中对我给予帮助的同学和张超老师表示感谢。 答辩记录： 第一章 1.1.(5) x (t)=sin(t)/t -10t10 x=linspace(-10,10); y=sinc(x); plot(x,y); ylabel(x(t)); xlabel(t); title(Sa函数曲线); 1.2.(3) 已知LTI离散系统，x(n)=[1 1 1],h(n)=[0 1 2 3],求y(n) x=[1,1,1,]; h=[0,1,2,3,]; y=conv(x,h); subplot(2,2,1);stem([0:length(x)-1],x); ylabel(x(n));xlabel(Time index n); subplot(2,2,2);stem([0:length(h)-1],h); ylabel(h(n));xlabel(Time index n) subplot(2,2,3);stem([0:length(y)-1],y); ylabel(y(n)=x(n)*h(n));xlabel(Time index n); 第二章 2.1.2．用DFT计算下列信号的频谱： （1） N=30; %数据的长度 L=1024; %DFT的点数 f=1/16;fs=600; T=1/fs; ws=2*pi*fs; t=(0:N-1)*T; x=5*cos(2*pi*f*t+pi/4); X=fftshift(fft(x,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(X)); ylabel(幅度谱) 2.1.(3) N=30; L=1024; f1=0.5;f2=4;fs=600; T=1/fs; ws=2*pi*fs; t=(0:N-1)*T; x=2*sin(2*pi*f1*t)+sin(2*pi*f2*t); X=fftshift(fft(x,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(X)); ylabel(幅度谱) 第三章 5．采用脉冲响应不变法和双线性变换法设计巴特沃斯数字低通滤波器，满足下列指标： 通带边缘频率：0.4，通带衰减：0.5dB ; 阻带边缘频率：06，阻带衰减：50dB Wp=04*pi;Ws=0.6*pi;Ap=0.5;As=50; Fs=1; wp=Wp*Fs;ws=Ws*Fs; N=buttord(wp,ws,Ap,As,s); wc=wp/(10^(0.1*Ap)-1)^(1/2/N); [numa,dena]=butter(N,wc,s); [numd,dend]=impinvar(numa,dena,Fs); w=linspace(0,pi,512); h=freqz(numd,dend,w); norm=max(abs(h)); numd=numd/norm; plot(w/pi,20*log10(abs(h)/norm)) w=[Wp,Ws]; h=freqz(numd,dend,w); fprintf(Ap=%.4\n,-20*log10(abs(h(1)))); fprintf(As=%.4\n