数字信号处理技术研究与应用.pptxVIP

数字信号处理技术研究与应用1：卷积计算项目要求：编写序列起点不在零点时的卷积程序。用所学列表法或作图法计算以下两序列卷积，并用程序验证（以 x=[1,2,3]; h=[3,4,5,6,7]为例，需要计算结果及作图），并与例3-4序列起点在零点时的结果对比。下划线代表n=0的点。要求做良好的人机界面 x=[1,2,3]; h=[3,4,5,6,7]; y=conv(x,h) a1=-1;a2=1; b1=-1;b2=3; A1=a1+b1;B2=a2+b2; C=A1:B2; stem(C,y);通过程序设定x、h的起点在-1处，通过运行程序得到其输出图形因为conv函数仅在序列起点为零点时适用，因此若不设定x、h的起点为非零点时，其输出如下验证时域采样定理（1）书上的例子：分别以500HZ和200HZ对f0进行采样。 n=0:1000;f0=150;fs1=500;fs2=200; x1=sin(2*pi*f0*n/fs1); x2=sin(2*pi*f0*n/fs2); X1=fft(x1,512); X2=fft(x2,512); subplot(2,1,1),stem(X1); subplot(2,1,2),stem(X2);运行程序得输出结果为（2）再分别以100HZ和50HZ对f0进行128点采样。N=0:40;f0=40;fs1=100;fs2=50;x1=sin(2*pi*f0*N/fs1);x2=sin(2*pi*f0*N/fs2);X1=fft(x1,128);X2=fft(x2,128);subplot(2,1,1),stem(X1);subplot(2,1,2),stem(X2);结论 由输出图像可知采样1在频率f0之外的频率处频谱几乎均为0，而采样2在频率f0之外的频率处频谱并不为0，这是因为fs12f0，满足采样定理，而fs22f0，不满足采样定理，发生了频率混叠现象。可见只有在满足采样定理的条件下进行采样才能避免频谱混叠现象。?3：矩形序列频谱分析（1）设L=10，画出有限长序列RL(n)的幅度谱；（2）在上述矩形序列后面补10个0，画出序列的幅度谱，要求（1）（2）的幅度谱画在同一坐标系下，横坐标统一以pi为单位n=0:1:9;x=[ones(1,10),zeros(1,0)];k=-200:0;w1=k*pi/100;X=x* exp(-j).^(n*w1);nagX1=abs(X);subplot (2,2,1);stem(n,x);title(10点);xlabel(n);ylabel(x);subplot (2,2,2);plot(nagX1);title(非周期有限长序列幅度普L=10 N=10);n1=0:1:19;x1=[ones(1,10),zeros(1,10)];k=-200:0;w1=k*pi/100;X1=x1* exp(-j).^(n1*w1);nagX2=abs(X1);subplot (2,2,3);stem(n1,x1);title(20点);xlabel(n);ylabel(x);subplot (2,2,4);plot(nagX2);title(非周期有限长序列幅度普L=10 N=20);4：栅栏效应项目要求：以教材例4-11为例，分析栅栏效应与谱分辨率的不同。（1）教材上程序 N=128;fs=100;n=0:N-1;f1=20;f2=20.5;f3=30; x=sin(2*pi*n*f1/fs)+sin(2*pi*n*f2/fs)+sin(2*pi*n*f3/fs); X=fft(x); AX=abs(X(1:N/2)); stem(AX); 对x(n)截取128点，其幅度频谱如下：对x(n)截取512点，其程序如下： N=512;fs=100;n=0:N-1;f1=20;f2=20.5;f3=30; x=sin(2*pi*n*f1/fs)+sin(2*pi*n*f2/fs)+sin(2*pi*n*f3/fs); X=fft(x); AX=abs(X(1:N/2)); stem(AX);分辨率： 由于128点不满足最少采样点数的要求，所以DFT结果图中无法分辨五序列中所有频率成分，及分辨率太小。当截取512点，满足分辨率最小采样点数要求，因此能区分序列的所有频率成分。最小分辨率为F=f2-f1=0.5,则公式：F=Fs/Nmin 式中Fs=100HZ，求出最少截取点数Nmin=200。结论栅栏效应： 当截取点数大于等于最小采样点数时不会出现栅栏效应，及N＞=Nmin时，512＞200，此时不会出现栅栏效应。相反，当截取点数小于最小采样点数时，128＜200，就会产生栅栏效应，产生误差。 由以上现象可以得出，无法区分序列的所有频率成分的原因可能是由于截取