数字信号处理实验二,三资料.docxVIP

《数字信号处理》上机实验二实验报告 班级1302013 学号 姓名：贾昊禛 给定模拟信号，对其进行采样，用DFT（FFT）进行信号频谱分析。 确定最小采样频率和最小采样点数。 若以秒进行采样，至少需要取多少采样点？ 用DFT的点数画出信号的点DFT的幅度谱，讨论幅度谱结果。 分别为和，能否分辨出信号的所有频率分量。 在（３）和（４）的条件下做补0 FFT，分析结果。 在不满足最小采样点数的情况下做补0DFT，观察是否可以分辨出两个频率分量。 (1) w1=4*pi;w2=8*pi; f1=w1/2/pi;f2=w2/2/pi; disp(最小采样频率：) fs1=2*max(f1,f2) f=f2-f1; disp(最小采样点数：) N= ceil(fs1/abs(f)); 运行结果： 最小采样频率： fs1 = 8 最小采样点数： N = 4 (2) T=0.01; fs2=1/T; disp(以T=0.01n采样，最少的采样点个数为：) N0=ceil(fs2/abs(f)) 运行结果： 以T=0.01n采样，最少的采样点个数为： N0 = 50 （3）（4）N1=50;N2=100;N3=64;N4=60; n1=0:N1-1; n2=0:N2-1; n3=0:N3-1; n4=0:N4-1; x1=2*cos(w1*n1*T)+5*cos(w2*n1*T); x2=2*cos(w1*n2*T)+5*cos(w2*n2*T); x3=2*cos(w1*n3*T)+5*cos(w2*n3*T); x4=2*cos(w1*n4*T)+5*cos(w2*n4*T); X1=abs(fft(x1,N1)); X2=abs(fft(x2,N2)); X3=abs(fft(x3,N3)); X4=abs(fft(x4,N4)); figure(1) subplot(2,2,1),stem(n1,X1,.); title(N=50);xlabel(n1);ylabel(X1); subplot(2,2,2),stem(n2,X2,.); title(N=100);xlabel(n2);ylabel(X2); subplot(2,2,3),stem(n3,X3,.); title(N=64);xlabel(n3);ylabel(X3); subplot(2,2,4),stem(n4,X4,.); title(N=60);xlabel(n4);ylabel(X4); 运行结果： （5） Nb=200; nb=0:Nb-1; X5=abs(fft(x1,Nb)); X6=abs(fft(x2,Nb)); X7=abs(fft(x3,Nb)); X8=abs(fft(x4,Nb)); figure(2) subplot(2,2,1),stem(nb,X5,.); title(N=50补零到200后的幅度值);xlabel(nb);ylabel(X5); subplot(2,2,2),stem(nb,X6,.); title(N=100补零到200后的幅度值);xlabel(nb);ylabel(X6); subplot(2,2,3),stem(nb,X7,.); title(N=64补零到200后的幅度值);xlabel(nb);ylabel(X7); subplot(2,2,4),stem(nb,X8,.); title(N=60补零到200后的???度值);xlabel(nb);ylabel(X8); 运行结果： （6） N9=3; n9=0:N9-1; x9=2*cos(w1*n9*T)+5*cos(w2*n9*T); X9=abs(fft(x9,Nb)); figure(3) stem(nb,X9,.); title(N=3不满足最小采样点时补零到200后的幅度值);xlabel(nb);ylabel(X9); 设雷达发射线性调频信号，，采样率，采样点数。回波信号，，。 画出的频谱。 利用DFT的时延性质产生，比较直接在时域产生和在频域产生（再变换到时域）的结果是否相同。 匹配滤波的结果是，（“”表示线性卷积）。分别用直接线性卷积和DFT的卷积定理求解。比较二者结果，并记录两种方法的运行时间（用tic，toc指令）。 画出的频谱。 （1） mu=5e13; fs=2e9; Ts=1/fs; N=20000; tao1=1e-6; tao2=1.1e-6; delay1=ceil(tao1/Ts); delay2=ceil(tao2/Ts); n=0:N-1; t=n*Ts; h=exp(j*2*pi*mu*t.^2); H=abs(fft(h,N)); figure(