数字信号处理实验二,三.docx

《数字信号处理》上机实验二实验报告班级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(1);stem(n,H,.);运行结果：（2）Ns=N+max(delay1,delay2);s1=zeros(1,Ns);s2=zeros(1,Ns);s1(delay1+1:delay1+N)=exp(2*j*pi*mu*t.^2);s2