![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?5a7ff87)
電磁波の数値解析手法の1つである有限積分法において
係数行列を求めるのは何故なのでしょうか。
ちなみに有限積分法はマクスウェル方程式の積分形とYee格子を組み合わせた数値解析手法です。
http://www.akita-nct.jp/yamamoto/study/thesis/20 …
を参考に、直交座標での式(2.17)が
Ey [i + 1][j][k] - Ey [i + 1][j][k + 1] - Ez [i + 1][j][k] + Ez [i + 1][j + 1][k] = Bx [i][j][k]
- Ex [i][j + 1][k] + Ex [i][j + 1][k + 1] + Ez [i][j + 1][k] - Ez [i + 1][j + 1][k] = By [i][j][k]
Ex [i][j][k+1] - Ex [i][j + 1][k + 1] - Ey [i][j][k + 1] + Ey [i + 1][j][k + 1] = Bz [i][j][k]
※ただし[]内の添え字は格子点の座標を表す
となることまでは理解したつもりです。
このときの係数行列の求め方も分かりませんが、
係数行列を求める必要性が理解できません。
有限積分法についての情報も少なく困っています。
皆様のご指導のほどよろしくお願いいたします。
No.2ベストアンサー
- 回答日時:
No.1の方は質問の意図を全く理解されていないようです。
したがって無視してください。
係数行列は、書かれた式をそのまま行列形式で書き下すだけです。
コーディングする上で、係数行列を求める必然性は必ずしもないように思います。
この卒業研究以外に、周波数領域で解く場合がありますが、周波数領域で解く場合には、行列を用いる必然性があります。
それ以外は、定式化上、整然とするという点がメリットに感じます。
この回答への補足
「周波数領域で解く場合」とはどのような場合ですか。
放射された電磁波のスペクトルを見る場合などでしょうか。
浅学で申し訳ございませんが、
「周波数領域で解く場合の行列を用いる必然性」を教えて頂けませんか。
ご回答ありがとうございます。
>定式化上、整然とするという点がメリットに感じます。
おっしゃる通りです。
>係数行列は、書かれた式をそのまま行列形式で書き下すだけです。
各軸の計算格子点を独立に(空間を立方体でなく直方体に)変更したり、
グリッド形状を変更した際に係数行列が生成できたらと思いまして、
コーディングも苦手ですが頑張ってみます。
No.3
- 回答日時:
時間領域で解くのは、逐次的に斬化式を解くことになります。
周波数領域で解くのは、逆行列を求めることになります。
時間領域で解く場合は、周波数スペクトラムを一気に求められるメリットがありますが、低い周波数の解析が出来ない問題があります。
No.1
- 回答日時:
そこの資料をよく読まずに返答しちゃいます。
Ex,Ey,Ezの空間での各値が求めたいということでしょうか?でBx,By,Bzが既知であるということでしょうか?
これは境界値問題を呼ばれる計算手法を取る必要があります。反復計算で近似解を求める方法が多く使用されます。モノによってはFourier変換法で一発で解を導き出すこともできますが…境界条件に依存しますので、ここでは前者をお勧めします。前者は実はいろいろな方法があります。簡単な方法ではSOR法、安定して高速に解きたい場合にはCG法という方法があります。私はCG法を勧めますが理解がもしかしたら大変であれば、SOR法を駆使してゴリゴリ解いていって下さい。
各方法は数値計算の本に何がしら資料があると思いますので、そちらを参考にしてください。もし、何かありましたら再度返答ください。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 テンソル ひずみのマトリクス表記 3 2022/04/23 21:22
- 物理学 物理の問題 2 2022/12/22 22:11
- 物理学 ベクトルと座標系につきまして 1 2022/04/03 06:23
- 数学 球面と接する直線の軌跡が表す領域 4 2023/07/30 12:37
- その他(プログラミング・Web制作) パイソンのプログラミングについての質問です 2 2023/05/22 12:39
- 高校 数3 面積 4 2022/05/11 12:37
- 物理学 電磁波に関する問題 2 2023/01/31 13:52
- 数学 【完全微分方程式⠀】 分数で分母が0になり定義できない場合、分母を仮にtと置いてそれを極限t→0とし 1 2022/05/06 14:43
- 統計学 どの統計を使えばいいのか教えてください(EZ-Rを使用) 5 2022/10/11 13:28
- 数学 この問題が分かりません! 右図の直線①②の式は、y=-x+4①、 y=3/4x+1② である。2つの 3 2022/05/04 22:29
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
物理の有効数字2桁
-
この問題の有効数字の合わせ方...
-
物理
-
音波の振動数の求め方を教えて...
-
配向テンソルについて
-
反エルミート行列の性質について
-
力学についての質問です。 偏角...
-
2端子対回路網の問題について
-
定常速度偏差
-
高校の理科の答えの書き方
-
問3(3)で答えが『-nCpT0』だっ...
-
bcc(面心立方体)の四面体位置で...
-
単振動の運動方程式の解
-
電磁気学の面積分
-
工業力学についてです。 質量80...
-
大学物理 波動
-
レポートの字数に参考文献は含む?
-
大学生です。7000~10000字のレ...
-
レポートを過去レポ(3年以上前)...
-
大学生のレポート4000字って多...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
物理の有効数字2桁
-
bcc(面心立方体)の四面体位置で...
-
熱力学の問題で、(6)について、...
-
高校の理科の答えの書き方
-
1リットルは1辺が1cmの立方体が...
-
単振動の運動方程式の解
-
算数の問題がわかりません。
-
保存力→ポテンシャル→力 の計...
-
FEMのモード法と直接法の違...
-
整数で答えよ、のように書かれ...
-
pHについてです
-
水素原子の波動関数の直交性に...
-
局所ヌセルト数から平均ヌセル...
-
慣性モーメント
-
この問題の(2)が 衝突直後の速...
-
全微分のdx,dyの意味
-
有効数字の桁に注意して、次の...
-
高さHのビルの屋上から初速v₀で...
-
(4)の解き方が分からないので...
-
変数分離が成功したからといっ...
おすすめ情報