前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >【数字信号处理】相关函数 ( 自相关函数示例 )

【数字信号处理】相关函数 ( 自相关函数示例 )

作者头像
韩曙亮
发布2023-03-30 11:54:59
5270
发布2023-03-30 11:54:59
举报
文章被收录于专栏:韩曙亮的移动开发专栏

文章目录

一、自相关函数 示例


给定一个 " 周期函数 " :

x(n) = A \sin (\omega n)

其中

\omega = \cfrac{2\pi}{N}

, 求该 " 周期函数 " 的 " 自相关函数 "

r_x(m)

" 周期信号 " 的 自相关函数 公式 :

r_x(m) = \cfrac{1}{N}\sum_{n = 0}^{N-1}x^*(n)x(n+m)

参考 【数字信号处理】相关函数 ( 周期信号 | 周期信号的自相关函数 ) 博客 ;

该信号是 " 实信号 " , 不是 " 复信号 " , 不需要使用共轭

^*

;

r_x(m) = \cfrac{1}{N}\sum_{n = 0}^{N-1}x(n)x(n+m)

x(n) = A \sin (\omega n)

代入到上面的式子中 ;

r_x(m) = \cfrac{1}{N}\sum_{n = 0}^{N-1} [ A \sin (\omega n) ] [ A \sin (\omega ( n + m )) ]

展开式子 , 计算得到 :

r_x(m) = \cfrac{1}{N}\sum_{n = 0}^{N-1} A^2 \sin (\omega n) \sin ( \omega n + \omega m )

使用 三角函数 和差化积 公式 , 参考 百度百科 https://baike.baidu.com/item/和差化积/6973039 ;

r_x(m) = \cfrac{A^2}{N} \cos \omega m \sum_{n = 0}^{N-1} \sin^2 \omega n + \cfrac{A^2}{N} \sin \omega m \sum_{n = 0}^{N-1} \sin \omega n \cos \omega n

下面的式子

\sum_{n = 0}^{N-1} \sin \omega n \cos \omega n = 0

值为

0

,

n = 0

时 ,

\sin \omega n \cos \omega n = 0

;

n = 1

时 , 与

n = N-1

时 , 抵消了 ;

n = 2

时 , 与

n = N-2

时 , 抵消了 ;

则最终结果为 0 , 则有 :

\cfrac{A^2}{N} \sin \omega m \sum_{n = 0}^{N-1} \sin \omega n \cos \omega n = 0

当前的推导相关函数为 :

r_x(m) = \cfrac{A^2}{N} \cos \omega m \sum_{n = 0}^{N-1} \sin^2 \omega n

根据 三角函数公式 :

\sin^2 \alpha=\cfrac{(1-\cos2\alpha)}{2}

可得 :

\sin^2 \omega n = \cfrac{(1- \cos 2 \omega n)}{2}

带入到相关函数中 , 可得 :

r_x(m) = \cfrac{A^2}{N} \cos \omega m \sum_{n = 0}^{N-1} \cfrac{1}{2} (1 - \cos 2 \omega n)

下面的式子

\sum_{n = 0}^{N-1} \cos 2 \omega n = 0

值为

0

,

则最终结果为 :

r_x(m) = \cfrac{A^2}{2} \cos \omega m
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2022-03-01,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 文章目录
  • 一、自相关函数 示例
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档