• 検索結果がありません。

最適施設配置問題の社会情報学的考察

N/A
N/A
Protected

Academic year: 2021

シェア "最適施設配置問題の社会情報学的考察"

Copied!
20
0
0

読み込み中.... (全文を見る)

全文

(1)

最適施設配置問題の社会情報学的考察

若 林 信 夫

目 次 1.はじめに

2.

単純化 された最適施設配置

3.

よ り一般的な定式化

4.

施設利用のための利用圏

5.Voronoi

図 と

Delaunay

三角形図

6.Voronoi

図を求めるソフ トウェア

7.

結語的覚書

1 は じめに

公共施設であれ,民間施設であれ,利用者 と設置者の観点か ら,施設を最適 に配置す ることの問題 は,古代 よ り長 い歴史がある

最近 にな って,オペ レー ションズ ・リサーチや計算幾何学の分野で再 び,脚光を浴 びている。問題 自体 は,社会的重要性 もさることなが ら,その問題解決面 はコンピュータ処理,あ るいは,情報科学の飛躍的な発展が基礎 にある。

公共施設,例 えば,公園,学校,病院,郵便局 な どは,「よい もの」で,刺 用者 にとって最適 とは,「 近 い」 こと一 に同義で あ り,他方, ごみ処理場や汚水 処理施設 は,一般住民 にとってで きるだけ遠 くにあることを望む。

民間施設の例 として,店舗, レス トラ ン,配送セ ンターなどは,公共施設 よ りも制約条件が厳 しいので,整数型の線形計画問題 として古 くか ら定式化 され て きた。

計画者 は,既設の施設を第一の優先度で考え,新設の施設を第二義的に考え る傾向がある。投票型の最適施設配置や,ゲーム論的な最適施設配置 も考え ら

49〕

(2)

50

討 究 第 46巻 第

1

号 れ る。

本稿では,施設の最適配置の基礎的な考察を通 して,問題 に対す る社会情報 学的な分析を行 うことを 目的にす る。すなわち,管理科学的,情報科学的に問 J 題解決 しよ うとす る

1)0

2 単純化 された最適施設配置

公共施設 の 目的 は平面上 の座標位置,Xl

‑(x

l ,

y

.

)

,..

.,X n‑ (xn,yn)

にある需要点 ( 住民)をサー ビスす ることである.一般住民の合意が得 られる ような場所 に,公共施設を配置 しようとす るとき,単純化 して, 「 便利 さ」の 概念 として,距離 ( 近 さ)を考える。最適性の基準 としては,以下で説明す る 三種鞍がある。

2. 1 3

種類の最適性基準

望 ま しい公共施設を需要点 までの距離の平方和を最小 にす るよ うにす るに は, どこに配置すべ きであろ うか。 これは微分の教科書 にある標準的な問題で ある。 しかし ,なぜ,距離の平方和を最小化 しなければな らないかの理由は存 在 しない。 これは,む しろ,数学的な便法 による。

この問題を数学的に解 くことは

n n

f( X)

∑ [(

i

‑1 x‑

xi)2+b,yi)2]1

∑ l ‑ 1 x‑Xi 1 2

を最小化す る点 P s q , を見出す ことである。 ここで,X

‑(

x , y) である。また, f X‑Xi l はXか ら Xi までのユーク リッ ド距離である。 これを 「 平方和問題」

1

)最適施設配置の数理を考え る上で,重要 な文献 は,以下の

3

つであろ う。

岡部篤行 ・鈴木敦夫 , 「 最適配置の数理」( 朝倉書店)

,1992.

A.

Okade,B.Bootsand K.Sugihara,SpatialTessellations:Concepts andApplicationsofVoroTWiDiagraTn

S

,(JohnWiley

&

Sons),1992.

伊理正夫 ( 監修)・腰塚武志 ( 編集) , 「 計算幾何学 と地理情報処理 ( 第

2

)

」( 共

立出版)

,1993.

(3)

最適施設配置問題の社会情報学的考察

51

( P s q , ) と呼ぼう

.

距離の平方和は数学的に取 り扱いが容易であるが,距離の合計を最小に した 方が直観 にかなっている。距離の和を最小化 した場合,答えは変わって くるだ ろうかを次に考えてみる。距離の和を最小化する問題 は,発電所か ら町までの 電線の総延長を最小化 したいとき,あるいは,情報処理セ ンターか ら各利用サ イ トまで構内回線を最小ケーブル長で引き回 したいとき, 自然な問題設定 とな る。目的関数,

72

f

( X

)‑ ∑ lX‑XiE

(I‑I

がその最小値に達する点

Psum

を求める問題のことを 「 距離和問題

」 (Ps

u m) と 呼ぼ う ?)。 距離和問題の数学は微分の教科書 には出て こな

い 。

む しろ,幾何 学の教科書に出ている。

問題の別の自然な変形は,消防施設や, ごみ収集施設のような,施設か ら最 も遠 い地域 までの距離をできるだけ小 さくす るよ うな配置を求め ることであ る

この問題を 「 最大距離問題

」 (Pmax)

と呼ぼ う 3)。

f( X)‑maxiL X ‑X i F

の最小値を与える点

Pmax

にその解がある

。Pmax

もまた興味深い. 幾何的な特徴 を持 っている。

最小化すべき距離の関数が,距離の平方和,和,あるいは,最大値のどれを 選択すべきか とは別に,距離の定義それ自身を変更す ることができる。

Xl‑(x

l

,yl)

か ら

x2‑(x2,y2)

までの 「タクシー距離

」 4)

は,

2)1927

年の

A.Weber

の業績に肖って

,「Weber

問題」 とも呼ばれ る

。A.Weber

,

UberdenStanportde7・Industrie

' L

,Tubingen,1927.(TheTheory ofthe LocationofIndustries

,に再録。)

3

) この基準 は,弱者優先の社会正義のロールズ ( ㌔.

Rawls)基準 に相当す る。筆者

は,総当た り戦のゲームにおいて,待ち時間を最小 にす る試合順序が作成で きるこ とを最大距離問題の解を通 じて,知 ることができた。

4)

タクシー距離 は,マ ンハ ッタン距離,矩形距離などと呼ばれることがある.

(4)

52

商 学 討 究 第

46

巻 第

1

JXl‑x2lt

Fx l‑ x2r+

l y

l‑y2E

によって与え られ る

都市の街路 に見 られ るよ うに,横 ( 水平)方向か,たて ( 垂直)方 向にのみ,移動す ることがで きる場合, このよ うに計算で きる。

距離和 問題

(Psum)

と最大距離 問題

(Pma,)

は,ユー ク リッ ド距離が タク シー距離 によ って置 き換え られ るとき,それぞれ,「タクシー距離和問題」 と

「タクシー最大距離問題」 と呼び,それぞれの解を,

Qsum

Qmax

と表わそ う

.

/ ー

2. 2 1

次元問題 とその解

X l ,‥. ,

Xn

はすべて

x

事 由上 に位置 してい るとしよ う

その とき,平方和, 距離和,および,最大距離問題 は,それぞれ,

E l xi‑Cf 2 ,

zl

=l

n

∑ I xt ・‑cl ,

H

F e

および,

maxtlxi‑ CI

を最小化す るよ うな数値

C

を求めることである。

平方和 問題 に対す る解

,Psq

,は,標準 的な方法 ( 凸関数 ゆえ,微分

‑ 0)

に よって

n (1/a)∑xll i

であるとわか る。 これ は又,標本

x l, X2

,.‥,Xn の標本平均x ‑と解釈で きる。

距離和問題 については,最初,次の関係式が成立す ることに注意 しよ う

TZ

f(C)∑ I1xr CIx∑ (.<cC‑xi)+ ∑ (I.>cx i‑C)

もしすべての i について,

C

≠xi であれば,

f'(C)‑L‑R

である。 ここで,

L

は, Cの左方 向にある点の個数であ り,

R

は, Cの右方向にある点の個数であ

(5)

最 適施設配 置 問題 の社 会情 報学 的考 察

53

る。 もし,

L‑R

となる点

x i

があるな らば,すなわち,左側 と右側 に等 しい 個数の点があれば,

f

, xi

の左側で減少 し,右側で増加す る。最小点

Psum

は,

x

z で達成 され るだろう

すなわち

,x.,x 2

,...

,Xn

の中央値 ( 中位数, メディア ン,Medi

an)である。 さもない ときには,解 は一意ではな く,最小

値 は,開区間

(x i,Xi)

のどれかの点で達成 され る。そ こでは,開区間の任意 の x に対 し

Tf'(x)‑L‑R‑0

である。

最大距離問題 に対 しては, m ‑T ni n ( x i ) , M ‑m

d.x (x2L

) とお こう

その とき,

E T n‑Cローl M ‑cl ≧M 一m が成立す る。

したが って,m

a x

zlx i‑

C

l≧(M ‑m)/2

である

等号 は

,C‑(m+M)/2

のとき達成 され るので

pmax

‑( T n+M )

/2

である.

3

つの問題に対す る解 は,すべて統計学でよ く知 られている

。P

sq ,は,デー タ集合の平均,

Psum

は,中央値,

Pmax

は,中間範囲

(Mid‑range)である。

n ‑3

の場合には, もしも,

x1‑‑

i

,x2‑ 1,x3‑a for‑1<

a

< l

であるな らば

,P

s q ,( 標本平均)は ,

a/3

,

Ps

u m ( 中央値)は, a,

Pma

, ( 中間 範 囲) は,

0

で あ る。点 が x 軸上 にあ るとき,

Psum‑Qsum

お よび,

Pmax

Qmax

である

2. 3 2

次元の場合の解

距離の平方和問題の解 は,偏微分を考え ると,各座標の平均値

(x

‑, y l であ

り , 「セ ン トロイ ド」 と呼ばれ る。

距離和問題の解 は,フェルマー(

Fermat)

・シュタイクー(

Steiner)問題 :

「 与え られた鋭角三角形

ABC

内に

1

P

を とり,

P

か ら各頂点

A

,

B

,

C

にいたる距離の和を最小 にせよ。 」

の解である。

(6)

54

商 学 討 究 第

46

巻 第 1号

3

点が鋭角三角形を作 るとき

,120

度を仰 ぐ内部の点が解である。鈍角三角 形な らば鈍角の頂点に一致す る 5)0

上述 した ものは,標本点が 3個の場合であり, 4個以上になると幾何的な解 法は困難になる

作図 して容易に確かめ られるように,解 も一意には定まらな

い。

解析的には,f( X)‑ ∑写 ‑ 1 l X ‑X l l を最小化す る問題であるか ら,

∇f‑

0 となる点p‑X ‑ があれば,

1 p‑X i

7 =l l p‑X z ・ E

‑ 0.

を満たさなければな らない。 これは,

2

次元の中央値

(2‑Median)

と解釈 できる。

最後 に,最大距離問題の解を特長付けよう

点 f X . ,

x 2

,‥.

;x n)

の集合を

S

で表わ し,f( X)‑m

axilX‑Xi

l とす る。

もし

,mini(X)f(Pmax)

‑ r ,であれば,その時,すべての tについて, I X

i‑Pmaxl

≦ r であ り,特定の iについて等号が成立す る。 したが って,

S

は,

Pmax

に中心を置 く半径 r の円の内部又は周辺 にあるO このよ うな円を,

S

の 「 最小スバニ ング円」 といい,一意に定まる

Pmax

を求めるには,

2

つのケ‑スを考え る。第一 は, Sを含む円の直径の 反 対 側 に,

2

点X壬とX,が あ る とき, r‑ l X i ‑X,E /2 , か っ,

Pmax

( X

z

+Xj )/2 であるO

第 2は,上のよ うな 2点が見 出せないとき,最小 スバニ ング円上 に 3点が

あって, これ ら

3

点を含む弧が

180

度以上である。 したが って, これ ら

3

点を

5)

フェルマー ・シュタイナー点を求める幾何学的な方法は

2

種類ある。罪‑ は,任意

の三角形の各辺. の上 に,正三角形を外側 に作 ると,対頂の線分の長 さが等 しく,

1

点で交わる。第

2

は, トリチェリの作図法で,任意の三角形の

1

辺を選び,その辺

を 1辺 とす る正三角形を外側に作 る。外側の頂点 と,任意の三角形の残 りの頂点 と

を結ぶ線分 と外側正三角形の外接円の交点が求めるフェルマー ・シュタイナー点 と

なる。

(7)

最適施設配置問題 の社会情報学的考察

55

頂点 とす る三角形 は,鋭角三角形で あ り,

P max

は,外接 円の中心 ( 外心)で 三角形の内部 にある。

例題 と して,需要者 が,

3

点 の位 置,Xl ‑ (

0

,

2

),

x 2‑(‑ 1

,

0)

,

X3‑ (

1

,

0)にいるとき,公共施設を どこに置いた らよいかを,上記の 3 つ の基準 によ り,求めてみると以下のよ うになる6

) 。

図 1.

3

つの異 なる基準 に対応 した最適配置

さまざまな配置が得 られ るが, この ことは,施設配置が,立案者の価値判断 によって どのよ うにで も決 まることを見 る,社会計画論の絶好 の教材 となる。

6)B.Eisenberg andS.Khabbaz

,

"OptimalLocations",CollegeMathe‑

TnaticsJozL,rnal,Vol.22,no.4 (1992),pp.278‑289.参照O

(8)

56

46

巻 第

1

なお,上記の問題の拡張 として,以下の ことが考え られる。

1.上の問題は,公共施設を 1個求める問題であるが,複数個の施設を求め るように一般化すること。 これは,次節以降で述べる。

2.

タク. シー距離を採用 したとき,最適配置を求めことも同様に,可能であ る。

3.

直線 ( ユークリッ ド)距離 とタクシー距離を混在 した場合の計算 も今後 の課題である。

4

.点 と点の間の距離ではな く,点 と直線の距離を もとに直線の最適配置を 求めること。 (これについては計算幾何学の文献に詳 しい

7)。)

ここで,点 と直線の距離の公式は,

DE(p

i

,H)‑ bo+∑;1bjX"・l

V∑∴

である

ここで,点 は

,p

t

‑ (xli

, .‥

, Xdi)

,i

‑ 1

, ...

n

であり,直線 ( 超 平面)は,

H :b。+blX.+…

+bd X

d‑0

である。

3

よ り一般的な定式化

施設の配置の計画立案者は,何 らかの評価基準を設定 した上で,その目的に 照 らした最適解を求める。公共施設の建設の場合,何 もないところに新設する ことはほとんどな く,多 くは既存の施設 とのか らみにおいて,新設することが 多 い。

このような観点を取 り入れて,問題を数理計画問題 に定式化 しよう

既存の施設 iと新設予定の施設 j があって,それぞれ, 座標位置

(x i,

y f ) , ( aj , b J ) を もつ とす る。任意の 2点に対す る距離尺度 として,タクシー距

7)

例 え ば,

NikolaiM.Korneenko and HorstMartini

,

Hyperplane AP‑

Proximation and Related Topics

, J

anos Pach (ed.)New Trendsin

Discreteand CompzLtationalGeoTnetr

) ′

,Springer‑Verlag,1993,pp.135‑

161

.

(9)

最適施設配置問題の社会情報学的考察 57 離を用いる。

d

(X

z , P , I) ‑l

xi

‑ a j I +F y

i

‑b j I

新規予定の施設 は,

T

n個あるもの とす る

U‑1,2

,…

,m)

.既存の施設 は, n 個ある ( i ‑ 1

,2

, ‥. ,a) .多数の顧客 ( 既存の施設)をサー ビスす るため に,新施設の場所を決めたいが,その目的関数 は,新施設 と既存施設の間の距 離の合計を最小にす ることであるとす る。 しか し,既存施設 は既得権益を もつ ので,加重値を導入す るO加重値 w t ・ は,既得施設のサー ビス需要の相対的割 合 とす る。問題 は,次の型の数理計画問題 になる

8)。

Minimizef∑T1∑71Zi,Mid(Xi,P,)

subjectto

m

∑ z ij主 1

,(

i‑ 1,2

, .‥ ,a)

i‑I

ここで,既存施設 i は, どれかの新施設 j kよりサー ビスされ るので,上記 の制約条件がついている。すべての zt . , は,

0

又 は

1

である。

問題 は,前節の

Qsum

に属す るが,幾何 的な方法で解 くことは不可能で,数 理計画 ソフ ト (コンピュータ)によ らざるを得ない。

最近,有力な解法 として注 目されているのは,内点法 ( 内部経路法) と焼 き なま し法 (シ ミュレーテ ッ ド・アニー リング法)である。詳細 については,数 理計画法 の標準 的 なテキス トの新版 や

NetNews

のニ ュース グル ープ

sci.

op‑search

の議論 に譲 る。

8)

大山達雄 , 「 最適化モデル分析 (日科技連)

,1993.

は,種々の型の最適配置モデル を提示 している。

「 負の」施設の配置分析の比較については

,E.ErkutandS.Nellman,Com‑

parisonofFourModelsforDispersingFacilities,Infor,Vol.29,no.2 (199

1

)pp.6886

がある。

(10)

58

商 学 討 究 第

46

巻 第

1

4.

施設利用のための利用圏

施設配置の問題 は,施設の計画立案者や設置者の観点 と,利用者の観点 とは 異なった見方をす る。利用者 は,既存の施設を所与 として,で きるだけ効率的 な利用方法を考え,しか る後 に, 施設の置かれている状況を改善 しようとす る。

利用者が施設を効率的に利用 しよ うとす る場合,利用圏,あるいは,経済活 動圏を構成す ることがで きる。商業活動な らば,商圏である。

施設Pi の利用圏 とは,利用者 にとって,Pi の方が他の施設 P , . よ りもある評 価の下で,た とえば,距離的に近いとか,便利であるとか,価値が大 きい等の 理 由で選ばれ る点か らなる領域の ことである。

ここでは,評価 として,距離の最小化を踏襲す る。

ある利用者が,現在 いる地点を Cとして,利用す る施設をPi ,その距離を

d(C,Pi

)とす ると, 他 の施設 よ りも近 い とい う条件 は

,d(C,Pi)≧d(C,P

,

)

,

j

≠i と書 けるか ら,

v l

‑ i cF d( C, Pt

)≦d(C,P

,

)

,j

≠貢

となる領域を定めることがで きる。

v

i は,Pt の

Voronoi(

ポロノイ)領域 といい,すべての点 についての図 を

Voronoi

図 といい,次節で詳 しく述べ る。

す べ て の施 設 は, 同質 で はな く,差 別化 され て い るので,重 み付 きの

voronoi

領域,重み付 き

Voronoi

図を考えるはうがより一般的である。

重み付 き

Voronoi

領域の定義式 は,以下のようになる

vi

‑(

pIPi+td

(

p,pi)≦P]+td

(

p,p]))

ここでの解釈 は, 可園の店があって,それぞれは品ぞろえや価格 に他の店 と は差別化を している

消費者 は,その店 までの距離 と価格要素を考慮 して店を 訪れ るとす る。

Pi+td

b

,pi

) は, i店で商品を買 う場合,単価

P

tと間接

コス トがかか ることを表わ している

亡は,距離当 りの コス トを表わす。

(11)

最適施設配置問題 の社会情報学 的考察

5 Vor onoi 図と DeJ aunay 三角形図

本節では,前節で紹介 した

Voronoi

図についてより厳密な解説を行 う

平面における

2

点 x

,yに対 して,その距離を

d( x, y) で表わす.

59

5. 1

距離の定義

2

点の距離の一般的な定義 としては,第

1

節でみたような,ユークリッ ド距 離,タクシー距離を考えるのが普通であるが,ここでは, 一般座標における 「ミ

ンコフスキー距離」を導入 し,特殊ケースを誘導す る。 ミン' 37スキ‑距離 と は,

Tl

dL

p ( p

,Pi)‑ [∑ lxi‑X

j lF p ] 1/p

I1

で表 される。 ここで

, (x l,x 2

, ‥.

,Xi

,‥.

,X

n ) と

(

x, 1

,X

j 2 , …

,X ,t,‥ ., xj

n ) はそれぞれ,p と

p i

のデカル ト座標である。 pは,ベキ乗の次数であ り,その取 りうる範囲は,

1≦β<

‑とする。

β‑ 1

な らば,

n

d

L l (

p,Pi)∑ k i ‑Xjil

である。 これは,タクシー距離,マ ン‑ ツタン距離,矩形距離 と呼ばれるもの である

β‑ 2

な らば,

dL

l (

p,Pt)

‑ ∫ i ∑ n

=

1 ( x

l‑xji)2

である。 これは,ユークリッ ド距離にはかな らない。

p

‑‑な らば,

dL

(

p,Pi)‑TnaXi=1.2, .

n

(lxl‑Xji

l )

であ り, 最大距離, 優越距離 と呼ばれ る。 ( 極限操作 (ロピタルの定理)を使 っ

(12)

60

46

第 1号

て,得 られ る。)

xl‑ x2

平面 にこれ らを図示す ると下図が得 られる。

£2

3

p= 1 2 C()

0

2. p‑

1,

2

,

3

,‑に対応する ミシコフスキー距離 の等高線

以上 は,距離の数学的な概念であるが,社会学,心理学では, 「 心理距離」

がある。実験上,指数型,対数型,ベキ乗型の組み合わせが考え られる。登山 や,スポーツの距離では, S字型距離がよ くあてはま りそ うである。

5. 2 Voronoi

図 とその構成法

x‑ (x

l

,x2

,…

,Xn)⊂R2

を平面上の nl 国の互 いに異な る点 の集合 とす るO 領域,

R(xi)‑fxLxeR2

,d(

x,xi)<d(x,x

, ・ )

,j≠i)

を母点xi の

Voronoi

領域 とい う.

2

つの

Vor9nOi

領域が共有す る辺を

Ⅴ。ronoi

辺,

3

つ以上の

Voronoi

領域の境界が共有す る点を

Voronoi

点 と

(13)

最適施設配置問題 の社会情報学 的考察

61

い う

Voronoi

図 とは, この よ うに して定 ま る平面 の

Voronoi

領域,

Voronoi

辺,

Voronoi

点に分割 した全体の図をい

う 。

Voronoi

図 という名前は,数学者

GeorgeeFedos?evichVoronoi(1868

‑1908)

に由来す る。彼 は

,1907‑1908

年 に

2

次形式論について

2

篇の論文 をフランス語で書 き,それが後世,先駆的業績 と評価されているが,

2

次形式 が平行多面体の頂点を母点 とする分割 とみな し解析可能であることを示 したこ

とによる 9) 0

Voronoi

図の作図法を述べよう

3

点,Xl ,

x 2

,X

3

が与え られた とす る。

3

点 の

Voronoi

図 は,線分 X

Ix 2

,

x 2

X

3

,X

3

X l の垂直

2

等分線を作図すればよい。

Voronoi

点 は, Xl ,

x 2

,X

3

を通 る外接円の中心,すなわち,外心にはかな らない。

4

点以上の

Voronoi

図は,どのようになるか。 この作図法 としては,逐次 添加法が実用上( 処理速度, 構造の単純 さ) , もっとも優れているといわれる。

この方法は,既に

3

点の

Voronoi

図が作 られているので これを基 に

1

点を追 加 して更新 していくものである。いわば, 「ブー トス トラップ」法である。

今,図

3‑ 1

のように,Xl ,

x 2

, X

3

か ら

Voronoi

図が作 られていて,新た に,X

4

が添加 されたとする。X

4

と一番近いX

3

( X

3

とX

4

は同一の

Voronoi

領域に属する)との垂直

2

等分線を引き,

Voronoi

辺 との交点の一つを

Y

l と するo

Y

l の属する別の

Voronoi

領域の母点XlとX

4

の垂直

2

等分線を引き,

Voronoi

辺 との交点の一つを

Y2

とする。

Y2

の属する別の

Voronoi

領域の母 点

x 2

とX。 の

Voronoi

辺 との交点をY

3

とすると,Y

3

は最初の

Voronoi

辺へ の交点

y

l の別の交点に一致す る。よって,

y

l

y 2y 3

の三角形の内部が,

方 。

Voronoi

領域で,

Voronoi

図が再構成された ( 図

3‑ 2)

9)Voronoi,G.,Nouvellesapplicationsdesparam昌trescontinusalath60‑

riedesformesquadratiques

,

PremierM6moire:Surquelquespropri6t6sdesformesquadratiques positivesparfaites

,J.

ReineAngew.Math.133(1907)pp.97‑178.

Deuxi8meMemoire:Recherchessurlespara1161108dresprlmitifs

,

J.ReineAngew.Mach.134(1908),pp.198‑287.

(14)

62

46

巻 第

1

V ( X l )

3.逐次添加法 による Voronoi

図の作成

Voronoi

領域 は, 通常,多角形 とな るので

,Voronoi

多角形 と呼ぶ。いま,

Voronoi

多角形 の隣接関係 について注 目しよ う

.

例えば,公共施設 x に出掛 けて行 ったが,たまたま,事情があって,利用で きなか った とす る。その′ / とき には, x の周辺を見渡 して,一番便利のよい ( 近 い)所を探すであろうo x の 利用圏に隣接す るすべての利用圏が候補であることは間違 いないが, どれが最 適であろ うか。 方と隣接す る

Voronoi

領域 の母点を次 々に結んでい くと,辛 直線が得 られ るが,すべての領域 について同様の操作を続 けると,三角形図が 得 られ る。 ( 退化が起 こらない と仮定 して い る。) これが,Del

aunay

三角形 図あ るいは,Del

aunay

三角形分割図 と呼ばれ る。

Delaunay

三角形 図 は, ロシアの数学者

BorisDelaunay

(ドローネ と読 む。Del

one

とも書 く。 )の

1934

年の 「 空 白円」に関す る論文 に由来す る

10)

0

Voronoi

点で は, 少な くとも

3

本の

Voronoi

辺が集 まることと

Delaunay

三角形 は,

3

つの

Voronoi

領域 の母点 を結んでで きた ものであることか らそ れぞれの

Delaunay

三角形 には,一つずつ

Voronoi

点が対応 している。

10)Surlasph占resvide,Izv.Akad.NaukSSSR,OtdelenieMateTnaticheskiii EstestueTmyhaNauh7 (1934),793‑800.

(15)

最適施設配置問題 の社会情報学 的考察

63 Voronoi

点か ら

3

つの母点 までの距離が等 しいことか ら

Voronoi

点 を中 心 として

Delaunay

三角形の

3

つの頂点を通 る外接円を措 くことがで きる。

外接 円 は,他 に頂 点 が釆 な い とい う意 味で,空 白円 とい う

その とき,

Voronoi

点は,

Delaunay

三角形の外心である。

逆に,

Delaunay

三角形図の隣接関係を見 ると,それは

Voronoi

辺 にはか ならない。 このことか ら,

Voronoi

図 と

Delaunay

三角形図は,双対 ( そ う つい)関係 にあるとい う

この関係か ら,

Delaunay

三角形図は,

.Voronoi

図を基 に作 り出せ,

Voronoi

図は,

Delaunay

三角形図か ら構成 されるので ある。

実際的な応用 として,小樽市街にある

,28

個の郵便局の所在 と,それ らの

Voronoi

図,

Delaunay

三角形図を求めたので,以下に示す。 これは,利用 者の利用圏や,適正な配置かどうかを知 る上で,有効な社会的情報 となる

11)0

4. 1

:小樽市街 に実在す る郵便局の配置 ( 次図参照)

ll)

若林ゼ ミナールの学生であった田中 麻友 さんの計算に負 う所が大 きい。

(16)

商 学 討 究 第

46

巻 第

1

4. 2 :

小棒市街の郵便局 を母点 とするポ ロノイ図

4. 3 :

ドローネ三角形図

(17)

最適施設配置問題 の社会情報学 的考察

65

6 Voronoi

図を求めるソフ トウェア

Voronoi

図を作成す るコンピュータプログラムには,い くつかの ものが知 られている。イ ンターネ ッ ト上の

archie

を実行 し,Vor

onoi

を探せば,

1.Voronoi.tar.Z40929Apri1301993

2.new‑version‑convex‑hullDelaunay‑Voronoi2294May261994

が 国内サイ トで見つか る

また,Del

aunay

で探せば,Mat

hematica

で実行で きる

3.Delaunay.m 5706Aug161994

や,

4.Delaunay.SUN

が見つか る。さ らに,外国まで,範囲を広げると,相当数の公開されたソフ ト ウェアが様 々な分野に応用 されて存在す ることがわか る

12)

以下では,上記以外の

3

つを紹介す る。

第一 は,杉原 ・伊理 プログラム

(1989)13)

である。 これは,Fort

ranプログ

ラム2,300 行で提供 され,Sun . 7‑クステー ションで コンパイル,実行 して使

ユーザーは,母点の座標 さえ入力すれば,算法の詳細を知 らな くて も,杏 易に,Voronoi 図や

Delaunay

三角形分解図を得 ることがで きる。算法の詳 細を知れば,Fort

ranプログラムを修正 して, 自分 の好 む ものが得 られ る。

実際,ユーザーイ ンタフェースが貧弱であるので,プログラムを読んで使い易 くしなければな らなか った。 しか し,数値計算の精度,頑健性 に定評がある。

第 2は,qhu日と呼ばれ る, C言語で書かれた,登録商標付 きのフ リーソフ トの凸包 プログラムである

14)0

12)netscapehttp://opentext.uunet.ca:8080/

を実行 し

,Voronoi

をサーチ すると,

̀̀voronoi"‑375hitsin158pages

158

か所に

375

件が表示され,そこか ら関連情報を得る ( サーフィンする)こと ができる。

(1995

4

28

日現在)

13)Sugihara,K.andlri

,

M.,Voronoi2ReferenceManual,RMI89‑04

, 東京大学工学部計数工学科

14)rtpgeom.umn.edu

( または

,ftp 128.101.25.35)とすれば,入手すること

ができる。

(18)

66

46

1

現在,

2.02(1995.01.25)

版であ り,米国 ミネソタ大学の ジオメ トリ ・セ ンターが提供 している

.

凸包の計算,

Delaunay

三角形,

Voronoi

図を

3

次 元以上 まで,最 も高速に計算す る。グラフ出力まで行 うので,多数のオプショ

ンを持つ大 きなパ ッケー ジで奉るo

Rbox

と呼ばれ る

Qhul

lのための入力を 生成す るツール も備えている。

出力 として,

Mathematica

用の入力データを生成す ることができることも 特徴である。

qhul

lの使用法を簡単に示す と以下のようである。

使用法 :

qhull

入力 ( 標準入力) :次元,点の個数,点の座標 主なオプション :

d Delaunay

三角形

V Voronoi

(Delaunay

三角形図か ら)

Tv

結果の検証

主な出力オプション ( デフォル トでは,要約情報を表示する。):

f すべての場面を印刷す る。

G Geomview

出力

1 (Voronoi

の中心)各々の面に接続する頂点

m Mathematica

出力

p (Voronoi

の中心)点の座標

3

は,数式処理言語 として名高い,

Mathematica

Package

利用であ る。周知のように,

Mathematica

は,数値計算,数式処理,グラフィックス 処理,マルチメディア機能を統一的な環境の下においているが,第

2.

0 版より,

DiscreteMath ̀ComputationalGeometry

̀とい う

Package

が提供 され た。ユーザーは,

くく

Folder

名 ( す なわ ち,

DiscreteMath)

̀ P

ackage

名 ( す なわ ち,

(19)

最適施設配置問題 の社会情報学 的考察

67

Comput at i onal Geom et r y) ̀ によ り,直接読 み込んで,以下 の

13

個 の シン ボル名のい くつかを,実行 させ る。

1.DelaunyTriangulation 2.VoronoiDiagram

3.ConvexHull

4.NearestNeighbor

5.TriangularSurfacePlot 6.PlanarGraphPlot

7.DiagramPlot‑

8.DelaunayTriangulationQ 9.LabelPoints

10.Ray

ll.AllPoints

12.Hull

13.TrimPoints

関数

VoronoiDiagram [

Hx l , y l) ,( x

2

,y

2

) ,‥. ,i xn,ynH ] は,辛 面の多数の点の

Voronoi

図を作成す る。結果 は,

2

つの リス トによって表わ され る

Voronoi

図の頂点の座標 リス トと,頂点 の隣接 リス トである。関数 の引数 に, さ らに

,val

hul

lを加えることによって,

Voronoi

図をよ り効 率よ く計算す ることがで きる。

Mathematica

を用いて

Voronoi

図を求めるには,オ ンライ ンマニ ュアル ともいうべ きパ ッケージの中身を読むのがよい。パ ッケージには,

題なども 詳 しく盛 り込んである

15)0

15)

この点に関しては,平成

6

年度特定研究経費

(

「 数式処理言語

Mathematica

を利 用 した数理モデルの解析

」(2

1

) )に基き

,PowerMacintos

b上で実施 した。

記して感謝 したい。

(20)

68

46

1

7

結語的覚書

最適施設配置の問題 は,古 くて新 しい問題である。また,応用範囲 も広 く, 単に,地理的な施設配置にとどま らない。社会的な人員配置や ロケーション問 題であるところの雇用の調整,組織の編成の問題 などにも考え方 には共通の要 素がある

最適施設配置 は, 本稿のよ うな数学的な処理や, 情報科学的な処理 によって, 解を提案で きるが,その実施 は,政策者の決断にゆだね られ る。社会情報アプ

ローチとしては,データの収集,処理,整理を適切 に行い,問題の枠組みを設

定 し,定式化 し,問題の解を得 るなどして政策者の価値判断を助 けるために情

報を提供す ることにある

そのさい, どのような 「 仮定」が設定 されているか

注意す る必要がある

図 4. 2 : 小棒市街の郵便局 を母点 とするポ ロノイ図

参照

関連したドキュメント

日本オペレーションズ・リサーチ学会 2004年秋季研究発表会 1−C−8 施設配置問題に対する近似アルゴリズムの実験評価

アブストラクト: オンライン最適化問題とは、

ラウド・サービスの利用が一般的になっている.大規模 なソフトウェア・システムをクラウド上に構築する場 合,必要なソフトウェア要素 (CMS サーバ,

ほかによく遭遇する最適化の例として施 設配置の問題がある。ある地域に郵便局を

またこの問題を考えるために筆者らは[1]において近似アルゴリズムを提案 している。 それは, 議席配分問題[2]のアルゴリズムを適用した割り当て決定 の部分

α 〈wl,…,W〝) (6) 絶対値基準に対し,施設利用者は最も便利な場所 に施設が立地された場合に対して現在の立地点を評

競合性を考慮した好ましくない施設の配置問題 阪大 工宇野 剛史 (Takeshi Uno) 石井 博昭 (Hiroaki Ishii) 斎藤 誠慈 (Seiji Saito) Osaka

競合環境下における施設配置問題 大阪府立大学総合科学部 大角盛広 (OSUMI Shigehiro) 大阪大学 I 学部 石井博昭 (ISHII Hiroaki)