发明名称 基于时间延迟的离散频谱低频成分的校正方法
摘要 本发明涉及一种基于时间延迟的离散频谱低频成分的校正方法,其特征在于包括以下步骤:(1).将离散时间信号按照一定时间延迟分成等长度的三段,对每段时间信号分别求和,并获得三个求和值;(2).利用离散时间信号频谱的对称性,并利用傅立叶变换公式获得离散频谱低频成分的有关频谱校正的表达形式;(3).将三个求和值结合离散频谱低频成分的有关频谱校正的表达形式获得离散频谱中低频成分的频率和相位;根据所获得的频率和相位结合(2)频谱校正的表达形式并获得低频成分的幅值;(4).对相位进行修正,获得相位的校正值。其计算量小,计算精度高,简单易行,适用于任何窗函数,提高了离散频谱低频成分校正的效率和准确性。
申请公布号 CN102269803B 申请公布日期 2014.07.09
申请号 CN201010191774.9 申请日期 2010.06.04
申请人 北京化工大学 发明人 江志农;冯坤;秦强;何威;马波
分类号 G01R35/00(2006.01)I;G01R23/16(2006.01)I 主分类号 G01R35/00(2006.01)I
代理机构 北京思创毕升专利事务所 11218 代理人 刘明华
主权项 一种基于时间延迟的离散频谱低频成分的校正方法,其特征在于包括以下步骤:(1).将离散时间信号按照一定时间延迟分成等长度的三段,对每段时间信号分别求和,并获得三个求和值;(2).利用离散时间信号频谱的对称性,使用傅立叶变换公式进行频谱分析,以便事先得到离散频谱低频成分的有关频谱校正的表达形式;(3).将(1)所述的三个求和值结合(2)所述的离散频谱低频成分的有关频谱校正的表达形式中涉及低频成分的频率和相位的表达式,获得离散频谱中低频成分的频率和相位,再根据所得到频率和相位结合(2)所述的离散频谱低频成分的有关频谱校正的表达形式中涉及低频成分的幅度的表达形式,并获得低频成分的幅值;(4).根据(3)所述的频率、幅值和相位结合(2)所述的有关频谱校正的表达形式,对相位进行修正,获得相位的校正值;获得步骤(1)所述的三个求和值,包括以下步骤:(1.1)将待分析的离散时间信号设为x(t),采样长度为2N,采样时间为2T;(1.2)截取该信号[T/2,3T/2)区间上的部分作为x<sub>1</sub>(t),[0,T)区间上的部分作为x<sub>2</sub>(t),[T,2T)区间上的部分作为x<sub>3</sub>(t);(1.3)则每段信号时长为T,长度为N,相邻信号的时间延迟为T/2;(1.4)对每段时间信号分别求和,获得三个求和值:<img file="FDA0000486337700000011.GIF" wi="310" he="141" />其中X<sub>i</sub>表示编号为i的三段时间信号的求和值,i=1,2,3;获得步骤(2)所述的离散频谱低频成分的有关频谱校正的表达形式的方法如下:(2.1).设三段信号的傅立叶变换为X<sub>i</sub>(f),i=1,2,3,f表示频率;(2.2).利用离散时间信号频谱的对称性,使用傅立叶变换公式,可获得离散频谱低频成分的频谱校正的表达形式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>cos</mi><mrow><mo>(</mo><mi>&pi;</mi><msub><mi>Tf</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>X</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>X</mi><mn>3</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow><mrow><msub><mrow><mn>2</mn><mi>X</mi></mrow><mn>1</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000486337700000021.GIF" wi="552" he="142" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>tan</mi><mi>&theta;</mi><mo>=</mo><mfrac><mrow><msub><mi>X</mi><mn>3</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>&pi;</mi><msub><mi>Tf</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>X</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>X</mi><mn>3</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>&pi;</mi><msub><mi>Tf</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000486337700000022.GIF" wi="651" he="142" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfrac><mrow><msub><mi>X</mi><mn>3</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow><mrow><mi>W</mi><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000486337700000023.GIF" wi="370" he="142" /></maths>式中,f<sub>0</sub>表示离散时间信号中低频成分的频率校正值,X<sub>i</sub>(0)分别表示三段信号频谱0Hz处的谱线值,其中i=1,2,3,θ表示低频成分的相位,A表示低频成分的幅值校正值,W(f<sub>0</sub>)表示窗函数频谱在频率f<sub>0</sub>处的幅值;(2.3).利用离散傅立叶变换公式的性质,可获得X<sub>i</sub>(0)与X<sub>i</sub>之间的关系为:X<sub>i</sub>(0)=X<sub>i</sub>,其中i=1,2,3。
地址 100029 北京市朝阳区北三环东路15号