メイン　power.html　　　Last modified: Sep 01, 2009
<html>
<head>
<meta http-equiv="Content-Type" content="text/html;CHARSET=EUC-JP">
<link rel="shortcut icon" href="../favicon.ico">
<title>JavaScript</title>
<script src="io.js">document.write("io.js ファイルが見つかりません??<br>")</script>

<script language="JavaScript">
<!--

function check(e_value, evec, a, nv)
{
  var i, j, c
  for (i = 0; i < nv; i++) {
    c = 0
    for (j = 0; j < nv; j++) {
      c += a[i][j]*evec[j]
    }
    if (Math.abs((c-e_value*evec[i])/c) > 1e-5) {
      return 1
    }
  }
  return 0
}

function power(r, nv)
{
  var MAXLOOP = 500
  var loop, l2, i, j, ok, temp
  var e_value, eva0, evec, eve0, orig
  orig = new Array(nv)
  for (i = 0; i < nv; i++) {
    orig[i] = new Array(nv)
    for (j = 0; j < nv; j++) {
      orig[i][j] = r[i][j]
    }
  }
  evec = new Array(nv)
  eve0 = new Array(nv)
  for (loop = 0; loop < nv; loop++) {
    for (i = 0; i < nv; i++) {
      evec[i] = 0
    }
    evec[loop] = 1
    eva0 = 1
    for (l2 = 0; l2 < MAXLOOP; l2++) {
      for (i = 0; i < nv; i++) {
        temp = 0
        for (j = 0; j < nv; j++) {
          temp += r[i][j]*evec[j]
        }
        eve0[i] = temp
      }
      e_value = 0
      for (i = 0; i < nv; i++) {
        e_value += eve0[i]*eve0[i]
      }
      e_value = Math.sqrt(e_value)
      ok = 1
      for (i = 0; i < nv; i++) {
        temp = evec[i]
        evec[i] = eve0[i]/e_value
        if (Math.abs((evec[i]-temp)/temp) > 1e-7) {
          ok = 0
        }
      }
      if (ok == 1) {
        break
      }
      eva0 = e_value
    }
    if (check(e_value, evec, orig, nv) == 1) {
      printf("これ以上の解は求まりません\n")
      return
    }
    printf("第 %i 固有値 = %.6g    反復回数 %i\n", loop+1, e_value, l2)
    printf("対応する固有値ベクトル = (")
    for (i = 0; i < nv; i++) {
      printf(" %.6g", evec[i])
    }
    printf(" )\n\n")
    for (i = 0; i < nv; i++) {
      for (j = 0; j < nv; j++) {
        r[i][j] -= e_value*evec[i]*evec[j]
      }
    }
  }
  printf("★ 残差行列 ★\n\n")
  for (i = 0; i < nv; i++) {
    for (j = 0; j < nv; j++) {
      printf("%12.6g", r[i][j])
    }
    printf("\n")
  }
}

function calc(data_string)
{
  var i, j, nc, rank, r
  if ((r = getdata(data_string, 0)) != false) {
    nc = r.length
    rank = r[0].length
    if (nc != rank) {
      printf("正方行列でなくてはなりません\n")
    }
    else {
      printf("★ 入力された行列 ★\n\n")
      for (i = 0; i < rank; i++) {
        for (j = 0; j < rank; j++) {
          printf("%12.6g", r[i][j])
        }
        printf("\n")
      }

      printf("\n★ 固有値・固有ベクトル ★\n\n")
      power(r, rank)
    }
  }
  sep2(50)
  printf("\n")
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>パワー法による固有値・固有ベクトル</b></font>　<a  href="src/power.html"><img src="png/src.png" width=35 height=11 alt="src" align=top></a>　　　　　Last modified: Jun 01, 2006<hr noshade><p>
<font color="#ff0000" size="+2">以下のプログラムのサポートは終了しました。自己責任でお使い下さい。</font>

<form name=Result>
<script language="JavaScript">
<!--
//-->JavaScript がサポートされていないブラウザですか？
</script>
<input type="button" name="calcurate" value="計算開始" onClick="calc(this.form.data.value)">　　
<input type="button" name="clear" value="入力欄クリア" onClick="this.form.data.value=''">　　
<input type="button" name="clear" value="出力欄クリア" onClick="this.form.result.value=''"><p>
入力欄（<a href="exa/kugirimoji.html">区切り文字</a>で区切って，行単位に入力）<br><textarea name="data" ROWS=10 COLS=80></textarea><p>
出力欄<br><textarea name="result" ROWS=30 COLS=80></textarea>
</form>

<p><hr noshade>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/power-method/power.html">手法の解説</a><br>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <A HREF="javascript:history.go(-1)">直前のページへ戻る</A>　　<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../mail.html">E-mail to Shigenobu AOKI</a>
<p><center><IMG SRC="../gra/ume5.png" width=121 height=37 ALT="Made with Macintosh"></center>
</body>
</html>

サブ　io.js　　　Last modified: Mar 25, 2004
・ 直前のページへ戻る　　・ E-mail to Shigenobu AOKI

Made with Macintosh