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

实验1 常见离散信号产生和实现 一、实验目的 1、加深对常用离散信号的理解； 2、熟悉使用MATLAB在时域中产生一些基本的离散时间信号。 二、实验原理 1、单位抽样序列 在MATLAB中可以利用函数实现。 2、单位阶越序列 在MATLAB中可以利用函数实现: 3、正弦序列 在MATLAB中实现过程如下： 4、复指数序列 在MATLAB中实现过程如下： 5、指数序列 在MATLAB中实现过程如下： 三、预习要求 1、预先阅读实验讲义（MATLAB基础介绍）； 2、讨论正弦序列、复指数序列的性质。 A．绘出信号，当、时、、时的信号实部和虚部图；当时呢？此时信号周期为多少？ 程序如下： titlez1=-1/12+j*pi/6; titlez2=1/12+j*pi/6; z3=1/12; z4=2+j*pi/6; z5=j*pi/6; n=0:20; x1=exp(titlez1*n); x2=exp(titlez2*n); x3=exp(z3*n); x4=exp(z4*n); x5=exp(z5*n); subplot(5,2,1); stem(n,real(x1)); xlabel(n);ylabel(real(x1)); title(z1=-1/12+j*pi/6时) subplot(5,2,2); stem(n,imag(x1)); xlabel(n);ylabel(imag(x1)); title(z1=-1/12+j*pi/6时) subplot(5,2,3); stem(n,real(x2)); xlabel(n);ylabel(real(x2)); title(z2=1/12+j*pi/6时) subplot(5,2,4); stem(n,imag(x2)); xlabel(n);ylabel(image(x2)); title(z2=1/12+j*pi/6时) subplot(5,2,5); stem(n,real(x3)); xlabel(n);ylabel(real(x3)); title(z3=1/12时) subplot(5,2,6); stem(n,imag(x3)); xlabel(n);ylabel(image(x3)); title(z3=1/12时) subplot(5,2,7); stem(n,real(x4)); xlabel(n);ylabel(real(x4)); title(z4=2+j*pi/6时) subplot(5,2,8); stem(n,imag(x4)); xlabel(n);ylabel(image(x4)); title(z4=2+j*pi/6时) subplot(5,2,9); stem(n,real(x5)); xlabel(n);ylabel(real(x5)); title(z5=j*pi/6时) subplot(5,2,10); stem(n,imag(x5)); xlabel(n);ylabel(image(x5)); title(z5=j*pi/6时) 运行结果如下： 当Z=pi/6时，序列周期为12。 B．绘出信号的频率是多少？周期是多少？产生一个数字频率为0.9的正弦序列，并显示该信号，说明其周期。 当x1=1.5*sin(2*pi*0.1*n)时 程序如下： n=-20:20； x=1.5*sin(2*pi*0.1*n); stem(n,y) xlabel(时间序号); ylabel(幅度); title(x1=1.5*sin(2*pi*0.1*n)) 程序运行结果如下： 由图可知，x1的周期为10，频率为0.1. 当x2=1.5*sin(0.9*pi*0.1*n)时 程序如下： n=-20:20; x2=1.5*sin(0.9*n); stem(n,x2) xlabel(时间序号); ylabel(幅度); title(x2=1.5*sin(0.9*pi*0.1*n)) 程序运行结果如下： 由图可知，x2为非周期函数，由理论分析可得，0.9不是π的倍数，所以不是周期函数。 3、使用帮助功能学习square(方波),sawtooth(锯齿波)和sinc函数，并绘图。 程序如下： （1）方波程序： %%%%%%%%%%%%%%%%%%%% n=-2*pi:0.0001:2*pi y=square(n); plot(n,y) xlabel(时间序号); ylabel(幅度); title(square) 方波结果： %%%%%%%%%%%%%%%%%%% （2）锯齿波程序： n=-10:0.0001:10; y=sawtooth(n); plot(n,y) xlabel(时间序号); ylabel(幅度); t