忍者ブログ
20080511~ 13と7と11の倍数の論理積は13と7と11の積の倍数である。 和ァ・・・
[4314] [4313] [4312] [4311] [4310] [4309] [4308] [4307] [4306] [4305] [4303]
y''+ky=0

こいつの初期値の、初期位置y(0)と初速度y'(0)に適当な数を入れると、
実数から複素数への体の拡大みたいのが起こせます。


たとえば差分方程式にして

{y(n+1)-2y(n)-y(n-1)}/dt^2=-k*y(n)

ってして、
y0=i、y1=1ってしたら、複素平面にはみ出した振り子
つまりほぼ円になったオイラーの公式が出てきます。

初期値次第では真円の螺旋だけじゃなく楕円の螺旋にもできるみたいですね。

もちろん実数だけでやれば直線偏光みたいな振り子的解が出てきますが、
y0=±i、y1=±i(複合同順でない)とかだと、いわゆる振り子の円弧を1次元とみなすみたいな線形従属っぽい感じで
実数あるいは純虚数の振動解だけが有限になります。


今はあまり余裕がないので詳しくは書きませんが
y(n+1)=IMSUB(IMPRODUCT(IMSUB(2,k*dt^2),y(n)),y(n-1))

の式で、Excelでも出せますよ。

拍手[0回]

PR

コメント


コメントフォーム
お名前
タイトル
文字色
メールアドレス
URL
コメント
パスワード
  Vodafone絵文字 i-mode絵文字 Ezweb絵文字


忍者ブログ [PR]
カレンダー
02 2024/03 04
S M T W T F S
1 2
3 4 5 6 7 8 9
10 11 12 13 14 15 16
17 18 19 20 21 22 23
24 25 26 27 28 29 30
31
ブログランキング
ブログランキング参戦中
にほんブログ村 アニメブログ 深夜アニメへ
にほんブログ村 漫画ブログ SF・ファンタジー漫画へ
にほんブログ村 科学ブログ 自然科学へ
よかったらポチッとお願いします^^
最新CM
[12/30 buy steroids credit card]
[09/26 Rositawok]
[03/24 hydraTep]
[03/18 Thomaniveigo]
[03/17 Robertaverm]
最新TB
プロフィール
HN:
量子きのこ
年齢:
42
性別:
男性
誕生日:
1981/04/04
職業:
WinDOS.N臣T
趣味:
妄想・計算・測定・アニメ
自己紹介:
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます
例:H→7番目、XP→15(P)×26+23(X)=413番目。
A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいズコー
バーコード
ブログ内検索
アクセス解析