发明名称 一种模拟通道的自动校准方法及系统
摘要 本发明公开了一种模拟通道的自动校准方法,通过获取各个频率下调理电路的幅值线性关系,根据模拟信号的主频率,得到各个幅值线性关系的权值,采用加权最小二乘法得到调理电路的增益系数、零漂校准参数,发送到被测系统的处理器中,对模拟通道采集的信号进行补偿。同时,获取各个频率下调理电路的相差,依据权值,利用加权最小二乘法求得最优相位差,并转换为延迟时间,被测系统的处理器控制时,将延迟时间τ考虑到控制过程中。在本发明中,考虑到了信号主频率实际分布情况,采用了加权最小二乘法,使处于主频率附近的拟合曲线精度更高,更加客观地反映了调理电路的实际情况,从而提高了模拟信号主频率附近的校准精度。
申请公布号 CN102654771A 申请公布日期 2012.09.05
申请号 CN201210150919.X 申请日期 2012.05.16
申请人 电子科技大学 发明人 凡时财;邹见效;徐红兵;李星晔
分类号 G05B23/02(2006.01)I 主分类号 G05B23/02(2006.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 温利平
主权项 1.一种模拟通道的自动校准方法,其特征在于,包括以下步骤:(1)、幅值校准1.1)、获取调理电路的幅频数据输入幅值x=A<sub>1</sub>的模拟信号S<sub>1</sub>到调理电路,依次改变模拟信号S<sub>1</sub>频率ω<sub>1</sub>,ω<sub>2</sub>,...,ω<sub>i</sub>,...,ω<sub>m</sub>,分别测得调理电路输出端幅值y=A′<sub>11</sub>,A′<sub>12</sub>,...,A′<sub>1i</sub>,...,A′<sub>1m</sub>;输入幅值x=A<sub>2</sub>的模拟信号S<sub>2</sub>到调理电路,依次改变模拟信号S<sub>2</sub>频率ω<sub>1</sub>,ω<sub>2</sub>,...,ω<sub>i</sub>,...,ω<sub>m</sub>,分别测得调理电路输出端幅值y=A′<sub>21</sub>,A′<sub>22</sub>,...,A′<sub>2i</sub>,...,A′<sub>2m</sub>;其中,m是频率数,频率ω<sub>1</sub>,ω<sub>2</sub>,...,ω<sub>i</sub>,...,ω<sub>m</sub>覆盖了调理电路工作的频率范围;1.2)、确定调理电路在频率ω<sub>i</sub>,i=1,2,...,m下前后端幅值线性关系根据两点法,有:<maths num="0001"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>A</mi><mrow><mn>1</mn><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>a</mi><mi>i</mi></msub><msub><mi>A</mi><mn>1</mn></msub><mo>+</mo><msub><mi>b</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msubsup><mi>A</mi><mrow><mn>2</mn><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>a</mi><mi>i</mi></msub><msub><mi>A</mi><mn>2</mn></msub><mo>+</mo><msub><mi>b</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>则:<maths num="0002"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>a</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>A</mi><mrow><mn>1</mn><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mo>-</mo><msubsup><mi>A</mi><mrow><mn>2</mn><mi>i</mi></mrow><mo>&prime;</mo></msubsup></mrow><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>-</mo><msub><mi>A</mi><mn>2</mn></msub></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>b</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>A</mi><mrow><mn>1</mn><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mo>-</mo><mfrac><mrow><msubsup><mi>A</mi><mrow><mn>1</mn><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mo>-</mo><msubsup><mi>A</mi><mrow><mn>2</mn><mi>i</mi></mrow><mo>&prime;</mo></msubsup></mrow><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>-</mo><msub><mi>A</mi><mn>2</mn></msub></mrow></mfrac><msub><mi>A</mi><mn>1</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>从而得到调理电路在频率ω<sub>i</sub>下的前后端幅值线性关系:y=a<sub>i</sub>x+b<sub>i</sub>,这样得到对应于频率ω<sub>1</sub>,ω<sub>2</sub>,...,ω<sub>i</sub>,...,ω<sub>m</sub>下的m条前后端幅值线性关系:<maths num="0003"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>y</mi><mo>=</mo><msub><mi>a</mi><mn>1</mn></msub><mi>x</mi><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>&omega;</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mi>y</mi><mo>=</mo><msub><mi>a</mi><mn>2</mn></msub><mi>x</mi><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>&omega;</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo>=</mo><msub><mi>a</mi><mi>i</mi></msub><mi>x</mi><mo>+</mo><msub><mi>b</mi><mi>i</mi></msub><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>&omega;</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo>=</mo><msub><mi>a</mi><mi>m</mi></msub><mi>x</mi><mo>+</mo><msub><mi>b</mi><mi>m</mi></msub><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>&omega;</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>1.3)、拟合调理电路前后端最优幅值线性关系y=ax+b公式(3)每条直线对应的权值分别为<img file="FDA00001644638700014.GIF" wi="400" he="46" />并且<img file="FDA00001644638700015.GIF" wi="364" he="77" />ω'为模拟信号的主频率,则有:<img file="FDA00001644638700021.GIF" wi="986" he="448" />以输入幅值x<sub>1</sub>,x<sub>2</sub>,...,x<sub>j</sub>,...,x<sub>n</sub>在每条直线上取n个点,即y<sub>ij</sub>,其中i=1,2,...,m,表示对应频率ω<sub>i</sub>下直线,j=1,2,...,n,表示取点序号;这样,这些点与最优直线y=ax+b的加权离差平方和为:<img file="FDA00001644638700022.GIF" wi="1099" he="127" />将φ分别对a,b求偏导数<img file="FDA00001644638700023.GIF" wi="108" he="116" /><img file="FDA00001644638700024.GIF" wi="99" he="116" />并令其等于零,即:<img file="FDA00001644638700025.GIF" wi="1123" he="295" />通过(6)式可解得a,b分别为:<img file="FDA00001644638700026.GIF" wi="1136" he="290" /><img file="FDA00001644638700027.GIF" wi="1149" he="290" />a、b分别作为增益系数、零漂校准参数;将增益系数、零漂校准参数发送到被测系统的处理器,对模拟通道采集的数字信号进行补偿:将数字信号的幅值减去零漂b,然后除以增益系数a,得到幅值校准后的数字信号;(2)、相位校准输入模拟信号S,分别测量调理电路前后端相位θ与θ',可得相位差Δθ=θ′-θ,依次改变模拟信号的频率为ω<sub>1</sub>,ω<sub>2</sub>,...,ω<sub>i</sub>,...,ω<sub>m</sub>可分别得到相位差Δθ<sub>1</sub>,Δθ<sub>2</sub>,...,Δθ<sub>i</sub>,...,Δθ<sub>m</sub>;利用加权最小二乘法求得最优相位差Δθ:加权离差平方和为:<img file="FDA00001644638700031.GIF" wi="1039" he="122" />将φ<sub>Δθ</sub>对Δθ求偏导数<img file="FDA00001644638700032.GIF" wi="136" he="117" /><img file="FDA00001644638700033.GIF" wi="1070" he="125" />可解得最优相位差Δθ为:<img file="FDA00001644638700034.GIF" wi="1045" he="254" />将最优相位差Δθ作为相位差校准参数;将相位差校准参数发送到被测系统的处理器,对步骤(1)得到的幅值校准后的数字信号进行相位补偿:根据主频率ω′将最优相位差Δθ转换为延迟时间τ,被测系统的处理器控制时,将延迟时间τ考虑到控制过程中,从而完成调理电路的相位校准。
地址 611731 四川省成都市高新区(西区)西源大道2006号