如何计算由差分方程定义的信号的自相关?

信息处理 离散信号 自相关
2022-02-18 04:38:00

我没有差分方程的经验,我想学习如何计算以下内容,但我在网上找不到资源。任何帮助将不胜感激。

寻找:

E[d[n]d[n+k]]

在哪里:

  • d[n]是形式为的离散时间信号:

d[n]=αd[n1]+v[n](α[0;1])

  • v[n]是一系列不相关的零均值单位方差随机变量。
2个回答

有不同的方法来计算请求的自相关,但我想提供最简单的方法来满足你的情况。

  • 差分方程是表示 LTI 系统的 LCCDE 类型。
  • 系统的输入v[n]是 WSS(广义平稳)随机过程。

我假设在您的问题中都可以观察到这两者。然后,我将使用由(实际)WSS 输入驱动的(实际)LTI 系统的输入和输出自相关之间的以下众所周知的关系v[n]

rdd[k]=h[k]h[k]rvv[k]

WSS的自相关v[n]rvv[k]=σv2δ[k]=δ[k]对于不相关的、零均值的单位方差过程。h[n]是由 LCCDE 表示的 LTI 系统的脉冲响应。

然后是 WSS 输出的自相关序列d[n]

rdd[k]=h[k]h[k]=m=h[m]h[(km)]   ,    for k=0,±1,±2,...

所以假设你的 LCCDE :

d[n]α d[n1]=v[n]   ,    with |α|<1

表示因果 LTI 系统,则其脉冲响应为:

h[n]=(α)nu[n]

并产生自相关为:

rdd[k]=m=(α)mu[m](α)mku[mk]   ,    for k=0,±1,±2,...

我希望你能继续剩下的。

我会通过尝试表达来解决这个问题d(n)不同,因为差分方程可能会令人困惑。请注意:

d(n)=αd(n1)+v(n)=α2d(n2)+αv(n1)+v(n)=α3d(n3)+α2v(n2)+αv(n1)+v(n)=...

你可以从你可以表达的模式中看到d(n)作为总结:

d(n)=αnd(0)+i=0n1αiv(ni)

我假设过程开始于n=0. 为简单起见,让我采取d(0)=0这样我们就不必在整个计算过程中携带该常数。

然后我们剩下:

E[d(n)d(n+k)]=E[(i=0n1αiv(ni))(i=0n+k1αiv(n+ki))]

知道噪声与零均值不相关,这个表达式可以真正简化。

请注意,表达式取决于n因为这个过程不是 WSS(因为我假设它开始于n=0,即它从一开始就没有运行过)。