誤差論と最小二乗法

第8回付録 線形(線型)代数の基礎2

A2. 行列の演算

一次(線形)写像

image98image99ベクトル、image100image101ベクトル、image102image103行列としたとき

 

image104

 

は、image105次元ベクトルimage106image107次元ベクトルimage108へ写すということを示しています。これを一次(線形)写像といいます。image109の時に一次変換という場合もあります。

 

二次形式
image102image110対称行列、image9image101ベクトルとしたとき積

 

image111

 

を二次形式といいます。値はスカラー(数)です。image102が対称でない正方行列だとしても

 

image112

ですから、

 

image113

 

となり、二次形式は常に対称行列で書けることがわかります。

 

正定値行列

対称行列image102がすべてのimage114に対して、image115となるとき、正定値といいます。また、image116となるとき半正定値といいます。

 

対称行列image102が正定値なら、正則行列image117があって、

 

image118

となります。従って、image102も正則です。

 

以下の定理があります。
image102が正定値なら、その対角成分image120 
image102が半正定値なら、その対角成分 image121

 

また、


image102image122正定値行列、image123image124行列としたとき、
image125ならば、image126は正定値        (5a)
image127またはimage128ならば、image126は半正定値   (5b)
です。

(参考) 共分散行列は正定値または半正定値ですが、半正定値となるのは変数間に線形関係があるときに限られます。

 

例. 二つの確率変数に常にimage129という関係があるとしましょう。image130と書けるので、image131の共分散行列はimage132の分散image14から誤差伝播により、

 

image133

 

となります。この共分散行列は、二次形式を作ると

 

image134

であり、image135の時は0となるので半正定値です(上記(5b)の場合)。

 

連立一次方程式

未知数image136個のimage105個の連立方程式

 

image137

 

は行列の形に書けます。

 

image138

 

image102image139行列、image100image140image141image101です。 image142image102が正則ならば(6)の解がただ一つ、

 

image143

 

と求まります。

 

一般に、image144なら方程式の数が未知数の数より多いので解は存在しません。image145なら解は無限に存在します。次の定理があります。

 

image146に解が存在する必要十分条件は、image147。    (8)

ここで、image148は、image102に列image141を付け加えた行列です。


image149と書くと

image150

 

ですから、image151となるimage100があるimage152(image153の列空間)
image154image155、となります。

 

例1.

image156
image157のランクは、第1列と第2列を足すと第3列なので、

image158です。従って解をもちます。第3式×2+第2式=第1式×3ですから、実は1式は不要で解は一意に決まります(図1)。

 

image159図1.一意の解

 

 

例2.

image160
image161の列は一次独立で、image162となり、解を持ちません(図2)。このような場合、近似的に最適な解は、最小二乗法により

 

image163

 

を最小にすれば求めることができます。

 

image164

図2.解なし

 

 

例3.

image165

image166で第1式と第2式を足すと第3式になるので、実質2式であり
image158となり解をもちます。第1,2式より

 

image167

 

と表せるので、解は、

 

image168

 

となります。image169は任意なので解は無限にあります。

 

行列式

image110行列image170の行列式は、各行と列から1つずつ行列要素をとって掛け合わせ足した関数で次のように定義され、image171またはimage172と書きます、

 

image173

 

ここでimage174image175の並べ替え(置換)で、和はすべての置換について取り、符号は置換が偶なら正、奇なら負となります(置換はもとのimage175から2つずつ何回か入れ替えると得られますが、その回数が偶数なら偶置換、奇数なら奇置換といいます)。例えば、

 

image176

image177

 

 

行列式には次のような性質があります。
image178が正則でない⇒image179
image178が正則⇒image180

image181
image182
image178が正定値ならimage183
image178が正則ならimage184

 

また、正則な行列image153の逆行列は行列式を使って次のように書けます。

 

image185

 

ここで、 image186
例えば、(2 × 2) 行列なら

 

image187

B. 多次元正規分布の確率密度関数

多次元の正規分布は以前にも紹介しましたが、ベクトルと行列で表示すると便利であるとともに変数の関係などが見えやすいのでここにまとめておきます。

 

標準化(平均=0、分散=1)されたimage188個の独立な正規分布確率変数image189があるとします。その時の同時確率密度関数は、image190として、

 

image191

image192

 

となります。image193である一般の確率変数に戻すには

 

image194

 

と変換すればよく、image9の確率密度関数は、

 

image195

 

となります。

C. テイラー展開

image196image197)を滑らかな(何回でも微分可能な)実数値関数とし、image198image100の近傍の点とすると

 

image199

 

と表されます。image200は滑らかな関数です。第3項はimage201の二次になっていますから、image201が小さいとすると無視することができ

 

image202

image203

 

となりimage204の線形式で書けることがわかります。測量では、未知数と観測量の関係が線形でないとき、テイラー展開を使って線形化するのが一般的です。

 

Share on FacebookTweet about this on TwitterShare on LinkedInEmail this to someonePrint this page

パラメータ情報

【地殻変動補正提供サービス】パラメータの生成と評価:2020年4月1日~4月30日

2020年4月1日から2020年4月30日を有効期間とする

地殻変動補正提供サービス パラメータを生成しましたので、お知らせいたします。

 

パラメータの精度評価結果は以下の通りです。

水平方向の推定誤差の平均値は、 0.010 mです。

水平方向の推定誤差の標準偏差は、0.006 mです。

水平方向の推定誤差が4cm以下の地域の割合は、99.44%です。

 

なお、今回生成したパラメータは、すべての電子基準点のデータを使用し、

パラメータの生成および精度評価を行っています。

今後、パラメータと評価指標の精度向上のため、

解析条件を変更する可能性がありますので、予めご了承ください。

Share on FacebookTweet about this on TwitterShare on LinkedInEmail this to someonePrint this page

座標系講座 生活編

第5回 まとめ

平成29年第3期地理空間情報活用推進基本計画によれば、平成30年度より我が国の準天頂衛星4機体制を確立し、衛星測位サービス、サブメータ級測位補強サービス等を開始します。更に、平成35年度を目途とする持続測位が可能な7機体制の構築に向けた取り組みが行われています。こうした衛星測位の高精度化及び高度化に伴って、国民は正確な位置情報を得ることができます。しかし、これらの正確な位置情報は、測量法上の位置ではなく、地上における国民生活上の位置情報としては、測量法及び測量作業規程の整備が要求されてきていることは、既に述べてきたところです。以下にそれらの整備に関する内容を列挙します。

 

1. フレームの統一

日本の測地座標系(フレーム)は、西日本及び北海道はITRF94(元期1997.0)及び2011年東北地方太平洋沖地震地域はITRF2008(元期2011.4)の2つによって構成されています(連載第1回 図16)。従いまして、これら二つの境界付近での処理が複雑になっています(連載第1回  8.フレーム境界における座標の不整合参照)。統一したフレームの構築は、衛星測位をより効果的なものにするためにも、早期の実現が望まれます。

 

2. 4次元の位置情報

第3期地理空間情報活用推進基本計画は、“全国の電子基準点による地殻変動の即時把握の技術開発”が定められています。現在の測量法は、定常的な地殻変動のような時間変化に対応した規定は存在しません。地理空間情報活用推進基本法と同様に、位置情報として4次元対応の法律的根拠が必要になります。技術的には、リアルタイム変動に対応可能な当社が開発した「セミ・ダイナミック リダクション」の活用があります。

 

3. 2+1次元処理から3次元処理への移行

衛星測位の高度化に伴って既に述べたような測量法の整備が必要になります。更に、現在の測量法は、水平位置と高さは別々な処理となっています。測地学でいう「2+1次元」で3次元処理ではありません。具体的には、日本経緯度原点の測地原子として楕円体高を加えることです。また、トータルステーションにより得られた3次元データは、わざわざ平面と高さに分離し、手間のかかる計算精度の落ちる処理をしています。3次元データは、素直に3次元処理を行うような測量作業規程の改善が求められています。

 

4. 相対測位から絶対測位への対応

標観測が中心になり、あらゆる位置測量は、座標観測を中心とした細部測量へと移り代わることになります。こうした技術の進歩に見合った測量作業規程の改正が必要になって来るので、測量関係者はその取り組みを行う必要があります。

 

5. GRS80楕円体の物理定数の規定

以上考察してきましたように、測量法施行令に楕円体高の規定を行うこと及び水準測量と衛星測位から得られる2種類の標高の統一があります。加えて、現在の標高計算にあたって、高さを決めるとき用いる正標高補正及び正規正標高補正の根拠となるGRS80楕円体の物理定数である地心引力定数及び地球の自転角速度の規定が必要になります。

Share on FacebookTweet about this on TwitterShare on LinkedInEmail this to someonePrint this page