DSP实验指导书.doc

DSP 实验指导书 实验1 抽样定理 一．实验目的 了解电信号的采样与恢复原理。 验证采样定理。 二．实验仪器 实验室：信号与系统实验室 设备： THSCC-1信号与系统实验箱一台。 示波器一台。 三．实验内容 将连续的时间信号进行采样，观察采样后的离散时间信号波形。 将离散时间信号恢复，观察恢复后的信号，并与原信号进行比较。 改变采样频率，观察采样后的离散时间信号波形。观察恢复后的信号，并与原信号进行比较。 连续改变采样频率，观察采样后的离散时间信号波形。观察恢复后的信号，并与原信号进行比较，从失真到开始不失真，确定最低的采样频率f s (t)。 四. 实验步骤： （1）按下正弦波和f2波段开关，从实验箱输出一个100Hz~300Hz、幅度为1V的正弦波f (t)，加到采样器的输入端f (t)。 （2）用双踪示波器分别观察输入f (t)与采样信号S（t）的波形。并测量出频率（或周期）。 （3）用双踪示波器分别观察离散信号f s (t)与恢复信号f’（t）的波形，比较异同之处。 （4）用双踪示波器分别观察输入信号f (t)与恢复信号f’（t）的波形，比较异同之处。 （5）调节采样频率旋钮改变采样频率f s (t)，直到f’（t）失真。反向调节采样频率旋钮改变采样频率f s (t)，直到f’（t）开始不失真（即输入信号f (t)与恢复信号f’（t）相同），使恢复的正弦波信号f’(t)处于刚好不失真的临界点。测量出此时离散信号S（t）’（t）波形。分析各波形结果。 3. 测出最低的采样频率f s (t)。 4.思考题：按抽样定理：fs = 2 f max(t)，即可不失真还原，实际测出的最低的采样频率f s (t)远大于理论值2 f max(t)，为什么？ 5. 恢复信号时失真产生的原因是什么？ 实验2 离散傅里叶变换及其快速算法 采样，截取长度N分别选N =20和N =16，观察其DFT结果的幅度谱。 解 此时离散序列 ，即k=8。用MATLAB计算并作图，函数fft用于计算离散傅里叶变换DFT，程序如下： k=8; N1=20； n1=[0:1:N1-1]; xa1=sin(2*pi*n1/k); subplot(2,2,1) plot(n1,xa1) xlabel(t/T);ylabel(x(n)); xk1=fft(xa1);xk1=abs(xk1); subplot(2,2,2) stem(n1,xk1) xlabel(k);ylabel(X(k)); N2=16； n2=[0:1:N2-1]; xa2=sin(2*pi*n2/k); subplot(2,2,3) plot(n2,xa2) xlabel(t/T);ylabel(x(n)); xk2=fft(xa2);xk2=abs(xk2); ??????????????????????? subplot(2,2,4) stem(n2,xk2) xlabel(k);ylabel(X(k)); 计算结果如图，(a)和(b)分别是N=20时的截取信号和DFT结果，由于截取了两个半周期，频谱出现泄漏；(c) 和(d) 分别是N=16时的截取信号和DFT结果，由于截取了两个整周期，得到单一谱线的频谱。上述频谱的误差主要是由于时域中对信号的非整周期截断产生的频谱泄漏。 实验3 离散卷积 （1）线性卷积 例：y=x*h x=[1 2 3 2 ]; h=[1 1 1]; %线性卷积 y=conv(x,h); N=length(y)-1; n=0:1:N; stem(n,y); title(线性卷积结果); xlabel(( n ));ylabel(线性卷积 y(n)); (2)FFT卷积 x1=[1 2 3 2 ]; x2=[1 1 1]; N1=length(x1); N2=length(x2); N=N1+N2-1; %卷积 X1=fft(x1,N); X2=fft(x2,N); Y=X1.*X2; %点乘：频谱对应元素相乘 y=ifft(Y,N); n=0:1:N-1; stem(n,y); title(FFT 线性卷积结果); ylabel(卷积); (3)圆周卷积 %定义函数，单独保存cirshiftd.m %位移函数 function xm= cirshiftd(x,m,N); n=[0:1:N-1]; xm=x(mod(n-m,N)+1); %圆周卷积圆周卷积调用 x1=[1 2 3 2 ]; x2=[1 1 1]; %步骤1定义长度 N=length(x1)+ length(x2)-1; %N=8; %步骤2. 配置序列周期、长度