发明名称 基于DFT插值的色载波起始相位恢复方法
摘要 基于DFT插值的色载波起始相位恢复方法,(1)、提取没有波形畸变的色同步信号共64个采样点;(2)、使用DFT变换计算第6和第7个离散频点的DFT系数;(3)、使用线性插值计算起始相位:上述频谱具有线性相位,色同步采样点个数应为2的整数次方,64为最优值;计算两个DFT系数,直接采用DFT变换。本发明应用于对视频信号进行数字彩色解码的情况,无需使用色载波恢复锁相环路或者矩阵运算,使用乘累加器计算所需的两个DFT系数,接着通过线性插值恢复出本地色载波信号的起始相位,以进行高精度的色调解调,具有结构简单,实现方便,运算快速,适用面广等优点。
申请公布号 CN101370147B 申请公布日期 2011.02.09
申请号 CN200810156554.5 申请日期 2008.09.28
申请人 南京大学 发明人 沈庆宏;肖波;都思丹;韩涛;黄勇才
分类号 H04N9/45(2006.01)I 主分类号 H04N9/45(2006.01)I
代理机构 南京天翼专利代理有限责任公司 32112 代理人 汤志武;王鹏翔
主权项 基于DFT插值的色载波起始相位恢复方法,其特征是:(1)、提取没有波形畸变的色同步信号共64个采样点;(2)、使用DFT变换计算第6和第7个离散频点的DFT系数;设经采样后的色同步信号的表达式为x(n)=asin(2π(fc+Δfc)n/fs+θ0)        n=1~N其中a,fc,Δfc,fs,θ0,N,n分别为色同步信号的幅度,标称频率,偏移频率,采样频率,初始相角,采样点数,采样点的序数;为了得到色载波的起始相位,对其作傅立叶变换,根据频谱即可得到相关信息,但在进行离散频谱分析时存在两个主要误差;x(n)的DFT变换X(eiω): <mrow> <mi>X</mi> <mrow> <mo>(</mo> <msup> <mi>e</mi> <mi>j&omega;</mi> </msup> <mo>)</mo> </mrow> <msub> <mo>|</mo> <mrow> <mi>&omega;</mi> <mo>=</mo> <mi>k</mi> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> <mi>N</mi> </mfrac> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mi>N</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>&omega;</mi> <mi>c</mi> </msub> <mo>+</mo> <mi>&Delta;</mi> <msub> <mi>&omega;</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>N</mi> <mi>sin</mi> <mrow> <mo>(</mo> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>&omega;</mi> <mi>c</mi> </msub> <mo>+</mo> <mi>&Delta;</mi> <msub> <mi>&omega;</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mn>0</mn> </msub> <mo>-</mo> <mfrac> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> <mn>2</mn> </mfrac> <mi>&omega;</mi> <mo>)</mo> </mrow> </mrow> </msup> </mrow>(3)、使用线性插值计算起始相位:上述频谱具有线性相位,色同步信号采样点个数应为64;当采样率为50MHz,频谱分辨率为0.78125MHz,而色载波频率为4.43361875MH,所以频谱的真实谱峰在第6个和第7个离散频点之间,并更接近第7个离散频点,设其与第6个频点的归一化间隔为δ,显然0<δ<1;得第6和第7个谱点的幅度为: <mrow> <mo>|</mo> <mi>X</mi> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> <mo>|</mo> <mo>=</mo> <mfrac> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&pi;&delta;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>64</mn> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&pi;&delta;</mi> <mo>/</mo> <mn>64</mn> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> <mrow> <mo>|</mo> <mi>X</mi> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> <mo>|</mo> <mo>=</mo> <mfrac> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&pi;&delta;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>64</mn> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&pi;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&delta;</mi> <mo>)</mo> </mrow> <mo>/</mo> <mn>64</mn> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>对分母进行泰勒展开得到近似表达式: <mrow> <mo>|</mo> <mi>X</mi> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> <mo>|</mo> <mo>&ap;</mo> <mfrac> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&pi;&delta;</mi> <mo>)</mo> </mrow> </mrow> <mi>&pi;&delta;</mi> </mfrac> </mrow>和 <mrow> <mo>|</mo> <mi>X</mi> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> <mo>|</mo> <mo>&ap;</mo> <mfrac> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&pi;&delta;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&pi;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&delta;</mi> <mo>)</mo> </mrow> <mi></mi> </mrow> </mfrac> </mrow>两者相除得: <mrow> <mi>&delta;</mi> <mo>=</mo> <mfrac> <mi>&beta;</mi> <mrow> <mn>1</mn> <mo>+</mo> <mi>&beta;</mi> </mrow> </mfrac> </mrow>其中 <mrow> <mi>&beta;</mi> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <mi>X</mi> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <mi>X</mi> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> <mo>|</mo> </mrow> </mfrac> </mrow>于是线性插值得到的起始相位估计值为: <mrow> <msub> <mover> <mi>&theta;</mi> <mo>^</mo> </mover> <mn>0</mn> </msub> <mo>=</mo> <mi>angle</mi> <mrow> <mo>(</mo> <mi>X</mi> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&delta;</mi> <mo>)</mo> </mrow> <mo>*</mo> <mi>&pi;</mi> <mo>*</mo> <mn>63</mn> <mo>/</mo> <mn>64</mn> <mo>-</mo> <mi>&pi;</mi> <mo>/</mo> <mn>2</mn> </mrow>只需得到指定频率位置的相位,故只计算两个DFT系数,直接采用DFT变换。
地址 210093 江苏省南京市汉口路22号