数字信号处理实验讲义(1-4)资料

1 2sx(n)rejn3clf;(2)单位阶越序列clf;4(3)正弦函数clf;clf;x=2*exp(j*3*n);x=2*exp(j*3*n)');5(5)指数序列clf;z66z=2+j、z=j时的信号实部和虚部图；6clf;6i7clf;非周期的。理论计算中对第一个，N＝2*pi/(0.1*pi)＝10，第二个(1)方波绘图程序如下：8)sinc函数绘图程序如下：9求2．在实验报告中附上在实验过程中记录的各个信号序列的时域和幅频特性离散正弦序列就是一个连续的正弦信号被一系列冲激函数采样后的结的，但采样后的离散序列却不一定是周期的。就是说当w是pi的倍数时，此离散序列才是周期的。所以在本实验中不是pi的倍数。2．离散复指数序列性质：对于离散复指数函数x＝a＊exp(z*n)，只有当z是纯虚数，且纯虚数的系数是pi的倍数时，才是周期的。其它情况下均不是。这个性实验二：信号、系统及系统响应时域特性。③利用卷积方法观察分析系统的时域特性。④掌握序列傅里叶变换的计算机实验方法，利用序列的傅里叶变换对连续号、离散信号及系统响应进行频域分析。B以了解抽样前后信号时域和频域特性发生的变化以及信号信息不丢失的条aa(t)=x(t)δT(t)——(aaTt变换Xˆ(j业)为aaXˆ(j业)=1xwX[j(业-k业)]——(1.3)aTas(1.3)式表明Xˆ(j业)为X(j业)的周期延拓，其延拓周期为抽样角频率aa在计算机上用高级语言编程直接按(1.3)式计算理想抽样(t)的频谱aXˆ(j业)很不方便。下面导出用序列的傅里叶变换来计算Xˆ(j业)的公式。aa将(1.2)式代入(1.1)式并进行傅里叶变换，a=(t)e-j业tdta-www1数字信号处理教程(程佩青)axx(nT)e-j业nTan=-的nn比较(1.5)和(1.4)可知TkM/输出关系为m=-的m(1)认真复习抽样理论、离散信号与系统、线性卷积、序列的傅里叶变换及编制实验用主程序及相应子程序。②系统单位脉冲响应序列产生子程序。本实验要用到两种FIR系统。③有限长序列线性卷积子程序，用于完成两个给定长度的序列的卷积。v流程框图，供实验者参考。(3)(3)调通并运行实验程序，完成下述实验内始序列x(n)调用序列傅氏变换数值计算子程序，求X(ejωk)kYN列产生子程序，求h(n)Hejk)YN调用傅氏变换数值计算子程序，求Y(ejωk)①①分析抽样序列的特性。产生抽样信号序应采用的抽样频率因，并记录(打印)这时的|X(ejω)|曲线。结果一致。调用序列博里叶变换数值计算子程序，求得Y(ejωk),观察(1)在分析理想抽样序列特性的实验中，采样频率不同时，相应理想抽样序列的傅里叶变换频谱的数字频率度量是否都相同？它们所对应的模拟频列的傅里叶变换频谱的数字频率度量都相同，它们都是以2π为周期，只观察分析各种序列的频域特性，通常对X(ejω)在[0，2π]上进x实验三：应用FFT对信号进行频谱分析B(n+1,0n3n0n3始调用信号产生子程序产生实验信号调用绘图子程序(函数)绘制时间序列波形图调用绘图子程序(函数)绘制|X(k)|曲线0,0,0,0,(n+1,(n+1,其他n其他n0,0,0,(n+1,其他n(n+1,其他n0,0,0,0,0,其他n其他n其他n其他nn在实验四离散系统的变换域分析B域方程为Ndy(nk)=px(nk)kkk=0k=0变换域分析方法如下：频域y[n]=x[n]*h[n]=x[m]h[nm]一Y()=X()H()统的频率响应为H()=p()=p0+p1ej+...+pMejM01NZynxnhnxmhnmYzX(z)H(z)系统的转移函数为H(z)=p(z)=D(z)d+dz1+..