三次方程式の解

今日はアルゴリズムのお話を。 FEMのポスト処理ソフトを作っていると三次方程式の解き方が必要な時が多々でてきます。三次方程式の出現元は2階のテンソルの主値を求める処理。例えば - 応力テンソルから最大主応力を求めたり - 剛体の慣性テンソル（らしき…

2008-06-15

ヤコビ法

Program Technique アルゴリズム

ラプラシアン行列を勉強するために、とりあえずヤコビ法のプログラムを作ってみました。標準入力から最初に次元の数をもらって、その後各行を受け取るプログラムです。反復回数は適当に１０回とかにしてみました（本当は非対角要素の大きさを見た方がいいか…

2008-03-27

曲率半径

アルゴリズム

梁理論を人前で説明する必要があり…曲率半径が必要になりました。で、曲率半径がと確認するためにかかった時間が…１時間…あれぇ、大学時代なら３０分ぐらいで思い出せるはずなのに…。とにかく、を納得するまで…ながかった気がします…（反省）大変参考になっ…

2008-02-14

常微分方程式を復習

アルゴリズム

フーリエ変換関連の本を読んでいたら、ピアノとチェンバロの話がでてきた。その本には、ピアノは叩く事で弦を振動させる、つまり厳密なことをいわなければ初期速度で振動する。チェンバロは引っ張る事で、つまり初期変位で振動する。同じ微分方程式をとい…

2008-02-01

離散フーリエ変換

アルゴリズム

周波数、周期、波長などなど、大学時代（本当は高校時代か！）やっていた内容を一通り復習しました。ずいぶん思い出してきました。波長が短くなると、すなわち一秒間あたりに一回りする量は増えるので周波数が大きくなるとか。で、sinとcosの級数和で関数を…

2008-01-30

離散フーリエ変換の周波数

アルゴリズム

あいかわらず復習中です。離散フーリエ変換。今日復習した内容は周波数です。積分する範囲が0からの時、周波数は整数でどんどん分割されていきます。具体的には、、、…この周波数に関して離散化した場合のイメージは、入力されたデータの個数がN個で幅がdx…

2008-01-29

離散化フーリエ変換その３

アルゴリズム

昨日書いたプログラムもなんか自信ないのでもうすこしなやもうと思います。周波数の求め方とか特に。それと係数のほうも最後にデータの個数で割るのか長さでわるのかいまいちちゃんと説明できない…。

2008-01-28

離散フーリエ変換

アルゴリズム

前回のフーリエ変換のプログラムだと入力されるデータの間隔が1.0でないとうまく動かないという問題がありました。それと、周波数を求める箇所が間違えてました。ってことで改良したものを #include #include using namespace std; class Dft { public: int …

2008-01-26

昨日の離散フーリエ変換のプログラム間違ってました

アルゴリズム

よく考えると、sinとcosの係数を求める箇所、間違ってました。データの総数で割るところ…しっぱいしました…。

2008-01-25

離散フーリエ変換

アルゴリズム

なんとなく離散フーリエ変換をやってみたくてこんなプログラムを作ってみました。 #include #include using namespace std; class Dft { public: int m_sum; double* m_cos_coefficient; double* m_sin_coefficient; public: Dft(); ~Dft(); void computing(…