数字信号处理西电.docxVIP

数字信号处理上机第一次实验 实验一： 设给定模拟信号，的单位是ms。 利用MATLAB绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 用两个不同的采样频率对给定的进行采样。 eq \o\ac(○,1) 。 eq \o\ac(○,2) 。 比较两种采样率下的信号频谱，并解释。 实验一MATLAB程序： （1） eq \o\ac(○,1) clc; fs=5000; ts=1/fs; N=1000; t=(-N:N)*ts; s=exp(-abs(t)); plot(t,s,linewidth,1.5) xlabel(时间) ylabel(幅度) set(gca,fontweight,b,fontsize,12) SPL=N*100; figure sp=fftshift(fft(s,SPL)); sp=sp/max(sp)*100; freqb=-fs/2:fs/SPL:fs/2-fs/SPL; plot(freqb,abs(sp)) xlabel(频率) ylabel(频谱幅度) set(gca,fontweight,b,fontsize,12) yy=abs(abs(sp)-3); [aa,freqind]=min(yy); (freqind-SPL/2)*fs/SPL eq \o\ac(○,2) clc; fs=1000; ts=1/fs; N=1000; t=(-N:N)*ts; s=exp(-abs(t)); plot(t,s,linewidth,1.5) xlabel(时间) ylabel(幅度) set(gca,fontweight,b,fontsize,12) SPL=N*100; figure sp=fftshift(fft(s,SPL)); sp=sp/max(sp)*100; freqb=-fs/2:fs/SPL:fs/2-fs/SPL; plot(freqb,abs(sp)) xlabel(频率) ylabel(频谱幅度) set(gca,fontweight,b,fontsize,12) yy=abs(abs(sp)-3); [aa,freqind]=min(yy); (freqind-SPL/2)*fs/SPL 实验三： 设，，编写MATLAB程序，计算： ５点圆周卷积； ６点圆周卷积； 线性卷积； 画出的，和时间轴对齐。 a = [1,2,2]; b = [1,2,3,4]; y1 = cconv(a,b,5) y2 = cconv(a,b,6) y3 = conv(a,b) figure(1); subplot(311) stem(y1); grid on title(五点圆周卷积y1(n)); xlabel(n),ylabel(y1(n));axis([0 6 0 15]) subplot(312) stem(y2); grid on title(六点圆周卷积y2(n)); xlabel(n),ylabel(y2(n));axis([0 6 0 15]) subplot(313) stem(y3); grid on title(线性卷积y3(n)); xlabel(n),ylabel(y3(n));axis([0 6 0 15]) 实验四： 给定因果系统： 求系统函数并画出零极点示意图。 画出系统的幅频特性和相频特性。 求脉冲响应并画序列图。 提示：在ＭＡＴＬＡＢ中，zplane(b,a) 函数可画零极点图；Freqz(b,a,N)可给出范围内均匀间隔的点频率响应的复振幅；Impz(b,a,N)可求的逆变换（即脉冲响应）。 clc a = [1,0] b = [1,-0.9] figure(1) zplane(b,a); title(零极点分布图) w=[-3*pi:0.01:3*pi]; [h,phi]=freqz(b,a,w); figure(2); subplot(2,1,1); plot(w, abs(h)); grid on; title(幅频特性); xlabel(f/Hz),ylabel(H(f)); subplot(2,1,2); plot(w, phi); grid on; title(相频特性); xlabel(f/Hz),ylabel(W(f)); 数字信号处理第二次实验 给定模拟信号，对其进行采样，用DFT（FFT）进行信号频谱分析。 确定最小采样频率和最小采样点数。 若以秒进行采样，至少需要取多少采样点？ 用DFT的点数画出信号的点DFT的幅度谱，讨论幅度谱结果。 分别为和，能否分辨出信号的所有频率分量。 在（３）和（４）的条