发明名称 基于变长-迭代运算的单音信号频率估计方法
摘要 本发明涉及单音信号频率估计技术,特别是基于变长-迭代运算的单音信号频率估计方法,对于待估计单音信号的正弦信号<img file="988514dest_path_image001.GIF" wi="105" he="25" />,其中<img file="768251dest_path_image002.GIF" wi="21" he="25" />是真实频率,数据长度为<img file="101143dest_path_image003.GIF" wi="20" he="18" />,采样率为<img file="915515dest_path_image004.GIF" wi="22" he="18" />,计算分辨率为<img file="675661dest_path_image005.GIF" wi="60" he="42" />,对该真实频率<img file="942694dest_path_image002.GIF" wi="21" he="25" />利用FFT定义得到估计频率<img file="79278dest_path_image006.GIF" wi="21" he="25" />,并取估计频率<img file="748156dest_path_image006.GIF" wi="21" he="25" />左右两边各取一个频率,对比估计频率<img file="177738dest_path_image006.GIF" wi="21" he="25" />和左右两边频率的幅度谱大小,从而进行迭代算法,最终达到三个频率中,使得估计频率<img file="932068dest_path_image006.GIF" wi="21" he="25" />的幅度谱最大,从而得到频率估计值:<img file="872342dest_path_image007.GIF" wi="114" he="46" />;本发明原理简单、精度高、运算速度快、易于工程实现,通过不断改变长度,利用迭代运算,逐步逼近频率真实值,寻找最优的频率估计值,广泛适用于接收端载波估计及信号分析等众多领域;并且估计精度可以通过改变迭代步进调整,方便灵活,适合于硬件平台处理。
申请公布号 CN102420792A 申请公布日期 2012.04.18
申请号 CN201110380709.5 申请日期 2011.11.25
申请人 中国工程物理研究院电子工程研究所 发明人 权友波;王甲峰;岳旸;李兵;任亚博;李蕾
分类号 H04L27/00(2006.01)I 主分类号 H04L27/00(2006.01)I
代理机构 成都天嘉专利事务所(普通合伙) 51211 代理人 方强
主权项 1.基于变长-迭代运算的单音信号频率估计方法,对于待估计单音信号的正弦信号<img file="545629DEST_PATH_IMAGE001.GIF" wi="104" he="25" />,其中<img file="647578DEST_PATH_IMAGE002.GIF" wi="21" he="25" />是真实频率,数据长度为<img file="826886DEST_PATH_IMAGE003.GIF" wi="20" he="18" />,采样率为<img file="367589DEST_PATH_IMAGE004.GIF" wi="22" he="18" />,计算分辨率为<img file="181699DEST_PATH_IMAGE005.GIF" wi="60" he="42" />,其特征在于:首先对所述的待估计单音信号的正弦信号<img file="782445DEST_PATH_IMAGE006.GIF" wi="33" he="24" />做<img file="183470DEST_PATH_IMAGE003.GIF" wi="20" he="18" />点傅里叶变换,<img file="527864DEST_PATH_IMAGE007.GIF" wi="178" he="25" />,得到幅度谱最大值<img file="697945DEST_PATH_IMAGE008.GIF" wi="152" he="53" />,<img file="469592DEST_PATH_IMAGE009.GIF" wi="24" he="25" />对应的序号<img file="357914DEST_PATH_IMAGE010.GIF" wi="18" he="25" />,估计频率为<img file="240419DEST_PATH_IMAGE011.GIF" wi="105" he="42" />;其中,<img file="592903DEST_PATH_IMAGE010.GIF" wi="18" he="25" />为幅度谱最大值<img file="971669DEST_PATH_IMAGE009.GIF" wi="24" he="25" />对应序号,也为频率估计值对应的序号,真实频率<img file="409604DEST_PATH_IMAGE012.GIF" wi="178" he="46" />;这时,真实频率<img file="33483DEST_PATH_IMAGE002.GIF" wi="21" he="25" />对应的序号<img file="240474DEST_PATH_IMAGE013.GIF" wi="17" he="25" />位于<img file="229289DEST_PATH_IMAGE014.GIF" wi="57" he="25" />之间或者位于<img file="154520DEST_PATH_IMAGE015.GIF" wi="58" he="25" />之间,其中<img file="815046DEST_PATH_IMAGE016.GIF" wi="20" he="24" />是幅度谱中最左端频率对应的序号,<img file="142122DEST_PATH_IMAGE017.GIF" wi="21" he="24" />是幅度谱中最左端频率对应的序号;当<img file="160894DEST_PATH_IMAGE013.GIF" wi="17" he="25" />位于<img file="511104DEST_PATH_IMAGE018.GIF" wi="58" he="25" />之间时,迭代步进为<img file="539103DEST_PATH_IMAGE019.GIF" wi="82" he="26" />;由<img file="658369DEST_PATH_IMAGE020.GIF" wi="105" he="42" />处幅度谱<img file="113621DEST_PATH_IMAGE021.GIF" wi="137" he="54" />,可以得到估计频率<img file="685547DEST_PATH_IMAGE022.GIF" wi="21" he="25" />的左边某一处频率<img file="251658DEST_PATH_IMAGE023.GIF" wi="20" he="24" />和<img file="983685DEST_PATH_IMAGE023.GIF" wi="20" he="24" />处的幅度谱<img file="609839DEST_PATH_IMAGE024.GIF" wi="22" he="24" />,以及估计频率<img file="731378DEST_PATH_IMAGE022.GIF" wi="21" he="25" />的右边某一处频率<img file="38863DEST_PATH_IMAGE025.GIF" wi="22" he="24" />和<img file="195038DEST_PATH_IMAGE025.GIF" wi="22" he="24" />处的幅度谱<img file="929776DEST_PATH_IMAGE026.GIF" wi="25" he="24" />;那么由频率<img file="538611DEST_PATH_IMAGE027.GIF" wi="110" he="25" />之间的关系:<img file="384208DEST_PATH_IMAGE028.GIF" wi="154" he="25" />,可得频率<img file="394889DEST_PATH_IMAGE027.GIF" wi="110" he="25" />分别对应的序号<img file="799063DEST_PATH_IMAGE029.GIF" wi="101" he="25" />之间的关系:<img file="895195DEST_PATH_IMAGE030.GIF" wi="145" he="25" />;然后判断<img file="544482DEST_PATH_IMAGE031.GIF" wi="85" he="25" />之间的关系,当满足<img file="409670DEST_PATH_IMAGE032.GIF" wi="86" he="25" />单调时,即<img file="548528DEST_PATH_IMAGE009.GIF" wi="24" he="25" />不是最大时,将<img file="69639DEST_PATH_IMAGE033.GIF" wi="177" he="26" />,再次得到新的<img file="584934DEST_PATH_IMAGE031.GIF" wi="85" he="25" />,判断新的<img file="242311DEST_PATH_IMAGE031.GIF" wi="85" he="25" />三者幅度谱的大小,继续迭代运算,直到某次迭代后<img file="552070DEST_PATH_IMAGE031.GIF" wi="85" he="25" />关系满足<img file="59012DEST_PATH_IMAGE034.GIF" wi="125" he="25" />,则将长度<img file="112419DEST_PATH_IMAGE003.GIF" wi="20" he="18" />记为<img file="889882DEST_PATH_IMAGE035.GIF" wi="34" he="25" />,得到频率估计值:<img file="370542DEST_PATH_IMAGE036.GIF" wi="114" he="46" />当<img file="866245DEST_PATH_IMAGE013.GIF" wi="17" he="25" />位于<img file="457764DEST_PATH_IMAGE015.GIF" wi="58" he="25" />之间,采用同样的方式进行判断得到频率估计值;当进行迭代运算时,<img file="152050DEST_PATH_IMAGE037.GIF" wi="114" he="26" />长度变化,分辨率<img file="6874DEST_PATH_IMAGE038.GIF" wi="60" he="42" />随之变化。
地址 621900 四川省绵阳市游仙区919信箱523分箱