线性相关,卷积运算与toeplitz矩阵.docxVIP

线性相关,卷积运算与toeplitz矩阵（仅供学习参考）------------例题--------------设x=[1 2 3] , y=[4 5 6], 求r_yx,c_yx 用竖式法对应竖式法，构造toeplitz矩阵用matlab完成运算---------------求线性相关运算--------------原理：表达式为，简洁表达式为： ，表达式中，m表示位移量，m0表示x序列左移，m0表示x序列右移，不同的m得到不同的值，如， 的值大于0表示有同相成分的存在，小于0表示有反相成分的存在，等于0表示两序列正交或者相互独立。注意：线性相关运算的长度为length(x) + length(y) -1 ，与是不同的。竖式法：写开了就是：由此构造矩阵xM*y’: ,上面的矩阵xM就是一个toeplitz矩阵，故在matlab程序中可以用toeplitz矩阵乘法运算代表线性相关运算。注意：xM的行数是线性相关运算结果的长度3+3-1=5，列数对应y的行数。Matlab程序：%------------------------------------------x=[1,2,3]; x1=fliplr(x);y=[4,5,6];xM=toeplitz([x1 zeros(1,length(x)-1)],[x1(1) zeros(1,length(x)-1)]);r_yx=xM*y’;%---------------------------------------------------------------程序注释：1. fliplr使得x左右翻转变成[3,2,1],这样做的目的是为了构造下一步的toeplitz矩阵；2. 依照竖式法的形式构造toeplitz矩阵的时候，一般来说toeplitz矩阵的格式为—toeplitz（c,r）,C代表列，r代表行，以c为第一列，r为第一行就可以构造出一个toeplitz矩阵，如（toeplitz矩阵可以不是方阵），其中，。注意，当c与r的第一个元素不相同时，矩阵的第一个元素是c的第一个元素。3. toeplitz矩阵补零的问题，由于x，y的长度都是3，则二者线性相关运算之后的长度是7，由x1生成的toeplitz矩阵的行应该是7，故以x1中的元素作为toeplitz矩阵第一列的前3个元素，然后再在后面补上2个0，列数为了与y的行数保持一致，设定为3，以x1的第一个元素为toeplitz矩阵第一行的第一个元素，之后再补2个0。---------------求线性卷积运算--------------原理：表达式为，简洁表达式为： ，计算卷积时，首先反折得到，m0表示右移，m0表示序列左移，不同的m得到不同的值。线性卷积运算的长度为length(x) + length(y) -1 ，与是不同的，也就是说线性卷积满足交换律，这一点与线性相关运算不同。---------------比较线性相关与线性卷积的关系--------------线性相关：，，线性卷积：，，两种运算有如下关系：-----------------------------------------------------------------------竖式法：写开了就是：由此构造矩阵xM*y’: ,上面的矩阵xM就是一个toeplitz矩阵，故在matlab程序中可以用toeplitz矩阵乘法运算代表线性卷积运算。注意：xM的行数是线性卷积运算结果的长度3+3-1=5，列数对应y的行数。Matlab程序：%------------------------------------------x=[1,2,3]; y=[4,5,6];xM=toeplitz([x zeros(1,length(x)-1)],[x(1) zeros(1,length(x)-1)]);r_yx=xM*y’;%---------------------------------------------------------------程序注释：1. toeplitz矩阵补零的问题，由于x，y的长度都是3，则二者线性卷积运算之后的长度是7，由x1生成的toeplitz矩阵的行应该是7，故以x1中的元素作为toeplitz矩阵第一列的前3个元素，然后再在后面补上2个0，列数为了与y的行数保持一致，设定为3，以x1的第一个元素为toeplitz矩阵第一行的第一个元素，之后再补2个0。2.比较前面线性相关的matlab程序，两种运算通过矩阵运算来表示时，唯一不同之处就是x生成的toeplitz矩阵不同，线性相关运算首先需要将x左右反转，进而再构成toeplitz矩阵。