两个首尾对称序列的自相关

设确定性离散序列x1(n)、x2(n)的长度均为N,若序列x1(n)=[ x(0), x(1), …, x(N-2), x(N-1) ],序列x2(n)=[ x(N-1), x(N-2), …, x(1), x(0) ],则两个序列的自相关相等。

  1. 结论证明
  • 序列x1(n)的自相关
    在这里插入图片描述
  • 序列x2(n)的自相关
    在这里插入图片描述
  • n=0时
    Rx1(n)=Rx2(n)=x(0)⋅x(0)+x(1)⋅x(1)+...+x(N−2)⋅x(N−2)+x(N−1)⋅x(N−1)R_{x1}(n)=R_{x2}(n)=x(0)\cdot x(0)+x(1)\cdot x(1)+...+x(N-2)\cdot x(N-2)+x(N-1)\cdot x(N-1)Rx1(n)=Rx2(n)=x(0)x(0)+x(1)x(1)+...+x(N2)x(N2)+x(N1)x(N1)
  • n=1时
    Rx(n)=x(0)⋅x(1)+x(1)⋅x(2)+...+x(N−2)⋅x(N−1)R_{x}(n)=x(0)\cdot x(1)+x(1)\cdot x(2)+...+x(N-2)\cdot x(N-1)Rx(n)=x(0)x(1)+x(1)x(2)+...+x(N2)x(N1)
    Rx2(n)=x(N−1)⋅x(N−2)+x(N−2)⋅x(N−3)+...+x(1)⋅x(0)R_{x2}(n)=x(N-1)\cdot x(N-2)+x(N-2)\cdot x(N-3)+...+x(1)\cdot x(0)Rx2(n)=x(N1)x(N2)+x(N2)x(N3)+...+x(1)x(0)
    所以,n=1时,Rx1(1)=Rx2(1)R_{x1}(1)=R_{x2}(1)Rx1(1)=Rx2(1)
  • n=2时
    Rx1(n)=x(0)⋅x(2)+x(1)⋅x(3)+...+x(N−3)⋅x(N−1)R_{x1}(n)=x(0)\cdot x(2)+x(1)\cdot x(3)+...+x(N-3)\cdot x(N-1)Rx1(n)=x(0)x(2)+x(1)x(3)+...+x(N3)x(N1)
    Rx2(n)=x(N−1)⋅x(N−3)+...+x(3)⋅x(1)+x(2)⋅x(0)R_{x2}(n)=x(N-1)\cdot x(N-3)+...+x(3)\cdot x(1)+x(2)\cdot x(0)Rx2(n)=x(N1)x(N3)+...+x(3)x(1)+x(2)x(0)
    所以,n=2时,Rx1(2)=Rx2(2)R_{x1}(2)=R_{x2}(2)Rx1(2)=Rx2(2)
  • n=i ( i=0, 1, …, N-1 )时
  • n=N-1时
    Rx1(N−1)=Rx2(N−1)=x(0)⋅x(N−1)R_{x1}(N-1)=R_{x2}(N-1)=x(0)\cdot x(N-1)Rx1(N1)=Rx2(N1)=x(0)x(N1)
    综上所述,Rx1(n)=Rx2(n)R_{x1}(n)=R_{x2}(n)Rx1(n)=Rx2(n),即两个序列的自相关相等。
  1. 程序代码
clc
clear all
close all

N=3;
x=[1 2 3];
y=[3 2 1];
rxx=xcorr(x);              %x的自相关序列
ryy=xcorr(y);              %y的自相关序列
Rxx=zeros(N);
Ryy=zeros(N);
for i=1:N                  %x的自相关矩阵
    for j=1:N
        Rxx(i,j)=rxx(N+i-j);
    end
end
for i=1:N                  %y的自相关矩阵
    for j=1:N
        Ryy(i,j)=ryy(N+i-j);
    end
end
Rxx
Ryy
  1. 仿真结果
    在这里插入图片描述
    仿真结果与理论证明一致,验证了该理论的正确性。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值