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

本科生课程设计论文题 目：数字信号处理学生姓名：学 号：专 业：通信工程班 级： 指导教师： 2013年 12 月 27日 内蒙古科技大学课程设计答辩书姓名班级通信2班学号指导教师分数设计总结：通过本次的课程设计，不仅对课本所学知识有了进一步的理解，而且对matlab软件的使用也熟练了许多。在对题目做预期结果分析的同时，巩固了在课堂上所学的诸如“信号fourier变换”、“滤波器设计”、“功率谱估计”等重要内容。用matlab软件对所做课题实现编程的同时，深切的体会到了计算机辅助软件对于信号处理的快捷与方便，matlab贴近数学公式化的编程思想另复杂的数学计算简单了许多，尤其对于fft运算的方法。 尽管开始对于程序的设计不是很明了，但通过查阅课本及相关参考数目，对于matlab在信号处理方面的基本语法有所掌握，让我能够编写出一些简单的语句，并且进行调试和纠错。本次设计最大的收获摸过于亲手对课题进行分析并编程实现后，真正的将理论与实际联系起来，正所谓学以至用，数字信号处理的在通信工程系统中的重要意义在设计中得到了体现，以及这门的在实际处理领域的用处。 对于在设计过程中对我给予帮助的同学和张超老师表示感谢。答辩记录：第一章1.1.(5) x (t)=sin(t)/t -10t10x=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=%.4n,-20*log10(abs(h(1);fprintf(as=%.4n,-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,s);numd,dend=bilinear(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=%.4n,-20*log10(abs(h(1);fprintf(as=%.4n,-20*log10(abs(h(1);第四章3已知一含有平稳高斯白噪声的序列xk= sin (0.8k) + s k，试分别用ld算法和burg算法实现该序列的功率谱估计，并估计其ar模型参数。n=512;nfft=1024;fs=2*pi;n=0:n-1; xn=sin(0.8*pi*n)+randn(size(n); order=50; figure(1) pburg(xn,order,nfft,fs);xlabel (50阶); r=xcorr(xn);order1=80; figure(2)pburg(xn,order1,nfft,fs);xlabel (80阶); l=levinson(r,order)l = columns 1 through 9 1.0000 0.0742 -0.2620 -0.5832 -0.2657 0.1715 0.3907 -0.3734 0.0043 columns 10 through 18 -0.7425 -0.2535 0.8478 -1.0250 -0.1673 1.7374 -1.0348 0.1818 0.1069 columns 19 through 27 -2.3389 0.8189 1.1640 -3.3557 1.0312 1.7547 -2.3801 0.7517 0.7413 columns 28 through 36 -3.3721 1.3212 1.3154 -2.7126 0.3671 1.4414 -0.9785 0.7407 -0.3178 columns 37 through 45 -1.0306 0.9474 0.1621 -0.9682 0.3659 0.3619 -0.1908 -0.2529 -0.5377 columns 46 through 51 -0.0414 0.3843 0.5893 0.1865 -0.3803 -0.5296n=512;nfft=1024;fs=2*pi;n=0:n-1; xn=sin(0.8*pi*n)+randn(size(n); order=50; figure(1) pyulear(xn,order,nfft,fs);xlabel (50阶); r=xcorr(xn);order1=80; figure(2) pyulear(xn,order1,nfft,fs);xlabel (80阶); l=levinson(r,order)l = columns 1 through 9 1.0000 -0.3463 -0.2629 0.4612 -0.4435 -0.0316 0.2078 -0.7380 0.9364 columns 10 through 18 0.7226 -1.1189 0.7047 -0.3691 -0.7053 0.9265 -0.2019 -0.2709 0.8969 columns 19 through 27 -0.9410 0.6964 1.0837 -0.7729 0.4398 -0.8230 -0.0867 0.0718 0.2445 columns 28 through 36 -0.3921 0.5204 -1.2555 0.5562 0.5678 -0.7