keisan.casio.jpに4x4行列の逆行列をUP!
です。どうやって解こうかと思ったんだが、直接Maximaで解いたものを入れるのが簡単かな、と思ってやってみた。こんな感じで...ちなみに5x5からはもう表示するのが複雑すぎてエラーがでてしまう...
matrix([((a22*a33-a23*a32)*a44+(a24*a32-a22*a34)*a43+(a23*a34-a24*a33)*a42)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a12*a33-a13*a32)*a44+(a14*a32-a12*a34)*a43+(a13*a34-a14*a33)*a42)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a12*a23-a13*a22)*a44+(a14*a22-a12*a24)*a43+(a13*a24-a14*a23)*a42)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a12*a23-a13*a22)*a34+(a14*a22-a12*a24)*a33+(a13*a24-a14*a23)*a32)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41)],[-((a21*a33-a23*a31)*a44+(a24*a31-a21*a34)*a43+(a23*a34-a24*a33)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a33-a13*a31)*a44+(a14*a31-a11*a34)*a43+(a13*a34-a14*a33)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a11*a23-a13*a21)*a44+(a14*a21-a11*a24)*a43+(a13*a24-a14*a23)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41)],[((a21*a32-a22*a31)*a44+(a24*a31-a21*a34)*a42+(a22*a34-a24*a32)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a11*a32-a12*a31)*a44+(a14*a31-a11*a34)*a42+(a12*a34-a14*a32)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a22-a12*a21)*a44+(a14*a21-a11*a24)*a42+(a12*a24-a14*a22)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a11*a22-a12*a21)*a34+(a14*a21-a11*a24)*a32+(a12*a24-a14*a22)*a31)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41)],[-((a21*a32-a22*a31)*a43+(a23*a31-a21*a33)*a42+(a22*a33-a23*a32)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a32-a12*a31)*a43+(a13*a31-a11*a33)*a42+(a12*a33-a13*a32)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a11*a22-a12*a21)*a43+(a13*a21-a11*a23)*a42+(a12*a23-a13*a22)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41)])
« 「タイタンの妖女」を読んだ。 | トップページ | 箱箱箱 »
「学問・資格」カテゴリの記事
- RF Weekly Digest (Gemini 3.1 Pro・Google AI Studio BuildによるAIで高周波・RF情報の週刊まとめアプリ)2026/4/5-4/12(2026.04.12)
- GLM-5.1(Ollamaから利用)でPythonのscikit-rfを使ってTouchstoneフォーマットのSパラメータファイルを読んでdB, 位相, スミスチャート, TDRを表示するGUIアプリを作ってもらった。5分など長く考えた後、Gemma 4:31bよりさらに出来が良く、思った通りのものができた。(2026.04.09)
- 高周波・RFニュース 2026年4月9日 QorvoがXバンドフェイズドアレイレーダ解説、atisがニューロシンボリックコグニティブアシスタントNESYのホワイトペーパー発行、Mini-CircuitsがADSを使ったアンプシミュレーション解説、Vishayがウィルキンソンデバイダ・コンバイナ発表(2026.04.09)
- 高周波・RFニュース 2026年4月8日 BroadcomがOFC2026での出展(200T AI時代)を解説、QualcommがMWCバルセロナ2026でのAIネイティブ6Gについて解説、iFixitがApple AirPods Max 2を分解、TDKが5GHz対応ノイズサプレッションフィルタ発表(2026.04.08)


コメント