配列・ベクトル・行列(11)

前の配列・ベクトル・行列(10)の余談８にて「４．尾固定小行列循環シフト型(仮)」の順列に基づいた行列式の作り方を模索してみた。それよりも因数分解表示はそれほど悪くはないと思っていたが、展開した形がこれほどまでにシンプルで最適だとは予想外だった。ネストが全然増えないのが素晴らしい。四則演算に演算子がほぼ使えない状況だから演算をまとめた方が綺麗なのかも。

 

まずは展開した行列式　(複素数版を修正して実数版も作ってみた。実数版は１行の文字数増加率が高いので１行にしたかったものが１行半になってしまったり。　r070506)

 

2x2　元行列=　AD16:AE17

実　AC23=　=AD16*AE17-AD17*AE16

複　AC23=　=IMSUM(IMPRODUCT(+1,AD16,AE17),IMPRODUCT(-1,AD17,AE16))

 

3x3　元行列=　AD16:AF18

実　AC25=　=AD16*AE17*AF18-AD17*AE16*AF18+AD18*AE16*AF17-AD16*AE18*AF17+AD17*AE18*AF16-AD18*AE17*AF16

複　AC25=　=IMSUM(
IMPRODUCT(+1,AD16,AE17,AF18),IMPRODUCT(-1,AD17,AE16,AF18),
IMPRODUCT(+1,AD18,AE16,AF17),IMPRODUCT(-1,AD16,AE18,AF17),
IMPRODUCT(+1,AD17,AE18,AF16),IMPRODUCT(-1,AD18,AE17,AF16))

 

4x4　元行列=　AD16:AG19

実　AC27=　=AD16*AE17*AF18*AG19-AD17*AE16*AF18*AG19+AD18*AE16*AF17*AG19-AD16*AE18*AF17*AG19
+AD17*AE18*AF16*AG19-AD18*AE17*AF16*AG19-AD19*AE16*AF17*AG18+AD16*AE19*AF17*AG18
-AD17*AE19*AF16*AG18+AD19*AE17*AF16*AG18-AD16*AE17*AF19*AG18+AD17*AE16*AF19*AG18
+AD18*AE19*AF16*AG17-AD19*AE18*AF16*AG17+AD16*AE18*AF19*AG17-AD18*AE16*AF19*AG17
+AD19*AE16*AF18*AG17-AD16*AE19*AF18*AG17-AD17*AE18*AF19*AG16+AD18*AE17*AF19*AG16
-AD19*AE17*AF18*AG16+AD17*AE19*AF18*AG16-AD18*AE19*AF17*AG16+AD19*AE18*AF17*AG16

複　AC27=　=IMSUM(
IMPRODUCT(+1,AD16,AE17,AF18,AG19),IMPRODUCT(-1,AD17,AE16,AF18,AG19),
IMPRODUCT(+1,AD18,AE16,AF17,AG19),IMPRODUCT(-1,AD16,AE18,AF17,AG19),
IMPRODUCT(+1,AD17,AE18,AF16,AG19),IMPRODUCT(-1,AD18,AE17,AF16,AG19),
IMPRODUCT(-1,AD19,AE16,AF17,AG18),IMPRODUCT(+1,AD16,AE19,AF17,AG18),
IMPRODUCT(-1,AD17,AE19,AF16,AG18),IMPRODUCT(+1,AD19,AE17,AF16,AG18),
IMPRODUCT(-1,AD16,AE17,AF19,AG18),IMPRODUCT(+1,AD17,AE16,AF19,AG18),
IMPRODUCT(+1,AD18,AE19,AF16,AG17),IMPRODUCT(-1,AD19,AE18,AF16,AG17),
IMPRODUCT(+1,AD16,AE18,AF19,AG17),IMPRODUCT(-1,AD18,AE16,AF19,AG17),
IMPRODUCT(+1,AD19,AE16,AF18,AG17),IMPRODUCT(-1,AD16,AE19,AF18,AG17),
IMPRODUCT(-1,AD17,AE18,AF19,AG16),IMPRODUCT(+1,AD18,AE17,AF19,AG16),
IMPRODUCT(-1,AD19,AE17,AF18,AG16),IMPRODUCT(+1,AD17,AE19,AF18,AG16),
IMPRODUCT(-1,AD18,AE19,AF17,AG16),IMPRODUCT(+1,AD19,AE18,AF17,AG16))

 

5x5　元行列=　AD16:AH20

実　AC29=　=AD16*AE17*AF18*AG19*AH20-AD17*AE16*AF18*AG19*AH20+AD18*AE16*AF17*AG19*AH20-AD16*AE18*AF17*AG19*AH20
+AD17*AE18*AF16*AG19*AH20-AD18*AE17*AF16*AG19*AH20-AD19*AE16*AF17*AG18*AH20+AD16*AE19*AF17*AG18*AH20
-AD17*AE19*AF16*AG18*AH20+AD19*AE17*AF16*AG18*AH20-AD16*AE17*AF19*AG18*AH20+AD17*AE16*AF19*AG18*AH20
+AD18*AE19*AF16*AG17*AH20-AD19*AE18*AF16*AG17*AH20+AD16*AE18*AF19*AG17*AH20-AD18*AE16*AF19*AG17*AH20
+AD19*AE16*AF18*AG17*AH20-AD16*AE19*AF18*AG17*AH20-AD17*AE18*AF19*AG16*AH20+AD18*AE17*AF19*AG16*AH20
-AD19*AE17*AF18*AG16*AH20+AD17*AE19*AF18*AG16*AH20-AD18*AE19*AF17*AG16*AH20+AD19*AE18*AF17*AG16*AH20
+AD20*AE16*AF17*AG18*AH19-AD16*AE20*AF17*AG18*AH19+AD17*AE20*AF16*AG18*AH19-AD20*AE17*AF16*AG18*AH19
+AD16*AE17*AF20*AG18*AH19-AD17*AE16*AF20*AG18*AH19-AD18*AE20*AF16*AG17*AH19+AD20*AE18*AF16*AG17*AH19
-AD16*AE18*AF20*AG17*AH19+AD18*AE16*AF20*AG17*AH19-AD20*AE16*AF18*AG17*AH19+AD16*AE20*AF18*AG17*AH19
+AD17*AE18*AF20*AG16*AH19-AD18*AE17*AF20*AG16*AH19+AD20*AE17*AF18*AG16*AH19-AD17*AE20*AF18*AG16*AH19
+AD18*AE20*AF17*AG16*AH19-AD20*AE18*AF17*AG16*AH19-AD16*AE17*AF18*AG20*AH19+AD17*AE16*AF18*AG20*AH19
-AD18*AE16*AF17*AG20*AH19+AD16*AE18*AF17*AG20*AH19-AD17*AE18*AF16*AG20*AH19+AD18*AE17*AF16*AG20*AH19
+AD19*AE20*AF16*AG17*AH18-AD20*AE19*AF16*AG17*AH18+AD16*AE19*AF20*AG17*AH18-AD19*AE16*AF20*AG17*AH18
+AD20*AE16*AF19*AG17*AH18-AD16*AE20*AF19*AG17*AH18-AD17*AE19*AF20*AG16*AH18+AD19*AE17*AF20*AG16*AH18
-AD20*AE17*AF19*AG16*AH18+AD17*AE20*AF19*AG16*AH18-AD19*AE20*AF17*AG16*AH18+AD20*AE19*AF17*AG16*AH18
+AD16*AE17*AF19*AG20*AH18-AD17*AE16*AF19*AG20*AH18+AD19*AE16*AF17*AG20*AH18-AD16*AE19*AF17*AG20*AH18
+AD17*AE19*AF16*AG20*AH18-AD19*AE17*AF16*AG20*AH18-AD20*AE16*AF17*AG19*AH18+AD16*AE20*AF17*AG19*AH18
-AD17*AE20*AF16*AG19*AH18+AD20*AE17*AF16*AG19*AH18-AD16*AE17*AF20*AG19*AH18+AD17*AE16*AF20*AG19*AH18
+AD18*AE19*AF20*AG16*AH17-AD19*AE18*AF20*AG16*AH17+AD20*AE18*AF19*AG16*AH17-AD18*AE20*AF19*AG16*AH17
+AD19*AE20*AF18*AG16*AH17-AD20*AE19*AF18*AG16*AH17-AD16*AE18*AF19*AG20*AH17+AD18*AE16*AF19*AG20*AH17
-AD19*AE16*AF18*AG20*AH17+AD16*AE19*AF18*AG20*AH17-AD18*AE19*AF16*AG20*AH17+AD19*AE18*AF16*AG20*AH17
+AD20*AE16*AF18*AG19*AH17-AD16*AE20*AF18*AG19*AH17+AD18*AE20*AF16*AG19*AH17-AD20*AE18*AF16*AG19*AH17
+AD16*AE18*AF20*AG19*AH17-AD18*AE16*AF20*AG19*AH17-AD19*AE20*AF16*AG18*AH17+AD20*AE19*AF16*AG18*AH17
-AD16*AE19*AF20*AG18*AH17+AD19*AE16*AF20*AG18*AH17-AD20*AE16*AF19*AG18*AH17+AD16*AE20*AF19*AG18*AH17
+AD17*AE18*AF19*AG20*AH16-AD18*AE17*AF19*AG20*AH16+AD19*AE17*AF18*AG20*AH16-AD17*AE19*AF18*AG20*AH16
+AD18*AE19*AF17*AG20*AH16-AD19*AE18*AF17*AG20*AH16-AD20*AE17*AF18*AG19*AH16+AD17*AE20*AF18*AG19*AH16
-AD18*AE20*AF17*AG19*AH16+AD20*AE18*AF17*AG19*AH16-AD17*AE18*AF20*AG19*AH16+AD18*AE17*AF20*AG19*AH16
+AD19*AE20*AF17*AG18*AH16-AD20*AE19*AF17*AG18*AH16+AD17*AE19*AF20*AG18*AH16-AD19*AE17*AF20*AG18*AH16
+AD20*AE17*AF19*AG18*AH16-AD17*AE20*AF19*AG18*AH16-AD18*AE19*AF20*AG17*AH16+AD19*AE18*AF20*AG17*AH16
-AD20*AE18*AF19*AG17*AH16+AD18*AE20*AF19*AG17*AH16-AD19*AE20*AF18*AG17*AH16+AD20*AE19*AF18*AG17*AH16

複　AC29=　=IMSUM(
IMPRODUCT(+1,AD16,AE17,AF18,AG19,AH20),IMPRODUCT(-1,AD17,AE16,AF18,AG19,AH20),
IMPRODUCT(+1,AD18,AE16,AF17,AG19,AH20),IMPRODUCT(-1,AD16,AE18,AF17,AG19,AH20),
IMPRODUCT(+1,AD17,AE18,AF16,AG19,AH20),IMPRODUCT(-1,AD18,AE17,AF16,AG19,AH20),
IMPRODUCT(-1,AD19,AE16,AF17,AG18,AH20),IMPRODUCT(+1,AD16,AE19,AF17,AG18,AH20),
IMPRODUCT(-1,AD17,AE19,AF16,AG18,AH20),IMPRODUCT(+1,AD19,AE17,AF16,AG18,AH20),
IMPRODUCT(-1,AD16,AE17,AF19,AG18,AH20),IMPRODUCT(+1,AD17,AE16,AF19,AG18,AH20),
IMPRODUCT(+1,AD18,AE19,AF16,AG17,AH20),IMPRODUCT(-1,AD19,AE18,AF16,AG17,AH20),
IMPRODUCT(+1,AD16,AE18,AF19,AG17,AH20),IMPRODUCT(-1,AD18,AE16,AF19,AG17,AH20),
IMPRODUCT(+1,AD19,AE16,AF18,AG17,AH20),IMPRODUCT(-1,AD16,AE19,AF18,AG17,AH20),
IMPRODUCT(-1,AD17,AE18,AF19,AG16,AH20),IMPRODUCT(+1,AD18,AE17,AF19,AG16,AH20),
IMPRODUCT(-1,AD19,AE17,AF18,AG16,AH20),IMPRODUCT(+1,AD17,AE19,AF18,AG16,AH20),
IMPRODUCT(-1,AD18,AE19,AF17,AG16,AH20),IMPRODUCT(+1,AD19,AE18,AF17,AG16,AH20),
IMPRODUCT(+1,AD20,AE16,AF17,AG18,AH19),IMPRODUCT(-1,AD16,AE20,AF17,AG18,AH19),
IMPRODUCT(+1,AD17,AE20,AF16,AG18,AH19),IMPRODUCT(-1,AD20,AE17,AF16,AG18,AH19),
IMPRODUCT(+1,AD16,AE17,AF20,AG18,AH19),IMPRODUCT(-1,AD17,AE16,AF20,AG18,AH19),
IMPRODUCT(-1,AD18,AE20,AF16,AG17,AH19),IMPRODUCT(+1,AD20,AE18,AF16,AG17,AH19),
IMPRODUCT(-1,AD16,AE18,AF20,AG17,AH19),IMPRODUCT(+1,AD18,AE16,AF20,AG17,AH19),
IMPRODUCT(-1,AD20,AE16,AF18,AG17,AH19),IMPRODUCT(+1,AD16,AE20,AF18,AG17,AH19),
IMPRODUCT(+1,AD17,AE18,AF20,AG16,AH19),IMPRODUCT(-1,AD18,AE17,AF20,AG16,AH19),
IMPRODUCT(+1,AD20,AE17,AF18,AG16,AH19),IMPRODUCT(-1,AD17,AE20,AF18,AG16,AH19),
IMPRODUCT(+1,AD18,AE20,AF17,AG16,AH19),IMPRODUCT(-1,AD20,AE18,AF17,AG16,AH19),
IMPRODUCT(-1,AD16,AE17,AF18,AG20,AH19),IMPRODUCT(+1,AD17,AE16,AF18,AG20,AH19),
IMPRODUCT(-1,AD18,AE16,AF17,AG20,AH19),IMPRODUCT(+1,AD16,AE18,AF17,AG20,AH19),
IMPRODUCT(-1,AD17,AE18,AF16,AG20,AH19),IMPRODUCT(+1,AD18,AE17,AF16,AG20,AH19),
IMPRODUCT(+1,AD19,AE20,AF16,AG17,AH18),IMPRODUCT(-1,AD20,AE19,AF16,AG17,AH18),
IMPRODUCT(+1,AD16,AE19,AF20,AG17,AH18),IMPRODUCT(-1,AD19,AE16,AF20,AG17,AH18),
IMPRODUCT(+1,AD20,AE16,AF19,AG17,AH18),IMPRODUCT(-1,AD16,AE20,AF19,AG17,AH18),
IMPRODUCT(-1,AD17,AE19,AF20,AG16,AH18),IMPRODUCT(+1,AD19,AE17,AF20,AG16,AH18),
IMPRODUCT(-1,AD20,AE17,AF19,AG16,AH18),IMPRODUCT(+1,AD17,AE20,AF19,AG16,AH18),
IMPRODUCT(-1,AD19,AE20,AF17,AG16,AH18),IMPRODUCT(+1,AD20,AE19,AF17,AG16,AH18),
IMPRODUCT(+1,AD16,AE17,AF19,AG20,AH18),IMPRODUCT(-1,AD17,AE16,AF19,AG20,AH18),
IMPRODUCT(+1,AD19,AE16,AF17,AG20,AH18),IMPRODUCT(-1,AD16,AE19,AF17,AG20,AH18),
IMPRODUCT(+1,AD17,AE19,AF16,AG20,AH18),IMPRODUCT(-1,AD19,AE17,AF16,AG20,AH18),
IMPRODUCT(-1,AD20,AE16,AF17,AG19,AH18),IMPRODUCT(+1,AD16,AE20,AF17,AG19,AH18),
IMPRODUCT(-1,AD17,AE20,AF16,AG19,AH18),IMPRODUCT(+1,AD20,AE17,AF16,AG19,AH18),
IMPRODUCT(-1,AD16,AE17,AF20,AG19,AH18),IMPRODUCT(+1,AD17,AE16,AF20,AG19,AH18),
IMPRODUCT(+1,AD18,AE19,AF20,AG16,AH17),IMPRODUCT(-1,AD19,AE18,AF20,AG16,AH17),
IMPRODUCT(+1,AD20,AE18,AF19,AG16,AH17),IMPRODUCT(-1,AD18,AE20,AF19,AG16,AH17),
IMPRODUCT(+1,AD19,AE20,AF18,AG16,AH17),IMPRODUCT(-1,AD20,AE19,AF18,AG16,AH17),
IMPRODUCT(-1,AD16,AE18,AF19,AG20,AH17),IMPRODUCT(+1,AD18,AE16,AF19,AG20,AH17),
IMPRODUCT(-1,AD19,AE16,AF18,AG20,AH17),IMPRODUCT(+1,AD16,AE19,AF18,AG20,AH17),
IMPRODUCT(-1,AD18,AE19,AF16,AG20,AH17),IMPRODUCT(+1,AD19,AE18,AF16,AG20,AH17),
IMPRODUCT(+1,AD20,AE16,AF18,AG19,AH17),IMPRODUCT(-1,AD16,AE20,AF18,AG19,AH17),
IMPRODUCT(+1,AD18,AE20,AF16,AG19,AH17),IMPRODUCT(-1,AD20,AE18,AF16,AG19,AH17),
IMPRODUCT(+1,AD16,AE18,AF20,AG19,AH17),IMPRODUCT(-1,AD18,AE16,AF20,AG19,AH17),
IMPRODUCT(-1,AD19,AE20,AF16,AG18,AH17),IMPRODUCT(+1,AD20,AE19,AF16,AG18,AH17),
IMPRODUCT(-1,AD16,AE19,AF20,AG18,AH17),IMPRODUCT(+1,AD19,AE16,AF20,AG18,AH17),
IMPRODUCT(-1,AD20,AE16,AF19,AG18,AH17),IMPRODUCT(+1,AD16,AE20,AF19,AG18,AH17),
IMPRODUCT(+1,AD17,AE18,AF19,AG20,AH16),IMPRODUCT(-1,AD18,AE17,AF19,AG20,AH16),
IMPRODUCT(+1,AD19,AE17,AF18,AG20,AH16),IMPRODUCT(-1,AD17,AE19,AF18,AG20,AH16),
IMPRODUCT(+1,AD18,AE19,AF17,AG20,AH16),IMPRODUCT(-1,AD19,AE18,AF17,AG20,AH16),
IMPRODUCT(-1,AD20,AE17,AF18,AG19,AH16),IMPRODUCT(+1,AD17,AE20,AF18,AG19,AH16),
IMPRODUCT(-1,AD18,AE20,AF17,AG19,AH16),IMPRODUCT(+1,AD20,AE18,AF17,AG19,AH16),
IMPRODUCT(-1,AD17,AE18,AF20,AG19,AH16),IMPRODUCT(+1,AD18,AE17,AF20,AG19,AH16),
IMPRODUCT(+1,AD19,AE20,AF17,AG18,AH16),IMPRODUCT(-1,AD20,AE19,AF17,AG18,AH16),
IMPRODUCT(+1,AD17,AE19,AF20,AG18,AH16),IMPRODUCT(-1,AD19,AE17,AF20,AG18,AH16),
IMPRODUCT(+1,AD20,AE17,AF19,AG18,AH16),IMPRODUCT(-1,AD17,AE20,AF19,AG18,AH16),
IMPRODUCT(-1,AD18,AE19,AF20,AG17,AH16),IMPRODUCT(+1,AD19,AE18,AF20,AG17,AH16),
IMPRODUCT(-1,AD20,AE18,AF19,AG17,AH16),IMPRODUCT(+1,AD18,AE20,AF19,AG17,AH16),
IMPRODUCT(-1,AD19,AE20,AF18,AG17,AH16),IMPRODUCT(+1,AD20,AE19,AF18,AG17,AH16))

 

6x6以降は省略

 

2x2

各項目についてはほぼ前と同じだが違う点もある。

　元行列の起点をAC15からAD16に変更。

　行列式は元行列の１列目の真下ではなくその一つ左側。

　sag1は元行列と小行列式はみ出しシフト(仮)の１列目ではなく最後の列(最右列)同士の積。作業セルの取り込みも下から順に。

　imsub関数は使わなくなった。

　罫線の囲みが揃っているので行挿入列挿入がしやすくなったか。

　展開作業では検索=「),」置換=「,追加でかけるセル),」等も駆使して何とか乗り切った。

　検索=「+」置換=「_」、検索=「-」置換=「+」、検索=「_」置換=「-」により「improduct(-1,～)」を減らしたり。

 

逆行列の分母　(ここでは省略したが約分作業が別途必要だ。)

　AM22=　=TEXT(IMABS(AC23)^2,"#,##0")　(カンマ区切りテキストにしてみた)

 

行列式の共役複素数

　AM24=　=IMCONJUGATE(AC23)

 

逆行列の分子

　AM19=　=IMPRODUCT(AJ19,$AM$24)

　AM19:AN20=　AM19をコピペ

 

検算２

　AM16=　=IMSUM(IMPRODUCT($AD16,AM$19),IMPRODUCT($AE16,AM$20))

　AM16:AN17=　AM16をコピペ

 

3x3

 

7x7

 

余談１

　リブレオフィスの6x6の行列式では作業セルが２つから３つに増えてしまった。リブレオフィスの違うバージョンを使ったからかも知れないが。要調査。

 

配列・ベクトル・行列(10)

複素数行列の行列式と余因子行列と逆行列を求める数式の作り方

エクセル、グーグルドライブ、リブレオフィスにて行列式と余因子行列と逆行列をmdeterm関数やminverse関数に頼らずに組み立ててみる。

 

最初は実数版を作ってみたが複素数版も手間は大差なく手順もほぼ同じなので前者は省略。(虚部を0にすれば実数の結果を見ることは可能にした。)

 

最初に複素数対応の行列式の数式

 

2x2　元行列=　AC15:AD16　

AC22=　=IMSUB(improduct(AC15,AD16),improduct(AC16,AD15))

 

3x3　元行列=　AC15:AE17

AC24=　=imsum(
IMPRODUCT(AC15,IMSUB(improduct(AD16,AE17),improduct(AD17,AE16))),
IMPRODUCT(AC16,IMSUB(improduct(AD17,AE15),improduct(AD15,AE17))),
IMPRODUCT(AC17,IMSUB(improduct(AD15,AE16),improduct(AD16,AE15))))

 

4x4　元行列=　AC15:AF18　(「+1,」はエクセルでは「1,」にされる)

AC26=　=imsum(
IMPRODUCT(+1,AC15,imsum(
IMPRODUCT(AD16,IMSUB(improduct(AE17,AF18),improduct(AE18,AF17))),
IMPRODUCT(AD17,IMSUB(improduct(AE18,AF16),improduct(AE16,AF18))),
IMPRODUCT(AD18,IMSUB(improduct(AE16,AF17),improduct(AE17,AF16))))),
IMPRODUCT(-1,AC16,imsum(
IMPRODUCT(AD17,IMSUB(improduct(AE18,AF15),improduct(AE15,AF18))),
IMPRODUCT(AD18,IMSUB(improduct(AE15,AF17),improduct(AE17,AF15))),
IMPRODUCT(AD15,IMSUB(improduct(AE17,AF18),improduct(AE18,AF17))))),
IMPRODUCT(+1,AC17,imsum(
IMPRODUCT(AD18,IMSUB(improduct(AE15,AF16),improduct(AE16,AF15))),
IMPRODUCT(AD15,IMSUB(improduct(AE16,AF18),improduct(AE18,AF16))),
IMPRODUCT(AD16,IMSUB(improduct(AE18,AF15),improduct(AE15,AF18))))),
IMPRODUCT(-1,AC18,imsum(
IMPRODUCT(AD15,IMSUB(improduct(AE16,AF17),improduct(AE17,AF16))),
IMPRODUCT(AD16,IMSUB(improduct(AE17,AF15),improduct(AE15,AF17))),
IMPRODUCT(AD17,IMSUB(improduct(AE15,AF16),improduct(AE16,AF15))))))

 

5x5　元行列=　AC15:AG19

AC28=　=imsum(
IMPRODUCT(AC15,imsum(
IMPRODUCT(+1,AD16,imsum(
IMPRODUCT(AE17,IMSUB(improduct(AF18,AG19),improduct(AF19,AG18))),
IMPRODUCT(AE18,IMSUB(improduct(AF19,AG17),improduct(AF17,AG19))),
IMPRODUCT(AE19,IMSUB(improduct(AF17,AG18),improduct(AF18,AG17))))),
IMPRODUCT(-1,AD17,imsum(
IMPRODUCT(AE18,IMSUB(improduct(AF19,AG16),improduct(AF16,AG19))),
IMPRODUCT(AE19,IMSUB(improduct(AF16,AG18),improduct(AF18,AG16))),
IMPRODUCT(AE16,IMSUB(improduct(AF18,AG19),improduct(AF19,AG18))))),
IMPRODUCT(+1,AD18,imsum(
IMPRODUCT(AE19,IMSUB(improduct(AF16,AG17),improduct(AF17,AG16))),
IMPRODUCT(AE16,IMSUB(improduct(AF17,AG19),improduct(AF19,AG17))),
IMPRODUCT(AE17,IMSUB(improduct(AF19,AG16),improduct(AF16,AG19))))),
IMPRODUCT(-1,AD19,imsum(
IMPRODUCT(AE16,IMSUB(improduct(AF17,AG18),improduct(AF18,AG17))),
IMPRODUCT(AE17,IMSUB(improduct(AF18,AG16),improduct(AF16,AG18))),
IMPRODUCT(AE18,IMSUB(improduct(AF16,AG17),improduct(AF17,AG16))))))),
IMPRODUCT(AC16,imsum(
IMPRODUCT(+1,AD17,imsum(
IMPRODUCT(AE18,IMSUB(improduct(AF19,AG15),improduct(AF15,AG19))),
IMPRODUCT(AE19,IMSUB(improduct(AF15,AG18),improduct(AF18,AG15))),
IMPRODUCT(AE15,IMSUB(improduct(AF18,AG19),improduct(AF19,AG18))))),
IMPRODUCT(-1,AD18,imsum(
IMPRODUCT(AE19,IMSUB(improduct(AF15,AG17),improduct(AF17,AG15))),
IMPRODUCT(AE15,IMSUB(improduct(AF17,AG19),improduct(AF19,AG17))),
IMPRODUCT(AE17,IMSUB(improduct(AF19,AG15),improduct(AF15,AG19))))),
IMPRODUCT(+1,AD19,imsum(
IMPRODUCT(AE15,IMSUB(improduct(AF17,AG18),improduct(AF18,AG17))),
IMPRODUCT(AE17,IMSUB(improduct(AF18,AG15),improduct(AF15,AG18))),
IMPRODUCT(AE18,IMSUB(improduct(AF15,AG17),improduct(AF17,AG15))))),
IMPRODUCT(-1,AD15,imsum(
IMPRODUCT(AE17,IMSUB(improduct(AF18,AG19),improduct(AF19,AG18))),
IMPRODUCT(AE18,IMSUB(improduct(AF19,AG17),improduct(AF17,AG19))),
IMPRODUCT(AE19,IMSUB(improduct(AF17,AG18),improduct(AF18,AG17))))))),
IMPRODUCT(AC17,imsum(
IMPRODUCT(+1,AD18,imsum(
IMPRODUCT(AE19,IMSUB(improduct(AF15,AG16),improduct(AF16,AG15))),
IMPRODUCT(AE15,IMSUB(improduct(AF16,AG19),improduct(AF19,AG16))),
IMPRODUCT(AE16,IMSUB(improduct(AF19,AG15),improduct(AF15,AG19))))),
IMPRODUCT(-1,AD19,imsum(
IMPRODUCT(AE15,IMSUB(improduct(AF16,AG18),improduct(AF18,AG16))),
IMPRODUCT(AE16,IMSUB(improduct(AF18,AG15),improduct(AF15,AG18))),
IMPRODUCT(AE18,IMSUB(improduct(AF15,AG16),improduct(AF16,AG15))))),
IMPRODUCT(+1,AD15,imsum(
IMPRODUCT(AE16,IMSUB(improduct(AF18,AG19),improduct(AF19,AG18))),
IMPRODUCT(AE18,IMSUB(improduct(AF19,AG16),improduct(AF16,AG19))),
IMPRODUCT(AE19,IMSUB(improduct(AF16,AG18),improduct(AF18,AG16))))),
IMPRODUCT(-1,AD16,imsum(
IMPRODUCT(AE18,IMSUB(improduct(AF19,AG15),improduct(AF15,AG19))),
IMPRODUCT(AE19,IMSUB(improduct(AF15,AG18),improduct(AF18,AG15))),
IMPRODUCT(AE15,IMSUB(improduct(AF18,AG19),improduct(AF19,AG18))))))),
IMPRODUCT(AC18,imsum(
IMPRODUCT(+1,AD19,imsum(
IMPRODUCT(AE15,IMSUB(improduct(AF16,AG17),improduct(AF17,AG16))),
IMPRODUCT(AE16,IMSUB(improduct(AF17,AG15),improduct(AF15,AG17))),
IMPRODUCT(AE17,IMSUB(improduct(AF15,AG16),improduct(AF16,AG15))))),
IMPRODUCT(-1,AD15,imsum(
IMPRODUCT(AE16,IMSUB(improduct(AF17,AG19),improduct(AF19,AG17))),
IMPRODUCT(AE17,IMSUB(improduct(AF19,AG16),improduct(AF16,AG19))),
IMPRODUCT(AE19,IMSUB(improduct(AF16,AG17),improduct(AF17,AG16))))),
IMPRODUCT(+1,AD16,imsum(
IMPRODUCT(AE17,IMSUB(improduct(AF19,AG15),improduct(AF15,AG19))),
IMPRODUCT(AE19,IMSUB(improduct(AF15,AG17),improduct(AF17,AG15))),
IMPRODUCT(AE15,IMSUB(improduct(AF17,AG19),improduct(AF19,AG17))))),
IMPRODUCT(-1,AD17,imsum(
IMPRODUCT(AE19,IMSUB(improduct(AF15,AG16),improduct(AF16,AG15))),
IMPRODUCT(AE15,IMSUB(improduct(AF16,AG19),improduct(AF19,AG16))),
IMPRODUCT(AE16,IMSUB(improduct(AF19,AG15),improduct(AF15,AG19))))))),
IMPRODUCT(AC19,imsum(
IMPRODUCT(+1,AD15,imsum(
IMPRODUCT(AE16,IMSUB(improduct(AF17,AG18),improduct(AF18,AG17))),
IMPRODUCT(AE17,IMSUB(improduct(AF18,AG16),improduct(AF16,AG18))),
IMPRODUCT(AE18,IMSUB(improduct(AF16,AG17),improduct(AF17,AG16))))),
IMPRODUCT(-1,AD16,imsum(
IMPRODUCT(AE17,IMSUB(improduct(AF18,AG15),improduct(AF15,AG18))),
IMPRODUCT(AE18,IMSUB(improduct(AF15,AG17),improduct(AF17,AG15))),
IMPRODUCT(AE15,IMSUB(improduct(AF17,AG18),improduct(AF18,AG17))))),
IMPRODUCT(+1,AD17,imsum(
IMPRODUCT(AE18,IMSUB(improduct(AF15,AG16),improduct(AF16,AG15))),
IMPRODUCT(AE15,IMSUB(improduct(AF16,AG18),improduct(AF18,AG16))),
IMPRODUCT(AE16,IMSUB(improduct(AF18,AG15),improduct(AF15,AG18))))),
IMPRODUCT(-1,AD18,imsum(
IMPRODUCT(AE15,IMSUB(improduct(AF16,AG17),improduct(AF17,AG16))),
IMPRODUCT(AE16,IMSUB(improduct(AF17,AG15),improduct(AF15,AG17))),
IMPRODUCT(AE17,IMSUB(improduct(AF15,AG16),improduct(AF16,AG15))))))))

6x6以降は省略

 

リブレオフィスカルクにて

2x2

(リブレオフィスでも複素数のmdeterm関数はエラーのはずだが0になっているのは自動計算をオフにしてあるからかも知れない)

 

下準備　元行列は乱数

　AB11:AD13は乱数用。元行列は乱数で変化する。

 

下準備　はみ出し部分

　AC17=　=AC15　　AC17:AD18=　AC17をコピペ　　(以前は配列数式で元行列の範囲丸ごと参照していたが(AC17:AD18=　{元行列}等)行挿入列挿入時に不便だったので不採用。配列数式大好きなのに今回は珍しく１か所も使ってない。)

　AE15=　=AC15　　AE15:AF18=　AE15をコピペ

 

小行列式はみ出しシフト(仮)

　AD21=　=AD16　(これは元行列の１行１列の小行列式)

　AG21=　AD21をカトペ　(「カトペ」はカット＆ペーストの意味)

　再度　AD21=　=AD16　(AD21を数式コピペでも良い。ただ作業セルが増えると面倒になる)

　AG21:AH22=　AG21をコピペ

 

余因子行列転置前(仮)

　小行列式はみ出しシフト(仮)を符号調整する

　AJ21=　=IMPRODUCT(AG21,-1^(ROW(A1)+COLUMN(A1)))

　AJ21:AK22=　AJ21をコピペ

 

余因子行列

　余因子行列転置前(仮)を転置する

　AJ18=　=OFFSET($AJ$21,COLUMN(A1)-1,ROW(A1)-1)

　AJ18:AK19=　AJ18をコピペ　(元々はtranspose関数を使っていたが行挿入列挿入時に不便だったので代替手段を考えてみた。)

 

検算

　元行列と余因子行列の行列積を求める

　AJ15=　=IMSUM(IMPRODUCT($AC15,AJ$18),IMPRODUCT($AD15,AJ$19))

　AJ15:AK16=　AJ15をコピペ

　左上から右下の対角線が行列式になって他が0になるととりあえずひと安心。

 

小行列式循環シフト(仮)

　AD22=　AD21をコピペ

　AD22を選択して「メニュー」「編集」「検索と置換」

　検索=17　置換=15　現在の選択範囲のみにチェック　検索場所は「数式」　で「すべて置換」

　(検索の17は下はみ出し部分の１行目の行番号)

　AE21:AE22=　AD21:AD22をコピペ

　AE21:AE22を選択して「メニュー」「編集」「検索と置換」

　検索=AE　置換=AC　現在の選択範囲のみにチェック　検索場所は「数式」　で「すべて置換」

　(検索のAEは右はみ出し部分の１列目の列名)

 

sag1　行列式前準備１

　符号調整と元行列の１列目と小行列式はみ出しシフト(仮)の１列目の積の和(内積)を求める

　AA20=　=IMSUM(AA21:AA22)

　AA21=　=IMPRODUCT(+1,AC15,AG21)

　AA22=　=IMPRODUCT(-1,AC16,AG22)

 

sag2　行列式前準備２

　sag1の数式の最後の項を元行列のセルだけの数式に変更する

　AB20=　AA20をコピペ

　AB21=　AA21の数式をコピペ(セルのコピペではなくセル内容のコピペ)

　AD21の数式をコピー(「=」を除外)、AB21「=IMPRODUCT(AC15,AG21)」の「AG21」を選択してペースト

　AD22の数式をコピー(「=」を除外)、AB22「=IMPRODUCT(AC16,AG22)」の「AG22」を選択してペースト

 

行列式

　sag2の全体を１つのセルにまとめる

　AC22=　=IMSUM(0,0)

　AB21の数式をコピー(「=」を除外)、AC22の１つ目の0を選択してペースト

　AB22の数式をコピー(「=」を除外)、AC22の２つ目の0を選択してペースト

　AC22=　=IMSUM(IMPRODUCT(+1,AC15,AD16),IMPRODUCT(-1,AC16,AD15))

この説明では-1倍を使ったが実際には2x2のみ「行列式」ではIMSUBを使った。4x4以上では前者の方が()のネストが少ないし順番変更も不要なのでシンプルかと思う。

　AC22=　=IMSUB(IMPRODUCT(AC15,AD16),IMPRODUCT(AC16,AD15))

　作業で乱れた罫線は必要に応じて書き直す。

 

3x3

 

下準備　シートのコピーと3x3化

　現シート「2x2」のタブを右クリックして「シートを移動またはコピー」

　「コピー」「- 末尾に挿入 -」新しい名前を「3x3」等としてコピー

　22行を選択して行番号の所で右クリック「上に行を挿入」

　19行の上と16行の上にも行挿入　AK列、AH列、AD列の左に列挿入

 

下準備　その他

　元行列は乱数が必要なら同じ数式をコピペすればよい。はみ出し部分も空白部分を埋めるだけ。

　小行列式循環シフト(仮)と小行列式はみ出しシフト(仮)の内容とsag2の内容を消去する

 

小行列式はみ出しシフト(仮)

　「3x3」シートのAD23=　「2x2」シートのAC22をコピペ

　「3x3」シートのAH23=　AD23をカトペ

　「3x3」シートのAD23=　「2x2」シートのAC22を再度コピペ

　(以下「3x3」シート)

　AH23:AJ25=　AH23をコピペ

 

余因子行列転置前(仮)

　小行列式はみ出しシフト(仮)を符号調整する(奇数行列では調整不要)

　AL23=　=IMPRODUCT(AH23,1^(ROW(A1)+COLUMN(A1)))　(奇数行列の場合「=AH23」だけで良いがこの形で残しておけば偶数行列での手間が軽減する)

　AL23:AN25=　AL23をコピペ

 

余因子行列

　余因子行列転置前(仮)を転置する

　AL19:AN21=　AL19をコピペ

 

検算

　元行列と余因子行列の行列積を求める

　AL15=　=IMSUM(IMPRODUCT($AC15,AL$19),IMPRODUCT($AD15,AL$20),IMPRODUCT($AE15,AL$21))

　(元の数式の最後の「,improduct(～,～)」を数式コピー、文末の「)」の前にカーソルを置いてペーストしてから修正するのが楽)

　AL15:AN17=　AL15をコピペ

 

小行列式循環シフト(仮)

　AD24=　AD23をコピペ

　AD24を選択して「メニュー」「編集」「検索と置換」

　検索=18　置換=15　現在の選択範囲のみにチェック　検索場所は「数式」　で「すべて置換」

　AD25=　AD24をコピペ

　AD25を選択して「メニュー」「編集」「検索と置換」

　検索=18　置換=15　現在の選択範囲のみにチェック　検索場所は「数式」　で「すべて置換」

　(4x4以降もこの作業を１行ずつ一番下まで行う)

　AE23:AE25=　AD23:AD25をコピペ

　AE23:AE25を選択して「メニュー」「編集」「検索と置換」

　検索=AF　置換=AC　現在の選択範囲のみにチェック　検索場所は「数式」　で「すべて置換」

　AF23:AF25=　AE23:AE25をコピペ

　AF23:AF25を選択して「メニュー」「編集」「検索と置換」

　検索=AF　置換=AC　現在の選択範囲のみにチェック　検索場所は「数式」　で「すべて置換」

　(4x4以降もこの作業を１列ずつ一番右まで行う)

 

sag1　行列式前準備１

　元行列の１列目と小行列式はみ出しシフト(仮)の１列目の積の和(内積)を求める(奇数行列では符号調整なし)

　AA22=　=IMSUM(AA23:AA25)　(変更なしで使えると思う)

　AA23=　=IMPRODUCT(AC15,AH23)　(変更なしで使えると思う)

　AA24:AA25=　AA23をコピペ　(AA24を埋めるだけで使えると思う)

 

sag2　行列式前準備２

　sag1の数式の最後の項を元行列のセルだけの数式に変更する

　AB22=　AA22をコピペ

　AB23=　AA23の数式をコピペ(セルのコピペではなくセル内容のコピペ)

　AD23の数式をコピー(「=」を除外)、AB23「=IMPRODUCT(AC15,AH23)」の「AH23」を選択してペースト

　AD24の数式をコピー(「=」を除外)、AB24「=IMPRODUCT(AC16,AH24)」の「AH24」を選択してペースト

　AD25の数式をコピー(「=」を除外)、AB25「=IMPRODUCT(AC17,AH25)」の「AH25」を選択してペースト

 

行列式

　sag2の全体を１つのセルにまとめる

　AC24=　=IMSUM(0,0,0)

　AB23の数式をコピー(「=」を除外)、AC24の１つ目の0を選択してペースト

　AB24の数式をコピー(「=」を除外)、AC24の２つ目の0を選択してペースト

　AB25の数式をコピー(「=」を除外)、AC24の３つ目の0を選択してペースト

　AC24=

　=IMSUM(
IMPRODUCT(AC15,IMSUB(IMPRODUCT(AD16,AE17),IMPRODUCT(AD17,AE16))),
IMPRODUCT(AC16,IMSUB(IMPRODUCT(AD17,AE15),IMPRODUCT(AD15,AE17))),
IMPRODUCT(AC17,IMSUB(IMPRODUCT(AD15,AE16),IMPRODUCT(AD16,AE15))))

 

4x4以降

　2x2を元に3x3を作ったのと同様、3x3を元に4x4を作る。偶数の場合は行列式の直前のsag1作成で符号が「+」「-」交互になるが奇数の場合は全て「+」となる。

 

「8x8」まで作ってみた。行列式を求めるために集計セルの他に作業セルが2*7*8=112セル必要だった。

小行列式はみ出しシフト(仮)　(作業セルがある場合)

　集計セル群と各作業セル間は１セル開けるようにした。

　下方向には作業セルがあるので簡単にコピペできない。

　まず1つ小さい行列(5x5なら4x4)の行列式を持ってきたら作業セルを枠外下方向に分離(カトペ)。集計セルを右１下１にコピペしてから作業セルも右１下１にコピペ。そのまま真下に隣接するようにカトペ。そのまま左１にコピペ。集計セルも左１にコピペ。この作業を一番下までやる。

　右方向へは集計セルも作業セルも最初の列をそのままコピペ。

 

小行列式循環シフト(仮)　(作業セルがある場合)

　はみ出しシフトとほぼ同様だが下に１つコピペするごとに検索と置換で１マスはみ出した行を「15」行に変えていく。

　右に１つコピペするごとに検索と置換で１マスはみ出した列を「AC」列に変えていく。

　次で使う訳ではないので２列目以降は作らなくても支障はない。

 

行列式　(作業セルがある場合)

　集計セルと各作業セル間は間をあけないように移動した。次で利用するときのために。

 

エクセルにて　余談２や余談４の後に配置を変更して逆行列も求めてからスクショした。逆行列は手作業で約分する必要があるようだ。

2x2

 

3x3

 

9x9

 

グーグルドライブにて　余談２や余談４の後に配置を変更して逆行列も求めてからスクショした。逆行列は手作業で約分する必要があるようだ。

2x2

 

3x3

 

9x9

 

 

余談１　IMSUBを使う？

　IMSUBかマイナス1倍かの問題でどちらかに統一した方が良いかと悩んだが、2x2のみIMSUBで4x4・6x6・8x8～は-1倍にするのが良さそうだ。

 

余談２　実数行列用のmdeterm関数等

　2x2のAI14、3x3のAK14には「=mdeterm(value(元行列))」が入っている。行列が大きくなると隠れてしまうので「AF11=　＝AI14」等と入れたがそこにmdeterm式を入れても良い。複素数の場合はエラー値が表示される。

　AB14は元々はsumproductで元行列の１列目と余因子行列転置前(仮)の１列目の積の合計を求めていたが、複素数では機能しないので個別の積を合計した。

　3x3の例　AB14=　=IMSUM(IMPRODUCT(AC15,AL23),IMPRODUCT(AC16,AL24),IMPRODUCT(AC17,AL25))

 

余談３　乱数作成

　AC15=　=complex(round((rand()*($AC$12-$AC$13)),0)+$AC$13,round((rand()*($AD$12-$AD$13)),0)+$AD$13)

 

余談４　逆行列

　imdiv関数を使って余因子行列を行列式で割れば良い。自分は小数になるのが嫌なので次のようにする。

　2x2の場合の逆行列の分子

　　AJ24=　=IMPRODUCT(AJ18,IMCONJUGATE($AC$22))　(IMCONJUGATE関数は共役複素数を求める関数)

　　AJ24:AK25=　AJ24をコピペ

　逆行列の分母

　　AG24=　=IMABS(AC22)^2　(IMABS関数は複素数の絶対値を求める関数)

 

余談５

　カトペ(移動)やコピペ(コピー)を多用したため、たとえば小行列式循環シフト(仮)なのに小行列式はみ出しシフト(仮)の数式が残ってしまう等のミスがあるかもしれない。しかし結果の行列式や検算が乱数変化で異常なしなので大きな問題はないと思う。ただ「行列式」にミスが混入すると次で確実に困る事になる。

 

余談６

　はみ出し部分は最後の１行１列は無駄なので削っても良い。

　小行列式はみ出しシフト(仮)は残さなくても良いがそうすると小行列式循環シフト(仮)作成中は「検算」等が乱れるためミスが分かりにくくなる。

 

余談７

　エクセルの9x9の行列式を求める際、集計セルが１つに入りきらなかったためその集計セルを作業セルとして更に上位の集計セルが必要になった。「作業セル　3033=6*7*8*9+9」以前作った時は作業セルに数式を追加していったからそんな事態にはなりにくかったと思うが、今回は簡略化のため集計セルに数式を追加していった。あとたとえば作業セル３つを２つに詰め込む事も可能かも知れないがあまりにも手間が大きいのでそれは考えない事にした。

 

作業セル数	エクセル	リブレオフィス	グーグルドライブ
6x6	6	2	0
7x7	6*7=42	2*7=14	7
8x8	6*7*8=336	2*7*8=112	7*8=56
9x9	6*7*8*9+9 =3033	未作成	7*8*9 =504

6x6の行列式を求める際、5x5の行列式が１セルに

　エクセルでは１つしか入らなかった。

　リブレオフィスでは３つ入った。４つ入ったかどうかは忘れたが６つを２セルに分けた。

　グーグルドライブでは全部(６つ)入った。

 

余談８

　今回の行列式は配列・ベクトル・行列(6)の順列列挙で言うと「２．頭固定小行列循環シフト型(仮)」で作った形になる。「123 132」ではなく因数分解されて「1(23 32)」のようになっているから同じではないが。「４．尾固定小行列循環シフト型(仮)」で作るとどうなるのか、ちょっと気になる。

 

余談９

　試運転版(数日間は凡ミス等を修正する予定。)　r070509試運転版解除。だが以後修正しないという意味ではない。

 

配列・ベクトル・行列(9)

小行列の「はみ出しシフト(仮)」や「循環シフト(仮)」を説明してみた。

元行列 = 黒枠内

１．小行列について

　C6=　=J2　　F6=　=M2

　C8=　=J4　　F8=　=M4

　他略

２．小行列はみ出しシフト(仮)について

　H19=　=O7　　F19=　=M7

　H16=　=O4　　F16=　=M4

　他略

３．小行列循環シフト(仮)について

　Q19=　=J2　　O19=　=M2

　Q16=　=J4　　O16=　=M4

　他略

４．小行列式について

　A_23　難易度がやや高いので今回は省略

５．小行列式はみ出しシフト(仮)について

　A_23'　=mdeterm(f16:i19)

６．小行列式循環シフト(仮)について

　A_23'　=mdeterm(o16:i19)

７．余因子

　a_32　今回は省略

８．余因子はみ出しシフト(仮)

　a_32　=mdeterm(f16:i19) * 符号調整

９．余因子循環シフト(仮)

　a_32　=mdeterm(o16:i19) * 符号調整

１０．余因子の７～９は符号調整により３つとも同じ値にする

　符号調整　元行列が奇数*奇数では無調整、偶数*偶数では-1^(i+j)倍となる

 

エクセル的には「～はみ出しシフト(仮)」と「～循環シフト(仮)」は異なる式だが数学的には同じだから両方を区別せず「～シフト(仮)」と言うこともあるかも知れない。

 

seriessum関数でべき級数を求めてみる

 seriessum関数の正しい使い道を模索してみた。

r070414　主にグーグルドライブで作成したがダウンロードしてエクセルとリブレオフィスの分も追加した。都合により数学的な式の画像とセル式を分離した。

 

(r070416　式１～５の画像が少々ぼけていたので変更)

「べき級数　円周率」で検索

　https://www.procrasist.com/entry/pi_1　(PROCRASISTさん)

(式１)

 

(式２)

 

(式３)　ラマヌジャン式

 

「ラマヌジャン」で検索

　https://mathlog.info/articles/3583　(Mathlogさん)

(式４)　ラマヌジャンの円周率公式

 

(式５)　Chudnovskyの公式

 

エクセルにて

(表)

(式１)

　c5=　{=(SERIESSUM(1,1,1,ROW(INDIRECT("A1:A"&B5))^-2)*6)^0.5}　(配列数式)

　d5=　{=SUM((ROW(INDIRECT("A1:A"&B5))^-2)*6)^0.5}　(配列数式)

(式２)

　c6=　{=SERIESSUM(-1,0,1,1/(2*(ROW(INDIRECT("A1:A"&B6))-1)+1))*4}　(配列数式)

　d6=　{=SUM((-1)^(ROW(INDIRECT("A1:A"&B6))-1)/(2*(ROW(INDIRECT("A1:A"&B6))-1)+1))*4}　(配列数式)

(式３)

　c7=　{=SERIESSUM(-1*882^-2*4^-4,0,1,FACT((ROW(INDIRECT("A1:A"&B7))-1)*4)*(1123+21460*(ROW(INDIRECT("A1:A"&B7))-1))*FACT(ROW(INDIRECT("A1:A"&B7))-1)^-4/882)^-1*4}　(配列数式)

　d7=　{=SUM((-1)^(ROW(INDIRECT("A1:A"&B7))-1)*FACT(4*(ROW(INDIRECT("A1:A"&B7))-1))*(1123+21460*(ROW(INDIRECT("A1:A"&B7))-1))*882^-(2*(ROW(INDIRECT("A1:A"&B7))-1)+1)*(4^(ROW(INDIRECT("A1:A"&B7))-1)*FACT(ROW(INDIRECT("A1:A"&B7))-1))^-4)^-1*4}　(配列数式)

(式４)

　c8=　{=(SERIESSUM(396^-4,0,1,FACT(4*(ROW(INDIRECT("A1:A"&B8))-1))*(FACT(ROW(INDIRECT("A1:A"&B8))-1))^-4*(26390*(ROW(INDIRECT("A1:A"&B8))-1)+1103))*2^(3/2)*99^-2)^-1}　(配列数式)

　d8=　{=(SUM(FACT(4*(ROW(INDIRECT("A1:A"&B8))-1))*FACT(ROW(INDIRECT("A1:A"&B8))-1)^-4*(26390*(ROW(INDIRECT("A1:A"&B8))-1)+1103)*396^(-4*(ROW(INDIRECT("A1:A"&B8))-1)))*2^(3/2)*99^-2)^-1}　(配列数式)

(式５)

　c9=　{=(12*SERIESSUM(640320^-3,0,1,-1^(ROW(INDIRECT("A1:A"&B9))-1)*FACT(6*(ROW(INDIRECT("A1:A"&B9))-1))*FACT(3*(ROW(INDIRECT("A1:A"&B9))-1))^-1*FACT(ROW(INDIRECT("A1:A"&B9))-1)^-3*(545140134*(ROW(INDIRECT("A1:A"&B9))-1)+13591409)*640320^-(3/2)))^-1}　(配列数式)

　d9=　{=(12*SUM(-1^(ROW(INDIRECT("A1:A"&B9))-1)*FACT(6*(ROW(INDIRECT("A1:A"&B9))-1))*FACT(3*(ROW(INDIRECT("A1:A"&B9))-1))^-1*FACT(ROW(INDIRECT("A1:A"&B9))-1)^-3*(545140134*(ROW(INDIRECT("A1:A"&B9))-1)+13591409)*640320^-(3*(ROW(INDIRECT("A1:A"&B9))-1)+3/2)))^-1}　(配列数式)

　f2=　=PI()

　e5:e9=　{=C5:C9-D5:D9}　(配列数式)

　f5:g9=　{=C5:D9-$F$2}　(配列数式)

 

グーグルドライブのスプレッドシートにて

(表)

(式１)

　c5=　=ArrayFormula((SERIESSUM(1,1,1,ROW(indirect("A1:A"&B5))^-2)*6)^0.5)

　d5=　=ArrayFormula(SUM((ROW(indirect("A1:A"&B5))^-2)*6)^0.5)

(式２)

　c6=　=ArrayFormula(SERIESSUM(-1,0,1,1/(2*(ROW(indirect("A1:A"&B6))-1)+1))*4)

　d6=　=ArrayFormula(sum((-1)^(ROW(indirect("A1:A"&B6))-1)/(2*(ROW(indirect("A1:A"&B6))-1)+1))*4)

(式３)

　c7=　=ArrayFormula(SERIESSUM(-1*882^-2*4^-4,0,1,FACT((ROW(indirect("A1:A"&B7))-1)*4)*(1123+21460*(ROW(indirect("A1:A"&B7))-1))*FACT(ROW(indirect("A1:A"&B7))-1)^-4/882)^-1*4)

　d7=　=ArrayFormula(SUM((-1)^(row(indirect("A1:A"&B7))-1)*fact(4*(row(indirect("A1:A"&B7))-1))*(1123+21460*(row(indirect("A1:A"&B7))-1))*882^-(2*(row(indirect("A1:A"&B7))-1)+1)*(4^(row(indirect("A1:A"&B7))-1)*fact(row(INDIRECT("A1:A"&B7))-1))^-4)^-1)*4

(式４)

　c8=　=ArrayFormula(seriessum(396^-4,0,1,fact(4*(row(indirect("A1:A"&B8))-1))*(fact(row(indirect("A1:A"&B8))-1))^-4*(26390*(row(INDIRECT("A1:A"&B8))-1)+1103))*2^(3/2)*99^-2)^-1

　d8=　=ArrayFormula(sum(fact(4*(row(indirect("A1:A"&B8))-1))*fact(row(indirect("A1:A"&B8))-1)^-4*(26390*(row(indirect("A1:A"&B8))-1)+1103)*396^(-4*(row(indirect("A1:A"&B8))-1)))*2^(3/2)*99^-2)^-1

(式５)

　c9=　=ArrayFormula(12*seriessum(640320^-3,0,1,-1^(row(indirect("A1:A"&B9))-1)*fact(6*(row(indirect("A1:A"&B9))-1))*fact(3*(row(indirect("A1:A"&B9))-1))^-1*fact(row(indirect("A1:A"&B9))-1)^-3*(545140134*(row(indirect("A1:A"&B9))-1)+13591409)*640320^-(3/2)))^-1

　d9=　=ArrayFormula(12*sum(-1^(row(indirect("A1:A"&B9))-1)*fact(6*(row(indirect("A1:A"&B9))-1))*fact(3*(row(indirect("A1:A"&B9))-1))^-1*fact(row(indirect("A1:A"&B9))-1)^-3*(545140134*(row(indirect("A1:A"&B9))-1)+13591409)*640320^-(3*(row(indirect("A1:A"&B9))-1)+3/2)))^-1

　f2=　=pi()

　e5:e9=　=ArrayFormula(C5:C9-D5:D9)

　f5:g9=　=ArrayFormula(C5:D9-$F$2)

 

リブレオフィスにて

(表)

(式１)

　c5=　{=(SERIESSUM(1,1,1,ROW(INDIRECT("A1:A"&B5))^-2)*6)^0.5}　(配列数式)

　d5=　{=SUM((ROW(INDIRECT("A1:A"&B5))^-2)*6)^0.5}　(配列数式)

(式２)

　c6=　{=SERIESSUM(-1,0,1,1/(2*(ROW(INDIRECT("A1:A"&B6))-1)+1))*4}　(配列数式)

　d6=　{=SUM((-1)^(ROW(INDIRECT("A1:A"&B6))-1)/(2*(ROW(INDIRECT("A1:A"&B6))-1)+1))*4}　(配列数式)

(式３)

　c7=　{=SERIESSUM(-1*882^-2*4^-4,0,1,FACT((ROW(INDIRECT("A1:A"&B7))-1)*4)*(1123+21460*(ROW(INDIRECT("A1:A"&B7))-1))*FACT(ROW(INDIRECT("A1:A"&B7))-1)^-4/882)^-1*4}　(配列数式)

　d7=　{=SUM((-1)^(ROW(INDIRECT("A1:A"&B7))-1)*FACT(4*(ROW(INDIRECT("A1:A"&B7))-1))*(1123+21460*(ROW(INDIRECT("A1:A"&B7))-1))*882^-(2*(ROW(INDIRECT("A1:A"&B7))-1)+1)*(4^(ROW(INDIRECT("A1:A"&B7))-1)*FACT(ROW(INDIRECT("A1:A"&B7))-1))^-4)^-1*4}　(配列数式)

(式４)

　c8=　{=(SERIESSUM(396^-4,0,1,FACT(4*(ROW(INDIRECT("A1:A"&B8))-1))*(FACT(ROW(INDIRECT("A1:A"&B8))-1))^-4*(26390*(ROW(INDIRECT("A1:A"&B8))-1)+1103))*2^(3/2)*99^-2)^-1}　(配列数式)

　d8=　{=(SUM(FACT(4*(ROW(INDIRECT("A1:A"&B8))-1))*FACT(ROW(INDIRECT("A1:A"&B8))-1)^-4*(26390*(ROW(INDIRECT("A1:A"&B8))-1)+1103)*396^(-4*(ROW(INDIRECT("A1:A"&B8))-1)))*2^(3/2)*99^-2)^-1}　(配列数式)

(式５)

　c9=　{=(12*SERIESSUM(640320^-3,0,1,-1^(ROW(INDIRECT("A1:A"&B9))-1)*FACT(6*(ROW(INDIRECT("A1:A"&B9))-1))*FACT(3*(ROW(INDIRECT("A1:A"&B9))-1))^-1*FACT(ROW(INDIRECT("A1:A"&B9))-1)^-3*(545140134*(ROW(INDIRECT("A1:A"&B9))-1)+13591409)*640320^-(3/2)))^-1}　(配列数式)

　d9=　{=(12*SUM(-1^(ROW(INDIRECT("A1:A"&B9))-1)*FACT(6*(ROW(INDIRECT("A1:A"&B9))-1))*FACT(3*(ROW(INDIRECT("A1:A"&B9))-1))^-1*FACT(ROW(INDIRECT("A1:A"&B9))-1)^-3*(545140134*(ROW(INDIRECT("A1:A"&B9))-1)+13591409)*640320^-(3*(ROW(INDIRECT("A1:A"&B9))-1)+3/2)))^-1}　(配列数式)

　f2=　=PI()

　e5:e9=　{=C5:C9-D5:D9}　(配列数式)

　f5:g9=　{=C5:D9-$F$2}　(配列数式)

 

元頁(blogger)は2023.04.28 19:07公開

 

余談１
　n乗の部分が多いとseriessum関数が効率化に役立つ感じかな。

余談２　r070414　余談３にあわせて文頭を微修正
　表の最初の２行、個数が多くなると誤差も大きくなるのは仕方ない。200万にしても多分もっと増やしても正解には近づかない。あと関数で組み立てた式が本当に合っているのか検証できないのも困りもの(特に下の３つ)ｗ　今回はよく知られた定数だから良いが、多分本来の使い方はもっと未知の数字を求めるのだと思う。色々と大変そうだ。

余談３　r070414追記
　表の下の３行はもう一つプラスするとpi()との誤差がゼロになる(この表の精度の限界)。リブレオフィスだけは近似が良くない。と思ってよく見たら定数pi()の精度が悪いだけだった。

 

余談４　r070414追記

　グーグルドライブで作った次の式を

　　=arrayformula(<中略>)^-1

　エクセルやリブレオフィスのファイルとしてダウンロードすると

　　={<中略>^-1}　(配列数式)

　となってしまいエラーだったり異なる結果になったりする場合がある。

　　={(<中略>)^-1}　(配列数式)

　に修正することで正常化できた。arrayformula関数の()を頼ると失敗の元。次からは作成時に気を付けようと思う。

 

余談５　r070414追記

　また引っ越し先を考えなくてはならないようだ。

 

余談６　没画像　r070416

　　

４点を通る球面(３次元空間内)

エクセル等の表計算ソフトの逆行列関数を使って４点を通る球面を求める。(３次元空間内)

「３点を通る円(２次元平面内)」と全く同じ方法で求める。

まず求める球面の式を

　x^2+y^2+z^2+l*x+m*y+n*z+o=0　(一般形　r070719追記)

とすると

　l*x+m*y+n*z+o=-(x^2+y^2+z^2)

なので４点を

　p1(x1,y1,z1)、p2(x2,y2,z2)、p3(x3,y3,z3)、p4(x4,y4,z4)　(表のD4:F7)

とすると

　l*x1+m*y1+n*z1+o=-(x1^2+y1^2+z1^2)

　l*x2+m*y2+n*z2+o=-(x2^2+y2^2+z2^2)

　l*x3+m*y3+n*z3+o=-(x3^2+y3^2+z3^2)

　l*x4+m*y4+n*z4+o=-(x4^2+y4^2+z4^2)

となる。行列の式で書けば(libreofficecalcでmathをオブジェクトとして使用)(r070416画像を変更)

となり、逆行列を使えば一発で答えが出る。(r070416画像を変更)

逆行列と配列定数を使用してl,m,n,oを求める。配列定数を使って前半は４行３列の行列の１～３列目の順を逆にして４列目を追加して４列目を１とした。後半は２乗して横方向の和。

　D8:D11=　=ArrayFormula(mmult(MINVERSE(mmult(D4:F7,{0,0,1,0;0,1,0,0;1,0,0,0})+{0,0,0,1;0,0,0,1;0,0,0,1;0,0,0,1}),-MMULT(D4:F7^2,{1;1;1})))

次に中心座標と半径を求める。

　中心z D12=　=-D10/2

　中心y E12=　=-D9/2

　中心x F12=　=-D8/2

　半径^2 D13=　=((D8^2/4+D9^2/4+D10^2/4)-D11)

　半径 D14=　=((D8^2/4+D9^2/4+D10^2/4)-D11)^(1/2)

 

今は3Dグラフを描画する力はないのでzを固定した輪切りグラフ複数個で状況確認

z=-0.6

z=-1

z=-2

z=-4

z=-5.4

 

元頁(blogger)は2023.10.07 20:49公開

 

余談１　没画像　r070416