dsp数字信号处理实验报告华工.doc

实验1常见离散信号产生与实现 一、 实验目得 加深对常用离散信号得理解； 2、 熟悉使用MATLAB在时域中产生一些基本得离散时间信号。 二、 实验原理 单位抽样序列 在MATLAB中可以利用Z纣兀0函数实现。 x = 1 zeros(l, N—1) 2、单位阶越序列 0 n 0 在MATLAB中可以利用⑷佗SO函数实现: x = onss( ; 3、正弦序列 x(w) = ^4sin( 2^fhFs + 0) 在MATLAB中实现过程如下： n = O:N?1 x = A* sin( 2*pz */* rdFs + jai) 4、复指数序列 x(n) = r -穴 在MATLAB中实现过程如下： 丹=0：时一1 x= r exp(j*m*li) 5、指数序列 x(ri) = tz: 在MATLAB中实现过程如下： n = 0:N-l 无=a. n 三、预习要求 预先阅读实验讲义(MATLAB基础介绍); 2、 讨论正弦序列、复指数序列得性质。 6时得信号实部与虚部图；当 6时呢？此时信号周期为多少？ 程序如下： titlezl=-l/12+j *pi/6; titlez2 = l/12+j *pi/6; z3=l/12; z4=2+j *pi/6; z5=j*pi/6; n=0:20; xl=exp(titlezl*n); x2=exp(titlez2*n); x3=exp(z3*n); x4=exp(z4*n); x5=exp(z5*n); subplot (5,2,1); stem(nz real(xl)); xlabel(n1);ylabel(1 real(xl) 1); title(1zl=-l/12+j *pi/6 时■) subplot(5,2,2); stem(n,imag(xl)); xlabel(1n *);ylabel(* imag(xl) 1); title ( * zl=-l/12+j *pi/6 II寸) subplot(5,2,3); stem(nz real(x2)); xla bel ( ) ; yla b ( 1 real (x2) 1 ); title ( * z2 = l/12+j*pi/6 (I寸?) subplot(5,2,4); stem(nfimag(x2)); xlabel(* n *);ylabel(1 image(x2) 1); title(* z2 = l/12+j *pi/6 时?) subplot(5,2,5); stem(nz real(x3)); xlabel(W);ylabJ(* real(x3) 1); title (1 z3=l/12 (14* ) subplot(5,2,6); stem(nz imag(x3)); xlabel(1n1);ylabel(1 image(x3) 1); title(* z3=l/12 时，) subplot(5,2,7); stem(nr real(x4)); xlabel();ylabel(1 real(x4) 1 ); title (,z4=2-Fj*pi/6 Ihf ) subplot(5,2,8); stem(n,imag(x4)); xlabel(1n *);ylabel(1 image(x4) *); title (fz4=2 + j*pi/6 111) subplot(5,2,9); stem(nz real(x5)); xlabel(* n *);ylabel(1 real(x5) 1); title ( f z5=j*pi/6 llt) subplot(5,2,10); stem(nz imag(x5)); xlabel(1n *);ylabel(1 image(x5) 1); title ( * z5=j *pi/6 时) 运行结果如下: 当Z二pi/6时，序列周期为12。 B.绘出信号x(?) =1.5sin(2^*0」町得频率就是多少？周期就是多少？产生一个数字频率 为0、9得正弦序列，并显示该信号,说明其周期。 当 x1=K 5*sin(2*pi*0, 1*n)时 程序如下： n=-20:20; x=l5*sin(2*pi*0l*n); stem(nz y) xlabel ( 1 时间!Z； 1 ); yla)0「幅度?)； title(*xl = l. 5*sin(2*pi*0. l*n) 1 ) 程序运行结果如下: x1=1.5*sin(2*pi*0.1*n) 时间序号 由图可知,xl得周期为10,频率为0、1、 当 x2=h 5*sin(0. 9*pi*0% 1*n)时 程序如下： n=-20:20; x2=l5*sin (09*n); stem(nz x2) xlabel (1时间序号