数字信号处理实验指导书15级.docxVIP

数字信号处理实验指导书实验一常见序列的产生及基本运算的MATLAB实现实验目的掌握离散信号的产生及基本运算的MATLAB实现方法。实验原理 利用MATLAB编写程序实现常用离散信号的产生及其基本运算。实验设备 电脑（windows xp 以上系统），MATLAB7.1以上版本软件。骤与内容实验步骤（一）、离散信号及其MATLAB实现单位冲激序列例1：用MATLAB产生64点的单位冲激序列clear all;N=64;x=zeros(1,N);x(1)=1;xn=0:N-1;stem(xn,x)axis([-1 65 0 1.1])任意序列例2：用MATLAB画出如下表达式的脉冲序列clear all;N=8;x=zeros(1,N);x(1)=8.0;x(2)=3.4x(3)=1.8;x(4)=5.6;x(5)=2.9;x(6)=0.7;xn=0:N-1;stem(xn,x)axis([-1 8 0 8.2])单位阶跃序列例3：用MATLAB实现单位阶跃函数clear all;N=32;x=ones(1,N);xn=0:N-1;stem(xn,x)axis([-1 32 0 1.1])斜坡序列例4：用MATLAB实现g(n)=3(n-4)点数为32的斜坡序列clear all;N=32;k=4B=3;t0=1;x=[zeros(1,k) ones(1,N-k)];for i=1:N x(i)=B*x(i)*(i-k);endxn=0:N-1;stem(xn,x)axis([-1 32 0 90])正弦序列例5：用MATLAB实现幅度A=3，频率f=100，初始相位Φ=1.2，点数为32的正弦信号clear all;N=32;A=3;f=100;fai=1.2;xn=0:N-1;x=A*sin(2*pi*f*(xn/N)+fai);stem(xn,x)axis([-1 32 -3.2 3.2])实指数序列例6：用MATLAB实现，点数为32的实指数序列clear all;N=32;A=3;a=0.7;xn=0:N-1;x=A*a.^xn;stem(xn,x)复指数序列例7：用MATLAB实现幅度A=3，a=0.7，角频率ω=314，点数为32的实指数序列clear all;N=32;A=3;a=0.7;w=314;xn=0:N-1;x=A*exp((a+j*w)*xn);stem(xn,x)随机序列利用MATLAB产生两种随机信号：rand(1,N)在区间上产生N点均匀分布的随机序列randn(1,N)产生均值为0，方差为1的高斯随机序列，即白噪声序列例8：用MATLAB产生点数为32的均匀分布的随机序列与高斯随机序列clear all;N=32;x_rand=rand(1,N);x_randn=randn(1,N);xn=0:N-1;figure(1);stem(xn,x_rand)figure(2);stem(xn,x_randn)（二）离散信号的基本运算信号的延迟给定离散信号x(n)，若信号y(n)定义为：y(n)=x(n-k)，那么y(n)是信号x(n)在时间轴上右移k个抽样周期得到的新序列。例1：正弦序列y(n)=sin(100n)右移3个抽样周期后所得的序列，MATLAB程序如下：clear all;N=32;w=100;k=3;x1=zeros(1,k);xn=0:N-1;x2=sin(100*xn);figure(1)stem(xn,x2)x=[x1 x2];axis([-1 N -1.1 1.1])N=N+k;xn=0:N-1;figure(2)stem(xn,x)axis([-1 N -1.1 1.1])信号相加若信号，值得注意的是当序列和的长度不相等或者位置不对应时，首先应该使两者的位置对齐，然后通过zeros函数左右补零使其长度相等后再相加例2：用MATLAB实现两序列相加clear all;n1=0:3x1=[2 0.5 0.9 1];figure(1)stem(n1,x1)axis([-1 8 0 2.1] )n2=0:7x2=[ 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7];figure(2)stem(n2,x2)axis([-1 8 0 0.8] )n=0:7;x1=[x1 zeros(1,8-length(n1))];x2=[ zeros(1,8-length(n2)),x2];x=x1+x2;figure(3)stem(n,x)axis([-1 8 0 2.1])信号相乘信号序列和相乘所得信号的表达式为：这是样本与样本之间的点乘运算，在MATLAB中可采用“.*”来实现，但是在信号序列相乘之前，应对其做与相加运算一样的操作。例3：用MATLAB实现上例中两序列相