メイン asyr.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="inv.js">document.write("inv.js ファイルが見つかりません??<br>")</script>
<script src="io.js">document.write("io.js ファイルが見つかりません??<br>")</script>

<script language="JavaScript">
<!--
var alpha, beta, gamma
var epsinv = 1e-8

function mreg(data, nc)
{
  var mean, variance, sd, ss, r, nv, i, j, k, b0, b, nv1, temp2
  nv = 3
  mean = new Array(nv)
  variance = new Array(nv)
  sd = new Array(nv)
  ss = new Array(nv)
  r = makeMatrix(nv, nv)
  b = new Array(nv)

// Netscape communicator のバグ?配列ベクトルは初期化されていない
  for (i = 0; i < nv; i++) {
    mean[i] = 0
    for (j = 0; j < nv; j++) {
      r[i][j] = 0
    }
  }
  for (i = 0; i < nc; i++) {
    for (j = 0; j < nv; j++) {
      mean[j] += data[i][j]
    }
  }
  for (j = 0; j < nv; j++) {
    mean[j] /= nc
  }

  for (i = 0; i < nc; i++) {
    for (j = 0; j < nv; j++) {
      for (k = 0; k < nv; k++) {
        r[j][k] += (mean[j]-data[i][j])*(mean[k]-data[i][k])
      }
    }
  }
  for (j = 0; j < nv; j++) {
    ss[j] = r[nv-1][j]
    variance[j] = r[j][j]/(nc-1)
    sd[j] = Math.sqrt(variance[j])
  }
  for (j = 0; j < nv; j++) {
    for (k = 0; k < nv; k++) {
      r[j][k] /= (nc-1)*sd[j]*sd[k]
    }
  }

  nv1 = nv-1
  if (inv(r, nv1, epsinv) == 1) {
    printf("逆行列を求められません\n")
    return 1
  }
  b0 = mean[nv1]
  for (i = 0; i < nv1; i++) {
    temp2 = 0
    for (j = 0; j < nv1; j++) {
      temp2 += r[i][j]*r[j][nv1]
    }
    b[i] = temp2*sd[nv1]/sd[i]
    b0 -= b[i]*mean[i]
  }
  alpha = b[0]
  beta = b[1]
  gamma = b0
  return 0
}

function init(data, nc)
{
  var k
  if (nc == 4) {
    return (4*data[3]+data[2]-5*data[1])/(4*data[2]+data[1]-5*data[0])
  }
  else if (nc == 5) {
    return (4*data[4]+3*data[3]-data[2]-6*data[1])/(4*data[3]+3*data[2]-data[1]-6*data[0])
  }
  else if (nc == 6) {
    return (4*data[5]+4*data[4]+2*data[3]-3*data[2]-7*data[1])/(4*data[4]+4*data[3]+2*data[2]-3*data[1]-7*data[0])
  }
  else if (nc == 7) {
    return (data[6]+data[5]+data[4]-data[2]-2*data[1])/(data[5]+data[4]+data[3]-data[1]-2*data[0])
  }
  else {
    k = nc/7
    return (data[Math.floor(6*k)]+data[Math.floor(5*k)]+data[Math.floor(4*k)]-data[Math.floor(2*k)]-2*data[Math.floor(k)])/(data[Math.floor(5*k)]+data[Math.floor(4*k)]+data[Math.floor(3*k)]-data[Math.floor(k)]-2*data[0])
  }
}

function calc(data_string)
{
  var data, nc, nv, i, j, d, u0, u1, loop
  var datax = new Array()
  var datay = new Array()
  if ((data = getdata(data_string, 2)) == false) return
  nc = data.length
  if (nc <= 3) {
    printf("有効データ組数が %i 個しかないので,分析できません\n", nc)
    return
  }
  for (i = 0; i < nc; i++) {
    datax[i] = data[i][0]
    datay[i] = data[i][1]
  }

  d = new Array(nc)
  for (i = 0; i < nc; i++) {
    d[i] = new Array(3)
  }

  printf("★ 漸近指数曲線( Y = a・b^X+c)へのあてはめ ★\n\n")
  u1 = u0 = init(datay, nc)
  printf("初期値 b = %g\n", u0)
  for (loop = 1; ; loop++) {
    if (loop % 20 == 0 && confirm("反復回数が "+loop+" 回になりましたが,まだ収束しません。続けますか?") == false) break
    for (i = 0; i < nc; i++) {
      d[i][0] = Math.pow(u1, datax[i])
      d[i][1] = datax[i]*Math.pow(u1, datax[i]-1)
      d[i][2] = datay[i]
    }
    if (mreg(d, nc) == 1) {
      return
    }
    u0 = u1
    u1 += beta/alpha
    if (Math.abs((u0-u1)/u0) < 0.000001) break
  }
  beta = u1
  printf("予測値 = %g * %g ^ x + %g\n\n", alpha, beta, gamma)
  printf("%12s%12s%12s\n", "X", "Y", "予測値")
  for (i = 0; i < nc; i++) {
    printf("%12.6g%12.6g%12.6g\n", datax[i], datay[i], alpha*Math.pow(beta,  datax[i])+gamma)
  }
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>漸近指数曲線へのあてはめ</b></font> <a  href="src/asyr.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>
独立変数と従属変数のデータを入力して,漸近指数曲線(Y = a・b^X+c)にあてはめを行う。<br>
本プログラムは線形近似によるものなので,<a href="fit-asymptotic.html">非線形回帰</a>を行う方が望ましい。<br>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Regression/expreg/type4.html">手法の解説</a> <img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/asyr.html">使用法</a>
<hr>
<table><tr>
<td><input type="button" name="calcurate" value="計算開始" onClick="calc(this.form.data.value)">  </td>
<td><input type="button" name="clear" value="入力欄クリア" onClick="this.form.data.value=''">  </td>
<td><input type="button" name="clear" value="出力欄クリア" onClick="this.form.result.value=''"></td>
</tr></table>

<p>
<table><tr>
<td>入力欄<br><textarea name="data" ROWS=15 COLS=15></textarea></td>
<td>出力欄<br><textarea name="result" ROWS=15 COLS=80></textarea></td>
</tr></table>
</form>

<p><hr noshade>
<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><IMG SRC="../gra/ume5.png" width=121 height=37 ALT="Made with Macintosh">
</body>
</html>

サブ inv.js   Last modified: Mar 25, 2004
サブ io.js   Last modified: Mar 25, 2004

・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI

Made with Macintosh