数字信号处理实验1.doc

《数字信号处理实验1.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、桂林理工大学实验报告1.实验目的（1）掌握信号表示方法。（2）掌握离散信号的延迟、相加、相乘及平移、反折变换2.实验原理与方法首先，在用MATLAB表示离散序列并将其可视化时，我们还应注意以下几点：第一、离散时间序列无法用符号运算来表示，要用矩阵的形式；第二、由于在MATLAB中矩阵的元素个数有限，因此，MATLAB无法表示无限长的序列；第三、在绘制离散信号波形的函数stem命令，而不是plot命令。习题1-2计算卷积y[k]=x[k]*g[k](2)x[k=g[k]={k+1；0<=k<=3}实验程序为x=[

2、1,2,3,4];kx=0:3;g=[1,2,3,4];kg=0:3;y1=conv(x,g);xf=fliplr(x)y2=conv(xf,g);k1=kx(1)+kg(1):kx(end)+kg(end);stem(k1,y1);xlabel('k1');ylabel('y1');习题1-3（1）已知序列想x[k]={1,2,3,4;k=0,1,2,3},y[k]={-1,1,-2,3;k=0,1,2,3},试计算序列x[k]的自相关函数rx[n],以及序列x[k]与y[k]的互相关函数rxy[k]和ryx

3、[k]。（2）若某两个序列的自相关函数rx[k]相同，则该两个序列之间存在何种关系？（3）若信号s[k]=x[k]+n[k],x[k]与n[k]相关性很小，如何判断s[k]中存在x[k]?x=[1,2,3,4];kx=0:3;y=[-1,1,-2,3];ky=0:3;xf=fliplr(x);rx=conv(x,x+1);rxy=conv(xf,y);ryx=conv(y,xf);k1=kx(1)+ky(1):kx(end)+ky(end);kxf=-fliplr(kx);k2=kxf(1)+ky(1):kxf

4、(end)+ky(end);subplot(1,3,1)stem(k1,rx);xlabel('k1');ylabel('rx');subplot(1,3,2)stem(k2,rxy);xlabel('k2');ylabel('rxy');subplot(1,3,3);stem(k2,ryx);xlabel('k2');ylabel('ryx');（2）因为rx[n]=x[n]*x[-n]，所以这两个序列相等或互为翻转序列。（3）两个相互无关的函数（即对于所有k，两函数的互相关均为0）之和的自相关函数等于各自自

5、相关函数之和。因为x[k]与n[k]相关性很小，可看做两个相互无关的函数。所以，若rs[n]=rx[n]+rn[n],则s[k]中存在x[k]。M1-1已知g1(t)=cos(6*pi*t),g2(t)=cos(14*pi*t),g3=(26*pi*t),以抽样频率fasm=10Hz对上述三个信号进行抽样。在同一张图上画出g1(t),g2(t),g3(t)及其抽样点，对所得结果进行讨论。程序为fs=10;t=0:0.01:2;g1=cos(6*pi*t);g2=cos(14*pi*t);g3=cos(26*pi

6、*t);T=1/10;k=0:T:2;g11=cos(6*pi*k);g22=cos(14*pi*k);g33=cos(26*pi*k);subplot(3,2,1);plot(t,g1);title('g1')subplot(3,2,2);stem(k,g11);title('g11');subplot(3,2,3);plot(t,g2);title('g2');subplot(3,2,4);stem(k,g22);title('g22');subplot(3,2,5);plot(t,g3);title('g

7、3');subplot(3,2,6);stem(k,g33);title('g33');