试验三时域采样与频域采样

1、%求出序列长度%位置为左上%时域波形%标题%求出频率%位置为右上%幅频特性曲线%标题%观察时间Tp=64ms%采样率300hz%采样间隔%fft给定点数%序列从0开始，至少要达到fft的点数实验三程序代码及实验结果图：(1)时域采样理论的验证。 给定模拟信号，？?？= ?产gin a。？?,现用DFT(FFT) 求该模拟信号的幅频特性，以验证时域采样理论。按照？?？勺幅频特性曲线，选取三种采样频率，即?=1kHz , 300Hz , 200Hz。观测时间选？?= 50?要求：编写实验程序，计算 ? ?、? ?和？ ？的幅度特性，并绘图显示，观察分析频谱混叠失真。实验程序代码及 结果如下：Tp=

2、64/1000;%观察时间 Tp=64msFs=1000;% 采样率 1khzT=1/Fs;%采样间隔M=Tp*Fs;%fft给定点数n=0:M-1;%序列从0开始，至少要达到 fft的点数%产生模拟信号对应的离散序列x1(n)A=444.128;alph=pi*50*2A0.5;omega=pi*50*2A0.5;xnt=A*exp(-alph*n*T).*sin(omega*n*T);% 无需乘 u(t),因为序列从 0 开始Xk=T*fft(xnt,M);%频谱函数n1 =0: (length(xnt)-1);n2 =0: (length(Xk)-1);subplot(3,2,1);st

3、em(n1,xnt);title(时域1000hz采样波形,);fk=n2/Tp;subplot(3,2,2);stem(fk,abs(Xk);title(频域 1000hz 采样)；%300hzTp=64/300;Fs=300;T=1/Fs;M=Tp*Fs;n=0:M-1;%产生模拟信号对应的离散序列x1(n)A=444.128;alph=pi*50*2A0.5;omega=pi*50*2A0.5;xnt=A*exp(-alph*n*T).*sin(omega*n*T);% 无需乘 u(t),因为序列从 0 开始Xk=T*fft(xnt,M);%频谱函数n1 =0: (length(xnt)

4、-1);%求出序列长度fft的点数n2 =0: (length(Xk)-l);subplot(3,2,3);stem(n1,xnt);title(时域 300hz 采样)；fk=n2/Tp;subplot(3,2,4);stem(fk,abs(Xk);title(频域 300hz 采样)；%200hzTp=64/200;Fs=200;T=1/Fs;M=Tp*Fs;n=0:M-1;%位置为左中%时域波形标题%求出频率%位置为右中%幅频特性曲线标题%观察时间Tp=64ms%采样率200hz%采样间隔%fft给定点数%序列从0开始，至少要达到%产生模拟信号对应的离散序列x1(n)A=444.128;

5、alph=pi*50*2A0.5;omega=pi*50*2A0.5;0开始xnt=A*exp(-alph*n*T).*sin(omega*n*T);% 无需乘 u(t),因为序列从Xk=T*fft(xnt,M);%频谱函数n1 =0: (length(xnt)-l);% 求出序列长度n2 =0: (length(Xk)-1);subplot(3,2,5);stem(n1,xnt);%时域波形title(时域 200hz 采样)；%位置在右下%幅频特性曲线%标题fk=n2/Tp;subplot(3,2,6);stem(fk,abs(Xk);title(频域 200hz 采样);(2)频域采样理

6、论的验证。给定信号如下：n 10 _n _13x(n) =27 -n 14 n 26 0 其它 J编写程序分别对频谱函数 X ?/? = ?q?在区间0,2?上等间隔采样32和16点。分 别画出？?2 ?、？?6(?和X ?的幅度谱，并绘图显示 ？?？、?2 ?和？?6 ?的波形，进行 对比和分析，验证总结频域采样定理。实验程序代码及结果如下：M =27;N = 32;n = 0 : M;xa = 0 : floor(M / 2);xb = ceil(M / 2) - 1 : -1 : 0;xn = xa,xb;%三角波输入序列Xk = fft(xn,1024);X32k = fft(xn,3

7、2);x32n = ifft(X32k);X16k = X32k(1:2:N);x16n = ifft(X16k);%1024点fft,近似看作原信号的频谱%32点fft , 32点频域采样%32点频域采样还原至时域%对32点频域采样作隔点抽取，得到 16点频域采样%16点频域采样还原至时域离散subplot(3,2,2);stem(n,xn,.);title(三角波序列x(n);xlabel(n);ylabel(x(n);axis(0,32,0,20);%位置为右上%标题%横坐标轴名称%纵坐标轴名称%坐标轴范围k = 0 : 1023;wk = 2 * k / 1024;subplot(3,

8、2,1); plot(wk,abs(Xk);title(FTx(n);xlabel(omega / pi); ylabel(|X(e A j A omega)|); axis(0,1,0,200);%此处如取0, 2, 0, 200,k = 0 : N / 2 - 1;subplot(3,2,3);stem(k,abs(X16k),.);title(16点频率采样,);xlabel(k);ylabel(|X16k|); axis(0,8,0,200);n1 = 0 : N / 2 - 1;subplot(3,2,4);stem(n1,x16n,.);title(16 点 IDFTX16k);xlabel(n);ylabel(x16n); axis(0,32,0,20);k = 0 : N - 1;subplot(3,2,5);stem(k,abs(X32k),.);title(32点频域采样,);xlabel(k);ylabel(|X32k|); axis(0,16,0,200);n1 = 0 : N - 1;subplot(3,2,6);stem(n1,x32n,.);title(32 点 IDFTX32k);xlabel(n);ylabel(x32n); axis(0,32,0,20);%位置为左上%标题%横坐标轴名称%纵坐标轴名称%坐标轴范围则后面的N