next up previous
Next: 記述関数法による安定判別 Up: 記述関数法(Describing Function) Previous: 記述関数法(Describing Function)

記述関数(等価伝達関数)

5.3(a)に示すように、非線形要素に正弦波入力を加えた場合、 出力波形は正弦波ではなく形が歪むが、入力と同じ同期を持つ周期的な波形となる。 したがって基本波およびその高周波数成分から構成されていると見なされる。 通常制御系の各種要素は高い周波数ではゲインが減衰する特性を有するから、 非線形要素と結合して使用されるときは高調波成分は減衰され無視できると考える。 そこで出力の基本成分のみを考え、これを(b)図に示すように等価出力とする。 その場合入力と出力は共に同じ周期の正弦波となるから、両者の間の関係は 等価線形化された要素と見なされ、これを記述関数(又は等価伝達関数)という。
図 5.3:
\includegraphics[scale=0.60]{eps/5-2-1.eps}

いま非線形要素に加わる入力を

\begin{displaymath}
v(t)=V_0\sin \omega t
\end{displaymath} (5.18)

の正弦波とし、出力波形をフーリエ級数で表示すると、次のようになる。
\begin{displaymath}
u(t)=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\cos n\omega t+b_n\sin \omega
t)
\end{displaymath} (5.19)

\begin{eqnarray*}
\mbox{ただし} a_n & = & \frac{\displaystyle{2}}{\displaystyle...
...i u(t)\sin n\omega t\mathrm{d}\omega t \\
& & (n=0,1,2,\cdots)
\end{eqnarray*}

出力波形が正負対称であれば$a_0=0$となる。基本波形成分のみを考えると

$\displaystyle u_1(t)$ $\textstyle =$ $\displaystyle a_1\cos \omega t+b_1\sin \omega t$  
  $\textstyle =$ $\displaystyle \sqrt{a_1^2+b_1^2}\sin (\omega t+\varphi )$ (5.20)

ただし $\varphi =\tan ^{-1}(a_1/b_1)$
となる。これが等価出力であるから、式5.18の入力と比較すると等価な 伝達関数、すなわち記述関数$G_n(V_0)$のゲインと位相は
\begin{displaymath}
\left\{ \begin{array}{ccl}
ゲイン & : & \vert G_n(V_0)\vert...
...gle G_n(V_0)=\varphi =\tan ^{-1}(a_1/b_1)
\end{array} \right.
\end{displaymath} (5.21)

で定義される。

[例1]ヒステリシスにある2位置リレー要素

5.4に示すような$h$の幅のヒステリシスを持つ2位置リレー要素に、 正弦波の入力を加えると図5.5に示すように入力が$h$に達したところで オン・オフが切り換り、出力は方形波状になる。

図 5.4:
\includegraphics[scale=0.60]{eps/5-2-2.eps}
図 5.5:
\includegraphics[scale=0.60]{eps/5-2-3.eps}

この基本成分の$a_1,b_1$を(5.19)式より求めると

$\displaystyle a_1$ $\textstyle =$ $\displaystyle \frac{2}{\pi} \left\{ \int_{0}^{\varphi}(-A)\cos \omega t\mathrm{d}\omega t
+\int_{\varphi}^{\pi}A\cos \omega t\cdot \mathrm{d}\omega t \right\}$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi} \{ [-\sin \omega t]_0^{\varphi}
+[\sin \omega t]_{\varphi} ^{\pi} \}$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi} \{ -\sin \varphi -\sin \varphi \}$  
  $\textstyle =$ $\displaystyle -\frac{4A}{\pi}\sin \varphi$ (5.22)
$\displaystyle b_1$ $\textstyle =$ $\displaystyle \frac{2}{\pi} \left\{ \int_{0}^{\varphi}(-A)\sin \omega t\mathrm{d}\omega t
+\int_{\varphi}^{\pi}A\sin \omega t\cdot \mathrm{d}\omega t \right\}$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi} \{ [\cos \omega t]_0^{\varphi}
-[\cos \omega t]_{\varphi} ^{\pi} \}$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi} \{ \cos \varphi -1+1+\cos \varphi \}$  
  $\textstyle =$ $\displaystyle \frac{4A}{\pi}\cos \varphi$ (5.23)

となる。したがって記述関数のゲイン及び位相は次のようになる。
\begin{displaymath}
\vert G_n(V_0)\vert=\sqrt{ \left( \frac{4A}{\pi} \right) ^2\...
...ac{4A}{\pi} \right) ^2\cos ^2\varphi}/V_0
=\frac{4A}{\pi V_0}
\end{displaymath} (5.24)


\begin{displaymath}
\angle G_n(V_0)=\tan ^{-1}
\left\{ \left( \frac{-4A}{\pi}\s...
.../\left( \frac{4A}{\pi}\cos \varphi \right)
\right\} =-\varphi
\end{displaymath} (5.25)

この$\varphi$は図5.5より
\begin{displaymath}
\begin{array}{l}
V_0 \sin \varphi = h\\
\varphi =\sin ^{-1} (h/V_0)
\end{array}\end{displaymath} (5.26)

の関係があることがわかる。

この例からもわかるように記述関数のゲインも位相も入力正弦波の振幅$V_0$の 関数となっている。

[例2]ヒステリシスのある3位置リレー要素

5.6に示すようなヒステリシスのある3位置リレーの場合

$\displaystyle \theta = \cos ^{-1} \left( \frac{d+h}{V_0} \right)$     (5.27)
$\displaystyle \psi = \cos ^{-1} \left( \frac{d-h}{V_0} \right)$     (5.28)

とすると図5.7より
$\displaystyle a_1$ $\textstyle =$ $\displaystyle \frac{2}{\pi}\int_{\frac{\pi}{2}-\theta}^{\frac{\pi}{2}+\psi}
A\cos \omega t\mathrm{d}\omega t$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi} \left[ \sin \left( \frac{\pi}{2}+\psi \right)
-\sin \left( \frac{\pi}{2}-\theta \right) \right]$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi}(\cos \psi -\cos \theta)$ (5.29)
$\displaystyle b_1$ $\textstyle =$ $\displaystyle \frac{2}{\pi}\int_{\frac{\pi}{2}-\theta}^{\frac{\pi}{2}+\psi}
A\sin \omega t\mathrm{d}\omega t$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi} \left[ -\cos \left( \frac{\pi}{2}+\psi \right)
+\cos \left( \frac{\pi}{2}-\theta \right) \right]$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi}(\sin \psi +\sin \theta)$ (5.30)

となる。これより記述関数のゲインと位相を求めると
$\displaystyle \vert G_n(V_0)\vert$ $\textstyle =$ $\displaystyle \frac{2A}{\pi}
\sqrt{ (\cos \psi -\cos \theta)^2+(\sin \psi +\sin \theta)^2 }
/V_0$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi V_0}
\sqrt{ 2 \{ 1-\cos \psi \cos \theta
+\sin \psi \sin \theta \} }$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi V_0}\sqrt{ 2 \{ 1-\cos(\psi +\theta) \} }$  
  $\textstyle =$ $\displaystyle \frac{2A}{\pi V_0}
\sqrt{ 2 \left\{ 2\sin ^2 \left (\frac{\psi +\theta}{2} \right)
\right\} }$  
  $\textstyle =$ $\displaystyle \frac{4A}{\pi V_0}\sin \left( \frac{\psi +\theta}{2} \right)$ (5.31)
$\displaystyle \angle G_n(V_0)$ $\textstyle =$ $\displaystyle \tan ^{-1} \left( \frac{\cos \psi -\cos \theta}{\sin \psi +\sin \theta}
\right)$  
  $\textstyle =$ $\displaystyle \tan ^{-1}
\left( \frac{ -2 \sin \frac{\displaystyle{\psi +\theta...
...aystyle{2}} \cos \frac{\displaystyle{\psi -\theta}}{\displaystyle{2}} }
\right)$  
  $\textstyle =$ $\displaystyle \tan ^{-1} \left( -\tan \frac{\psi -\theta}{2} \right)$  
  $\textstyle =$ $\displaystyle -\frac{\psi -\theta}{2}$ (5.32)

となる。
図 5.6:
\includegraphics[scale=0.60]{eps/5-2-4.eps}
図 5.7:
\includegraphics[scale=0.60]{eps/5-2-5.eps} j

[例3]バックラッシュ(ヒステリシス)要素

5.8に示すようなバックラッシュ要素の場合、正弦波入力に対する 出力波形は図5.9のようになる。

図 5.8:
\includegraphics[scale=0.60]{eps/5-2-6.eps}
図 5.9:
\includegraphics[scale=0.60]{eps/5-2-7.eps}


\begin{displaymath}
\theta =\sin ^{-1} \left( \frac{V_0-2h}{V_0} \right)
\end{displaymath} (5.33)

とした場合、出力は

\begin{displaymath}
\left\{
\begin{array}{lcl}
K(V_0\sin \omega t-h) & \hbox{...
...& \hbox{ } & \pi -\theta<\omega t\leq \pi
\end{array} \right.
\end{displaymath}

となるので
$\displaystyle a_1$ $\textstyle =$ $\displaystyle \frac{2}{\pi}
\left[ \int_0^{\frac{pi}{2}}K(V_0\sin \omega t-h)\cos \omega t\mathrm{d}\omega t
\right.$  
    $\displaystyle +\int_{\frac{\pi }{2}}^{\pi -\theta}K(V_0-h)\cos \omega t\mathrm{d}\omega t$  
    $\displaystyle \left. +\int_{\pi -\theta}^{\pi}K(V_0\sin \omega t+h)\cos \omega t\mathrm{d}\omega t
\right]$  
  $\textstyle =$ $\displaystyle \frac{2K}{\pi}
\left\{ \left[ \frac{V_0}{2}\sin ^2\omega t-h\sin ...
...}}
+ \left[ (V_0-h)\sin \omega t \right] _{\frac{\pi}{2}}^{\pi -\theta}
\right.$  
    $\displaystyle \left. + \left[ \frac{V_0}{2}\sin ^2\omega t+h\sin \omega t \right]
_{\pi -\theta}^{\pi} \right\}$  
  $\textstyle =$ $\displaystyle \frac{2K}{\pi}
\left\{ \left( \frac{V_0}{2}-h \right) +(V_0-h)\sin (\pi -\theta)
-(V_0-h) \right.$  
    $\displaystyle \left.-\frac{V_0}{2}\sin ^2(\pi -\theta)-h\sin (\pi -\theta)
\right\}$  
  $\textstyle =$ $\displaystyle \frac{2K}{\pi}
\left\{ -\frac{V_0}{2}+(V_0-2h)\sin (\pi -\theta)
-\frac{V_0}{2}\sin ^2(\pi -\theta) \right\}$  
  $\textstyle =$ $\displaystyle \frac{2K}{\pi}
\left\{ -\frac{V_0}{2}+\frac{(V_0-2h)^2}{V_0}
-\frac{(V_0-2h)^2}{2V_0} \right\}$  
  $\textstyle =$ $\displaystyle \frac{K}{\pi V_0}(-4V_0h+4h^2)$  
  $\textstyle =$ $\displaystyle \frac{-4Kh}{\pi} \left( 1-\frac{h}{V_0} \right)$ (5.34)
$\displaystyle b_1$ $\textstyle =$ $\displaystyle \frac{2}{\pi}
\left[ \int_0^{\frac{pi}{2}}K(V_0\sin \omega t-h)\sin \omega t\mathrm{d}\omega t
\right.$  
    $\displaystyle +\int_{\frac{\pi}{2}}^{\pi -\theta}K(V_0-h)\sin \omega t\mathrm{d}\omega t$  
    $\displaystyle \left. +\int_{\pi -\theta}^{\pi}K(V_0\sin \omega t+h)\sin \omega td
\omega t \right]$  
  $\textstyle =$ $\displaystyle \frac{2K}{\pi}
\left\{ \left[ \frac{V_0\cdot \omega t}{2}-\frac{V...
...}
- \left[ (V_0-h)\cos \omega t \right] _{\frac{\pi}{2}}
^{\pi -\theta} \right.$  
    $\displaystyle \left. + \left[ \frac{V_0\cdot \omega t}{2}-\frac{V_0}{4}\sin 2\omega t
-h\cos \omega t \right] _{\pi -\theta}^{\pi}
\right\}$  
  $\textstyle =$ $\displaystyle \frac{22K}{\pi}
\left\{ \left( \frac{V_0\pi}{4}-h \right) +(V_0-h...
... -\theta)
\left( \frac{V_0\pi}{2}+h \right) -\frac{V_0(\pi -\theta)}{2}
\right.$  
    $\displaystyle \left. +\frac{V_0}{4}\sin (2\pi -2\theta)+h\cos (2\pi -\theta)
\right\}$  
  $\textstyle =$ $\displaystyle \frac{2K}{\pi}
\left\{ \frac{V_0\pi}{4}+\frac{V_0\theta}{2}+(V_0-h)\cos \theta
-\frac{V_0}{2}\sin \theta \cos \theta -h\cos \theta
\right\}$  
  $\textstyle =$ $\displaystyle \frac{2K}{\pi}
\left\{ \frac{V_0\pi}{4}+\frac{V_0\theta}{2}
+\frac{(2V_0-2h-V_0\sin \theta -2h)\cos \theta}{2}
\right\}$  
  $\textstyle =$ $\displaystyle \frac{2K}{\pi}
\left\{ \frac{V_0\pi}{4}+\frac{V_0\theta}{2}+\frac{(V_0-2h)
\left( \frac{2\sqrt{hV_0-h^2}}{V_0} \right)}{2}
\right\}$  
  $\textstyle =$ $\displaystyle \frac{KV_0}{\pi}
\left\{ \frac{\pi}{2} +\sin ^{-1} \left( 1-\frac...
...\left( 1-\frac{2h}{V_0} \right)
\sqrt{\frac{h}{V_0}-\frac{h^2}{V_0^2}} \right\}$ (5.35)

となり、これを(5.21)式に代入してゲインおよび位相が求められる。

飽和要素

\begin{displaymath}
\begin{array}{ccc}
a_1=0 & , & b_1=\frac{\displaystyle{KV_0}}{\displaystyle{\pi}}(2\alpha +\sin 2\alpha )
\end{array}\end{displaymath}

ただし $\alpha=\sin ^{-1}\frac{\displaystyle{A}}{\displaystyle{KV_0}}=\sin ^{-1}\frac{\displaystyle{S}}{\displaystyle{V_0}}$

不感帯要素

\begin{displaymath}
\begin{array}{ccc}
a_1=0 & , & b_1=\frac{\displaystyle{K}}{\displaystyle{\pi}}(\pi-2\alpha +\sin 2\alpha )
\end{array}
\end{displaymath}

ただし $\alpha=\sin ^{-1}\frac{\displaystyle{d}}{\displaystyle{V_0}}$


next up previous
Next: 記述関数法による安定判別 Up: 記述関数法(Describing Function) Previous: 記述関数法(Describing Function)
Yasunari SHIDAMA
平成15年7月28日