JPH04364477A - Phase computing apparatus - Google Patents
Phase computing apparatusInfo
- Publication number
- JPH04364477A JPH04364477A JP13986691A JP13986691A JPH04364477A JP H04364477 A JPH04364477 A JP H04364477A JP 13986691 A JP13986691 A JP 13986691A JP 13986691 A JP13986691 A JP 13986691A JP H04364477 A JPH04364477 A JP H04364477A
- Authority
- JP
- Japan
- Prior art keywords
- phase
- value
- degrees
- sin
- phase shift
- 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
- 230000010363 phase shift Effects 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims description 22
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Landscapes
- Measuring Phase Differences (AREA)
Abstract
Description
【0001】0001
【産業上の利用分野】本発明は演算処理装置を用いて位
相を算出する位相算出装置に関するものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a phase calculation device that calculates a phase using an arithmetic processing unit.
【0002】0002
【従来の技術】例えばディジタルリレーにおいては、電
力系統より電流を検出し、一定の同期した周期にてその
電流信号をサンプリングしてマイクにコンピュータに導
入し、所定のプログラムに基づいてSinθとCosθ
とを求めた後に位相を算出することが行われている。[Prior Art] For example, in a digital relay, a current is detected from a power system, the current signal is sampled at a certain synchronized period, and is introduced into a computer through a microphone, and based on a predetermined program, Sinθ and Cosθ
The phase is calculated after determining the .
【0003】図3は従来の位相算出方式を示したもので
、1は演算処理部で、ISinθとICosθの信号を
導入して演算部2でtanθを求める。変換部3では変
換部(ROM)4に蓄えられたtan−1用のデータを
参照しながら位相角に変換し、位相信号として出力する
。FIG. 3 shows a conventional phase calculation method. Reference numeral 1 denotes an arithmetic processing unit, in which ISin θ and ICos θ signals are introduced, and an arithmetic unit 2 calculates tan θ. The conversion unit 3 converts the data into a phase angle while referring to the data for tan-1 stored in the conversion unit (ROM) 4, and outputs it as a phase signal.
【0004】0004
【発明が解決しようとする課題】図3で示す従来の方法
では、位相を求めるためにtan−1用のデータを用意
しなければならない。このデータは通常ROM領域に置
かれるためプログラム容量が減少する。例えば、位相角
の刻みを0.1度とすると、データは90度分必要であ
り、且つ1データ2バイトであるので(90×10)×
2=1800バイト必要となる。これは、1命令を平均
3バイトとすると600命令分のデータが必要となり、
その分プログラム容量が減少する問題を有している。In the conventional method shown in FIG. 3, data for tan-1 must be prepared in order to determine the phase. Since this data is usually placed in the ROM area, the program capacity is reduced. For example, if the phase angle increments are 0.1 degrees, 90 degrees of data are required, and each data is 2 bytes, so (90×10)×
2=1800 bytes are required. This means that if one instruction is 3 bytes on average, data for 600 instructions is required.
The problem is that the program capacity decreases accordingly.
【0005】[0005]
【課題を解決するための手段】本発明は位相を求めるも
のにおいて、SinθとCosθの信号をもとに一定角
度の移相演算をする移相演算部と、求められた移相値に
対応したオフセット値を決定するオフセット値決定部と
、移相された角度の領域に分けて級数展開する級数展開
部と、展開された数値に数値に各領域のオフセット値を
加算して位相を算出する加算部とを設けたものである。[Means for Solving the Problems] The present invention provides a device for determining a phase, which includes a phase shift calculation section that performs a phase shift calculation of a constant angle based on signals of Sin θ and Cos θ, and a An offset value determination unit that determines the offset value, a series expansion unit that performs series expansion by dividing the phase-shifted angle into regions, and an addition that calculates the phase by adding the offset value of each region to the expanded value. It has a section.
【0006】[0006]
【作用】移相演算部はSinθ,CosθからSinθ
を移相してSin(θ−K′)を求め、オフセット値決
定部において、Sin(θ−K′)の値が0<Sin(
θ−K′)<α(但しα<1)とし、級数展開部および
加算部において、Sin−1θ=Sin−1(θ−K′
)+Kを求めることによって位相を算出する。[Operation] The phase shift calculation section changes from Sinθ, Cosθ to Sinθ
is phase-shifted to obtain Sin(θ-K'), and in the offset value determination section, the value of Sin(θ-K') is determined as 0<Sin(
θ-K') < α (however, α<1), and in the series expansion section and addition section, Sin-1θ=Sin-1(θ-K'
)+K to calculate the phase.
【0007】[0007]
【実施例】図1は本発明の一実施例を示したものである
。同図において、10は演算処理部で、この処理部10
には、SinθとCosθの交流信号から複数の異なる
移相角のSin(θ−K′)を求める移相演算部11と
、各移相角に対応したオフセット値Kを決定するオフセ
ット値決定部12と、移相演算部11にて求められたS
in(θ−K′)をもとに級数展開する級数展開部13
と、この級数展開部によって得られた値にオフセット値
を加算する加算部14とを有している。20はSinθ
とCosθとを記憶する記憶部、30は交流信号を導入
するインタフェースである。Embodiment FIG. 1 shows an embodiment of the present invention. In the same figure, 10 is an arithmetic processing section, and this processing section 10
includes a phase shift calculation unit 11 that calculates Sin(θ−K′) of a plurality of different phase shift angles from an AC signal of Sinθ and Cosθ, and an offset value determination unit that determines an offset value K corresponding to each phase shift angle. 12 and S obtained by the phase shift calculation unit 11
Series expansion unit 13 that performs series expansion based on in(θ−K′)
and an addition section 14 that adds an offset value to the value obtained by the series expansion section. 20 is Sinθ
and Cos θ, and 30 is an interface for introducing an AC signal.
【0008】以上のように構成された本発明において、
図2を用いてその動作を説明する。インタフェース30
に、例えば交流のアナログ信号が導入されたとすると、
このインタフェース30でディジタル信号に変換された
後、演算処理部10でSinθおよびCosθの演算が
なされて記憶部20に一旦記憶される。移相演算部11
は、ステップS1においてSinθとCosθを用い、
且つ移相角K′を60度と仮定して
x=SinθCosθ60°−CosθSin60°を
演算してSin(θ−60°)を求める。[0008] In the present invention configured as above,
The operation will be explained using FIG. 2. interface 30
For example, if an AC analog signal is introduced into
After being converted into a digital signal by this interface 30, the arithmetic processing unit 10 calculates Sin θ and Cos θ and temporarily stores it in the storage unit 20. Phase shift calculation unit 11
uses Sinθ and Cosθ in step S1,
Assuming that the phase shift angle K' is 60 degrees, x=SinθCosθ60°−CosθSin60° is calculated to obtain Sin(θ−60°).
【0009】オフセット値決定部12は、ステップS2
において、S1で求められたxが0より大きいか否かを
判断し、大の場合にはS3でSin(θ−60°)にお
けるオフセット値Kを60に決定して記憶部に保存する
。[0009] The offset value determining unit 12 performs step S2.
In S1, it is determined whether or not the obtained x is greater than 0. If it is, in S3, the offset value K at Sin (θ-60°) is determined to be 60 and is stored in the storage unit.
【0010】S2でxが0より小の場合には、S4で再
度移相演算部11において移相演算する。この場合の移
相角は30度と仮定し、x=SinθCosθ30°−
CosθSin30°を演算する。オフセット値決定部
12は、ステップS5においてS4で求められたxが0
より大きいか否かを判断し、大の場合にはS6でオフセ
ット値Kを30に決定して記憶部に記憶する。また、x
が0より小の場合にはステップS7においてx=Sin
θとしてオフセット値を0とする。そしてステップS8
で級数展開部13は、移相演算部11において求められ
たx=Sin(θ−K′)但しK′は0°,30°,6
0°,90°の30°間隔)を用いて(0°≦x≦30
°),(30°<x≦60°),(60°<x≦90°
)の領域に分けて0°≦x≦30°の範囲で級数展開を
し、Sin−1xを求める。ステップS9で加算部14
は、30°<x≦60°の範囲では30に、60°<x
≦90°の範囲では60に夫々領域毎に決定されている
各領域のオフセット値をSin−1xに加算することに
よって位相を算出する。If x is smaller than 0 in S2, the phase shift calculation unit 11 performs phase shift calculation again in S4. In this case, the phase shift angle is assumed to be 30 degrees, and x=SinθCosθ30°−
Calculate CosθSin30°. In step S5, the offset value determining unit 12 determines that x obtained in S4 is 0.
If it is larger, the offset value K is determined to be 30 in S6 and is stored in the storage unit. Also, x
is smaller than 0, x=Sin in step S7
The offset value is set to 0 as θ. And step S8
Then, the series expansion unit 13 calculates x=Sin(θ−K′) obtained in the phase shift calculation unit 11, where K′ is 0°, 30°, 6
(0°≦x≦30
°), (30°<x≦60°), (60°<x≦90°
) and perform series expansion in the range 0°≦x≦30° to find Sin-1x. In step S9, the adder 14
is 30 in the range of 30°<x≦60°, and 60°<x
In the range of ≦90°, the phase is calculated by adding the offset value of each region determined for each region to 60 to Sin-1x.
【0011】なお上記実施例においては、角度を30度
分割にしてオフセット値を決定しているが、この角度は
、例えば15度でもよく、その場合には実施例の倍であ
る6個のオフセット値決定のための演算が必要となる。
級数展開部11は、[0011] In the above embodiment, the offset value is determined by dividing the angle into 30 degrees, but this angle may also be, for example, 15 degrees. Calculations are required to determine the value. The series expansion unit 11 is
【0012】0012
【数1】[Math 1]
【0013】(但し│x│≦│)の演算を行うが、xが
0〜0.5程度であれば右辺3項までで充分近似できる
。
すなわち、x=0.5のとき、 真値
0.523598……
第3
項までの近似 0.523177……で、その誤
差が略0.025度となる。このため前述した3領域に
分けて(1)式第3項までの演算でSin−1xを求め
ることで実用的には充分である。また、角度の分割領域
を30度より小さく(18度)し、第3項までの近似で
x=0.309の場合には誤差は更に小さくなり約0.
0007度となる。The following calculation is performed (where |x|≦|); if x is about 0 to 0.5, it can be sufficiently approximated with up to three terms on the right side. That is, when x=0.5, the true value
0.523598...
Third
The approximation to the term is 0.523177..., and the error is approximately 0.025 degree. Therefore, it is practically sufficient to calculate Sin-1x by dividing into the three regions described above and calculating up to the third term of equation (1). Furthermore, if the angle division area is made smaller than 30 degrees (18 degrees) and x=0.309 in the approximation up to the third term, the error becomes even smaller, about 0.
It becomes 0007 degrees.
【0014】[0014]
【発明の効果】以上本発明によれば、Sinθを任意角
度に移相し、移相された領域に分けて或る範囲内で級数
展開してSin−1xを求め、更に各領域のオフセット
値を加算することによって比較的少ないプログラムによ
って位相を算出するようにしたものであるから、従来の
ようなtan−1θ用のデータが不要となり、その分の
ROM領域が確保できるものである。[Effects of the Invention] According to the present invention, the phase of Sin θ is shifted to an arbitrary angle, and Sin-1x is obtained by dividing the phase-shifted region into series expansion within a certain range, and further calculating the offset value of each region. Since the phase is calculated using a relatively small number of programs by adding , there is no need for conventional data for tan-1θ, and the ROM area can be secured accordingly.
【図1】本発明の一実施例を示す構成図。FIG. 1 is a configuration diagram showing an embodiment of the present invention.
【図2】本発明を説明するための位相算出演算のフロー
チャート。FIG. 2 is a flowchart of phase calculation calculation for explaining the present invention.
【図3】従来の位相検出装置の構成図。FIG. 3 is a configuration diagram of a conventional phase detection device.
10…演算処理部、11…移相演算部、12…オフセッ
ト値決定部、13…級数展開部、14…加算部。DESCRIPTION OF SYMBOLS 10... Arithmetic processing section, 11... Phase shift calculation section, 12... Offset value determination section, 13... Series expansion section, 14... Addition section.
Claims (1)
て位相を算出するものにおいて、前記SinθとCos
θの信号と任意位相角の値K′をもとに移相演算して移
相値Sin(θ−K′)を求める移相演算部と、この移
相演算部によって求められた値と所定比較値と比較し、
比較条件成立時に前記任意位相角時におけるオフセット
値Kと判断するオフセット値決定部と、前記移相演算部
とオフセット値決定部とによって複数位相角に対応した
移相値を求め、各位相値の領域に分けて級数展開する級
数展開演算部と、この演算部によって求められた値に前
記各オフセット値を加算する加算部を備えたことを特徴
とする位相算出装置。[Claim 1] In a device that calculates a phase by introducing signals of Sin θ and Cos θ, the signals of Sin θ and Cos θ are
A phase shift calculation section that calculates a phase shift value Sin (θ-K') by performing a phase shift calculation based on the θ signal and the arbitrary phase angle value K', and Compare with the comparison value,
An offset value determination unit that determines the offset value K at the arbitrary phase angle when the comparison condition is satisfied, and the phase shift calculation unit and offset value determination unit calculate phase shift values corresponding to a plurality of phase angles, and determine the phase shift values for each phase value. 1. A phase calculation device comprising: a series expansion calculation unit that performs series expansion in areas; and an addition unit that adds each of the offset values to the values obtained by the calculation unit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3139866A JP3042026B2 (en) | 1991-06-12 | 1991-06-12 | Phase calculator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3139866A JP3042026B2 (en) | 1991-06-12 | 1991-06-12 | Phase calculator |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH04364477A true JPH04364477A (en) | 1992-12-16 |
JP3042026B2 JP3042026B2 (en) | 2000-05-15 |
Family
ID=15255378
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP3139866A Expired - Fee Related JP3042026B2 (en) | 1991-06-12 | 1991-06-12 | Phase calculator |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3042026B2 (en) |
-
1991
- 1991-06-12 JP JP3139866A patent/JP3042026B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP3042026B2 (en) | 2000-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH0783964A (en) | Method for estimating amplitude and frequency of sine-wave signal, method for obtaining display of signal characteristic and system for monitoring plurality of analog-signal states in electronic circuit | |
JPH04364477A (en) | Phase computing apparatus | |
JP2001516544A (en) | Tone detection using discrete Fourier transform technology | |
JP3508643B2 (en) | Resolver signal processing device and signal processing method | |
JP2508036B2 (en) | Digital Directional Ground Fault Relay | |
JPH0627798B2 (en) | Automatic correction circuit of direction finding device | |
JPH08285933A (en) | Gps receiver | |
JPS59125014A (en) | Composite bow bearing measuring device | |
JPS6140332B2 (en) | ||
JPH08285932A (en) | Gps receiving apparatus | |
JPH01257274A (en) | Frequency estimating system | |
JPH033424A (en) | Receiver | |
KR100189937B1 (en) | Rotation angle estimation device and method using sinusoidal wave | |
JPS60205269A (en) | Under water sound detector | |
JP2002009856A (en) | Inverse tangent arithmetic circuit in digital signal processing | |
JP2616601B2 (en) | Disk servo controller | |
JPH07110349A (en) | Measuring device for phase angle | |
JPH0746745A (en) | Digital-type protective relay | |
JPS58216910A (en) | Position measuring method in radio navigation system | |
JPH0426066B2 (en) | ||
JPS60144878A (en) | Detector for constant of circle | |
JPH05172872A (en) | Digital phase detecting device | |
JP3609320B2 (en) | Angle demodulator, demodulation characteristic correction apparatus, and recording medium | |
JPH09212340A (en) | Trigonometric function calculating circuit | |
JPS63288503A (en) | Fm demodulator |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |