实验二：时域采样与频域采样

实验二：时域采样与频域采样

1. 实验目的

时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。

2. 实验原理与方法

对模拟信号以间隔T进行时域等间隔理想采样，形成的采样信号的频谱是原模拟信号频谱以采样角频率()为周期进行周期延拓。公式为：

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

实验内容及步骤

%物联一班 胡洪 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('图1 Fs=1000Hz');

subplot(3,2,2);plot(n/Tp,abs(Xk));title('图2 Fs=1000Hz幅度');

Fs=300;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,3);

stem(n,xnt,'.');axis([0,M,-10,150])

title('图3 Fs=300Hz');

subplot(3,2,4);plot(n/Tp,abs(Xk));title('图4 Fs=300Hz幅度');

Fs=200;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,5);

stem(n,xnt,'.');axis([0,M,-10,150])

title('图5 Fs=200Hz');

subplot(3,2,6);plot(n/Tp,abs(Xk));title('图6 Fs=200Hz幅度');

图1

%物联一班 胡洪 201313060110

%2015年10月24日

%实验二：程序2

n=0:13;xa=n+1;

n=14:26;xb=27-n;

xn=[xa,xb];n=0:26;

subplot(3,2,2);stem(n,xn,'.');

title('三角波序列x(n)');axis([0,32,0,15])

Xk=fft(xn,1024);

k=0:1023;wk=2*k/1024;

subplot(3,2,1);plot(wk,abs(Xk));ylabel('|X(e^(j*w))|');xlabel('w/ pi')

axis([0,1,0,200]);

Xk32=fft(xn,32);

subplot(3,2,3);k=0:31;

stem(k,abs(Xk32),'.');axis([0,16,0,200]);xlabel('k');ylabel('|Xk32|')

xn32=ifft(Xk32);

subplot(3,2,4);k=0:31;

stem(k,xn32,'.');axis([0,32,0,15]);

Xk16=Xk32(1:2:32);

subplot(3,2,5);k=0:15;

stem(k,abs(Xk16),'.');axis([0,8,0,200]);xlabel('k');ylabel('|Xk16')

xn16=ifft(Xk16);

subplot(3,2,6);k=0:15;

stem(k,xn16,'.');axis([0,32,0,15]);

图2

4. 思考题：

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

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

学生实验 心得

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