正弦信号的采样与恢复

参考文献

[1] 高西全、丁玉美编著.数字信号处理.西安:西安电子科技大学出版社，2008. [2] 郑君里等编.信号与系统.北京:高等教育出版社，2000.

[3] 刘树棠译.数字信号处理——使用MATLAB.西安:西安交通大学出版社，2002. [4] 刘敏.MATLAB通信仿真与应用.北京：国防工业出版社.

[5] 陈怀琛等编著.MATLAB及在电子信息课中的应用.北京：电子工业出版社，2002. [6]樊昌信.通信原理.北京：国防工业出版社，2002.

[7] 刘卫国主编.MATLAB程序设计与应用（第二版）. 北京：高等教育出版社，2006.

14

致谢

通过两周的课程设计，我们都学到了许多东西，体会到了从书本学习与实际应用中的不同，这种感同身受必将对我们今后的学习与生活带来很大的帮助。在两周的时间中，我们的指导老师时时在我们的身旁引导我们，帮助我们，倾注他们所有的才华，用心血让我们学会从理论走向实际这一目标十分艰辛，对我们来说是飞跃的过程。并且不会忘记他们为了我们紧缩的眉头，焦急的眼神，疲倦的笑容，忙碌的身影！谢谢我们的指导老师对我们的帮助！

15

附录

1、连续信号y=sin(t)和其对应的频谱

x1=0:pi/10:(10*pi);

w=linspace(0,2.5*pi,length(x1)); figure subplot(211)

plot(x1,sin(x1)); %原时域连续信号y=sin(t) xlabel('t');ylabel('x(t)'); title('原时域连续信号y=sin(t)'); grid

sin1=sin(x1); n=0:(length(x1)-1); subplot(212)

plot(w,fft1(w,sin1,n)); %其对应频域信号Y=FFT（sin(t)） xlabel('w');ylabel('x(w)');

title('其对应频域信号Y=FT（sin(t)）'); grid result=fft1(w,hanshu,n) a=cell(1,length(w)); for i=1:length(w);

m=hanshu.*((exp(-1i*(1i-1)*pi/100)).^n); a{i}=sum(m); end for i=1:length(w); end

2、对连续信号y=sin(t)进行抽样并产生其频谱 n1=input('请输入采样点数n:'); n=0:n1; zb=size(n); figure

sinf=sin(10*pi*n/zb(2)); subplot(211); stem(n,sinf,'.');

16

xlabel('n');ylabel('x(n)'); title('采样后的时域信号y=x(n)'); w=0:(pi/100):5*pi; subplot(212)

plot(w,fft1(w,sinf,n)); xlabel('w');ylabel('x(w)');

title('采样后的频域信号y=FT（sin(n)）'); grid 3、低通滤波器

[B,A]=butter(8,350/500); [H,w]=freqz(B,A,128,1500); figure;

plot(w*2000/(2*pi),abs(H)); xlabel('Hz');ylabel('频率响应幅度'); title('低通滤波器'); grid; 4、恢复

n1=input('请输入采样点数n:'); n=0:n1; zb=size(n); figure

sinf=sin(10*pi*n/zb(2)); subplot(211); stem(n,sinf,'.');

xlabel('n');ylabel('x(n)'); title('采样后的时域信号y=x(n)'); w=0:(pi/100):5*pi; subplot(212)

plot(w,fft1(w,sinf,n)); xlabel('w');ylabel('x(w)');

title('采样后的频域信号y=FT（sin(n)）'); grid

17