• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 15
    • 下载费用:30 金币  

    重庆时时彩可以玩吗: 一种低频弦波信号相位差测量方法.pdf

    关 键 词:
    一种 低频 信号 相位差 测量方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610845924.0

    申请日:

    2016.09.23

    公开号:

    CN106353596A

    公开日:

    2017.01.25

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G01R 25/00申请日:20160923|||公开
    IPC分类号: G01R25/00 主分类号: G01R25/00
    申请人: 西北大学
    发明人: 唐升; 王天翔; 雷浩丹; 邓智峰; 杨育婷; 王万恒
    地址: 710127 陕西省西安市长安区学府大道1号
    优先权:
    专利代理机构: 西安铭泽知识产权代理事务所(普通合伙) 61223 代理人: 俞晓明
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610845924.0

    授权公告号:

    |||

    法律状态公告日:

    2017.03.01|||2017.01.25

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明提供一种低频弦波信号相位差测量方法,属于相位差测量领域。解决现有相位差测量存在电路结构复杂,测量精度低的问题。包括:将同步采集到的两路同频弦波信号进行无符号化处理,获取FPGA中待处理的第一信号和第二信号;依次确定第一信号的自相关测量值,第二信号的自相关测量值和第一信号与第二信号的互相关测量值;根据噪声和信号不相关特性,获取第一信号的自相关理论值,第二信号的自相关理论值和第一信号与第二信号的互相关理论值;根据两路同频弦波信号的幅度,第一信号的自相关理论值,第二信号的自相关理论值和第一信号与第二信号的互相关理论值以及任意角度函数关系式,确定两路同频弦波信号之间的相位差。

    权利要求书

    1.一种低频弦波信号相位差测量方法,其特征在于,包括:
    将同步采集到的两路同频弦波信号进行数字化处理,将数字化处理的所述两路同频弦
    波信号进行无符号化处理,获取如公式(1)所示的FPGA中待处理的第一信号和第二信号;
    通过下列公式(2)依次确定所述第一信号的自相关测量值,所述第二信号的自相关测
    量值和所述第一信号与所述第二信号的互相关测量值;
    根据噪声和信号不相关特性,获取如公式(3)所示的第一信号的自相关理论值,所述第
    二信号的自相关理论值和所述第一信号与所述第二信号的互相关理论值;
    根据所述两路同频弦波信号的幅度,所述第一信号的自相关理论值,所述第二信号的
    自相关理论值和所述第一信号与所述第二信号的互相关理论值以及任意角度函数关系式,
    通过公式(4)确定所述两路同频弦波信号之间的相位差;
    公式(1)如下所示:
    x(n)=x0(n)+M
    y(n)=y0(n)+M
    公式(2)如下所示:
    <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>x</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>
    <mrow> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mi>y</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>
    <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>y</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow>
    公式(3)如下所示:
    <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>x</mi> <mo>_</mo> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>t</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>A</mi> <mn>2</mn> </msup> <mn>2</mn> </mfrac> <mo>+</mo> <msup> <mi>M</mi> <mn>2</mn> </msup> </mrow>
    <mrow> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mi>y</mi> <mo>_</mo> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>t</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>B</mi> <mn>2</mn> </msup> <mn>2</mn> </mfrac> <mo>+</mo> <msup> <mi>M</mi> <mn>2</mn> </msup> </mrow>

    公式(4)如下所示:

    其中,x(n)为FPGA中待处理的第一信号,y(n)为FPGA中待处理的第二信号,x0(n)为采集
    到第一弦波信号,y0(n)为采集到的第二弦波信号,M为正向偏移量,Rxx_measured(0)为第一信
    号的自相关测量值,Ryy_measured为第二信号的自相关测量值,Rxy_measured(0)为第一信号与第
    二信号的互相关测量值,Rxx_theoretical(0)为第一信号的自相关理论值,Ryy_theoretical(0)为第
    二信号的自相关理论值,Rxy_theoretical(0)为第一信号与第二信号的互相关理论值,N为数据
    采样点数,A和B分别两路同频弦波信号的幅度,为两路同频弦波信号之间的相位差。
    2.如权利要求1所述的方法,其特征在于,所述无符号化处理包括,将所述数字化处理
    的第一信号和所述数字化处理的第二信号同时叠加一个相同的正向偏移量,其中所述正向
    偏移量大于所述两路同频弦波信号的幅度。

    说明书

    一种低频弦波信号相位差测量方法

    技术领域

    本发明属于相位差测量技术领域,更具体的涉及一种低频弦波信号相位差测量方
    法。

    背景技术

    在测试计量技术领域,测量低频弦波信号的相位差在工程中有着重要地位。例如,
    电网电能计量中高精度功率参数的分析和计算就涉及50Hz工频电压弦波信号与相应电流
    信号间相位差的准确测量。在时频基准计量领域,目前被中国计量研究院、美国喷气推进实
    验室等国内外研究机构和GPS导航卫星等重大工程装置广泛采用的双混频时差精密频率测
    量技术方案中,也是将5MHz或10MHz等标准频率下变频成10Hz或100Hz的低频弦波信号进行
    相位差比对测量。

    测量弦波信号的相位差最惯用的方法是使用具备相位差测量功能的通用仪器,例
    如示波器、计数器等。此类通用仪器通常采用基于脉冲填充计数的相位差测量方法,需要将
    两路待测弦波信号整形转换为方波信号,然后利用方波信号上升沿或下降沿开关闸门,利
    用脉冲填充计数的方式实现两路信号的相位差测量。这种测量方式必须有非常高的填充频
    率才能获得好的测量精度,而且不可避免存在±个字的计数误差,测量精度一般较低。

    测量弦波信号的相位差也可以使用专用的比相仪,例如HP/Agilent公司的K34-
    59991(100KHz-10MHz)、Microsemi公司的MMS多通道测试仪(1MHz-20MHz)等。但是,这些专
    用的相位差测量仪器的可测量对象的频率范围一般较高,不能满足某些领域低频信号相位
    差的测量需求。例如,电网电能计量中待测信号频率为50Hz。

    综上所述,现有的相位差测量仪器存在电路结构复杂,体积较大,且测量精度低的
    问题。

    发明内容

    本发明实施例提供一种低频弦波信号相位差测量方法,用以解决现有相位差测量
    仪器存在电路结构复杂,体积较大,且测量精度低的问题。

    本发明实施例提供一种低频弦波信号相位差测量方法包括:

    将同步采集到的两路同频弦波信号进行数字化处理,将数字化处理的所述两路同
    频弦波信号进行无符号化处理,获取如公式(1)所述的FPGA中待处理的第一信号和第二信
    号;

    通过下列公式(2)依次确定所述第一信号的自相关测量值,所述第二信号的自相
    关测量值和所述第一信号与所述第二信号的互相关测量值;

    根据噪声和信号不相关特性,获取如公式(3)所示的第一信号的自相关理论值,所
    述第二信号的自相关理论值和所述第一信号与所述第二信号的互相关理论值;

    根据所述两路同频弦波信号的幅度,所述第一信号的自相关理论值,所述第二信
    号的自相关理论值和所述第一信号与所述第二信号的互相关理论值以及任意角度函数关
    系式,通过公式(4)确定所述两路同频弦波信号之间的相位差;

    公式(1)如下所示:

    x(n)=x0(n)+M

    y(n)=y0(n)+M

    公式(2)如下所示:

    <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>x</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>

    <mrow> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mi>y</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>

    <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>y</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow>

    公式(3)如下所示:

    <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>x</mi> <mo>_</mo> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>t</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>A</mi> <mn>2</mn> </msup> <mn>2</mn> </mfrac> <mo>+</mo> <msup> <mi>M</mi> <mn>2</mn> </msup> </mrow>

    <mrow> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mi>y</mi> <mo>_</mo> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>t</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>B</mi> <mn>2</mn> </msup> <mn>2</mn> </mfrac> <mo>+</mo> <msup> <mi>M</mi> <mn>2</mn> </msup> </mrow>


    公式(4)如下所示:


    其中,x(n)为FPGA中待处理的第一信号,y(n)为FPGA中待处理的第二信号,x0(n)
    为采集到第一弦波信号,y0(n)为采集到的第二弦波信号,M为正向偏移量,Rxx_measured(0)为
    第一信号的自相关测量值,Ryy_measured(0)为第二信号的自相关测量值,Rxy_measured(0)为第一
    信号与第二信号的互相关测量值,Rxx_theoretical(0)为第一信号的自相关理论值,
    Ryy_theoretical(0)为第二信号的自相关理论值,Rxy_theoretical(0)为第一信号与第二信号的互
    相关理论值,N为数据采样点数,A和B分别两路同频弦波信号的幅度,为两路同频弦波信
    号之间的相位差。

    优选地,所述无符号化处理包括,将所述数字化处理的第一信号和所述数字化处
    理的第二信号同时叠加一个相同的正向偏移量,其中所述正向偏移量大于所述两路同频弦
    波信号的幅度。

    本发明实施例中,提供的本发明实施例提供一种低频弦波信号相位差测量方法,
    包括:将同步采集到的两路同频弦波信号进行数字化处理,将数字化处理的所述两路同频
    弦波信号进行无符号化处理,获取如公式x(n)=x0(n)+M,y(n)=y0(n)+M所述的FPGA中待处
    理的第一信号和第二信号;通过下列公式
    依次
    确定所述第一信号的自相关测量值,所述第二信号的自相关测量值和所述第一信号与所述
    第二信号的互相关测量值;根据噪声和信号不相关特性,获取如公式
    所示的
    第一信号的自相关理论值,所述第二信号的自相关理论值和所述第一信号与所述第二信号的
    互相关理论值;根据所述两路同频弦波信号的幅度,所述第一信号的自相关理论值,所述第二
    信号的自相关理论值和所述第一信号与所述第二信号的互相关理论值以及任意角度函数关系
    式,通过公式
    确定所述两路同频弦波信号之间的相位差;其中,x(n)为FPGA中待处理的第一信号,y(n)为
    FPGA中待处理的第二信号,x0(n)为采集到第一弦波信号,y0(n)为采集到的第二弦波信号,M
    为正向偏移量,Rxx_measured(0)为第一信号的自相关测量值,Ryy_measured(0)为第二信号的自相
    关测量值,Rxy_measured(0)为第一信号与第二信号的互相关测量值,Rxx_theoretical(0)为第一信
    号的自相关理论值,Ryy_theoretical(0)为第二信号的自相关理论值,Rxy_theoretical(0)为第一信
    号与第二信号的互相关理论值,N为数据采样点数,A和B分别两路同频弦波信号的幅度,
    为两路同频弦波信号之间的相位差。在本发明实施例中,采用FPGA芯片,对采集到的两路同
    频弦波信号进行处理,解决了现有顺序指令执行CPU构架模式不可避免的速度瓶颈问题,确
    保了测量的实时性,在1Hz-500Hz的频率范围内,采用FPGA芯片的系统测量相对误差小于
    1.38%,从而提高的测量精度。进一步地,用于FPGA芯片体积比较小,便于携带,从而扩大了
    该测量工具的应用场合。

    附图说明

    为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现
    有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本
    发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以
    根据这些附图获得其他的附图。

    图1为本发明实施例提供一种低频弦波信号相位差测量测量方法流程图;

    图2为本发明实施例提供的测量电路结构图;

    图3为本发明实施例提供的信号调理电路和数据同步采集电路;

    图4为本发明实施例提供的FPGA与AD芯片之间的时序控制逻辑;

    图5为本发明实施例提供的FPGA电路;

    图6为本发明实施例提供的串口电平转换电路。

    具体实施方式

    下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完
    整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例?;?br />本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他
    实施例,都属于本发明?;さ姆段?。

    图1为本发明实施例提供的一种低频弦波信号相位差测量测量方法流程图,如图1
    所示,本发明实施例提供的一种低频弦波信号相位差测量测量方法包括以下步骤:

    步骤101,将同步采集到的两路同频弦波信号进行数字化处理,将数字化处理的所
    述两路同频弦波信号进行无符号化处理,获取如公式(1)所示的FPGA中待处理的第一信号
    和第二信号;

    步骤102,通过下列公式(2)依次确定所述第一信号的自相关测量值,所述第二信
    号的自相关测量值和所述第一信号与所述第二信号的互相关测量值;

    步骤103,根据噪声和信号不相关特性,获取如公式(3)所示的第一信号的自相关
    理论值,所述第二信号的自相关理论值和所述第一信号与所述第二信号的互相关理论值;

    步骤104,根据所述两路同频弦波信号的幅度,所述第一信号的自相关理论值,所
    述第二信号的自相关理论值和所述第一信号与所述第二信号的互相关理论值以及任意角
    度函数关系式,通过公式(4)确定所述两路同频弦波信号之间的相位差;

    公式(1)如下所示:

    x(n)=x0(n)+M

    y(n)=y0(n)+M

    公式(2)如下所示:

    <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>x</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>

    <mrow> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mi>y</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>

    <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>y</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow>

    公式(3)如下所示:

    <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>x</mi> <mo>_</mo> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>t</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>A</mi> <mn>2</mn> </msup> <mn>2</mn> </mfrac> <mo>+</mo> <msup> <mi>M</mi> <mn>2</mn> </msup> </mrow>

    <mrow> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mi>y</mi> <mo>_</mo> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>t</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>B</mi> <mn>2</mn> </msup> <mn>2</mn> </mfrac> <mo>+</mo> <msup> <mi>M</mi> <mn>2</mn> </msup> </mrow>


    公式(4)如下所示:


    其中,x(n)为FPGA中待处理的第一信号,y(n)为FPGA中待处理的第二信号,x0(n)
    为采集到第一弦波信号,y0(n)为采集到的第二弦波信号,M为正向偏移量,Rxx_measured(0)为
    第一信号的自相关测量值,Ryy_measured(0)为第二信号的自相关测量值,Rxy_measured(0)为第一
    信号与第二信号的互相关测量值,Rxx_theoretical(0)为第一信号的自相关理论值,
    Ryy_theoretical(0)为第二信号的自相关理论值,Rxy_theoretical(0)为第一信号与第二信号的互
    相关理论值,N为数据采样点数,A和B分别两路同频弦波信号的幅度,为两路同频弦波信
    号之间的相位差。

    需要说明的是,本发明实施例所提供的基于FPGA的低频弦波信号相位差测量方
    法,需要基于图2提供的测量电路结构图实施,以下结合图2所提供的测量电路,具体介绍本
    发明实施例提供的基于FPGA的低频弦波信号相位差测量方法。

    如图2所示,该测量电路主要包括:信号调理电路201、数据同步采集电路202、FPGA
    电路203和串口电平转换电路204。

    具体地,信号调理电路201分别对采集到的第一弦波信号和采集到的第二弦波信
    号进行调理,将具有不同幅度且同频的第一弦波信号和第二弦波信号进行适当幅度的放大
    或衰减,使通过信号调理电路201的第一弦波信号和第二弦波信号满足后级采样电路的输
    入要求。

    数据同步采集电路202的作用是保证对调理后的两路弦波信号的每一次数据采样
    严格在同一时刻完成,同时对采样数据进行数字化处理。

    在实际应用在,数据同步采集电路202为了避免两个通道的采样时间间隔产生Δ
    t,则会终累积到采集到的两路弦波信号的相位差中,具体如以下公式(5)所示:


    在公式(5)中,f为两路弦波信号的标称频率值,为两路弦波信号间的实际相
    位差,2πf·Δt是由于非同步采样所引入的测量误差。

    在步骤101中,将同步采集到的两路同频弦波信号进行数字化处理,其中两路同频
    弦波信号分别为第一弦波信号和第二弦波信号。

    在实际应用中,对第一弦波信号和第二弦波信号进行数字化处理,主要是将两路
    同步采集的弦波信号进行量化、编码转化为数字信号。

    进一步地,将数字化处理的两路同频弦波信号进行无符号化处理,其中,无符号化
    处理包括,将数字化处理的第一信号和数字化处理的第二信号同时叠加一个相同的正向偏
    移量,其中正向偏移量大于两路同频弦波信号的幅度。

    具体地,将数字化处理的两路同频弦波信号进行无符号化处理之后,得到如公式
    (1)所示的FPGA中待处理的第一信号和第二信号:

    <mrow> <mtable> <mtr> <mtd> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>M</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>y</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>M</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>

    在公式(1)中,x(n)为FPGA中待处理的第一信号,y(n)为FPGA中待处理的第二信
    号,x0(n)为采集到第一弦波信号,y0(n)为采集到的第二弦波信号,M为正向偏移量。

    在步骤102中,FPGA电路203利用硬件描述语言执行相位差测量算法,用于两路弦
    波信号的相位差的求解运算。

    具体地,根据FPGA的相关运算主要包括三部分内容:计算通道一对应数据x(n)的
    自相关值、计算通道二数据y(n)的自相关值、计算x(n)与y(n)的互相关值,在本发明实施例
    中,通过下列公式(2)依次确定第一信号的自相关测量值,第二信号的自相关测量值和第一
    信号与第二信号的互相关测量值:

    <mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>x</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mi>y</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>y</mi> <mo>_</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>

    在公式(2)中,Rxx_measured(0)为第一信号的自相关测量值,Ryy_measured(0)为第二信
    号的自相关测量值,Rxy_measured(0)为第一信号与第二信号的互相关测量值,N为数据采样点
    数。

    在信号处理领域,可以通过下列公式(6)表示两路待测低频弦波信号的数学模型
    可表示为:

    x(t)=Asin2πf·t+NX(t)+M


    公式(6)中,A和B分别是两路信号x(t)和y(t)的幅度,NX(t)和Ny(t)分别为叠加在
    两路信号上的噪声,为两信号的相位差,M为数据预处理时叠加的偏移量,M稍大于信号
    的幅度值A和B。

    对x(t)和y(t)在时间T(一般为待测信号周期的整数倍)内进行相关运算,则可以
    有:


    需要说明的是,在信号处理领域公式(7)所示的互相关(有时也称为“互协方差”)
    是用来表示两个信号之间相似性的一个度量,通常通过与已知信号比较用于寻找未知信号
    中的特性。它是两个信号之间相对于时间的一个函数,有时也称为滑动点积。叠加在两路弦
    波信号上噪声NX(t)和Ny(t)主要来源包括AD采样量化噪声、电路内部噪声(主要包括了热噪
    声、散弹噪声等)和其他外部噪声。这些噪声叠加在一起通常具有随机白噪声的特征,所以
    可认为噪声NX(t)、Ny(t)和信号不具备相似性,即噪声和信号不相关。又因为噪声是随机的,
    所以噪声与噪声也可认为不相关。那么,信号x(t)和y(t)的互相关函数值取决于两者的相
    位差和延迟量τ有关。当τ=0时,互相关函数值Rxy(0)就主要与x(t)和y(t)之间的相位差
    有关。

    在步骤103中,根据上述分析,可以通过公式(3),依次获取第一信号的自相关理论
    值,第二信号的自相关理论值和第一信号与第二信号的互相关理论值:


    在公式(3)中,Rxx_theoretical(0)为第一信号的自相关理论值,Ryy_theoretical(0)为第
    二信号的自相关理论值,Rxy_theoretical(0)为第一信号与第二信号的互相关理论值,为两
    路同频弦波信号之间的相位差。

    通过上述公式(3)可以得到下列公式(8):

    <mrow> <mtable> <mtr> <mtd> <mrow> <mi>A</mi> <mo>=</mo> <msqrt> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mi>x</mi> <mo>_</mo> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>t</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>M</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>B</mi> <mo>=</mo> <msqrt> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mi>y</mi> <mo>_</mo> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>t</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>M</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </msqrt> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>

    进一步地,通过公式(8)和公式(3),可以得到下列公式(9):


    在实际应用中,由三角函数理论可知,对于任意存在下列公式(10):


    进一步地,通过上述公式(9)和公式(10),可以通过公式(4)确定所述两路同频弦
    波信号之间的相位差,具体地,公式(4)如下:


    综上所述,在本发明实施例中,采用FPGA芯片,对采集到的两路同频弦波信号进行
    处理,解决了现有顺序指令执行CPU构架模式不可避免的速度瓶颈问题,确保了测量的实时
    性,在1Hz-500Hz的频率范围内,采用FPGA芯片的系统测量相对误差小于1.38%,从而提高
    的测量精度。进一步地,用于FPGA芯片体积比较小,便于携带,从而扩大了该测量工具的应
    用场合。

    为了更清楚的介绍本发明实施例提供的一种低频弦波信号相位差测量方法,以下
    结合附图3,图4,图5和图6对本发明的技术方案进行详细说明。

    图3为本发明实施例提供的信号调理电路和数据同步采集电路,图4为本发明实施
    例提供的FPGA与AD芯片之间的时序控制逻辑,图5为本发明实施例提供的FPGA电路,图6为
    本发明实施例提供的串口电平转换电路。

    如图3所示,本发明实施例中前端的信号调理电路基于TI公司的运算放大器芯片
    OPA4132(U1)设计,此外包括了12个匹配电阻R1、R2、R3、R4、R5、R6、R7、R8、R9、R10、R11和
    R12,主要作用是对输入的模拟弦波信号进行等比例电平转换,使之输出的两路差分信号符
    合16位AD芯片ADS8364(U2)的±2.5V输入电压范围。两路差分信号经由两组匹配电阻R5、R6
    和R11、R12输出至U1的CHA0采样通道和CHA1采样通道。U1的采样参考电压REFin由芯片内部
    参考电压REFout提供;同时REFout电压经过电容C1、C2稳压滤波后输出至信号调理电路作
    为差分电压的参考。U1的4条控制总线CLK、EOC、RD和HOLDA负责交互传递它与FPGA之间的逻
    辑控制指令。4条控制总线与16条数据总线D0~D15的时序逻辑如图4所示。U1的16条数据总
    线D0~D15将采样量化的数据传送至FPGA的并行16位数据总线。

    本发明实施例中的FPGA选用Altera公司的Cyclone EP2C70F672C8(芯片U3)。U3的
    16条数据总线与U1的16条数据总线相连,U3的串行数据输出DX_0连接至串行口驱动电路。
    同时,U3利用硬件描述语言在芯片内部构建5部分功能电路:时序逻辑电路、数据预处理电
    路、相关器、反正切变换器和串口驱动电路,如图5所示。

    使用硬件描述语言构建时序逻辑电路,产生时钟信号CLK、采样控制信号
    数据读取使能信号输送至AD芯片U1,并读取U1反馈回来的数据准备好信号EOC,如图3所
    示。时序逻辑电路中,设计分频器对FPGA外部50MHz晶振进行5000分频作为采样控制信号,
    故信号频率为10KHz,从而实现1s采集单路弦波信号10000个离散数据点。

    时序逻辑电路中设计计数器,对采集到的离散数据点进行计数。当计数器中数值
    数累加到10000时,时序逻辑电路产生一个标识信号,作为相关器中数据锁存器的使能信
    号,同时作为数据预处理??榈氖鼓苄藕?。当计数器中数值累加到10001时,时序逻辑???br />产生第二个标识信号,作为相关器中乘加器的清零信号,同时作为反正切变换器的使能信
    号。

    本发明实施例中的串口驱动电路由波特率发生器、UART接收器和UART发送器构
    成。波特率发生器实际上就是一个分频器??梢愿莞ǖ南低呈敝悠德?50MHz晶振)和要
    求的波特率算出波特率分频因子,算出的波特率分频因子作为分频器的分频数。计算公式
    为:50000000/(16*9600)=326。UART发送器的用途是将准备输出的并行数据按照基本UART
    帧格式转为TXD信号串行输出。UART接收器接收RXD串行信号,并将其转化为并行数据。

    本返实施例中相位差测量结果选用串口RS232标准(协议)输出。RS232对逻辑电平
    作了如下规定:逻辑1为-3V~-15V,逻辑0为+3~+15V。而U3的I/O电平是:逻辑1或无信号为
    3.3V,逻辑0为0V。为了使它们进行有效信息传递,必须进行电平转换。实施例中采用
    MAX232ESE(U4)实现二者之间的电平转换,具体的硬件接口电路如图6所示。其中C3和C4为
    匹配电容,C5、C6和C7为稳压滤波电容,DX_0为U3的输出信号,DX_1是DX_0经过电平转换后
    的RS232标准串行输出。

    本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序
    产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实
    施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机
    可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产
    品的形式。

    本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程
    图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流
    程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合??商峁┱庑┘扑慊绦?br />指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产
    生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实
    现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的系统。

    这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特
    定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指
    令系统的制造品,该指令系统实现在流程图一个流程或多个流程和/或方框图一个方框或
    多个方框中指定的功能。

    这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计
    算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或
    其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一
    个方框或多个方框中指定的功能的步骤。

    尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造
    性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优
    选实施例以及落入本发明范围的所有变更和修改。

    显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精
    神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围
    之内,则本发明也意图包含这些改动和变型在内。

    关于本文
    本文标题:一种低频弦波信号相位差测量方法.pdf
    链接地址://www.4mum.com.cn/p-6014145.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 辽宁35选7风采走势图 浙江快乐12任选五遗漏 山西泳坛夺金开奖结果 教师退休后赚钱 吉林吉祥棋牌官方下载 球中奖规则 中国福利彩票官网杭州 海南环岛赛开奖查询 足彩分析基础知识 篮球让分胜负怎么看 砂锅米线赚钱吗 安徽25选5开奖视频 贵州十一选五开奖结果查询结果 云南时时彩走势图 福建快3开奖结果今天 中国福利彩票开奖结果