Sub 連立一次方程式(係数(), 定数(), 解(), r)
Rem (c) 2014 Kaoru Fueda
Dim 係数行列() As Double, 最大値() As Double
Dim i As Integer, j As Integer, k As Integer, 最大行 As Integer, 最小行 As Integer
Dim temp As Double, 最大比 As Double, 最小値 As Double
ReDim 係数行列(r, r + 1), 最大値(r)
For i = 1 To r
For j = 1 To r
係数行列(i, j) = 係数(i, j)
Next j
係数行列(i, r + 1) = 定数(i)
Next i
Rem iはピボット
For i = 1 To r
Rem ピボット選択
If i < r Then
最大比 = 0
最大行 = i
For j = i To r
最大値(j) = 0
For k = i To r + 1
If Abs(係数行列(j, k)) > 最大値(j) Then
最大値(j) = Abs(係数行列(j, k))
End If
Next k
If 最大値(j) > 0 Then
If 最大比 < Abs(係数行列(j, i) / 最大値(j)) Then
最大比 = Abs(係数行列(j, i) / 最大値(j))
最大行 = j
End If
End If
Next j
If 最大行 > i Then
For k = 1 To r + 1
temp = 係数行列(i, k)
係数行列(i, k) = 係数行列(最大行, k)
係数行列(最大行, k) = temp
Next k
End If
Rem ピボットが0の場合は最も0=0に近い行を第i行に移動して、後に無視する
If 最大比 = 0 Then
最小値 = 最大値(i)
最小行 = i
For j = i + 1 To r
If 最小値 > 最大値(j) Then
最小値 = 最大値(j)
最小行 = j
End If
Next j
If 最小行 > i Then
For k = 1 To r + 1
temp = 係数行列(i, k)
係数行列(i, k) = 係数行列(最小行, k)
係数行列(最小行, k) = temp
Next k
End If
End If
End If
Rem ピボットを1に
If 係数行列(i, i) <> 0 Then
For k = i + 1 To r + 1
係数行列(i, k) = 係数行列(i, k) / 係数行列(i, i)
Next k
Else
MsgBox ("解けないので解の一つを0とします。")
Rem 本来は解ベクトルの長さが最小とすべき。
For k = i + 1 To r + 1
係数行列(i, k) = 0
Next k
End If
係数行列(i, i) = 1
Rem ピボット以外を0に
For j = 1 To r
If j <> i Then
For k = i + 1 To r + 1
係数行列(j, k) = 係数行列(j, k) - 係数行列(i, k) * 係数行列(j, i)
Next k
係数行列(j, i) = 0
End If
Next j
Next i
For i = 1 To r
解(i) = 係数行列(i, r + 1)
Next i
End Sub
中央部の2/3×√(1200/1450)≒1/1.65インチ部分しか使っていないのかな？と思い、実焦点距離と35mm判換算焦点距離の換算倍率を見てみました。
参考として1インチセンサーのソニーのRX100とも比べました。
| 機種 | 実焦点距離 | 換算焦点距離 | 換算倍率 |
|---|---|---|---|
| XQ1 | 6.4 | 25 | 3.91 |
| XQ1 | 25.6 | 100 | 3.91 |
| X20 | 7.1 | 28 | 3.94 |
| X20 | 28.4 | 112 | 3.94 |
| RX100 | 10.4 | 28 | 2.69 |
| RX100 | 37.1 | 100 | 2.70 |
色々試してみたところ
fanboy-annoyances
https://secure.fanboy.co.nz/fanboy-addon.txt
を使うとfacebookの広告も消すことができ、日本独自の広告を消すにはABP Japanese Filtersの代わりに
fanboy-Japanese
https://secure.fanboy.co.nz/fanboy-japanese.txt
を使えばgoogle検索窓が消えないことが分かりました。
