CN101865651B - 一种旋转变压器角度信号解码方法 - Google Patents
一种旋转变压器角度信号解码方法 Download PDFInfo
- Publication number
- CN101865651B CN101865651B CN2010101972318A CN201010197231A CN101865651B CN 101865651 B CN101865651 B CN 101865651B CN 2010101972318 A CN2010101972318 A CN 2010101972318A CN 201010197231 A CN201010197231 A CN 201010197231A CN 101865651 B CN101865651 B CN 101865651B
- Authority
- CN
- China
- Prior art keywords
- complex
- cos
- sin
- sequence
- rotary transformer
- 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.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
- Transmission And Conversion Of Sensor Element Output (AREA)
Abstract
Description
技术领域
本发明是属于传感器测量技术领域,涉及一种旋转变压器角度信号解码方法。
背景技术
旋转变压器作为一种角度检测元件,具有精度高、结构简单,抗干扰能力强,已经广泛应用于环境恶劣而精度要求高的场合下的角度测量。
由于旋转变压器的输出实际上是调制了角度信息的载波信号,因此从中提取角度信号需要设计旋转变压器角度信号解码方法。目前的解码方法都是在信号的时间域内设计算法,例如在峰值处采样,然后利用求反正切的方法来计算角度等等,存在着易于受干扰影响,鲁棒性差等问题。
发明内容
本发明所要解决的技术问题是针对现有技术的上述不足,提出一种易于实现鲁棒性好的基于旋转变压器的角度信号解码方法。本发明采用如下的技术方案:
一种旋转变压器角度信号解码方法,包括下列步骤:
(1)设旋转变压器的励磁信号为uin(t)=Uincos(2πfint),其两路输出信号分别为ucos(t)=KUincos(2πfint)cos(θ)和usin(t)=KUincos(2πfint)sin(θ),其中,t表示时间,θ为角度信号,fin为励磁信号uin(t)的频率,Uin为励磁信号的幅值,K为由旋转变压器的结构确定的变比常数;取采样频率fs=nfin,n为大于2的正整数,对ucos(t)和usin(t)两路输出信号连续采集N次,令N=n*m,m为大于2的正整数,得到如下的两路数字输出信号:ucos(k)=KUincos(2πk/n)cos(θ)和
usin(k)=KUincos(2πk/n)sin(θ);
(2)利用ucos(k)和usin(k)构造复数数字信号序列:ucomplex(k)=ucos(k)+i·usin(k),k=0,1,2,...,N-1;
(3)对复数数字信号序列ucomplex(k)进行离散傅里叶变换,得到一个长度为N的复数序列,当N为奇数时,该复数序列可表示为:
Ucomplex(f-L),......Ucomplex(f0),......,Ucomplex(fL),其中L=(N-1)/2,
当N为偶数时,该复数序列可表示为:
Ucomplex(f-L),......Ucomplex(f0),......,Ucomplex(fL-1),其中L=(N)/2;
设fl表示频率域的离散自变量,fl的代表的实际频率值由公式计算,则傅立叶变换复数序列的解析表达式为:
(4)当l=m时,Ucomplex(fl=fm)=KUin(cos(θ)+i·sin(θ)),在傅立叶变换复数序列中找到fl=fm处的取值Ucomplex(fm),然后计算Ucomplex(fm)的复数角,即得到位置信号θ值。
由于在实际应用过程中各种干扰往往对信号的时域幅度信息影响较大,例如脉冲干扰,载波信号畸变等等,而对信号的频率特性影响不明显。利用这个特点,本发明提出一种基于傅立叶变换的旋转变压器角度信号解码的方法,首先把旋转变压器的输出信号变换到频率域,然后通过计算在载波频率处的傅立叶变换结果恢复角度信号。由于本方法在旋转变压器输出信号的频率域设计角度信号解码算法,对各种干扰源的具有较高的鲁棒性和抗干扰能力。
具体实施方式
旋转变压器是一种输出信号电压随转子的转角变化的角度传感器,主要由一相励磁绕组和两相输出绕组组成。当励磁绕组输入一定频率的交流电压时,两相输出绕组分别调制了转角的余弦值和正弦值,如下式所示:
uin(t)=Uincos(2πfint) (1)
ucos(t)=KUincos(2πfint)cos(θ)(2)
usin(t)=KUincos(2πfint)sin(θ)(3)
其中t表示时间,θ为角度信号,uin(t)是(t)励磁信号,ucos(t)和usin(t)分别是旋转变压器的输出信号。fin为励磁信号的频率,Uin为励磁信号的幅值,K为变比常数(由旋转变压器的结构确定)。
对旋转变压器信号进行解码的目的就是利用uin(t),ucos(t)和usin(t),从中回复θ值。
1.信号采集
为了对信号进行离散傅里叶变换,需要把模拟信号ucos和usin转换为数字信号。取采样频率fs=nfin,n为大于2的正整数。假设在时间间隔Δt内,θ值的变化可以忽略不计,对ucos(t)和usin(t)连续采集N次。令N=n*m,m为大于2的正整数,得到数字信号如下:
ucos(k)=KUincos(2πk/n)cos(θ)(4)
usin(k)=KUincos(2πk/n)sin(θ)(5)
利用ucos(k)和usin(k)构造复数数字信号序列:
ucomplex(k)=ucos(k)+i·usin(k),k=0,1,2,...,N-1
2.离散傅里叶变换
对复数信号ucomplex(k)进行离散傅里叶变换,可以得到一个与ucomplex(k)等长的复数序列。在实际计算ucomplex(k)的傅里叶变换时,可采用快速傅里叶算法以提高计算速度。
对长度为N的复数数字信号序列ucomplex(k)进行离散傅里叶变换,得到一个长度为N的复数序列。当N为奇数时,该复数序列可表示为:
Ucomplex(f-L),......Ucomplex(f0),......,Ucomplex(fL),其中L=(N-1)/2(6)
当N为偶数时,该复数序列可表示为:
Ucomplex(f-L),......Ucomplex(f0),......,Ucomplex(fL-1),其中L=(N)/2;(7)
如果利用fl表示频率域的离散自变量,该复数序列可统一表示为:
fl代表的实际频率值可按由公式计算。
3.计算角度值θ
根据离散傅里叶变换的定义,Ucomplex(fl)可由下式计算:
在不考虑噪声及信号畸变的情况下,利用(4)、(5)计算上式可得:
(8)
由(8)式可知,当l=m时,
Ucomplex(fl=fm=fin)=KUin(cos(θ)+i·sin(θ)),
因此,在(6)或(7)的离散傅里叶变换结果中找到fl=fm处的取值Ucomplex(fm),然后计算Ucomplex(fm)的复数角,即得到位置信号θ值。
Claims (1)
1.一种旋转变压器角度信号解码方法,包括下列步骤:
(1)设旋转变压器的励磁信号为uin(t)=Uincos(2πfint),其两路输出信号分别为ucos(t)=KUincos(2πfint)cos(θ)和usin(t)=KUincos(2πfint)sin(θ),其中,t表示时间,θ为角度信号,fin为励磁信号uin(t)的频率,Uin为励磁信号的幅值,K为由旋转变压器的结构确定的变比常数;取采样频率fs=nfin,n为大于2的正整数,对ucos(t)和usin(t)两路输出信号连续采集N次,令N=n*m,m为大于2的正整数,得到如下的两路数字输出信号:
ucos(k)=KUincos(2πk/n)cos(θ)和usin(k)=KUincos(2πk/n)sin(θ);
(2)利用ucos(k)和usin(k)构造复数数字信号序列:ucomplex(k)=ucos(k)+i·usin(k),k=0,1,2,...,N-1;
(3)对复数数字信号序列ucomplex(k)进行离散傅里叶变换,得到一个长度为N的复数序列,当N为奇数时,该复数序列表示为:
Ucomplex(f-L),......Ucomplex(f0),......,Ucomplex(fL),其中L=(N-1)/2,
当N为偶数时,该复数序列表示为:
Ucomplex(f-L),......Ucomplex(f0),......,Ucomplex(fL-1),其中L=(N)/2;
(4)当l=m时,Ucomplex(fl=fm)=KUin(cos(θ)+i·sin(θ)),在傅立叶变换复数序列中找到fl=fm处的取值Ucomplex(fm),然后计算Ucomplex(fm)的复数角,即得到位置信号θ值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101972318A CN101865651B (zh) | 2010-06-11 | 2010-06-11 | 一种旋转变压器角度信号解码方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101972318A CN101865651B (zh) | 2010-06-11 | 2010-06-11 | 一种旋转变压器角度信号解码方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101865651A CN101865651A (zh) | 2010-10-20 |
CN101865651B true CN101865651B (zh) | 2011-08-10 |
Family
ID=42957469
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101972318A Expired - Fee Related CN101865651B (zh) | 2010-06-11 | 2010-06-11 | 一种旋转变压器角度信号解码方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101865651B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721362B (zh) * | 2012-06-27 | 2015-05-20 | 同济大学 | 旋转变压器位置测量***及方法 |
CN103762912B (zh) * | 2013-12-24 | 2017-11-03 | 联合汽车电子有限公司 | 旋转变压器软解码*** |
CN104796051B (zh) * | 2015-04-17 | 2017-08-01 | 天津大学 | 基于fpga的旋转变压器线性变换方法及变换器 |
CN108981758B (zh) * | 2018-08-15 | 2020-10-16 | 江门市蒙德电气股份有限公司 | 一种旋转变压器的运动误差补偿方法及装置 |
CN112382488B (zh) * | 2020-10-20 | 2021-12-28 | 上海赢双电机有限公司 | 极数大于齿槽数的磁阻式旋转变压器及其制备方法 |
CN113131800B (zh) * | 2021-03-25 | 2022-12-27 | 格至控智能动力科技(上海)有限公司 | 永磁同步电机旋变信号波动测量装置及其使用方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3693280B2 (ja) * | 1999-08-16 | 2005-09-07 | オークマ株式会社 | レゾルバの組立装置 |
US6754610B2 (en) * | 2001-05-16 | 2004-06-22 | Raytheon Company | Digital signal processing of resolver rotor angle signals |
JP5041419B2 (ja) * | 2007-12-28 | 2012-10-03 | 東芝機械株式会社 | レゾルバ装置およびレゾルバの角度検出装置とその方法 |
-
2010
- 2010-06-11 CN CN2010101972318A patent/CN101865651B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101865651A (zh) | 2010-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101865651B (zh) | 一种旋转变压器角度信号解码方法 | |
CN101825660B (zh) | 欠采样下的正弦信号频率的高效测量方法及实施装置 | |
CN100463371C (zh) | 正交频分复用***中的粗糙频率同步方法和装置 | |
CN103941089B (zh) | 基于dft的正弦信号频率估计方法 | |
CN105738696A (zh) | 全相位时移相位差频率估计方法及装置 | |
CN104796051B (zh) | 基于fpga的旋转变压器线性变换方法及变换器 | |
CN102487265A (zh) | 用于自适应补偿解算器的位置误差的方法 | |
CN103364783B (zh) | 基于单通道sar的运动目标径向速度无模糊估计方法 | |
CN104391331A (zh) | 一种用于距离探测的巨磁阻抗传感器装置及信号处理方法 | |
CN102269803B (zh) | 基于时间延迟的离散频谱低频成分的校正方法 | |
CN107490722A (zh) | 一种低信噪比实信号的频率估计方法 | |
CN108318737A (zh) | 一种高精度数字鉴相方法 | |
CN105137180A (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
CN101561501A (zh) | 基于dsp的雷达目标跟踪识别器 | |
CN103018555A (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN103776471A (zh) | 一种基于双同步旋转坐标系的磁编码器 | |
CN106534014A (zh) | 一种多分量lfm信号的精确检测与分离方法 | |
CN109342813A (zh) | 一种基于dft和二分法的正弦信号频率估计方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN101718816B (zh) | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN103162614A (zh) | 一种角位移传感器在线自标定方法 | |
CN102401664A (zh) | 位置检测器 | |
CN102809688B (zh) | 一种基于迭代傅里叶变换计算三相电压实时相位的方法 | |
CN103575979A (zh) | 一种数字化测量交流电频率的方法 | |
CN106027454A (zh) | 基于频率补偿的qam信号的频偏估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110810 Termination date: 20120611 |