Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JPH0830712B2 - Exponential average wind vane - Google Patents
[go: Go Back, main page]

JPH0830712B2 - Exponential average wind vane - Google Patents

Exponential average wind vane

Info

Publication number
JPH0830712B2
JPH0830712B2 JP2142508A JP14250890A JPH0830712B2 JP H0830712 B2 JPH0830712 B2 JP H0830712B2 JP 2142508 A JP2142508 A JP 2142508A JP 14250890 A JP14250890 A JP 14250890A JP H0830712 B2 JPH0830712 B2 JP H0830712B2
Authority
JP
Japan
Prior art keywords
value
angle
exponential average
wind direction
exponential
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 - Lifetime
Application number
JP2142508A
Other languages
Japanese (ja)
Other versions
JPH0434372A (en
Inventor
晴彦 平野
利光 横田
弘 鈴木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yokogawa Electric Corp
Original Assignee
Yokogawa Electric Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yokogawa Electric Corp filed Critical Yokogawa Electric Corp
Priority to JP2142508A priority Critical patent/JPH0830712B2/en
Publication of JPH0434372A publication Critical patent/JPH0434372A/en
Publication of JPH0830712B2 publication Critical patent/JPH0830712B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は、風向を測定する計器に利用する。特に、指
数平均風向に基づき風向方位を決定する演算手段に関す
る。
DETAILED DESCRIPTION OF THE INVENTION [Industrial field of application] The present invention is used for an instrument for measuring a wind direction. In particular, the present invention relates to a calculation means for determining the wind direction based on the exponential average wind direction.

〔概要〕〔Overview〕

本発明は、風向瞬間値を示す入力値に基づき風向方位
を決定する指数平均風向計において、 風向瞬間値を示す入力値に基づきテーブルから抽出し
た角度値を用いて指数平均演算を行うことにより、 指数平均演算の繰し返しを避けかつ計算アルゴリズム
が複雑な逆正接値の算出を省略して短い処理時間で風向
方位を決定するすることができるようにしたものであ
る。
The present invention, in the exponential average anemometer for determining the wind direction based on the input value indicating the wind direction instantaneous value, by performing the index averaging operation using the angle value extracted from the table based on the input value indicating the wind direction instantaneous value, It is possible to determine the wind direction in a short processing time by avoiding repetition of exponential averaging and omitting calculation of a complex arctangent value by a calculation algorithm.

〔従来の技術〕[Conventional technology]

従来例では、風向の瞬間入力値の正弦成分と余弦成分
との各々の指数平均値を求め、この値に基づき算出され
た正弦値の逆正接値をさらに指数平均して指数平均値風
向を決定している。
In the conventional example, the exponential average value of each of the sine component and the cosine component of the instantaneous wind direction input value is obtained, and the arctangent value of the sine value calculated based on this value is further exponentially averaged to determine the exponential average value wind direction. are doing.

〔発明が解決しようとする問題点〕[Problems to be solved by the invention]

このような従来例では、指数平均演算を繰り返し行
い、しかも、計算アルゴリズムが複雑な逆正接値の算出
を必要とするので、CPUの処理時間が長くなり、したが
って、サンプリング周期の短い入力処理するには高速度
演算能力を持つ計算機を必要とする欠点があった。
In such a conventional example, the exponential averaging operation is repeatedly performed, and moreover, the calculation algorithm needs to calculate a complex arctangent value, so that the processing time of the CPU becomes long, and therefore, the input processing with a short sampling cycle is performed. Has the drawback of requiring a computer with high speed computing power.

本発明はこのような欠点を除去するもので、簡単な計
算アルゴリズムで平均風向を決定することができる指数
平均風向計を提供することを目的とする。
The present invention eliminates such drawbacks, and an object of the present invention is to provide an exponential average wind vane capable of determining the average wind direction with a simple calculation algorithm.

〔課題を解決するための手段〕[Means for solving the problem]

本発明は、風向瞬間値を示す入力値に基づき風向方位
を決定する指数平均風向計において、角度0°ないし36
0°の範囲が上記入力値の分解能に相応の区分に分割さ
れ、この区分を代表する角度値がそれぞれの区分に割り
付けられたデータテーブルを格納し、この区分を代表す
る角度値が上記入力値でアクセスできるROMと、このROM
からアクセスした角度値と前回までの入力値に対する指
数平均角度値との角度差の値を参照してこの角度値に相
当の絶対角度差を演算する絶対角度演算手段と、この絶
対角度演算手段で得られた絶対角度に対する指数平均値
を演算する指数平均演算手段と、この指数平均演算手段
で演算された指数平均値を風向方位に変換する方位変換
手段とを備えたことを特徴とする。
The present invention is an exponential average wind vane that determines a wind direction based on an input value indicating an instantaneous wind direction, and an angle of 0 ° to 36 °
The range of 0 ° is divided into divisions corresponding to the resolution of the input value, and the data table in which the angle values representing this division are assigned to each division is stored. ROM that can be accessed with and this ROM
The absolute angle calculation means for calculating the absolute angle difference corresponding to this angle value by referring to the value of the angle difference between the angle value accessed from An exponential average calculating means for calculating an exponential average value with respect to the obtained absolute angle and an azimuth converting means for converting the exponential average value calculated by the exponential average calculating means into a wind direction.

〔作用〕[Action]

風向瞬間値を示す入力値に基づきデータテーブルに登
録された角度値を抽出する。今回値と前回までの平均値
との差が180°以上であるか否かを判定し、この判定結
果で絶対角度差を演算し、つづいて、指数平均演算を行
う。この演算値が360°以上であるときに360°を減算
し、この指数平均演算値を方位名に対応付ける。例え
ば、16方位に対応付けるには、11.25°を加算して22.5
°で除算した無名数が所定の範囲内にあれば、この範囲
に対応する方位名をもってこの指数平均演算値に対する
方位名とする。
The angle value registered in the data table is extracted based on the input value indicating the wind direction instantaneous value. It is determined whether or not the difference between the current value and the average value up to the previous time is 180 ° or more, the absolute angle difference is calculated based on the result of this determination, and then the exponential average calculation is performed. When this calculated value is 360 ° or more, 360 ° is subtracted and this exponential average calculated value is associated with the bearing name. For example, to associate with 16 bearings, add 11.25 ° and add 22.5
If the anonymous number divided by ° is within a predetermined range, the direction name corresponding to this range is used as the direction name for this exponential average calculation value.

この過程では、指数平均演算の繰り返しが避けられか
つ計算アルゴリズムが複雑な逆正接値の算出が省略され
るので、短い処理時間で風向方位名を決定することがで
きる。
In this process, repetition of exponential averaging is avoided and calculation of arctangent value, which is complicated in the calculation algorithm, is omitted, so that the wind direction name can be determined in a short processing time.

〔実施例〕〔Example〕

以下、本発明の一実施例を図面に基づき説明する。 An embodiment of the present invention will be described below with reference to the drawings.

まず、この実施例の構成を第1図および第2図に基づ
き説明する。この実施例は、第1図に示すように、0°
ないし360°の範囲を有し、風向瞬間値を示す入力値の
分解能に相応の区分にこの範囲が分割され、区分を代表
する角度値がそれぞれの区分に割り付けられ、この区分
を代表する角度値を入力値でアクセスできるデータテー
ブル11を格納するROM1と、このROM1からアクセスした角
度値と前回までの入力値に対する指数平均角度値との角
度差の値を参照してこの角度値に相当の絶対角度差を演
算する絶対角度演算手段2と、この手段で得られた絶対
角度に対する指数平均値を演算する指数平均演算手段3
と、この手段で演算された指数平均値を方位に変換する
方位変換手段4を備える。
First, the structure of this embodiment will be described with reference to FIGS. In this embodiment, as shown in FIG.
To 360 °, the range is divided into sections corresponding to the resolution of the input value indicating the wind direction instantaneous value, the angle values representing the sections are assigned to the respective sections, and the angle values representing this section are assigned. Refer to the value of the angle difference between the ROM1 that stores the data table 11 that can be accessed with the input value and the angle value accessed from this ROM1 and the exponential average angle value for the input value up to the previous time, and the absolute value equivalent to this angle value. Absolute angle calculation means 2 for calculating the angle difference, and exponential average calculation means 3 for calculating the exponential average value for the absolute angle obtained by this means.
And an azimuth conversion means 4 for converting the exponential average value calculated by this means into an azimuth.

次に、この実施例の動作を説明する。 Next, the operation of this embodiment will be described.

風向瞬間値を検出する感部からの入力値は、エンコー
ダとそれを挟むフォト・カプラとにより2n種類(nは、
フォト・カプラの組数)のコードにディジタル化されて
グレイコードになる。ROM1では、このグレイコードに基
づきデータテーブル11がアクセスされてこの入力値に対
応する角度値が抽出される。絶対角度演算手段2は、前
回までの入力値に対する指数平均角度値と今回の入力値
に対してデータテーブル11から抽出された角度値との角
度差を演算し(ステップS1)、この角度差が180°未満
であるときにプラスフラグをセットし(ステップS5)、
180°以上であるときにマイナスフラグをセットし(ス
テップS3)、プラスフラグがセットされた状態では、角
度差が−180°以下であるときにこの角度差に360°を加
算した値を絶対角度差とし(ステップS7)、また、角度
差が−180°超過であるときにこの角度差を絶対角度差
とし(ステップS8)、また、マイナスフラグをセットさ
れている状態では、360°から角度差を減算した値を絶
対角度差とする(ステップS4)。指数平均演算手段3
は、この絶対角度差を(1)式に代入して指数平均角度
値を算定する(ステップS9)。ここで、この式の右辺の
第1項と第2項とを結合する演算符号±は、プラスフラ
グがセットされているときは符号+が選択され、マイナ
スフラグがセットされているときは符号−が選択され
る。
The input value from the sensing unit that detects the instantaneous value of the wind direction is 2 n types (n is
It is digitized into a code of the number of pairs of photo couplers and becomes a gray code. In the ROM1, the data table 11 is accessed based on this Gray code and the angle value corresponding to this input value is extracted. The absolute angle calculation means 2 calculates the angle difference between the exponential average angle value for the previous input value and the angle value extracted from the data table 11 for the current input value (step S1). If it is less than 180 °, set the plus flag (step S5),
When the angle is 180 ° or more, the minus flag is set (step S3), and when the plus flag is set, when the angle difference is −180 ° or less, the value obtained by adding 360 ° to this angle difference is the absolute angle. Difference (step S7), and when the angle difference exceeds −180 °, this angle difference is used as the absolute angle difference (step S8). Also, with the minus flag set, the angle difference from 360 ° The value obtained by subtracting is the absolute angle difference (step S4). Exponential averaging means 3
Substitutes this absolute angle difference into the equation (1) to calculate the exponential average angle value (step S9). Here, as the operation sign ± that connects the first term and the second term on the right side of this expression, the sign + is selected when the plus flag is set, and the sign − when the minus flag is set. Is selected.

Yn=Yn-1±Yd(1−ε-t/T) …(1) ただし、 Yn:平均値〔°〕 Yn-1:前回までの平均値〔°〕 Yd:絶対角度差〔°〕 ε:自然対数の底(=2.71828…) T:時間定数 t:データ更新間隔 次に、方位変換手段4では、指数平均演算手段3で演
算された指数平均角度値が360°以上であるときにこの
指数平均角度値から360°を減算した値を出力値とし
(ステップS11)、360°未満のときにこの指数平均角度
値を出力値とし、さらに、指数平均角度値が負の値であ
るときにこの指数平均角度値の360°を加算した値を出
力値とし(ステップS12およびS13)、この指数平均角度
値が正の値または0°に等しいときにこの指数平均角度
値を出力値とし(ステップS12)、次に、この出力値に1
1.5°を加算して22.5°で除算して得られる数値に対応
する方位名を16方位から選択する(ステップS15)。
Y n = Y n-1 ± Y d (1-ε -t / T ) (1) where Y n : average value [°] Y n-1 : previous average value [°] Y d : absolute Angle difference [°] ε: Base of natural logarithm (= 2.71828 ...) T: Time constant t: Data update interval Next, in the azimuth conversion means 4, the exponential average angle value calculated by the exponential average calculation means 3 is 360 °. When it is above, the output value is the value obtained by subtracting 360 ° from this exponential average angle value (step S11). When it is less than 360 °, this exponential average angle value is used as the output value, and the exponential average angle value is negative. Value of this exponential average angle value is added as the output value (steps S12 and S13), and when this exponential average angle value is a positive value or 0 °, this exponential average angle value is As the output value (step S12), and then add 1 to this output value.
The azimuth name corresponding to the numerical value obtained by adding 1.5 ° and dividing by 22.5 ° is selected from 16 azimuths (step S15).

この実施例では、16方位名の場合を説明したが32方位
名を含む他の方位名系に適用しても本発明を実施するこ
とができる。
In this embodiment, the case of 16 azimuth names has been described, but the present invention can be implemented even when applied to another azimuth name system including 32 azimuth names.

〔発明の効果〕〔The invention's effect〕

本発明は、以上説明したように、入力値に対して正弦
成分と余弦成分とに分けることなく代表値としての角度
をテーブル化し、この値により指数平均演算を行うの
で、指数平均演算の回数が従来例に比べて1回減り、か
つ、従来例で行っていた時間のかかる逆正接の計算を行
わないので、平均処理時間を大幅に短縮することができ
る効果がある。
As described above, according to the present invention, the angle as a representative value is tabulated without dividing the input value into the sine component and the cosine component, and the exponential averaging operation is performed using this value. This is one time less than in the conventional example, and since the time-consuming arctangent calculation, which is performed in the conventional example, is not performed, the average processing time can be significantly shortened.

【図面の簡単な説明】[Brief description of drawings]

第1図は、本発明実施例の全体構成を示すブロック構成
図。 第2図は、本発明実施例のの処理を示すフローチャー
ト。 1……ROM、2……絶対角度演算手段、3……指数平均
演算手段、4……方位変換手段、11……データテーブ
ル。
FIG. 1 is a block diagram showing the overall configuration of the embodiment of the present invention. FIG. 2 is a flowchart showing the processing of the embodiment of the present invention. 1 ... ROM, 2 ... absolute angle calculation means, 3 ... exponential average calculation means, 4 ... direction conversion means, 11 ... data table.

Claims (1)

【特許請求の範囲】[Claims] 【請求項1】風向瞬間値を示す入力値に基づき風向方位
を決定する指数平均風向計において、 角度0°ないし360°の範囲が上記入力値の分解能に相
応の区分に分割され、この区分を代表する角度値がそれ
ぞれの区分に割り付けられたデータテーブルを格納し、
この区分を代表する角度値が上記入力値でアクセスでき
るROMと、 このROMからアクセスした角度値と前回までの入力値に
対する指数平均角度値との角度差の値を参照してこの角
度値に相当の絶対角度差を演算する絶対角度演算手段
と、 この絶対角度演算手段で得られた絶対角度に対する指数
平均値を演算する指数平均演算手段と、 この指数平均演算手段で演算された指数平均値を風向方
位に変換する方位変換手段と を備えたことを特徴とする指数平均風向計。
1. In an exponential average anemometer for determining a wind direction based on an input value indicating an instantaneous wind direction value, an angle range of 0 ° to 360 ° is divided into sections corresponding to the resolution of the input value. Stores a data table in which representative angle values are assigned to each category,
Corresponding to this angle value by referring to the ROM that can access the angle value that represents this category with the above input value, and the angle difference value between the angle value accessed from this ROM and the exponential average angle value for the previous input value The absolute angle calculation means for calculating the absolute angle difference of, the exponential average calculation means for calculating the exponential average value for the absolute angle obtained by this absolute angle calculation means, and the exponential average value calculated by this exponential average calculation means An exponential average wind vane, comprising: an azimuth converting means for converting the wind direction.
JP2142508A 1990-05-31 1990-05-31 Exponential average wind vane Expired - Lifetime JPH0830712B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2142508A JPH0830712B2 (en) 1990-05-31 1990-05-31 Exponential average wind vane

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2142508A JPH0830712B2 (en) 1990-05-31 1990-05-31 Exponential average wind vane

Publications (2)

Publication Number Publication Date
JPH0434372A JPH0434372A (en) 1992-02-05
JPH0830712B2 true JPH0830712B2 (en) 1996-03-27

Family

ID=15316982

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2142508A Expired - Lifetime JPH0830712B2 (en) 1990-05-31 1990-05-31 Exponential average wind vane

Country Status (1)

Country Link
JP (1) JPH0830712B2 (en)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5650316A (en) * 1979-09-29 1981-05-07 Ricoh Co Ltd Automatic exposure control circuit of camera
JPS60192266A (en) * 1984-03-14 1985-09-30 Fuji Facom Corp Mean bearing measuring system
JPH0233101B2 (en) * 1984-12-06 1990-07-25 Meisei Electric Co Ltd FUKOFUSOKUKEIOYOBISONOENZANSOCHI

Also Published As

Publication number Publication date
JPH0434372A (en) 1992-02-05

Similar Documents

Publication Publication Date Title
JPS5873231A (en) Analog-to-digital converter
CN115659662A (en) Vehicle weight estimation method and device, storage medium and engineering vehicle
JPH0830712B2 (en) Exponential average wind vane
JPH0644290A (en) Method and device for conducting high-speed fourier transform
CN114706342B (en) Oxygen concentration detection method based on nonlinear interpolation compensation algorithm
JPS62218813A (en) Pressure detector
JP3289961B2 (en) Circuit element constant measuring device
JP2880347B2 (en) Structured system processing time estimation device
JPH0434373A (en) Anemoscope
JPH0743267B2 (en) Position detection method by resolver / digital converter
JPS63222268A (en) Rotational speed detector
JP2003057073A (en) Interpolation error estimating method, its device, and interpolation error correcting apparatus using the same
JPH0821867B2 (en) Measuring value converter for measuring instruments
SU1619198A1 (en) Device for measuring modulus of gain ratio of four-terminal networks
JP3121609B2 (en) Digital signal processor
JPH074658Y2 (en) Logarithmic conversion circuit
SU729839A1 (en) Voltage-to-code converter
CN117781846A (en) Angle compensation method, device, equipment and storage medium
JPS58167971A (en) Digital type maximum/minimum detection circuit
CN117830403A (en) A lane diverging point and merging point detection method, device, equipment and medium
JP3078696B2 (en) Reciprocal arithmetic unit
JPH0743566B2 (en) Angle calculator
JPS62191932A (en) Character string processing system
JP3403650B2 (en) Displacement measuring device
KR100493010B1 (en) High accuracy division method for digital signal processing