信号的自然采样与恢复.doc

《信号的自然采样与恢复.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验九信号的自然采样与恢复一、实验目的：1、理解信号的采样及采样定理以及自然采样信号的频谱特征。2、掌握和理解信号自然采样以及信号重建的原理，并能用MATLAB实现。二、实验原理及方法：本实验主要涉及采样定理的相关内容以及低通滤波器恢复原连续信号的相关知识。信号的抽样与恢复示意图如图7-1所示。图7-1信号的抽样与恢复示意图信号抽样与恢复的原理框图如图7-2所示。图7－2信号抽样与恢复的原理框图由原理框图不难看出，A/D转换环节实现抽样、量化、编码过程；数字信号处理环节对得到的数字信号进行必要的处理；D/A转换环节实现数/模转换，得到连续时间信号；低通滤波器的作用是滤除截止频率以外的

2、信号，恢复出与原信号相比无失真的信号。原信号得以恢复的条件是，其中为采样频率，B为原信号占有的频带宽度。为最低采样频率，当时，采样信号的频率会发生混迭，所以无法用低通滤波器获得原信号频谱的全部内容。三、实验内容及步骤：给定带限信号f(t)，其频谱为1、画出此信号的频谱图（ω的取值：-0.5π<ω<0.5π，精度取0.01rad）。答：画出f(t)的频谱图即F(W)的图像程序代码如下：#include#include#definePI3.14doublef(doublew){if(w>=-0.5*PI&&w<=0.5*PI)returncos(w);e

3、lsereturn0;}main(){doublew,F;FILE*fp;for(w=-0.5*PI;w<=0.5*PI;w+=0.01){F=f(w);printf("w=%.2f,F(w)=%f",w,F);fp=fopen("d:\2.txt","w");fprintf(fp,"%ft",F);}system("pause");}③F(W)的图像2、对此频域信号进行傅里叶逆变换，得到相应的时域信号，画出此信号的时域波形f(t)（t的取值：-20s

4、：#include#include#definePI3.14doublef(doublet){doublew=-0.5*PI,f=0;doubleans;while(w<=0.5*PI){f=f+cos(w)*cos(w*t)*0.01;w=w+0.01;}ans=f/(2*PI);returnans;}main(){doublet,ans=0;for(t=-20;t<20;t+=0.1){ans=f(t);printf("t=%.1f,f(t)=%f",t,ans);}system("pause");}将实验数据导入matlab，对数据进行仿真

5、绘图。3分别用三种不同抽样频率f=0.2Hz,0.5Hz,1.0Hz的周期矩形脉冲信号（矩形脉冲的幅度E取1，宽度τ取0.01s）对f(t)进行抽样，画出抽样后的信号的频谱图（ω的取值：-10rad<ω<10rad，精度取0.01rad）。答：三种不同频率的抽样分别用三种不同抽样频率f=0.2Hz,0.5Hz,1.0Hz的周期矩形脉冲信号（矩形脉冲的幅度E取1，宽度τ取0.01s）对f(t)进行抽样，画出抽样后的信号的频谱图程序代码如下：#include#include#definePI3.14#defineINF2000doublesa(doubl

6、ex){if(x==0)return1;elsereturnsin(x)/x;}doublef1(doublew){if(w>=-0.5*PI&&w<=0.5*PI)returncos(w);elsereturn0;}doublef2(doublew,doublef){doubles=0,ans;intn;for(n=-INF;n<=INF;n+=1){s=s+f1(w-n*2*PI*f)*sa(n*2*PI*f*0.01/2);}ans=s*0.01*f;returnans;}main(){doublew,fs=0,f;f=0.25;/*0.2,0.5,1*/for(w=-10;w

7、<=10;w+=0.1){fs=f2(w,f);printf("w=%.2f,Fs(w)=%f",w,fs);}system("pause");}③不同抽样频率下的频谱图：f=0.2Hzf=0.5Hzf=1.0Hz4针对3中抽样所得的矩形抽样信号，用滤波器对所得信号进行滤波，所得恢复信号f(t)的频谱记为F‘(w)，与原信号的频谱F(w)进行比较（ω的取值：-2rad<ω<2rad，精度取0.01rad）。将恢复信号的频谱图与原信号的频谱图进行比较①程