实验二时域采样与频域采样Word文档下载推荐.doc

实验二时域采样与频域采样Word文档下载推荐.doc

《实验二时域采样与频域采样Word文档下载推荐.doc》由会员分享，可在线阅读，更多相关《实验二时域采样与频域采样Word文档下载推荐.doc（6页珍藏版）》请在冰豆网上搜索。

采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。

3.实验内容及步骤

%物联一班胡洪201313060110

%2015年10月24日

%实验二：

程序1

Tp=64/1000;

Fs=1000;

T=1/Fs;

M=ceil（Tp*Fs）;

n=0:

M-1;

A=444.128;

a=pi*50*2^0.5;

w=pi*50*2^0.5;

xnt=A*exp（-a*n*T）.*sin（w*n*T）;

Xk=fft（xnt,M）;

subplot（3,2,1）;

stem（n,xnt,'

.'

）;

axis（[1,65,-5,150]）;

title（'

图1Fs=1000Hz'

subplot（3,2,2）;

plot（n/Tp,abs（Xk））;

图2Fs=1000Hz幅度'

Fs=300;

subplot（3,2,3）;

axis（[0,M,-10,150]）

图3Fs=300Hz'

subplot（3,2,4）;

图4Fs=300Hz幅度'

Fs=200;

subplot（3,2,5）;

图5Fs=200Hz'

subplot（3,2,6）;

图6Fs=200Hz幅度'

图1

程序2

13;

xa=n+1;

n=14:

26;

xb=27-n;

xn=[xa,xb];

stem（n,xn,'

三角波序列x（n）'

axis（[0,32,0,15]）

Xk=fft（xn,1024）;

k=0:

1023;

wk=2*k/1024;

plot（wk,abs（Xk））;

ylabel（'

|X（e^（j*w））|'

xlabel（'

w/ pi'

）

axis（[0,1,0,200]）;

Xk32=fft（xn,32）;

31;

stem（k,abs（Xk32）,'

axis（[0,16,0,200]）;

k'

|Xk3 2|'

xn32=ifft（Xk32）;

stem（k,xn32,'

axis（[0,32,0,15]）;

Xk16=Xk32（1:

2:

32）;

15;

stem（k,abs（Xk16）,'

axis（[0,8,0,200]）;

|Xk16 '

xn16=ifft（Xk16）;

stem（k,xn16,'

图2

4.思考题：

1.如果序列x（n）的长度为M，希望得到其平铺X（e^（jw））在[0,2pi]上的N点等间隔采样，当N<

M,如何用一次最小的点数DFT得到该频谱采样？

答：

先对原序列x（n）以N为周期进行周期延拓后取主值区序列，再计算N点DFT则得到N点频域采样

学生实验心得

1.通过本次实验，了解到了MATLAB软件的开发环境、函数库、语言、图形处理系统、应用程序接口，最主要的是用MATLAB对信号做的相关处理。

同时熟悉了MATLAB语言的特点，知道MATLAB是一种面向科学与工程计算的高级语言，允许用数学形式的语言编写程序。

它自身有独立的编写文件格式（.M），可以像脚本一样在对应平台上运行，调试同时强大的图形处理功能是这个软件的一大特点。

2.同时，用MATLAB许多函数用于产生常用的基本信号，可对时域采样与疲于采样的理论进行验证。

通过对MATLAB的仿真掌握了模拟信号采样的频谱的变化，以及如何选择采样频率才能是采样后的信息不丢失信息

3.明白了频域采样点数的选择，对《数字信号处理》以后的学习打下了基础。

指导

教师

评语

成绩评定：

指导教师（签名）：

年月日