php中的矩阵乘法

时间:2014-08-09 16:28:13

标签: php matrix

虽然矩阵的顺序应该没问题,但以下代码会抛出异常。这可能是我无法注意到的一件小事,但无法弄清楚。

<?php
  $mat1 = array(5,1);
  $mat2 = array(1,5);
  function matrixmult($m1,$m2){
    $r=count($m1);
    $c=count($m2[0]);
    $p=count($m2);
    if(count($m1[0])!=$p){throw new Exception('Incompatible matrixes');}
      $m3=array();
      for ($i=0;$i< $r;$i++){
        for($j=0;$j<$c;$j++){
          $m3[$i][$j]=0;
          for($k=0;$k<$p;$k++){
            $m3[$i][$j]+=$m1[$i][$k]*$m2[$k][$j];
          }
        }
      }
    }
    return($m3);
  }
  matrixmult($mat1,$mat2);
?>

3 个答案:

答案 0 :(得分:1)

您以两种方式指定测试矩阵错误:

  1. 数组不是二维的(即数组数组)。
  2. 即使你在它们周围缠绕了另一个array( ),第一个矩阵的宽度等于第二个矩阵的高度的条件也不会与[5 1]和[1 5]成立,是2宽1高。
  3. 您需要的是

    $mat1 = array(array(5,1));
    $mat2 = array(array(1),array(5));
    

答案 1 :(得分:0)

只是为了解决这个问题。这是一个有效的解决方案:

function M_mult($_A,$_B) {
  // AxB outcome is C with A's rows and B'c cols
  $r = count($_A);
  $c = count($_B[0]);
  $in= count($_B); // or $_A[0]. $in is 'inner' count

  if ( $in != count($_A[0]) ) {
    print("ERROR: need to have inner size of matrices match.\n");
    print("     : trying to multiply a ".count($_A)."x".count($_A[0])." by a ".count($_B)."x".count($_B[0])." matrix.\n");
    print("\n");
    exit(1);
  }

  // allocate retval
  $retval = array();
  for($i=0;$i< $r; $i++) { $retval[$i] = array(); }
    // multiplication here
    for($ri=0;$ri<$r;$ri++) {
      for($ci=0;$ci<$c;$ci++) {
        $retval[$ri][$ci] = 0.0;
        for($j=0;$j<$in;$j++) {
          $retval[$ri][$ci] += $_A[$ri][$j] * $_B[$j][$ci];
        }
      }
    }
    return $retval;
  }
}

答案 2 :(得分:0)

返回前必须删除花括号。

<?php
  $mat1 = array(5,1);
  $mat2 = array(1,5);
  function matrixmult($m1,$m2){
    $r=count($m1);
    $c=count($m2[0]);
    $p=count($m2);
    if(count($m1[0])!=$p){throw new Exception('Incompatible matrixes');}
    $m3=array();
    for ($i=0;$i< $r;$i++){
      for($j=0;$j<$c;$j++){
        $m3[$i][$j]=0;
        for($k=0;$k<$p;$k++){
          $m3[$i][$j]+=$m1[$i][$k]*$m2[$k][$j];
        }
      }
    }
    return($m3);
  }
  matrixmult($mat1,$mat2);
?>