数字信号处理上级实验作业.docVIP

数字信号处理Matlab实验作业一 题目：数字信号采样，根据以下情况编写程序并说明。 1泄露； 2防止泄露； 3单一信号正常采样； 4两个信号正常采样； 5不满足采样定理； 6分辨率问题及解决办法； 7两个频率混迭及解决办法 一、泄露 clear all; clc;f1=2; fs=20;N=256;dt=1/fs;t=dt*(0:N-1);x=sin(f1*2*pi*t);y=fft(x,N);figure(1)p=abs(y)/N; f=fs*(0:N-1)/N;plot (f(1:end/2),p(1:end/2));grid axis([0 6 0 0.4]); clear all; clc; f1=2; fs=20; N=256; dt=1/fs; t=dt*(0:N-1); x=sin(f1*2*pi*t); y=fft(x,N); figure(1) p=abs(y)/N; f=fs*(0:N-1)/N; plot (f(1:end/2),p(1:end/2)); grid axis([0 6 0 0.4]); 本程序中Fs=fs=20，N=256，Ft=2，采样时间为12.75秒，因此M=12.75/0.5=25.5,故Fs/N=20/256=0.078125.而Ft/M=2/25.5=0.07843,因此两者不相等，所以造成了泄露问题。 二、解决办法 为了避免频谱泄漏一般有两种方法，一种就是采取相关采样，另外一种则需要用到窗函数进行补偿，减少频谱泄漏，但是不能完全消除，一般窗函数有不同种类，可根据实际情况选择适合的窗函数。6 Z??} k7 i* A% P关于相关采样，需要遵循一定的规则：1）需要将被采样信号划分周期；2）Fs/N=Ft/M, Fs为采样频率，N为采样点的数目，Ft为被采样信号频率，M为被采样信号所采样的周期的数目；3）此外，M和N需要互为素数且都需要为整数，这样可避免重复采样，提高效率。? ?? ? 换句话说，就是我们用FFT算法的时候，信号在时域的时候，应该为周期信号，这样进行采样的话，就不会发生诸如频谱泄漏的现象。方法一 采用相关采样 clear all; clc;f1=2; fs=32; N=256; dt=1/fs; t=dt*(0:N-1); x=sin(f1*2*pi*t); y=fft(x,N); figure(1) p=abs(y)/N; f=fs*(0:N-1)/N; plot (f(1:end/2),p(1:end/2)); grid axis([0 6 0 0.6]); 分析：本程序中Fs=fs=32，N=256，Ft=2，采样时间为7.9688秒，近似认为是8秒，因此M=8/0.5=16,故Fs/N=32/256=0. 125.而Ft/M=2/16=0. 125,因此两者相等，所以避免了泄露问题。 方法二： 加窗 clear all; clc; f1=2; fs=20; N=256; dt=1/fs; t=dt*(0:N-1); x=sin(f1*2*pi*t); han1=hanning(N); x=x.*han1; y=fft(x,N); figure(1) p=abs(y)/N; f=fs*(0:N-1)/N; plot (f(1:end/2),p(1:end/2)); grid axis([0 6 0 0.4]); 分析：同样是20赫兹的采样频率，加窗后的效果好了很多 Hamming窗比hanning窗的效果更好一些 三、单一信号正常采样 %single signal frequency normal clear all; clc; fs=32; f=4; N=256; t=0:1/fs:256/fs; a=10*cos(2*pi*f*t); figure(1); plot (t,a); i=1:256; fi(i)=(i-1)*fs/N; b=fft(a); figure(2); plot (fi(1:end/2),abs(b(1:end/2))); 四、两个信号叠加正常采样 %two normal frequency clear all clc fs=64; f=2; N=512; t=0:1/fs:512/fs; a=10*cos(2*pi*f*t)+10*cos(2*pi*4*f*t); figure(1); plot (t,a); i=1:512; fi(i)=(i-1)*fs/N; b=fft(a); figure(2); plot (fi(1:end/2),abs(b(1:end/2))); 五、不满足采样定理的采样 clear all; c