数字信号处理第三版用MATLAB上机实验VIP专享

实验二：时域采样与频域采样 一、 时域采样 1.用 MATLAB 编程如下: %1 时域采样序列分析fs=1000 A=444.128; a=222.144; w=222.144; ts=64*10^(-3); fs=1000;T=1/fs; n=0:ts/T-1; xn=A*exp((-a)*n/fs).*sin(w*n/fs); Xk=fft(xn); subplot(3,2,1);stem(n,xn);xlabel('n,fs=1000Hz'); ylabel('xn');title('xn'); subplot(3,2,2);plot(n,abs(Xk));xlabel('k,fs=1000Hz'); title('|X(k)|'); %1 时域采样序列分析fs=200 A=444.128; a=222.144; w=222.144; ts=64*10^(-3); fs=200;T=1/fs; n=0:ts/T-1; xn=A*exp((-a)*n/fs).*sin(w*n/fs); Xk=fft(xn);subplot(3,2,3);stem(n,xn);xlabel('n,fs=200Hz'); ylabel('xn');title('xn'); subplot(3,2,4);plot(n,abs(Xk));xlabel('k,fs=200Hz'); title('|X(k)|'); %1 时域采样序列分析fs=500 A=444.128; a=222.144; w=222.144; ts=64*10^(-3); fs=500;T=1/fs; n=0:ts/T-1; xn=A*exp((-a)*n/fs).*sin(w*n/fs); Xk=fft(xn); subplot(3,2,5);stem(n,xn);xlabel('n,fs=500Hz'); ylabel('xn');title('xn'); subplot(3,2,6);plot(n,abs(Xk));xlabel('k,fs=500Hz'); title('|X(k)|'); 2.经调试结果如下图: 020406080-2000200n,fs=1000Hzxnxn02040608005001000k,fs=1000Hz|X(k)|051015-2000200n,fs=200Hzxnxn0510150100200k,fs=200Hz|X(k)|010203040-2000200n,fs=500Hzxnxn0102030400500k,fs=500Hz|X(k)| 实验结果说明:对时域信号采样频率必须大于等于模拟信号频率的两倍以上,才 能使采样信号的频谱不产生混叠.fs=200Hz 时，采样信号的频谱产生了混叠，fs=500Hz 和 fs=1000Hz 时，大于模拟信号频率的两倍以上，采样信号的频谱不产生混叠。 二、 频域采样 1.用MATLAB 编程如下: %频域采样定理验证 M=26;N=32;n=0:M; n1=0:13;x1=n1+1; n2=14:26;x2=27-n2; x=[x1,x2];Xk=fft(x,512); X32k=fft(x,32); k=0:511;w=(pi/512)*k; subplot(321);stem(n,x);xlabel('n'); ylabel('xn');axis([0,31,0,15]); subplot(322);plot(w,abs(Xk));xlabel('k'); ylabel('|X(k)|');axis([0,1,0,200]) X16k=X32k(1:2:N); x32n=ifft(X32k);x16n=ifft(X16k,16); k1=0:31;k2=0:15; subplot(323);stem(k1,abs(X32k));xlabel('k'); ylabel('X32k');axis([0,31,0,200]); subplot(325);stem(k2,abs(X16k));xlabel('k'); ylabel('|X(k)|');axis([0,15,0,200]) n=0:31; subp...