资源描述：
《实验三 信号采样与重建(参考答案).ppt》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、实验三信号采样与重建LectureScript2013.5.5一、实验目的1、进一步学习MATLAB的函数及其表示。2、掌握及验证信号的SHANNON采样定理。3、由采样序列重构恢复原信号。二、实验内容1、对连续时间信号y(t)＝sin(24πt)+sin(40πt)，它有12Hz和20Hz两个等幅度分量。用MATLAB作图法求出Nyquist频率2fmax。%AnalogSignalDt=0.005;t=-1:Dt:1;%samplefrequency200Hzy=sin(24*pi*t)+s
2、in(40*pi*t);subplot(3,1,1);plot(t*4,y);xlabel('tin1/4sec.');ylabel('y(t)');title('AnalogSignal');subplot(3,1,2);plot(t*12,sin(24*pi*t));xlabel('tin1/12sec.');ylabel('sin(24*pi*t)');%f=12Hzsubplot(3,1,3);plot(t*20,sin(40*pi*t));xlabel('tin1/20sec.');y
3、label('sin(40*pi*t)');%f=20Hz复合时看不出频率分离成两个分量一般的分离成很多个正余弦，最高频率的2倍就是。另一解法频谱法（适合很多个正弦分量情况）t=0:0.01:1;%T＝0.01Sec，samplefrequency100Hzy=sin(24*pi*t)+sin(40*pi*t);Y=fft(y);f=(0:length(y)-1)/length(y)*100;％每点对应频率figure(1);plot(f,abs(Y));grid;%幅度figure(2);pl
4、ot(f,unwrap(angle(Y))*180/pi);grid;相位曲线不作折叠，而是全部展开的画法。两个频率分量采样频率2002：设连续信号x(t)=exp（－1000
5、t
6、）时A、求连续傅利叶变换X(jw)（用很密的点）所以：满足抽样定理，从而可以用xa[n]代替x(t),范围[－0.005到0.005之间每间隔5*10－5，20KHz，共200点]，从0到起点0.005时，e－5＝0.0067，已经足够小。0相位，相对幅度直流处最大0.0022：设连续信号x(t)=exp（－1000
7、
8、t
9、）时A、求连续时间傅利叶变换X(jw)t=-0.005:0.00005:0.005;%10mS长，201点密度x=exp(－1000*abs(t))；％共有201点Wmax=2*pi*2000;％只要算到频率2000Hz止M=500;m=0:1:M;W=m*Wmax/M;％得到W(1)~W(501)共501个频谱点值。Xa=x*exp(－j*t’*W)*0.00005;％计算谱501点Xa(1)~Xa(501).Xa=real(Xa);Xa=[fliplr(Xa),Xa(2:501)];%
10、左折但Xa(1)不动W=[－fliplr(W),W(2:501)];％频点也左折。figure(1);plot(t*1000,x,’r’);grid;%毫秒figure(2);plot(W/(2*pi*1000),Xa*1000,’r’);％幅度放大1000频谱波形幅度放大1000倍后2：设连续信号x(t)=exp（－1000
11、t
12、）时B、Fs=5000,和1000Hz先采样，在求离散DFT频谱。t=-0.005:0.00005:0.005;%原信号条件不变x1=exp(－1000*abs(t)
13、)；Ts=0.0002;n=－25:1:25;%Fs=5000,x=exp(－1000*abs(n*Ts))；M=500;m=0:1:M;W=pi*m/M;％频率点Xa=x*exp(－j*n’*W）;％计算离散谱Xa=real(Xa);Xa=[fliplr(Xa),Xa(2:M+1)];W=[－fliplr(W),W(2:M+1)];figure(1);plot(t*1000,x1,’r’);holdon;stem(n*Ts*1000,x);holdoff;figure(2);plot(W/pi
14、,Xa);频谱。序列C、减小Fs=1000Hz，先采样，再求离散DFT频谱。t=-0.005:0.00005:0.005;%连续的条件不变x1=exp(－1000*abs(t))；Ts=0.001;n=－5:1:5;%Fs=1000,Ts＝1mSx=exp(－1000*abs(n*Ts))；％n*Ts=－0.005～0.005M=500;m=0:1:M;W=pi*m/M;％频率点Xa=x*exp(－j*n’*W）;％计算离散谱Xa=real(Xa);Xa=[fliplr(Xa),Xa(2:M+1
