5 结论与问题讨论
5.1 课程设计总结
本次课程设计历时一个星期，最终各部分要求都得到了较好的实现，并且在此基础上设计了GUI界面，使得显示界面更加的清晰明了，让人一目了然。
5.2 心得体会
本次课程设计主要是运用到MATLAB软件编程，由于之前有过MATLAB的实验进行练习，对于MATLAB软件编程有一定的基础，但是并不是很熟练，在本次设计过程当中还存在一个摸索和熟练的过程。通过本次课程设计，经过我和同伴的共同研究和探索，我们对MATLAB软件的使用有了更多的了解，使用起来也更加的熟练自如；此外，我们通过查找资料，自学掌握了基于MATLAB的IIR和FIR滤波器的设计方法。在设计过程当中，我们也遇到不少的问题和困惑，但是通过组内讨论还有请教同学，最终经过多次的尝试，这些问题都得到了较好的解决。在为期三周的通信系统综合训练当中我们系统的学习了嵌入式开发的知识，通过本次综合训练让我了解了各种移动信道的知识，这次综合训练对将来从事这一工作的同学尤为重要，这等于让我们提前对这方面的知识有了进一步的认识。在这次综合训练过程中也当遇到了很多实际问题，比如说对有些关于嵌入式知 识的缺乏，不过在老师的进一步讲解之后，我们逐步认识和理解了其中的关键，而且又对所学的理论知识有了升华。
同时通过本次嵌入开发综训练，加深了我对所学过的各种理论数据的认识和理解，并在一定程度上掌握并会运用。我还学会了把学到的知识用于解决实际问题，培养、加强锻炼了我的动手实践能力。更为难得的是，在这次训练过程中，屡屡碰见一些问题，在解决这些问题的过程中，不断加强了我对的理解。嵌入式对于一些自己不清楚，不明白但平时又很难发现的知识点有了一次全面的巩固与复习。
在大学阶段，理论的学习和实践是密不可分的。离开了实践的理论是没有任何意义可言的。与此同时，理论是需要伴随着实践才能完善。同时，在与同学共同在解决一些问题的过程中，提高了我们的团队协作精神。
参考文献
1．《数字信号处理教程——MATLAB释义与实现》 陈怀琛编著 电子工业出版社
2.《辅助信号处理技术与应用 MATLAB 7》飞思科技产品研发中心编著 电子工业出版社
3.《自适应滤波器原理（第四版）》 西蒙·赫金编著 郑宝玉等译 电子工业出版社
附 录
close all clf
t=-1:0.001:1; n=1:256; N=512; fs=1000;
x1=cos(150*pi*t); f1=n*fs/N; figure(1); subplot(3,1,1) plot(x1);
title('x1时域波形');
xlabel('x1的时间');ylabel('x1的幅值');
x2=cos(320*pi*t); subplot(3,1,2) plot(x2);
title('x2的时域波形');
xlabel('x2的时间）');ylabel('x2的幅值'); axis([0,100,-1,1])
x3=cos(450*pi*t);
subplot(3,1,3) plot(x3);
title('x3的时域波形');
xlabel('x3的时间');ylabel('x3的幅值'); axis([0,100,-1,1])
y1=fft(x1,512); figure(2); subplot(3,1,1)
plot(f1,abs(y1(1:256))); title('x1的频域波形');
xlabel('x1的频率（Hz）');ylabel('x1的幅值'); axis([0,500,0,250])
y2=fft(x2,512); subplot(3,1,2)
plot(f1,abs(y2(1:256))); title('x2的频域波形');
xlabel('x2的频率（Hz）');ylabel('x2的幅值'); axis([0,500,0,250])
y3=fft(x3,512); subplot(3,1,3)
plot(f1,abs(y3(1:256))); title('x3的频域波形');
xlabel('x3的频率（Hz）');ylabel('x3的幅值'); axis([0,500,0,250]) %close all;
x=x1+x2+x3; figure(3) subplot(211); plot(x(1:100));
title('三个信号叠加的时域波形'); y=fft(x,512); subplot(212);
plot(f1,abs(y(1:256)));
title('三个信号叠加的频谱') axis([0,600,0,250])
%close all
Wp=2*pi*75;Ws=2*pi*100;Rp=1;As=30; [N,wc]=buttord(Wp,Ws,Rp,As,'s');
[b,a]=butter(N,wc,'s');
k=1:512;fk=0:1000/1024:1000; wk=2*pi*fk;
Hk=freqs(b,a,wk); figure(4);
plot(fk,abs(Hk));grid on;
xlabel('频率(Hz)');ylabel('幅度'); axis([0,200,0,1.1])
[Bz,Az]=bilinear(b,a,1000); y=filter(Bz,Az,x); Y1=fft(y,512); figure(5); subplot(211); plot(y);
title('滤波后x1的时域波形'); axis([0,200,-1,1]) subplot(212);
plot(f1,abs(Y1(1:256)));
title('滤波后x1的频谱波形');
%close all
wp=2*pi*[140,200]; ws=2*pi*[120,230]; Rp=1; As=30;
[N,wc]=buttord(wp,ws,Rp,As,'s'); [B1,A1]=butter(N,wc,'s'); k=0:511;
fk=0:1000/512:1000; wk=2*pi*fk;
Hk=freqs(B1,A1,wk);
figure(6)
plot(fk,20*log10(abs(Hk))); grid on
title('带通滤波器的频响2') xlabel('频率(Hz)'); ylabel('幅度(dB)') axis([100,400,-40,5])
[Bz1,Az1]=bilinear(B1,A1,1000); y2=filter(Bz1,Az1,x);
Y2=fft(y2,512); figure(7); subplot(211); plot(y2);
title('滤波后x2的时域波形'); axis([0,200,-1,1])
subplot(212);
plot(f1,abs(Y2(1:256)));
title('滤波后x2的频谱波形');
%%%%%%%%%%%%%%%%%%%%%% wp=2*pi*[260,320]; ws=2*pi*[240,350]; Rp=1; As=30;
[N,wc]=buttord(wp,ws,Rp,As,'s'); [B2,A2]=butter(N,wc,'s'); k=0:511;
fk=-1000:1000/512:1000; wk=2*pi*fk;
Hk=freqs(B2,A2,wk);
figure(8)
plot(fk,20*log10(abs(Hk))); grid on
title('带通滤波器的频响3') xlabel('频率(Hz)'); ylabel('幅度(dB)') axis([200,600,-40,5])
[Bz2,Az2]=bilinear(B2,A2,1000); y3=filter(Bz2,Az2,x); Y3=fft(y3,512); figure(9); subplot(211); plot(y3);
title('滤波后x3的时域波形'); axis([0,200,-1,1])
subplot(212);
plot(f1,abs(Y3(1:256)));
title('滤波后x3的频谱波形');