2

我必须在这里遗漏一些非常简单的东西。我有一些创建简单线性系统的 JS 代码（我正在尝试创建两条斜线之间的最短线）。我已经到了我有 Ax = b 的地步，需要求解 x。A 是 3 x 2 矩阵，b 是 3 x 1。

我有：

function build_equation_system(v1, v2, b) { 
    var a = [ [v1.x, v2.x], [v1.y, v2.y], [v1.z, v2.z] ];
    var b = [ [b.x], [b.y], [b.z]];
    return numeric.solve(a,b)
}

Numeric 返回一个 1 x 3 的 NaN 矩阵，即使存在解也是如此。

4

1 回答 1

2

使用数字，您可以执行以下操作：

创建一个计算 A 矩阵的伪逆的函数：

function pinv(A) {
    return numeric.dot(numeric.inv(numeric.dot(numeric.transpose(A),A)),numeric.transpose(A));
}

使用该函数求解线性最小二乘方程以获得系数。

var p = numeric.dot(pinv(a),b);

我尝试了您最初的使用方法，numeric.solve但也无法正常工作，所以我很想知道问题出在哪里。

一个简单的测试...

var x = new Array(10);
var y = new Array(10);
for (var i = 0; i < 10; ++i) { 
    x[i] = i;
    y[i] = i;
}
// Solve for the first order equation representing this data
var n = 1;
// Construct Vandermonde matrix.
var A = numeric.rep([x.length, n + 1], 1);
for (var i = 0; i < x.length; ++i) {
    for (var j = n-1; j >= 0; --j) {
        A[i][j] = x[i] * A[i][j+1];
    }
}

// Solves the system Ap = y
var p = numeric.dot(pinv(A),y);

p = [1, 2.55351295663786e-15]

我已经使用这种方法重新创建了 MATLAB 的 polyfit 以供 Javascript 使用。

于 2013-08-28T13:31:15.990 回答