付録 A3 調和解析器の機構と原理
調和解析器の作動原理とフーリエ係数作動原理、解析手法の詳細を記述するために、報告書の図 3.３を最初に再掲する。

レール

移動系Ａ（赤）歯車

移動系 Aa（青）図 1．報告書の図 3.3（再掲）
この機器の核心は、波形追跡ポインタＦ点の動きと、⑤歯車の回転（それは、⑦プラニメータのトレース点がどの様に動くか、という事だが）の関係を知る事にある。以下、この点に注目すると理解し易い。
図１に示した各指定点の記号を参考にして以下の原理の説明文をお読みください。
・①：移動系 A が動く固定レールで、解析すべき変動波形の OX 軸と直交して設定する。
・②：赤色台車、移動系 A ・③：水色台車、移動系 A 上を動く移動系 Aa
・④：梃子 FKS・⑤：半径 R の歯車で②移動系Ａに搭載
・計測波形の長さ a＝ｍ(=)となる様に m を調整した。・
1= = 20

１．波形の微少面積素と歯車⑤の動き（プラニメータのトレース点の動き）の関係
１）ポインタ F 点が上下（y 方向）だけに動いた場合、②移動系 A、③移動系 Aa は一緒に動くので歯車⑤（半径 R）の回転角 は変らない。プラニメータのトレース点P をセットする穴がある歯車⑤は、両移動系と共に一緒に上下するだけである。
歯車半径 R と歯車中心からトレース点P の穴までの距離 r（R>r）は少しだけ違う。
２）F 点が左右（x 方向）だけに動いた場合、K 点は x=a/2 の線上を動きながら回転する。
この K 点の回転が点 S を通じたスコッチヨーク機構により③移動系 Aa を動かす。その動きにより、歯車⑤が回転する。
即ち、F 点を x 方向に動かす場合，K 点は y 方向にも上下するが、歯車の回転 を考える場合は K 点を固定して回転運動のみを考えて良い。
３）F 点が x 方向に x 動いたときの歯車⑤の回転角 は、下記(9)式に示すが=となる．なお、 = 0で = 0となるように初期にセットする．

(1)４）プラニメータトレース点P は歯車の回転 により、下記の距離だけ動く。
上下方向：r ∙ sin； 左右方向：r(1 − cos )(2)５）計測波形を追跡するポインタ F 点の動きを次のようにした場合を考える。
イ 或る x に於いて関数値 y(x)だけ上方に動き： ○ロ 次に右に dx 動き○
ロ 次に下方に y(x)だけ動き○：

ニ 次に左に dx 動き、最初に戻る。
○

この時、プラニメータのトレース点 Ps と歯車の回転はどう動くか考える。図２を参照。
ロ○yイ○

y = −yハ○

ニ○

∆イ○dy

Pdx

図２．波形ポインタ F 点が微少面積素（y(x)dx）を追跡した時、歯車⑤の Ps 点の動き
歯車の中心を D とする。
2

イ の操作で点 Ps は、y 方向上に y○だけ上がる。しかし、歯車⑤は回転しない。
６）○イロ の操作で歯車の回転角 の増分は、(1)式より、∆○ロ○=で与えられる。

ハ の操作で点 Ps は、y 方向に y○○イだけ下がる。しかし、歯車⑤は回転しない。
二 で Ps は、上下dy○ニ○= −dy：

ロ○

左右dx

ニ
○

= −dx

ロ
○

だけ動いて元に戻る。

即ち、このポインタ F 点の微少部分の動きから歯車⑤の Ps の x、y 方向の変化は(2)式より次式で与えられえる。
a) ∆β による P の上下増分
dyロ○= r sin β + ∆− sin β ≒ r cos β sin ∆ロ○ロ○≒ r cos β ∙ ∆= r cos β ∙(3)b) ∆β による P の左右増分
dxロ○ロ○= −r cos β + ∆= r sin β ∙ロ○

− cos β ≒ +r sin β ∙ sin ∆ロ○≒ r sin β ∙ ∆ロ○(4)７）結局、ポインタ F 点が、 ( )

イ○
ロ○
ハ○
二のを一周した時、トレース点P は平行四辺形○領域を一巡する。この平行四辺形の面積は、下記である。図２を参照。
( ) ∙ ∙ sin∙=( )∙∙∙ sin ( ) ∙(5)

機器が与えられれば下線部 1,2 は決まる。歯車、梃子 FKS の設定が決まれば定数になるから、それを C とする。別言すると、この事は歯車の径を決める有効な情報となる。
結論を記すと、波形を追跡するポインタ F 点が x の位置で y(x)dx の微少面積素を一巡した時プラニメータのトレース点P は、(1)式を ( )に代入した式で表現すれば次式になる。
C∙ ( )( )∙= C ∙ ( ) sindx(6)即ち、プラニメータのトレース点P はこの図形を一巡する．トレース点 Ps と直交している
トレース点 Pc（DPs⊥DPc）の座標は、 →C∙ ( )( )∙+= C ∙ ( ) cos

とおいて次式になる。
dx(7)

両式を積分表示すれば、次節で詳説されるフーリエ級数を求める(12)、(13)式と同じである事が解る。この計測波形の微少面積素とプラニメータの動きの関係の考察を基にして、次節を読むと良く理解できるであろう。

3

２．作動原理の詳細
和 田 英 一 東 大 名 誉 教 授 （ IIJ 技 術 研 究 所 所 長 ） の ホ ー ム ペ ー ジ
http://parametron.blogspot.jp/2013/07/blog-post.html に「面積計を使う調和解析器」というタイトルで本計測器の作動原理が解説されている。以下ではその内容を抜粋し、詳細に説明を加えて紹介する。
面積計（プラニメータ）以外の部分（調和解析器と解析される変動波形）を以下の図３のように表す。プラニメータは報告書の３．１を参照してください。
１）この図は解析する変動波形を最初に設定した状況を示し、この図を使って解析する各種数値の説明をする。

・上下方向に動く物をガイドするレール（図１の①）
ポインタ F 点が上下移動すると台車②移動系 A 全体がこのレール上を上下移動する。
・F 点の x 方向の動きのみが台車③移動系 Aa（ZZ’）
を上下させ、歯車を回転させる
・F 点の y 方向の動きは台車②の先端にある K点を上下に動かす。その時、台車③移動系 Aaも一緒に動くので、歯車は回転しない、∙点

R( , )

= ∙

( 2)

==ℎyx
図３．変動波形 f(x)を解析する時に初期設定された状態図
図 1 を更に模式的に書いた図。ポインタ F 点が原点にある状況図である。

4

①上記波形曲線が解析すべき周期関数で、x=0 から x=a までが 1 周期である。その区間で一周する歯車を基本周期の歯車とする。n 倍周期成分はその区間で n 回転（n=1,2・・）する歯車に取り替えて解析する。
②上図 ZZ’と連結されている S 点はＫ点の回転により図１の台車②移動系 A で上下に動く。
③台車②移動系 A の右の腕の先に K という軸があり、梃子 FKS がその軸で回転できる。
ここで、= m である。この時の角度を とする。これは初期設定で決まる。

④最初、下記のように梃子の状態を設定する。
⊥の状態で X 軸の中点を F 点と一致させる。この時、図３に示す初期設定時、K の x 座標は=

すると、⑤梃子の下腕の長さ= /2、y 座標は、ℎ ==∥ ①レールである。
−2 となる。

2 / が求められる。
は、解析する波形の長さに応じてポインタ F 点が可動である。そのポインタ F 点で波形を解析者が追跡する。と直行する梃子の腕 KS の長さ の先にはローラーS があり、その上下で左の歯軌条⑥ZZ'（青色の移動系）が一緒に上下する。
⑦台車②移動系 A には点 D を中心とする半径 R の歯車⑤がある。それは、台車②に対する台車③の、即ち ZZ'の相対的上下移動に従って回転する。よって K が上下だけをしても歯車⑤は回転することはない。K 点が回転する時だけ歯車⑤は回転する。この事はこのシステムの重要な点である。出発位置での D の座標を(b,g)とする。
⑧この歯車には、中心 D から r だけ離れた左と上に穴 Ps と Pc があって、それに面積計⑦のトレース点（針）を接続する。直線 DPs と直線 DPc は直交している。これは、フーリエ係数の sin 成分と cos 成分を求めるために直交している。
⑩ポインタ F 点が波形 f(x)を追跡し、一周して来ると台車や梃子が動き、歯車も上下に移動しながら回転し、その歯車の 1 点の描く面積を面積計で計測する。
⑪半径 R の歯車はポインタ F 点が x=0 から x=a まで変動波形を追跡移動した時に、丁度n(Fourier 係数の次数)回転するように出来ている。
⑫次数 n に従って歯車の半径 R は変わり、別の歯車が用意されている。それぞれの歯車に対応して異なる歯車回転軸用の穴が用意されている。
（付録 A2 調和解析器写真集参照）
２）上述の操作によってフーリエ係数が求められる理由
ポインタ F 点が曲線 f(x)の(x,y)まで移動した時の様子を図４に示す。ｘだけ動くと K 点は上下とともに回転する。この時、x=a/2 の点を M とし、∠FKM =①この時の K 座標は(下記である。
∙とする。

, y + h)である。K 点から見て S 点が最初の点から下がった変位は− ∙=12−2−=(8)これが台車②に対する ZZ’の相対的上下移動である。故に、点 D を軸とする半径 R の歯車
5の回転角 β は、下記で与えられる。この式が(1)式で使われている。
=→=(9)

出発位置での歯車中心 D の座標を原点(0,0)にすると、原点に対して Ps の座標（X,Y）は、(− ∙, +ℎ−ℎ + ∙)= − ∙( ∙, +ℎ−ℎ + ∙=, +ℎ−ℎ + ∙となる。Ps と直行している Pc の座標（X,Y）は、 →となる。

∙+とおいて、−)=(10 ), +ℎ−ℎ + ∙(11 )・S の上下変位∙(K 点から見た上下変位

θ∙=点 K の座標: (θ, + ℎ)

初期セッティングにより常に=2= 0で m=a の様に KF の長さを調整した。

=h=−−M= /2

図４．変動波形追跡するポインタＦ点が、点(x,y)にある時の状況図

6

Pc の軌跡による（取り囲んだ）面積==+は、（10）式より以下となる。

|=

∙

(12 )

下線部で示した式は、解析するグラフ上でポインタ F 点を y=0 に沿って x=a から x=0 にもどす操作を表現している。最後にこの式に相当する操作が必要である。
同様に、Ps 点の軌跡による（取り囲んだ）面積=や=+|は、次式で得る事ができる。

=∙(13 )の値は面積計（インテグレーター）により計測される。

さて、区間[0〜a]で n 周期の波の波数を kn とすると、= (2 / ) である。n=1 の時は、区間[0〜a]内に一波長分が収まる場合である。n の場合は n 波長分が区間内に収まる場合である。
上式の角度部/が n(2 / )になる様に歯車の半径 R を選定すれば、n=1,2,3…の周期成分と取り出すための歯車の半径を次式で決められる。
1 11=22=1( 14 )即ち、ポインタ F 点が x=a まで動いたときに n 回転する歯車を用意する。即ち、解析する波形区間[0〜a]に含まれる n 周期成分を取り出す歯車を n 次用歯車と言うと、
n=1∶

=1;
2=2∶

=

1 1;
2 2で与えられる。ここで、本解析で使われた数値、データ長さ:= 20cm,=3 ∶m = 20cm,を代入してn = 1,2,3,∙∙∙ に対応する歯車半径=1 1∶ ∗∗∗∗∗∗
2 3= 20,(15 )を決める事ができる。 を決めた時、m とが比で与えられているから、長いデータ（ が大きい）であってもそれに対応して、m を大きくすれば、歯車の径は同じでも良いことになる。
（12）式、（13）式を、歯車径を n 次用歯車の径として表示した(14)式を使うと

==

∙

=
∙

=2∙

2∙722=( 16)=(17 )

となり、これらからフーリエ係数=,が下式で求められる。
：

面積計（インテグレータ）により計測された=,(18 )の値と、用いた 次歯車の中心点 D と穴迄の距離 によりフーリエ係数が上式で求められる。

ここで、歯車の中心点 D と Ps 点、Pc 点の間の距離 は、下記である。
１次用= 3.2：2 次用= 1.6：3 次用

8= 1.08