资源描述:
《数字信号处理实验二,三》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、《数字信号处理》上机实验二实验报告班级1302013学号:13020130082姓名:贾昊稹1.给定模拟信号x(0=2sin(4“+5cos(8“),对其进行采样,用DFT(FFT)进行信号•频谱分析。(1)确定最小采样频率和最小采样点数。(2)若以r=0.01n(^=0:/V-l)秒进行采样,至少需要取多少采样点?(3)用DFT的点数N=50,100®出信号的N点DFT的幅度谱,讨论幅度谱结果。(4)N分别为N=64和N=(O,能否分辨出信号的所有频率分量。(5)在(3)和(4)的条件下做补0FFT,分析结果。(6)
2、在不满足最小采样点数的情况下做补0DFT,观察是否可以分辨出两个频率分量。(1)wl=4*pi;w2=8*pi;fl=wl/2/pi;f2=w2/2/pi;dispf最小采样频率:Jfsl=2*max(flJ2)f=f2-fl;dispC最小釆样点数:’)N=ceil(fsl/abs(f));运行结果:最小采样频率:fsl=最小采样点数:4⑵T=0.01;fs2=l/T;disp(■以T=0.01n釆样,最少的釆样点个数为:’)N0=ceil(fs2/abs(f))运行结杲:以T=0.01n采样,最少的采样点个数为:N
3、0=50(3)(4)N1=50;N2=100;N3=64;N4=60;nl=0:Nl-l;n2=0:N2-l;n3=0:N3-l;n4=0:N4-l;xl=2*cos(wl*nl*T)+5*cos(w2*nl*T);x2=2*cos(wl*n2*T)+5*cos(w2*n2*T);x3=2*cos(wl*n3*T)+5*cos(w2*n3*T);x4=2*cos(wl*n4*T)+5*cos(w2*n4*T);Xl=abs(fft(xl,Nl));X2=abs(fft(x2,N2));X3=abs(fft(x3,N3)
4、);X4=abs(fft(x4,N4));figure(l)subplot(2,2,l),stem(nl,Xl/.');title(,N=50,);xlabel('nl,);ylabel(,Xl,);subplot(222),stem(n2,X2,T);title(,N=100,);xlabel(,n2');ylabel(,X2,);subplot(2/2,3),stem(n3/X3/'.');title('N=64');xlabel('n3');ylabel('X3');subplot(224),stem(n4,X4
5、,T);title('N=60,);xlabel('n4,);ylabel(,X4,);运行结果:(5)Nb=200;nb=O:Nb-l;X5=abs(fft(xl,Nb));X6=abs(fft(x2,Nb));X7=abs(fft(x3,Nb));X8=abs(fft(x4,Nb));figure(2)subplot(221),stem(nb,X5,T);title('N=50补零到200后的幅度值');xlabel('nb');ylabel('X5');subplot(222),stem(nb,X6,T);tit
6、le(*N=100补零到200后的幅度值'hxlabelCnb'bylabelCXe');subplot(223),stem(nb,X7,T);title('N=64补零到200后的幅度值');xlabel('nb');ylabel(X7J;subplot(224),stem(nb,X8,T);title('N=60补零到200后的幅度值,);xlabel(,nb,);ylabel(,X8,);运行结果:(6)N9=3;n9=0:N9-l;x9=2*cos(wl*n9*T)+5*cos(w2*n9*T);X9=abs(
7、fft(x9,Nb));figure(3)stem(nb,X9,T);title('N=3不满足最小采样点时补零到200后的幅度值,);xlabel(,nb,);ylabel(,X9');2.设雷达发射线性调频信号/i(r)=exp(j2^/r2),//=5xio,S采样率/,=2xio9,采样点数N=20000。冋波信号$(/)=〃(/_「)+/(_q),旺=10",丁2=11x10"。(1)画出/()的频谱。(2)利用DFT的时延性质产生W),比较直接在时域产生和在频域产生(再变换到时域)的结果是否相同。(3)匹配
8、滤波的结果是)()“(/)7・(t),(〔〃表示线性卷积)。分别用直接线性卷积和DFT的卷积定理求解⑴。比较二者结呆,并记录两种方法的运行时间(用tic,toe指令)。(4)画出$⑴的频谱。(1)mu=5el3;fs=2e9;Ts=l/fs;N=20000;taol=le-6;tao2=l.le-6;delayl=ceil(tao