数字信号处理系统设计—课程设计论文.doc
本科生课程设计论文题 目：数字信号处理学生姓名：学 号：专 业：通信工程班 级： 指导教师： 2013年 12 月 27日 内蒙古科技大学课程设计答辩书姓名班级通信2班学号指导教师分数设计总结：通过本次的课程设计，不仅对课本所学知识有了进一步的理解，而且对MATLAB软件的使用也熟练了许多。在对题目做预期结果分析的同时，巩固了在课堂上所学的诸如“信号Fourier变换”、“滤波器设计”、“功率谱估计”等重要内容。用MATLAB软件对所做课题实现编程的同时，深切的体会到了计算机辅助软件对于信号处理的快捷与方便，MATLAB贴近数学公式化的编程思想另复杂的数学计算简单了许多，尤其对于FFT运算的方法。 尽管开始对于程序的设计不是很明了，但通过查阅课本及相关参考数目，对于MATLAB在信号处理方面的基本语法有所掌握，让我能够编写出一些简单的语句，并且进行调试和纠错。本次设计最大的收获摸过于亲手对课题进行分析并编程实现后，真正的将理论与实际联系起来，正所谓学以至用，数字信号处理的在通信工程系统中的重要意义在设计中得到了体现，以及这门的在实际处理领域的用处。 对于在设计过程中对我给予帮助的同学和张超老师表示感谢。答辩记录：第一章1.1.(5) x (t)=sin(t)/t -10<t<10x=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.7690 0.8311 -0.3740 -0.3442 Columns 37 through 45 1.0592 -0.5897 -0.3410 0.7344 -1.4344 0.2595 0.7010 -0.5160 0.2881 Columns 46 through 51 -0.0099 -0.5054 0.5843 -0.3761 -0.6786 0.7528
