現在、C言語でプログラムを組んでいるのですが、複素数からなる行列の固有値を求めるプログラムがどうしてもつくることができません。
インターネットやプログラムのパッケージなどを調べてみたのですが、実数からなる行列の固有値を求めるプログラムしか載っておらず、勉強不足によりそれを複素数の場合に拡張することもできませんでした。
もしよければ、プログラムの組み方や実数のプログラムから複素数のプログラムへの替え方のこつ、もしくは「この本に載っていたよ」など、どんな情報でもかまいませんので教えて頂けないでしょうか?
すいませんがよろしくお願いします。
No.3
- 回答日時:
一部訂正。
typedef double[2] complex;
↓
typedef double complex[2];
ですね。(最近Javaばかりやってたもんでつい...。)
コンパイラによっては、
typedef double *complex;
でないと、関数定義が通らないかもしれません。
No.2
- 回答日時:
実数行列に関するアルゴリズムは分かっているわけですから、あとはそれを複素数に
拡張すればよいだけの話ですよね。(C++等の方がやりやすそうですが)
複素数をどのような形式で表現するかにもよりますが、例えば
typedef double[2] complex;
のようにしていたら、
void plus(complex c1, complex c2, complex result) {
result[0] = c1[0] + c2[0];
result[1] = c1[1] + c2[1];
}
void minus(complex c1, complex c2, complex result) {
result[0] = c1[0] - c2[0];
result[1] = c1[1] - c2[1];
}
void multiple(complex c1, complex c2, complex result) {
result[0] = c1[0] * c2[0] - c1[1] * c2[1];
result[1] = c1[0] * c2[1] + c1[1] * c2[0];
}
のように加算・減算・乗算を定義して(今回は除算はとりあえず不要でしょうから)
実数向けのアルゴリズムを置き換えればよいと思います。
No.1ベストアンサー
- 回答日時:
takashi5050 さんのお望みのアルゴリズムが載っているかどうか知りませんが、
技術評論社から出版されている Numerical Recipes in C に、
数値計算のプログラムが大量に載っています。
11章が固有値問題の数値計算法というタイトルなので、載ってるんじゃないかと思います
ちょっと高いですけど、それだけの価値はあると思います。
大き目の本屋さんなら置いていると思うので、覗いてみてください。
Numerical Recipes in C 日本語版
William H. Press・Saul A. Teukolsky・William T. Vetterling・Brian P. Flannery 共著
丹慶勝市・奥村晴彦・佐藤俊郎・小林誠 訳
B5変形 判 / 680 ページ / 本体価格 4757 円+税
ISBN4-87408-560-1
参考URL:http://bookweb.kinokuniya.co.jp/guest/cgi-bin/ws …
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
関連するカテゴリからQ&Aを探す
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
vba listviewにおけるtextのAli...
-
指定したフォルダーを最前面表...
-
c言語のサイコロを100回振って...
-
GitHubで拾ったプログラムが動...
-
世界一美しいソースコード
-
C言語primeについて
-
デジタルフィルタ
-
プロダクションコードとは?
-
初級者の質問――time_tについて
-
掲示板での出来事について
-
C言語でのソケット通信のclose
-
C言語で前面表示や、特定のファ...
-
フィボナッチヒープです
-
Access2013からADP廃止
-
SH2を使いたいです
-
パソコン CPU オーバークロッ...
-
Pythonで仮想通貨の自動売買が...
-
C言語について。
-
バックグラウンドで動作させる。
-
プログラミング、計算機科学、d...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プロダクションコードとは?
-
マイコンからプログラムを読み...
-
vba listviewにおけるtextのAli...
-
自分のHTMLサイトでRSSを読み込...
-
指定したフォルダーを最前面表...
-
Arduino nano(アルドゥイーノ ...
-
プログラミングの問題と解答例...
-
COBOLのS9タイプからXタイプへ...
-
2値化した画像の黒い部分のピク...
-
無限ループ中にある任意のキー...
-
初級者の質問――time_tについて
-
Excel VBA じゃんけん
-
どうやってパソコンなどの凄い...
-
三次元のグラフを書きたいので...
-
VBコンソールアプリ(.NET Frame...
-
NumLockの制御について
-
世界一美しいソースコード
-
expectを使用しない対話型処理
-
Excel VBA で、データベースフ...
-
C言語→MATLABへの変換で、困っ...
おすすめ情報