离散信号与系统分析(上)

瘦欲@ 提交于 2020-02-22 03:22:31

一、利用MATLAB产生离散信号

1.前言部分

stem(X,Y):在X的指定点处画出数据序列Y;

stem(X,Y,'filled'):以实心的方式画出茎秆;

axis([xmin,xmax,ymin,ymax]):设置当前图形的坐标范围,分别为x轴的最小、最大值,y轴的最小、最大值;

2.实验部分

茎状图

k=-2:1:4;
x=[2,1,1,-1,3,0,2];
stem(k,x);
axis([-3,5,-1.5,3.5])

实心茎状图

k=-2:1:4;
x=[2,1,1,-1,3,0,2];
stem(k,x,'filled');
axis([-3,5,-1.1,3.1])

实验结果图

3.结果分析部分

matlab只能表示一个有限长度的序列,而不能表示无穷大的序列。

二、利用MATLAB计算系列卷积与相关

1.前言部分

卷积函数:conv
w = conv(u,v): 返回向量 u 和 v 的卷积。如果 u 和 v 是多项式系数的向量,对其卷积与将这两个多项式相乘等效。

相关函数:xcorr
c = xcorr(x,y):返回矢量长度为2*N-1互相关函数序列,其中x和y的矢量长度均为N,如果x和y的长度不一样,则在短的序列后补零直到两者长度相等。

2.实验部分

x=[1,2,3,4];
h=[2,1,3];
y1=conv(x,h)
y2=xcorr(x,h)
y3=conv(h,x)
y4=xcorr(h,x)
y5=xcorr(x,x)
y6=xcorr(h,h)

输出结果

y1 =

 2     5    11    17    13    12

y2 =

-0.0000 3.0000 7.0000 13.0000 19.0000 10.0000 8.0000

y3 =

 2     5    11    17    13    12

y4 =

8.0000   10.0000   19.0000   13.0000    7.0000    3.0000   -0.0000

y5 =

4.0000   11.0000   20.0000   30.0000   20.0000   11.0000    4.0000

y6 =

 6     5    14     5     6

3.结果分析部分

从数字信号处理的角度方面来看,自相关运算可以用卷积运算来代替;在此我就不摆复杂公式了,简单的列举几个结论;

自相关函数:r[-n]=r[n] 偶对称序列,关于x=0对称;可以用xcorr[-n]=xcorr[n]表示;如实验结果y5,y6;

r[n]在n=0处的数值最大;如实验结果y5,y6中30.0000、14分别为其序列中的最大值;

互相关函数xcorr[X,Y]=-xcorr[Y,X],可见xcorr[X,Y]与xcorr[Y,X]互为其翻转序列;如实验结果y4,y6;

易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!