数字信号处理实验一信号、-系统及系统响应.docVIP

西 安 郵 電 學 院 数字信号处理课内实验 报 告 书 系部名称 ： 计算机系 学生姓名 ： 常成娟 专业名称 ： 电子信息科学与技术 班 级 ： 电科0603 学号 ：22号） 时间 ： 200 实验一： 信号、 系统及系统响应 一. 实验目的 (1) 熟悉连续信号经理想采样前后的频谱变化关系， 加深对时域采样定理的理解。 (2) 熟悉时域离散系统的时域特性。 (3) 利用卷积方法观察分析系统的时域特性。 (4) 掌握序列傅里叶变换的计算机实现方法， 利用序列的傅里叶变换对连续信号、 离散信号及系统响应进行频域分析。 二. 实验原理与方法 采样是连续信号数字处理的第一个关键环节。 对一个连续信号xa(t)进行理想采样的过程可用(10.3.1)式表示。 (10.3.1)其中 (t)为xa(t)的理想采样， p(t)为周期冲激脉冲， 即 (10.3.2) (t)的傅里叶变换 (jΩ)为 (10.3.3) 将(10.3.2)式代入(10.3.1)式并进行傅里叶变换， (10.3.4) 式中的xa(nT)就是采样后得到的序列x(n)， 即 x(n)的傅里叶变换为 (10.3.5) 比较(10.3.5)和(10.3.4)可知 （10.3.6) 在数字计算机上观察分析各种序列的频域特性，通常对X(ejω)在［0, 2π］上进行M点采样来观察分析。 对长度为N的有限长序列x(n), 有 (10.3.7) 其中 一个时域离散线性非移变系统的输入/输出关系为 (10.3.8) 上述卷积运算也可以在频域实现 图10.3.1 实验一的主程序框图 三. 实验内容及步骤 (1) 认真复习采样理论、 离散信号与系统、 线性卷积、 序列的傅里叶变换及性质等有关内容， 阅读本实验原理与方法。 (2) 编制实验用主程序及相应子程序。 ① 信号产生子程序， 用于产生实验中要用到的下列信号序列： xa(t)=Ae-at sin(Ω0t)u(t) 进行采样， 可得到采样序列 xa(n)=xa(nT)=Ae-anTsin(Ω0nT)u(n), 0≤n50 其中A为幅度因子， a为衰减因子， Ω0是模拟角频率， T为采样间隔。 这些参数都要在实验过程中由键盘输入， 图10.3.2 xa(t)的幅频特性曲线 产生不同的xa(t)和xa(n)。 b. 单位脉冲序列： xb(n)=δ(n) c. 矩形序列： xc(n)=RN(n), N=10 ② 系统单位脉冲响应序列产生子程序。 本实验要用到两种FIR系统。 a. ha(n)=R10(n); b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3) ③ 有限长序列线性卷积子程序， 用于完成两个给定长度的序列的卷积。 可以直接调用MATLAB语言中的卷积函数conv。 conv用于两个有限长度序列的卷积， 它假定两个序列都从n=0 开始。 调用格式如下： y=conv (x, h) (3) 调通并运行实验程序， 完成下述实验内容： ① 分析采样序列的特性。 a. 取采样频率fs=1 kHz, 即T=1 ms。 b. 改变采样频率， fs=300 Hz， 观察|X(ejω)|的变化， 并做记录(打印曲线)； 进一步降低采样频率， fs=200 Hz， 观察频谱混叠是否明显存在， 说明原因， 并记录(打印)这时的|X(ejω)|曲线。 源程序： A=444.128; a=50*sqrt(2)*pi; w=50*sqrt(2)*pi; n=0:49; fs=1000; x=A*exp((-a)*n/fs).*sin(w*n/fs); k=-200:200; w=(pi/100)*k; y=x*(exp(-j*pi/100)).^(n*k); %y=fft(x) subplot(1,2,1); stem(n,x); axis([0,50,-50,150]); xlabel(n); ylabel(Xa(n)); title(fs=1000); subplot(1,2,2); plot(w/pi,abs(y)) axis([-2,2,0,1000]); xlabel(w/pi); ylabel(/Xa(ejw)/); A=444.128; a=50*sqrt(2)*pi; w=50*sqrt(2)*pi; n=0:49; fs=500; x=A*exp((-a)*n/f