实验二时域采样与频域采样Word下载.docx
《实验二时域采样与频域采样Word下载.docx》由会员分享,可在线阅读,更多相关《实验二时域采样与频域采样Word下载.docx(7页珍藏版)》请在冰豆网上搜索。
Xa(t)d(t-nT)eJdt
+?
?
)?
a(jvy=蝌I[Xa(t)邋d(t-nT)e-jVVdt=
n=-?
n-?
在上式的积分号内只有当t=nT时,才有非零值,因此
*a(jvy二?
Xa(nT)e
n二-?
上式中,在数值上Xa(nT)=x(n),再将w=WT代入,得到:
X^jV)二?
Xa(n)e-jnww=wt=X(ejw)w=wtn二-?
上式说明采样信号的傅里叶变换可用相应序列的傅里叶变换得到,只要将自变量
用T代替即可。
2频域采样定理
对信号x(n)的频谱函数X(ej)在[0,2]上等间隔采样N点,得到
X(k)=X(ejw)2pk=0,1,2,L,N-1
w=——k
N
则有:
XN(n)=IDFT[X(k)h=[?
x(n+iN)]Rn(n)
i=-?
即N点IDFT[X(k)]得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值序列,
因此,频率域采样要使时域不发生混叠,则频域采样点数N必须大于等于时域
离散信号的长度M(即N3M)。
在满足频率域采样定理的条件下,XN(n)就是原序列x(n)。
如果N>
M,则Xn(n)比原序列x(n)尾部多N-M个零点,反之,时域发生混叠,Xn(n)与x(n)不等。
对比时域采样定理与频域采样定理,可以得到这样的结论:
两个定理具有对偶性,即“时域采样,频谱周期延拓;
频域采样,时域信号周期延拓”。
在数字
信号处理中,都必须服从这二个定理。
三实验内容
1.时域采样实验:
%时域采样实验
A=444.128;
a=50*sqrt
(2)*pi;
w0=50*sqrt
(2)*pi;
邛=64/1000;
F1=1000;
F2=300;
F3=200;
%观察时间,Tp=64ms
T1=1/F1;
T2=1/F2;
T3=1/F3;
%不同的米样频率
%产生不同的长度区间
n1=0:
Tp*F1-1;
n2=0:
Tp*F2-1;
n3=0:
Tp*F3-1;
stem(n1,x1,'
.'
)
n1,n2,n3
x1=A*exp(-a*n1*T1).*sin(w0*n1*T1);
x2=A*exp(-a*n2*T2).*sin(w0*n2*T2);
x3=A*exp(-a*n3*T3).*sin(w0*n3*T3);
f1=fft(x1,length(n1));
f2=fft(x2,length(n2));
f3=fft(x3,length(n3));
k1=0:
length(f1)-1;
fk1=k1/Tp;
k2=0:
length(f2)-1;
fk2=k2/Tp;
k3=0:
length(f3)-1;
fk3=k3/Tp;
subplot(3,2,1)
)%此处也可用
title('
(1)Fs=1000Hz'
);
xlabel('
n1'
ylabel('
x1(n)'
%产生采样序列x1(n)
%产生采样序列x2(n)
%产生采样序列x3(n)
%采样序列x1(n的FFT变换
%采样序列x2(n)的FFT变换
%采样序列x3(n)的FFT变换
%x1(n)的频谱的横坐标的取值
%x2(n)的频谱的横坐标的取值
%x3(n)的频谱的横坐标的取值
gribon;
%添加网络线
subplot(3,2,3)
stem(n2,x2,'
(3)Fs=300Hz'
n2'
x2(n)'
subplot(3,2,5)
stem(n3,x3,'
(5)Fs=200Hz'
n3'
x3(n)'
subplot(3,2,2)
plot(fk1,abs(f1))
(2)FT[xa(nT)],Fs=1000Hz'
f(Hz)'
幅度'
subplot(3,2,4)
plot(fk2,abs(f2))
(4)FT[xa(nT)],Fs=300Hz'
幅度'
subplot(3,2,6)
plot(fk3,abs(f3))
(6)FT[xa(nT)],Fs=200Hz'
)gribon;
时域采样波形:
2.频域采样实验
%频域采样实验
xa=0:
floor(M/2);
%floor
是向下取整例如floor(2.5)=2
xb=ceil(M/2)-1:
-1:
0;
%ceil(M/2)是取大于等于M/2的最小整数
xn=[xa,xb];
Xk=fft(xn,1024);
%1024
点FFT[x(n)],用于近似序列x(n)的TF
X32k=fft(xn,32);
%32
点FFT[x(n)]
x32n=ifft(X32k);
%32
点IFFT[X32(k)]得到x32(n)
X16k=X32k(1:
2:
N);
%
隔点抽取X32k得到X16(K)
x16n=ifft(X16k,N/2);
%16点IFFT[X16(k)]得到x16(n)subplot(3,2,2);
stem(n,xn,'
boxon
(2)三角波序列x(n)'
n'
x(n)'
axis([0,32,0,20])
k=0:
1023;
wk=2*k/1024;
subplot(3,2,1);
plot(wk,abs(Xk));
(1)FT[x(n)]'
\omega/\pi'
|X(eAjA\omega)|'
axis([0,1,0,200])k=0:
N/2-1;
subplot(3,2,3);
stem(k,abs(X16k),'
(3)16点频域采样'
k'
|X_1_6(k)|'
axis([0,8,0,200])
subplot(3,2,4);
stem(n1,x16n,'
(4)16点IDFT[X_1_6(k)]'
x_1_6(n)'
N-1;
subplot(3,2,5);
stem(k,abs(X32k),'
(5)32点频域采样'
|X_3_2(k)|'
axis([0,16,0,200])
subplot(3,2,6);
stem(n1,x32n,'
(6)32点IDFT[X_3_2(k)]'
x_3_2(n)'
频域采样波形:
四思考题
如果序列x(n)的长度为M,希望得到其频谱X(ej)在[0,2]上N点等间隔采样,当N<
M时,如何用一次最少点数的DFT得到该频谱采样?
答:
n<
m时,频域抽样不够密,(x)n以周期进行延拓,频域产生混叠,抽样信号不能还原原信号。
可将m分为n长度的k段,不足时域补零。
分段进行DFT。
此时DFT点数最少为N次。
五实验报告及要求
(1)由上图可得:
时域采样,对连续的信号进行等间隔采样形成采样信号,采样信号的频谱是原连续信号的频谱已采样信号为周期进行周期性的延拓形成的。
(2)由上图可得:
时域采样,采样频率越高,时域内信号分辨率就越高,采集到的信号就越接近原始信号,在频谱上的频带就越宽。
这有利于后期频域分析相位分量的相位改变是不影响该波的频率成分和幅值大
小,也就是说,在幅频内的本质是没有发生改变的,所以最终合成的波形幅值频谱是不会改变的