数字信号处理系统设计—课程设计论文

《数字信号处理系统设计—课程设计论文》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、本科生课程设计论文题目：数字信号处理学生姓名：学号：专业：通信工程班级：指导教师：2013年12月27日内蒙古科技大学课程设计答辩书姓名班级通信2班学号指导教师分数设计总结：通过本次的课程设计，不仅对课本所学知识有了进一步的理解，而且对MATLAB软件的使用也熟练了许多。在对题目做预期结果分析的同时，巩固了在课堂上所学的诸如“信号Fourier变换”、“滤波器设计”、“功率谱估计”等重要内容。用MATLAB软件对所做课题实现编程的同时，深切的体会到了计算机辅助软件对于信号处理的快捷与方便，MATLAB贴近数学公式化的编程思想另复杂的数学计算简单

2、了许多，尤其对于FFT运算的方法。尽管开始对于程序的设计不是很明了，但通过查阅课本及相关参考数目，对于MATLAB在信号处理方面的基本语法有所掌握，让我能够编写出一些简单的语句，并且进行调试和纠错。本次设计最大的收获摸过于亲手对课题进行分析并编程实现后，真正的将理论与实际联系起来，正所谓学以至用，数字信号处理的在通信工程系统中的重要意义在设计中得到了体现，以及这门的在实际处理领域的用处。对于在设计过程中对我给予帮助的同学和张超老师表示感谢。答辩记录：第一章1.1.(5)x(t)=sin(t)/t-10

3、);y=sinc(x);plot(x,y);ylabel('x(t)');xlabel('t');title('Sa函数曲线');1.2.(3)已知LTI离散系统，x(n)=[111],h(n)=[0123],求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('Timeindexn');subplot(2,2,2);stem([0:length(h)-1],h);ylabel('h(n)');

4、xlabel('Timeindexn')subplot(2,2,3);stem([0:length(y)-1],y);ylabel('y(n)=x(n)*h(n)');xlabel('Timeindexn');第二章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);p

5、lot(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;阻带边

6、缘频率：06，阻带衰减：50dBWp=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

7、,20*log10(abs(h)/norm))w=[Wp,Ws];h=freqz(numd,dend,w);fprintf('Ap=%.4',-20*log10(abs(h(1))));fprintf('As=%.4',-20*log10(abs(h(1))));Wp=04*pi;Ws=0.6*pi;Ap=0.5;As=50;Fs=0.5;wp=0.7265;ws=1.3764;N=buttord(wp,ws,Ap,As,'s');wc=wp/(10^(0.1*Ap)-1)^(1/2/N);[numa,dena]=butter(N,wc

8、,'s');[numd,dend]=bilinear(numa,dena,Fs);w=linspace(0,pi,512);h=freqz(numd,