phpで移動平均(SMA,EMA)、ボリンジャーバンドを計算する

$aData[] = 12;
$aData[] = 1;
$aData[] = 100;
$aData[] = 13;

// 平均値
$Average = array_sum($aData)/count($aArrayData);

trait traitTrader // traitという色々なクラスで useで取り込める部品として作る
{
  // 単純"移動"平均 Simple Moving Average
  public function SMA($position, $period, $aPrice)
  {
    $result = 0;
    if($position>=$period-1 && $period>0) {
      for($i=0; $i<$period;$i++) {
        $result+=$aPrice[$position-$i];
      }
      $result/=$period;
    }
    return $result;
  }

  // 上記関数では、指定した $position のところの移動平均値が計算される
  // これを全期間分やりたい場合
  // 下記のように  $period:移動平均を計算する期間 $aPrice:価格が入った配列
  // を渡して、$aSMA 配列を作成する
  public function getSmaArray($period, $aPrice)
  {
    $aSMA = array();
    foreach($aPrice as $position=>$price) {
      $aSMA[$position] = $this->SMA($position, $period, $aPrice);
    }
    return $aSMA;
  }
}

 // これも　trait traitTraderの中に作る


  // double ExponentialMA(const int position,const int period,const double prev_
  // value,const double &price[])
  // ↑ MQ4の中のサンプルソースの関数定義
  // ↓ それを基に作ったphp関数
  public function EMA($position, $period, $prev_value, $aPrice)
  {
    $result = 0;
    if($period>0) {
      $pr = 2 / ( $period + 1 );
      $result = $aPrice[ $position ]*$pr + $prev_value*(1-$pr);
    }
    return $result;
  }


  // 上記関数では、指定した $position のところの移動平均値が計算される
  // これを全期間分やりたい場合
  // 下記のように  $period:移動平均を計算する期間 $aPrice:価格が入った配列
  // を渡して、$aEMA 配列を作成する
  public function getEmaArray($period, $aPrice)
  {
    $aEma = array();
    $prev_value = $this->SMA(0, $period, $aPrice);
    foreach($aPrice as $position=>$price) {
      if($position<$period) {
        $ema = $this->SMA($position, $period, $aPrice);
      } else {
        $ema = $this->EMA($position, $period, $prev_value, $aPrice);
      }
      $aEma[$position] = $ema;
      $prev_value = $ema;
    }
    return $aEma;
  }

 // これも trait traitTraderの中に作る

  // 標準偏差を求める関数
  public function StdDev_Func($position, $aPrice, $aMAprice, $period)
  {
    $StdDev_dTmp = 0;
    if($position>=$period) {
      for($i=0; $i<$period; $i++) {
         $StdDev_dTmp += pow($aPrice[$position-$i]-$aMAprice[$position],2); 
         // 分散を計算 (平均と価格の差の２乗）
      }
      $StdDev_dTmp = sqrt($StdDev_dTmp/$period); 
      // 分散のスクエアルート（平方根)が標準偏差
    }
    return $StdDev_dTmp;
  }

// --------------------------
    // ボリンジャーバンドを計算するには
    // traitTraderをuseしているクラスの中で下記のように求める

    // 単純移動平均をまず求めて
    $aSMA = $this->getSmaArray($period, $aPrice);

    // $aPrice[0] = 100万円;
    // $aPrice[1] =  90万円; // のように価格が入っている配列
    foreach($aPrice as $i=>$price) {

      // 単純移動平均の標準偏差を計算する
      $StdDev = $this->StdDev_Func($i, $aPrice, $aSMA, $period);

      $DT['+3σ'][$i] = $aSMA[$i]+3*$StdDev; // 標準偏差σ の３倍を平均に加える
　　　　　　　　　　　　　　　　　　　　　　 // 99.73％区間の上限
      $DT['-3σ'][$i] = $aSMA[$i]-3*$StdDev; // 標準偏差σ の３倍を平均から引く
　　　　　　　　　　　　　　　　　　　　　　 // 99.73％区間の下限
      $DT['+2σ'][$i] = $aSMA[$i]+2*$StdDev; // 標準偏差σ の2倍を平均に加える
　　　　　　　　　　　　　　　　　　　　　　 // 95.45％区間の上限
      $DT['-2σ'][$i] = $aSMA[$i]-2*$StdDev; // 標準偏差σ の2倍を平均から引く
　　　　　　　　　　　　　　　　　　　　　　 // 95.45％区間の下限
    }