JP5989366B2 - 交流信号測定装置 - Google Patents
交流信号測定装置 Download PDFInfo
- Publication number
- JP5989366B2 JP5989366B2 JP2012062390A JP2012062390A JP5989366B2 JP 5989366 B2 JP5989366 B2 JP 5989366B2 JP 2012062390 A JP2012062390 A JP 2012062390A JP 2012062390 A JP2012062390 A JP 2012062390A JP 5989366 B2 JP5989366 B2 JP 5989366B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- voltage
- current
- rms
- effective value
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 claims description 82
- 238000006243 chemical reaction Methods 0.000 claims description 21
- 230000014509 gene expression Effects 0.000 claims description 21
- 238000000034 method Methods 0.000 claims description 14
- 230000009466 transformation Effects 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 description 45
- 238000010586 diagram Methods 0.000 description 6
- 239000002131 composite material Substances 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 229940050561 matrix product Drugs 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000003990 capacitor Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000001020 plasma etching Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
Images
Landscapes
- Measuring Phase Differences (AREA)
- Measurement Of Current Or Voltage (AREA)
- Drying Of Semiconductors (AREA)
Description
v’[k]=(v[k+1]−v[k-1])/(2・sin(ω))
i’[k]=(i[k+1]−i[k-1])/(2・sin(ω))
但し、v’[k]:v[k]の時間微分値
i’[k]:i[k]の時間微分値
を演算した後、
φ[k]=tan-1[(v[k]・i’[k]-v’[k]・i[k])/(v’[k]・i’[k]+v[k]・i[k]]
を演算して位相φを算出する方法である。
Vrms=√2・√(v1 2+v2 2)=V/√(2)
Irms=√2・√(i1 2+i2 2)=I/√(2)
φ=tan-1[(v1・i2−v2・i1)/(v2・i2+v1・i1)]
の簡単な演算式により交流電圧vと交流電流iの実効値Vrms,Irmsと、位相差φを高速かつ高精度で算出することができる。しかし、従来、これを実現する交流信号測定装置は全く提案されていない。
流の実効値の出力をIとすると、下記の(a5)式乃至(a7)式のいずれかを演算する
ことによって前記交流電圧と前記交流電流の位相差φを演算する位相差演算手段を更に備えるとよい(請求項5)。
Y1(s),Y2(s)は入力y1(t),y2(t)のラプラス変換
上記(6)式の演算により入力信号u1(t)=A・cos(ωo−ψ)+f(ωn・t)の基本波成分(角周波数ωoの成分)と同一の角周波数の余弦波A・cos(ωo−ψ)/2と正弦波A・sin(ωo−ψ)/2を算出するようにしている。
ya(t)=exp(−j・ωo・t)・F(s)・exp(j・ωo・t)・ua(t) …(11a)
yb(t)=exp(j・ωo・t)・F(s)・exp(−j・ωo・t)・ub(t) …(11b)
であるから、(11a)式と(11b)式はそれぞれ図7(a)の信号処理システム8と図7(b)の信号処理システム8'を表している。
[d(X(t))/dt]=[A]・[X(t)]+[B]・U(t)…(12a)
Y(t)=[C]・[X(t)]+D・U(t) …(12b)
但し、Y(t):出力 U(t):入力
[X(t)]:状態変数の行列 d(X(t))/dt:状態変数の時間微分値
[A]:状態係数の行列 [B]:入力係数の行列
[C]:出力係数の行列 D:直達係数
の状態微分方程式で表わされる。説明を簡単にするために、状態変数X(t)がスカラーの一次元システムで初期状態が零(零状態)の場合、(12a)式と(12b)式の出力Y(t)を与える解は、
ξ(t)’=exp(-A・t)・B・exp(j・ωo・t)・ua(t) …(15a)
y(t)=exp(-j・ωo・t)・C・exp(A・t)・ξ(t)+D・ua(t) …(15b)
但し、ξ(t)’=d(ξ(t))/dt
と表現することができる。
M(t)=exp(A・t)・exp(−j・ωo・t)
ζ(t)=M(t)・ξ(t)
の変数変換を行うと、(15a)式と(15b)式は、
ζ(t)’=M(t)’・ξ(t)+M(t)・ξ(t)’
=(A−j・ωo)・exp(A・t)・exp(−j・ωo・t)・ξ(t)+B・ua(t)
=(A−j・ωo)・ζ(t)+B・ua(t) …(16a)
y(t)=C・ζ(t)+D・ua(t) …(16b)
但し、ζ(t)’=d(ζ(t))/dt、M(t)’=d(M(t))/dt
と変形することができる。
ζ(t)’=(A+j・ωo)・ζ(t)+B・ua(t) …(17a)
y(t)=C・ζ(t)+D・ua(t) …(17b)
の状態微分方程式で表わされる数理モデルと考えることができる。
s・Xa(s)=Aa・Xa(s)+B・Ua(s)
Ya(s)=C・Xa(s)+D・Ua(s)
s・Xb(s)=Ab・Xb(s)+B・Ub(s)
Yb(s)=C・Xb(s)+D・Ub(s)
となり、伝達関数Fa(s),Fb(s)は、
Fa(s)=Ya(s)/Ua(s)=C・B/(s-Aa)+D…(19a)
Fb(s)=Yb(s)/Ub(s)=C・B/(s-Ab)+D…(19b)
で表わされる。
Fa(s)=C・B/(s’-A)+D
となり、この伝達関数Fa(s)は、信号処理ブロック6bの伝達関数F(s)のラプラス変数sをs’=s+j・ω0に変更したF(s+j・ω0)に相当している。(19b)式は、(19a)式に対して「ω0」を「−ω0」に変換することによって得られるから、(19b)式から得られる伝達関数Fb(s)は、伝達関数F(s+j・ω0)の「ω0」を「−ω0」に変換することによって得られ、F(s-j・ω0)となる。
Y1(s)=F11(s)・U1(s)+F12(s)・U2(s)
Y2(s)=F21(s)・U1(s)+F22(s)・U2(s)
で表わされ、出力y1(t)は、伝達関数F11(s)のフィルタに入力U1(s)を通した信号y11(t)と伝達関数F12(s)のフィルタに入力U2(s)を通した信号y12(t)の合成波で与えられ、出力y2(t)は、伝達関数F21(s)のフィルタに入力U1(s)を通した信号y21(t)と伝達関数F22(s)のフィルタに入力U2(s)を通した信号y22(t)の合成波で与えられる。
y1(t)=(A/2)・cos(ωo・t−ψ)
y2(t)=(A/2)・sin(ωo・t−ψ)
となり、(4)式に示した図5に示す信号処理部6の出力y1(t),y2と同じ結果が得られる。
Vrms=√(2)・√(ya1(t)2+ya2(t)2) …(20a)
の演算を行うことにより実効値Vrmsを算出する。
Irms=√(2)・√(yb1(t)2+yb2(t)2) …(20b)
の演算を行うことにより実効値Irmsを算出する。
Vrms=√(y1(t)2+y2(t)2) …(21a)
となる。
Irms=√(yb1(t)2+yb2(t)2) …(21b)
となる。
sin(φ)=sin(θ-(θ-φ))=sin(θ)・cos(θ-φ)-cos(θ)・sin(θ-φ)
=(a2・b1−a1・b2)/(A・B) …(22a)
cos(φ)=cos(θ-(θ-φ))=cos(θ)・cos(θ-φ)+ sin(θ)・sin(θ-φ)
=(a1・b1+a2・b2)/(A・B) …(22b)
tan(φ)=sin(φ)/cos(φ)
=(a2・b1−a1・b2)/(a1・b1+a2・b2) …(22c)
の演算式によって位相差φの三角関数値が得られるから、位相差演算部55は、(22a)式乃至(22c)式のいずれかを演算することによって位相差φを算出する。
2 インピーダンス整合装置
3 プラズマチャンバー
4 センサ(交流信号検出手段)
5 電圧/電流測定装置(交流信号測定装置)
51 電圧用のフィルタ部(フィルタ手段)
52 電流用のフィルタ部(フィルタ手段)
53 電圧用の実効値演算部(実効値演算手段)
54 電流用の実効値演算部(実効値演算手段)
55 位相差演算部(位相差演算手段)
6 信号処理部
6a αβ/dq変換部
6b フィルタ処理部
6c dq/αβ変換部
7 信号処理システム
L 伝送線路
Claims (6)
- 伝送線路の所定の測定点における交流信号を検出する交流信号検出手段と、
前記交流信号検出手段で検出された前記交流信号の基本波成分を抽出するフィルタ手段と、
前記交流信号の基本波成分の実効値を算出する実効値演算手段と、
を備えた交流信号測定装置であって、
前記フィルタ手段は、
前記交流信号検出手段で検出された前記交流信号を静止しているαβ座標系のα軸成分u1(t)とし、所定の初期値をβ軸成分u2(t)として両成分u1(t),u2(t)を前記交流信号の角周波数ωoで回転するdq座標系のd軸成分ud(t)とq軸成分uq(t)に変換し、所定のローパスフィルタ特性を有する伝達関数F(s)によってフィルタリング処理をして前記基本波成分のみを抽出した後、前記αβ座標系のα軸成分y1(t)とβ軸成分y2(t)に逆変換する信号処理と等価な下記の(a1)式の演算を行うことによって前記交流信号の基本波成分を抽出し、
前記実効値演算手段は、下記の(a2)式の演算を行うことによって実効値Armsを算出する、
ことを特徴とする交流信号測定装置。
- 前記交流信号検出手段が検出する交流信号は交流電圧若しくは交流電流であり、
前記実効値演算手段は、前記交流電圧若しくは前記交流電流の実効値を演算する、
請求項1に記載の交流信号測定装置。 - 伝送線路の所定の測定点における交流電圧と交流電流を検出する交流信号検出手段と、
前記交流信号検出手段で検出された前記交流電圧と前記交流電流の基本波成分をそれぞれ抽出する一対のフィルタ手段と、
前記交流電圧と前記交流電流の基本波成分の実効値をそれぞれ算出する一対の実効値演算手段と、
を備えた交流信号測定装置であって、
前記一対のフィルタ手段は、
前記交流信号検出手段で検出された前記交流電圧と前記交流電流を、静止しているαβ座標系のα軸成分u1(t)とし、所定の初期値をβ軸成分u2(t)として両成分u1(t),u2(t)を前記交流電圧若しくは前記交流電流の角周波数ωoで回転するdq座標系のd軸成分ud(t)とq軸成分uq(t)にそれぞれ変換し、所定のローパスフィルタ特性を有する伝達関数F(s)によってフィルタリング処理をして前記交流電圧と前記交流電流の基本波成分のみをそれぞれ抽出した後、前記αβ座標系のα軸成分y1(t)とβ軸成分y2(t)に逆変換する信号処理と等価な下記の(a1)式の演算を行うことによって前記交流電圧と前記交流電流の基本波成分をそれぞれ抽出し、
前記実効値演算手段は、下記の(a2)式の演算を行うことによって前記フィルタ手段から出力される前記交流電圧と前記交流電流の実効値Vrms,Irmsとを算出する
ことを特徴とする交流信号測定装置。
- 前記伝達関数F(s)は、1/(s・T+1)(T:時定数)で表わされる、請求項1乃至5のいずれかに記載の交流信号測定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012062390A JP5989366B2 (ja) | 2012-03-19 | 2012-03-19 | 交流信号測定装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012062390A JP5989366B2 (ja) | 2012-03-19 | 2012-03-19 | 交流信号測定装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2013195222A JP2013195222A (ja) | 2013-09-30 |
JP5989366B2 true JP5989366B2 (ja) | 2016-09-07 |
Family
ID=49394353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012062390A Active JP5989366B2 (ja) | 2012-03-19 | 2012-03-19 | 交流信号測定装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5989366B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106597079B (zh) * | 2016-12-29 | 2020-02-04 | 福建奥通迈胜电力科技有限公司 | 一种用于故障指示器正弦波振幅求值的方法 |
CN112710890B (zh) * | 2020-12-15 | 2024-02-27 | 特变电工西安电气科技有限公司 | 一种单相正弦交流相量实时计算方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3372177B2 (ja) * | 1996-12-04 | 2003-01-27 | 株式会社荏原製作所 | 直交2軸信号用フィルタ回路 |
JP5495172B2 (ja) * | 2009-09-30 | 2014-05-21 | 株式会社Ihi | 磁気軸受制御装置と方法 |
-
2012
- 2012-03-19 JP JP2012062390A patent/JP5989366B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2013195222A (ja) | 2013-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5762622B2 (ja) | 角度検出装置 | |
Duda et al. | Harmonic phasor estimation with flat-top FIR filter | |
US10584953B2 (en) | Angle detection device and electric power steering device | |
JP5989366B2 (ja) | 交流信号測定装置 | |
CN103353548B (zh) | 电网电压同步信号提取装置及方法 | |
CN104852662A (zh) | 永磁同步电机静态电感参数的测量方法及系统 | |
KR102596568B1 (ko) | 전동기의 회전자 각도 추정 장치 및 방법 | |
Li et al. | A method to remove odd harmonic interferences in square wave reference digital lock-in amplifier | |
Petrović et al. | Computational effective modified Newton–Raphson algorithm for power harmonics parameters estimation | |
JP2006300746A (ja) | 信号処理装置 | |
JP3805718B2 (ja) | 電力系統周波数変動発生検出方法、電力系統周波数偏差計測方法及び電力系統周波数計測方法 | |
JP4566906B2 (ja) | 回転機の定数測定装置 | |
JP6098827B2 (ja) | 永久磁石形同期電動機の制御装置 | |
US20200295693A1 (en) | Estimation of the amplitude of a periodic component in a measured signal through a delta-sigma modulator | |
JP2016201923A (ja) | モータの回転位置推定方法およびモータの制御装置 | |
Wu et al. | Effect of adding DC‐offset estimation integrators in there‐phase enhanced phase‐locked loop on dynamic performance and alternative scheme | |
Abramovitch | Low latency demodulation for atomic force microscopes, Part II: Efficient calculation of magnitude and phase | |
US10585126B2 (en) | Apparatus and method to measure AC power | |
CN113131820B (zh) | 一种永磁同步电机周期性误差的分析方法及装置 | |
JP5989327B2 (ja) | 不平衡率検出装置、および不平衡率検出方法 | |
CN105634362B (zh) | 基于三参数估计的异步电机参数辨识方法 | |
Chen et al. | Implementation of power measurement system with Fourier series and zero-crossing algorithm | |
KR101538738B1 (ko) | 삼단 이산 푸리에 변환을 이용한 주파수 측정 장치 및 방법 | |
Możdżyński | Simple digital integration algorithm with saturation and drift elimination based Second-Order Generalized Integrator | |
CN111884459A (zh) | 一种高精度高带宽的信号处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20150204 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20151215 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20151216 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160203 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20160203 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160628 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160713 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20160802 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160810 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5989366 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |