CN104391674A - 基于fpga的采样值线性插值运算器及运算方法 - Google Patents
基于fpga的采样值线性插值运算器及运算方法 Download PDFInfo
- Publication number
- CN104391674A CN104391674A CN201410564585.XA CN201410564585A CN104391674A CN 104391674 A CN104391674 A CN 104391674A CN 201410564585 A CN201410564585 A CN 201410564585A CN 104391674 A CN104391674 A CN 104391674A
- Authority
- CN
- China
- Prior art keywords
- value
- linear interpolation
- plus
- buddhist
- sent
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Electrophonic Musical Instruments (AREA)
Abstract
本发明涉及一种基于FPGA的采样值线性插值运算器及运算方法,该运算器包括:加/减法器、乘法器、运算数据调度器、查表单元和截位运算器,借助中间参数对线性插值运算公式逐级进行拆分,归类为只有加减法和乘法运算,采用时分复用方式按步加载中间参数进行运算,并根据当前的数据帧接收时刻间隔t2-t1和其预设下限值T查表读取比例系数k,用于与截位运算器一起代替除法运算用到线性插值运算中,获得要求的Vx值。本发明通过查表方式省去直接的除法运算,通过时分复用方式减少加/减法器、乘法器的数量,节省大量的资源。在完成既定功能的前提下使原本极为紧张的硬件资源占用情况得到有效缓解,并为后期不更换FPGA芯片的情况下进行产品升级留出了资源余量。
Description
技术领域
本发明涉及变电站采样值处理技术领域,尤其涉及一种基于FPGA的采样值线性插值运算器及运算方法。
背景技术
插值运算是智能变电站采样值处理装置的重要功能之一,用于将远方采集的数据流变换为等效于本地采集的数据流。线性插值运算是插值运算中最基本、最常用的一种,它利用在t1、t2时刻对变量Vt采集到的采样值V1、V2,同时假定在Vt-t坐标系中tx时刻的采样值Vx与V1、V2在一条直线上,从而以解线性方程的方式求得Vx的值,其运算公式为 。由于计算时需要用到除法,而除法运算占用的资源较大,导致FPGA芯片的硬件资源紧张;同时,运算过程中需要反复进行加法、减法、乘法运算,如果不采取措施则将进一步导致FPGA芯片的硬件资源紧张。 针对这些情况,若采用大容量的芯片会导致成本大幅增加,不便于该处理器的普及推广。
发明内容
为解决上述问题,本发明提供了一种基于FPGA的采样值线性插值运算器,通过对运算单元进行时分复用和查表方式降低采样值线性插值运算器的资源占用量。
本发明采用的技术方案是:
一种基于FPGA的采样值线性插值运算器,用于根据t1、t2时刻对变量Vt采集到的采样值V1、V2求得tx时刻的采样值Vx,包括:
一加/减法器,用于在运算数据调度器的指示下设置加法或减法模式,执行加法或减法运算;
一乘法器,用于在运算数据调度器的指示下执行乘法运算;
一运算数据调度器,用于按数据处理顺序加载数据至加/减法器和乘法器,向其发出运算指令并能对加/减法器进行加/减设置;
一查表单元,用于根据当前的数据帧接收时刻间隔t2-t1和其预设下限值T读取比例系数k,k用于与截位运算器一起代替除法运算用到线性插值运算中;
一截位运算器,用于对运算结果进行精确处理,获得要求的Vx值。
文中所述的帧间隔均指数据帧接收时刻间隔。
上述加/减法器所进行的运算包括t2-t1、t2-t1-T、tx-t1、tn=T-tm以及P=tm*V2+tn*V1中的加法运算。
上述乘法器所进行的运算包括tm*V2、tn*V1、Q=P*U及tm=(tx-t1)*k中的乘法运算,其中U的值通过对2q/T的计算结果四舍五入后取整获得。
上述截位运算器是对乘法器输出的Q值取绝对值后截去q位,即除以2q,之后进行四舍五入后恢复取绝对值之前的正负号,得到Vx值。
该运算器还包括一判越限器,用于将数据帧接收时刻间隔t2-t1与预设的帧间隔上限值、下限值进行比较,若越限则发出告警信号,外部电路收到此信号后对运算器各模块进行复位。
进一步,查表单元以t2-t1-T为地址查询一个FPGA芯片内置的存储区块,存储深度为未越限的帧间隔值的个数,每个帧间隔值对应于存储模块的一个存储单元的地址,查表单元读出存储单元的存储内容即为比例系数k。所述的存储区块在硬件结构上独立于周边逻辑电路,由RAM块构成。其中k值的计算存储方法是:离线计算得到每个处于帧间隔上下限值之间的t2-t1值对应的T/(t2-t1)的小数值之后,根据预定的k的实际取值范围以及系统设计要求达到的量化精度确定的k的位数来确定p的值,然后离线计算出全部的k值,存储在片外非易失存储器中,在FPGA上电时载入到FPGA片内的存储器(即表格)中以供查表操作直接读取。
上述加/减法器、乘法器、运算数据调度器、查表单元、截位运算器和判越限器在同一外部时钟信号驱动下同步工作,各时间间隔以时钟周期为单位进行计算。
该运算器在外部启动信号指示下启动,在外部复位信号指示下复位。
采用上述线性插值运算器进行采样值线性插值运算的方法,特点在于:借助中间参数对线性插值运算公式逐级进行拆分,归类为只有加减法和乘法运算,采用时分复用方式按步加载中间参数进行运算,并利用查表方式获得比例系数k,查表方式用于实现公式中除法的运算。所述的k值是针对每个处于帧间隔上下限值之间的帧间隔值t2-t1利用公式离线计算获得,并存储在外部存储器中,在FPGA上电时导入FPGA芯片内供查表操作直接读取,公式中T为帧间隔的下限值,p为由k的实际取值范围和量化精度确定的k的位数。
上述方法包括以下步骤:
(1)外部电路准备好输入数据t1、t2、V1、V2、tx后,发启动信号,启动插值运算器;
(2)运算数据调度器置加/减法器的运算模式为减法,将数据组t1、t2送入加/减法器,得到t2-t1并送到运算数据调度器;
(3)运算数据调度器将t2-t1及帧间隔下限值T送入加/减法器,得到t2-t1-T ;同时将t2-t1送入判越限器,将t2-t1与预设的帧间隔上限值、下限值进行比较,若越限则发出告警信号,外部电路收到此信号后对运算器各模块进行复位,停止所有正在执行步骤的操作,回到第(1)步;未越限,则进入下一步;
(4)运算数据调度器将数据组t1、tx送入加/减法器,得到tx-t1;同时,查表单元以t2-t1-T为地址查询一个存储区块,读出对应的存储单元中存储的比例系数k;
(5)运算数据调度器将数据组tx-t1、k送入乘法器,得到tm=(tx-t1)*k,将tm送入运算数据调度器,由其分配进入加/减法器、乘法器,分别与T相减得到tn=T-tm、与V2相乘得到tm*V2;
(6)将tn送入运算数据调度器,由其将tn、V1送入乘法器,得到tn*V1;
(7)运算数据调度器置加/减法器的运算模式为加法,将tm*V2、tn*V1送入其输入端,计算得到P=tm*V2+tn*V1;
(8)运算数据调度器将P与常数U送入乘法器,相乘得到Q=P*U,U的值通过对2q/T四舍五入后取整获得;
(9)将Q送入截位运算器,对Q取绝对值后除以2q,再四舍五入后恢复取绝对值之前的正负号,得到Vx值,然后整个运算器等待下一个启动信号。
本发明的有益效果:
1、本发明借助一查表单元调取比例系数k,提供给计算模块,实现除以t2-t1的运算,这样可以省去直接的除法运算,节省大量的资源占用量。在完成既定功能的前提下使原本极为紧张的硬件资源占用情况得到有效缓解,并为后期不更换FPGA芯片的情况下进行产品升级留出了资源余量。
2、本发明借助运算数据调度器实现对数据调用的管理,因此,仅通过一个加/减法器和一个乘法器就可实现所有的加/减法和乘法运算,节省了硬件成本。
3、本发明设有判越限器,可在数据帧接收时刻间隔超过界限值时告警并复位,确保运算器的正常运行。
4、本发明采用具有加/减法器和乘法器的任一型FPGA芯片即可实现。
附图说明
图1是本发明的电路原理框图,为便于画图,图中将运算数据调度器分为A、B两部分表示,在实际电路中属于同一个模块 ;
图2是运算数据调度器的工作原理图。
具体实施方式
一种基于FPGA的采样值线性插值运算器,用于根据t1、t2时刻对变量Vt采集到的采样值V1、V2求得tx时刻的采样值Vx,线性插值运算公式为:,如图1所示,包括加/减法器、乘法器、运算数据调度器(图中运算数据调度器A、B两部分实为同一个模块,后续描述中不予区分)、查表单元、截位运算器和判越限器,各模块通过数据传输线连接,进行数据传输。查表单元以t2-t1-T为地址查询FPGA芯片内置的存储区块,存储深度为未越限的帧间隔值的个数,每个帧间隔值对应于存储模块的一个存储单元的地址,查表单元读出存储单元的存储内容即为比例系数k。其中k值的计算存储方法是:离线计算得到每个处于帧间隔上下限值之间的t2-t1值对应的T/(t2-t1)的小数值之后,根据预定的k的实际取值范围以及系统设计要求达到的量化精度确定的k的位数来确定p的值,然后离线计算出全部的k值,存储在片外非易失存储器中,在FPGA上电时载入到FPGA片内的存储器(即表格)中以供查表操作直接读取。
所有模块在同一外部时钟信号驱动下同步工作,所有定时、时间间隔均以时钟周期为单位进行计算,整个运算器由唯一的外部启动信号启动,被唯一的外部复位信号复位。
采用上述线性插值运算器进行采样值线性插值运算的方法,包括以下步骤:
(1)外部电路准备好输入数据t1、t2、V1、V2、tx后,发启动信号,启动插值运算器;
(2)运算数据调度器置加/减法器的运算模式为减法,将数据组t1、t2送入减法器的减数、被减数端口 ,加/减法器输出计算结果t2-t1 并送到运算数据调度器;
(3)运算数据调度器将t2-t1及帧间隔下限值T送入加/减法器,得到t2-t1-T ;同时将t2-t1送入判越限器,将t2-t1与预设的帧间隔上限值、下限值进行比较,若越限则发出告警信号,外部电路收到此信号后对运算器各模块进行复位,停止所有正在执行步骤的操作,回到第(1)步;未越限,则进入下一步;
(4)运算数据调度器将数据组t1、tx送入加/减法器,得到tx-t1;同时,查表单元以t2-t1-T为地址查询一个存储区块,读出对应的存储单元中存储的比例系数k;
(5)运算数据调度器将数据组tx-t1、k送入乘法器,得到tm=(tx-t1)*k,将tm送入运算数据调度器,由其分配进入加/减法器、乘法器,分别与T相减得到tn=T-tm、与V2相乘得到tm*V2;
(6)将得到的tn送入运算数据调度器,由其将tn、V1送入乘法器,得到tn*V1;
(7)运算数据调度器置加/减法器的运算模式为加法,将tm*V2、tn*V1送入其输入端,计算得到P=tm*V2+tn*V1;
(8)运算数据调度器将P与常数U送入乘法器,相乘得到Q=P*U,U的值通过对2q/T四舍五入后取整获得;
(9)将Q送入截位运算器,对Q取绝对值后除以2q,去掉低q位,即回到“/T”的目的,再四舍五入后恢复取绝对值之前的正负号,得到Vx值,然后整个运算器等待下一个启动信号。这一步可将前述引入比例系数k运算时加入的一些参数再除去,并去除一些不精确的误差,以得到精确的Vx值。
Claims (10)
1. 一种基于FPGA的采样值线性插值运算器,用于根据t1、t2时刻对变量Vt采集到的采样值V1、V2求得tx时刻的采样值Vx,其特征在于:包括
一加/减法器,用于在运算数据调度器的指示下设置加法或减法模式,执行加法或减法运算;
一乘法器,用于在运算数据调度器的指示下执行乘法运算;
一运算数据调度器,用于按数据处理顺序加载数据至加/减法器和乘法器,向其发出运算指令并能对加/减法器进行加/减设置;
一查表单元,用于根据当前的数据帧接收时刻间隔t2-t1和其预设下限值T读取比例系数k,k用于与截位运算器一起代替除法运算用到线性插值运算中;
一截位运算器,用于对运算结果进行精确处理,获得要求的Vx值。
2.根据权利要求1所述的基于FPGA的采样值线性插值运算器,其特征在于:加/减法器所进行的运算包括t2-t1、t2-t1-T、tx-t1、tn=T-tm以及P=tm*V2+tn*V1中的加法运算。
3.根据权利要求1所述的基于FPGA的采样值线性插值运算器,其特征在于:乘法器所进行的运算包括tm*V2、tn*V1、Q=P*U及tm=(tx-t1)*k中的乘法运算,其中U的值通过对2q/T的计算结果四舍五入后取整获得。
4.根据权利要求3所述的基于FPGA的采样值线性插值运算器,其特征在于:截位运算器是对乘法器输出的Q值取绝对值后截去q位,即除以2q,之后进行四舍五入后恢复取绝对值之前的正负号,得到Vx值。
5.根据权利要求1至4任一项所述的基于FPGA的采样值线性插值运算器,其特征在于:还包括一判越限器,用于将数据帧接收时刻间隔t2-t1与预设的帧间隔上限值、下限值进行比较,若越限则发出告警信号,外部电路收到此信号后对运算器各模块进行复位。
6.根据权利要求5所述的基于FPGA的采样值线性插值运算器,其特征在于:查表单元以t2-t1-T为地址查询一个FPGA芯片内置的存储区块,存储深度为未越限的帧间隔值的个数,每个帧间隔值对应于存储模块的一个存储单元的地址,查表单元读出存储单元的存储内容即为比例系数k。
7.根据权利要求5所述的基于FPGA的采样值线性插值运算器,其特征在于:上述加/减法器、乘法器、运算数据调度器、查表单元、截位运算器和判越限器在同一外部时钟信号驱动下同步工作,各时间间隔以时钟周期为单位进行计算,取值为离散化的整数数字量。
8.应用权1所述线性插值运算器进行采样值线性插值运算的方法,其特征在于:借助中间参数对线性插值运算公式逐级进行拆分,归类为只有加减法和乘法运算,采用时分复用方式按步加载中间参数进行运算,并利用查表方式获得比例系数k,查表方式用于实现公式中除法的运算。
9.根据权利要求8所述的基于FPGA的采样值线性插值运算方法,其特征在于:所述的k值是针对每个处于帧间隔上下限值之间的帧间隔值t2-t1利用 公式离线计算获得,并存储在外部存储器中,在FPGA上电时导入FPGA芯片内供查表操作直接读取,公式中T为帧间隔的下限值,p为由k的实际取值范围和量化精度确定的k的位数。
10.根据权利要求9所述的基于FPGA的采样值线性插值运算方法,其特征在于:包括以下步骤:
(1)外部电路准备好输入数据t1、t2、V1、V2、tx后,发启动信号,启动插值运算器;
(2)运算数据调度器置加/减法器的运算模式为减法,将数据组t1、t2送入加/减法器,得到t2-t1并送到运算数据调度器;
(3)运算数据调度器将t2-t1及帧间隔下限值T送入加/减法器,得到t2-t1-T ;同时将t2-t1送入判越限器,将t2-t1与预设的帧间隔上限值、下限值进行比较,若越限则发出告警信号,外部电路收到此信号后对运算器各模块进行复位,停止所有正在执行步骤的操作,回到第(1)步;未越限,则进入下一步;
(4)运算数据调度器将数据组t1、tx送入加/减法器,得到tx-t1;同时,查表单元以t2-t1-T为地址查询一个存储区块,读出对应的存储单元中存储的比例系数k;
(5)运算数据调度器将数据组tx-t1、k送入乘法器,得到tm=(tx-t1)*k,将tm送入运算数据调度器,由其分配进入加/减法器、乘法器,分别与T相减得到tn=T-tm、与V2相乘得到tm*V2;
(6)将tn送入运算数据调度器,由其将tn、V1送入乘法器,得到tn*V1;
(7)运算数据调度器置加/减法器的运算模式为加法,将tm*V2、tn*V1送入其输入端,计算得到P=tm*V2+tn*V1;
(8)运算数据调度器将P与常数U送入乘法器,相乘得到Q=P*U,U的值通过对2q/T四舍五入后取整获得;
(9)将Q送入截位运算器,对Q取绝对值后除以2q,再四舍五入后恢复取绝对值之前的正负号,得到Vx值,然后整个运算器等待下一个启动信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410564585.XA CN104391674B (zh) | 2014-10-22 | 2014-10-22 | 基于fpga的采样值线性插值运算器及运算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410564585.XA CN104391674B (zh) | 2014-10-22 | 2014-10-22 | 基于fpga的采样值线性插值运算器及运算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104391674A true CN104391674A (zh) | 2015-03-04 |
CN104391674B CN104391674B (zh) | 2017-09-05 |
Family
ID=52609584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410564585.XA Active CN104391674B (zh) | 2014-10-22 | 2014-10-22 | 基于fpga的采样值线性插值运算器及运算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104391674B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104717508A (zh) * | 2015-03-17 | 2015-06-17 | 华中科技大学 | 一种适合硬件实现的jpeg-ls标准像素重建方法 |
CN107121581A (zh) * | 2017-04-25 | 2017-09-01 | 电子科技大学 | 一种数据采集系统的数据处理方法 |
RU182851U1 (ru) * | 2018-04-09 | 2018-09-04 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Михайловская военная артиллерийская академия" МО РФ | Устройство для комплексного оценивания уровня профессионально-должностной подготовки офицеров-артиллеристов тактического звена и определения их рейтинга в единой комплексной тактической задаче |
CN110502278A (zh) * | 2019-07-24 | 2019-11-26 | 福州瑞芯微电子股份有限公司 | 基于RiscV扩展指令的神经网络协处理器及其协处理方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1448838A (zh) * | 2003-04-24 | 2003-10-15 | 上海复旦微电子股份有限公司 | 固定时间间隔采样的插值计算方法及其装置 |
US7212139B1 (en) * | 2005-11-15 | 2007-05-01 | Texas Instruments Incorporated | System for suppressing aliasing interferers in decimating and sub-sampling systems |
CN101980416A (zh) * | 2010-10-15 | 2011-02-23 | 中国南方电网有限责任公司超高压输电公司柳州局 | 一种智能变电站系统中实现采样值同步的方法 |
CN102412806A (zh) * | 2011-10-24 | 2012-04-11 | 中兴通讯股份有限公司 | 基于逻辑电路的Farrow滤波器及其实现方法 |
CN103633741A (zh) * | 2013-12-05 | 2014-03-12 | 长园深瑞继保自动化有限公司 | 智能变电站采样值网数据交换方法和装置 |
CN103823996A (zh) * | 2014-03-18 | 2014-05-28 | 国家电网公司 | 适用于电力系统高速采样计算的插值方法 |
-
2014
- 2014-10-22 CN CN201410564585.XA patent/CN104391674B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1448838A (zh) * | 2003-04-24 | 2003-10-15 | 上海复旦微电子股份有限公司 | 固定时间间隔采样的插值计算方法及其装置 |
US7212139B1 (en) * | 2005-11-15 | 2007-05-01 | Texas Instruments Incorporated | System for suppressing aliasing interferers in decimating and sub-sampling systems |
CN101980416A (zh) * | 2010-10-15 | 2011-02-23 | 中国南方电网有限责任公司超高压输电公司柳州局 | 一种智能变电站系统中实现采样值同步的方法 |
CN102412806A (zh) * | 2011-10-24 | 2012-04-11 | 中兴通讯股份有限公司 | 基于逻辑电路的Farrow滤波器及其实现方法 |
CN103633741A (zh) * | 2013-12-05 | 2014-03-12 | 长园深瑞继保自动化有限公司 | 智能变电站采样值网数据交换方法和装置 |
CN103823996A (zh) * | 2014-03-18 | 2014-05-28 | 国家电网公司 | 适用于电力系统高速采样计算的插值方法 |
Non-Patent Citations (5)
Title |
---|
何成阶: "查表式线性插值方法在微机中的应用", 《火力与指挥控制》 * |
周亚凤 等: "基于查找表的高采样率FIR滤波器的设计", 《集成电路应用》 * |
周斌 等: "基于线性Lagrange插值法的变电站IED采样值接口方法", 《电力系统自动化》 * |
张芳玲 等: "一种数据保持的高效插值算法及其FPGA实现", 《现代雷达》 * |
李帅 等: "智能变电站采样值插值算法研究", 《2013第十四届全国保护和控制学术研讨会论文集》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104717508A (zh) * | 2015-03-17 | 2015-06-17 | 华中科技大学 | 一种适合硬件实现的jpeg-ls标准像素重建方法 |
CN107121581A (zh) * | 2017-04-25 | 2017-09-01 | 电子科技大学 | 一种数据采集系统的数据处理方法 |
CN107121581B (zh) * | 2017-04-25 | 2019-07-12 | 电子科技大学 | 一种数据采集系统的数据处理方法 |
RU182851U1 (ru) * | 2018-04-09 | 2018-09-04 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Михайловская военная артиллерийская академия" МО РФ | Устройство для комплексного оценивания уровня профессионально-должностной подготовки офицеров-артиллеристов тактического звена и определения их рейтинга в единой комплексной тактической задаче |
CN110502278A (zh) * | 2019-07-24 | 2019-11-26 | 福州瑞芯微电子股份有限公司 | 基于RiscV扩展指令的神经网络协处理器及其协处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104391674B (zh) | 2017-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104391674A (zh) | 基于fpga的采样值线性插值运算器及运算方法 | |
CN107798039B (zh) | 一种数据同步方法和装置 | |
CN106648839B (zh) | 数据处理的方法和装置 | |
CN109870172B (zh) | 计步检测方法、装置、设备及存储介质 | |
CN103389413B (zh) | 一种频谱直方图的实时统计方法 | |
CN103018557A (zh) | 归一化主从式谐波与间谐波实时分析方法 | |
US9882705B2 (en) | Communication apparatus, communication method, and computer readable medium using propagation delay for time synchronization | |
RU2537040C1 (ru) | Устройство для определения оптимального периода технического обслуживания изделия | |
CN103018547B (zh) | 归一化多处理器电功率计量方法 | |
CN105704173B (zh) | 一种集群系统数据分布方法及服务器 | |
CN109993286B (zh) | 稀疏神经网络的计算方法及相关产品 | |
CN101505173A (zh) | 一种随机接入检测方法和通信装置 | |
CN116544999A (zh) | 储能系统的装机容量确定方法、装置、设备和介质 | |
US10082956B2 (en) | Method and apparatus for downloading data including a progress bar indicating progress of downloading | |
CN102999311A (zh) | 基于Booth算法的48x30位乘法器 | |
RU2720382C1 (ru) | Устройство для определения значений характеристик готовности изделия к применению | |
CN102402270A (zh) | 中央处理器频率调节方法和系统 | |
CN204291185U (zh) | 一种时间同步图像采集记录仪 | |
RU2009144107A (ru) | Способы и устройство прогнозирования индикатора качества канала в системе связи | |
RU181260U1 (ru) | Вероятностный вычитатель | |
CN105469535A (zh) | 一种周围环境的反馈方法及终端 | |
CN102253924A (zh) | 开方运算的硬件实现方法以及开方运算器 | |
CN204808037U (zh) | 一种电子控制单元测试系统 | |
CN103854227A (zh) | 人际关系分析系统及方法 | |
CN116593768B (zh) | 一种电能计量方法、装置、电子设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |