资源描述:
《信号抽样实验报告材料》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实用标准文案大连理工大学实验报告学院(系):专业:班级:姓名:学号:组:___实验时间:实验室:实验台:指导教师签字:成绩:实验三信号抽样一、实验目的1学会运用MATLAB完成信号抽样及对抽样信号的频谱进行分析;2学会运用MATLAB改变抽样间隔,观察抽样后信号的频谱变化;3学会运用MATLAB对抽样后的信号进行重建。二、习题1.设有三个不同频率的正弦信号,频率分别为。现在用抽样频率对这三个正弦信号进行抽样,用MATLAB命令画出各抽样信号的波形及频谱,并分析频率混叠现象。解:分别写出三个频率正弦波的代码与图形:(f1=100HZ的正弦信号)代码如下:Ts=1/
2、3800;dt=0.0001;t1=-0.008:dt:0.008;ft=sin(2*pi*100*t1).*(uCT(t1+0.005)-uCT(t1-0.005));subplot(221);plot(t1,ft),gridon;axis([-0.0060.006-1.11.1]);xlabel('Time(sec)'),ylabel('f(t)')title('正弦信号波形');精彩文档实用标准文案N=5000;k=-N:N;W=2*pi*k/((2*N+1)*dt);Fw=dt*ft*exp(-j*t1'*W);subplot(222);plot(W,a
3、bs(Fw));gridon;axis([-300003000000.006]);xlabel('omega'),ylabel('F(w)');title('正弦信号的频谱');t2=-0.008:Ts:0.008;fst=sin(2*pi*100*t2).*(uCT(t2+0.005)-uCT(t2-0.005));subplot(223);plot(t1,ft,':'),holdon;stem(t2,fst),gridon;axis([-0.0050.005-1.11.1]);xlabel('Time(sec)'),ylabel('fs(t)');titl
4、e('抽样后的信号'),holdoff;Fsw=Ts*fst*exp(-j*t2'*W);subplot(224);plot(W,abs(Fsw)),gridon;axis([-300003000000.006]);xlabel('omega'),ylabel('Fs(w)');title('抽样信号的频谱');matlab波形如下:其中单个正弦信号(未经抽样)的频谱放大后如下:精彩文档实用标准文案(200HZ的正弦信号)代码如下:Ts=1/3800;dt=0.0001;t1=-0.003:dt:0.003;ft=sin(2*pi*200*t1).*(uCT(
5、t1+0.0025)-uCT(t1-0.0025));subplot(221);plot(t1,ft),gridon;axis([-0.0030.003-1.11.1]);xlabel('Time(sec)'),ylabel('f(t)')title('200HZ正弦信号波形');N=5000;k=-N:N;W=2*pi*k/((2*N+1)*dt);Fw=dt*ft*exp(-j*t1'*W);subplot(222);plot(W,abs(Fw));gridon;axis([-300003000000.003]);xlabel('omega'),ylabe
6、l('F(w)');title('200HZ正弦信号的频谱');t2=-0.003:Ts:0.003;fst=sin(2*pi*200*t2).*(uCT(t2+0.0025)-uCT(t2-0.0025));subplot(223);plot(t1,ft,':'),holdon;stem(t2,fst),gridon;axis([-0.0030.003-1.11.1]);xlabel('Time(sec)'),ylabel('fs(t)');title('200HZ抽样后的信号'),holdoff;精彩文档实用标准文案Fsw=Ts*fst*exp(-j*t2'
7、*W);subplot(224);plot(W,abs(Fsw)),gridon;axis([-300003000000.003]);xlabel('omega'),ylabel('Fs(w)');title('200HZ抽样信号的频谱');matlab波形如下:(3800HZ正弦信号)代码如下:Ts=1/3800;dt=0.00001;t1=-1/7600:dt:1/7600;ft=sin(2*pi*3800*t1).*(uCT(t1+1/7600)-uCT(t1-1/7600));subplot(221);plot(t1,ft),gridon;axis([
8、-1/76001/760