0

データセットがあり、そこから平面の方程式を取得しようとしています。つまり、a * x + b * y + c = z私の場合、dT = x、dTa = y、Constant = c、およびdV=zです。

私はこれをMatlabで非常に簡単に行うことができます、コード：

dT = [8.5; 3.5; .4; 12.9]

dT =

    8.5000
    3.5000
    0.4000
   12.9000

dTa = [8.5; 18; 22; 34.9]

dTa =

    8.5000
   18.0000
   22.0000
   34.9000

dV = [3; 1; .5; 3]

dV =

    3.0000
    1.0000
    0.5000
    3.0000

Constant = ones(size(dT))

Constant =

     1
     1
     1
     1

coefficients = [dT dTa Constant]\dV

coefficients =

    0.2535
   -0.0392
    1.0895

したがって、ここでは、係数=（a、b、c）です。

Pythonでこれを行う同等の方法はありますか？numpyモジュール（numpy.linalg）を使おうとしていますが、うまく機能していません。一つには、行列は正方形でなければなりません、そしてそれでも、それは非常に良い答えを与えません。例えば：

エラー：

>>> dT
[8.5, 3.5, 0.4, 12.9]
>>> dTa
[8.5, 18, 22, 34.9]
>>> dV
[3, 1, 0.5, 3]
>>> Constant
array([ 1.,  1.,  1.,  1.])
>>> numpy.linalg.solve([dT, dTa, Constant], dV)

Traceback (most recent call last):
  File "<pyshell#45>", line 1, in <module>
    numpy.linalg.solve([dT, dTa, Constant], dV)
  File "C:\Python27\lib\site-packages\numpy\linalg\linalg.py", line 312, in solve
    _assertSquareness(a)
  File "C:\Python27\lib\site-packages\numpy\linalg\linalg.py", line 160, in _assertSquareness
    raise LinAlgError, 'Array must be square'
LinAlgError: Array must be square

正方行列を使ったワーク：

>>> dT
array([  8.5,   3.5,  12.9])
>>> dTa
array([  8.5,  18. ,  34.9])
>>> dV
array([3, 1, 3])
>>> Constant
array([ 1.,  1.,  1.])
>>> numpy.linalg.solve([dT, dTa, Constant], dV)
array([ 2.1372267 ,  2.79746835, -1.93469505])

これらは私が以前に得た値にさえ近くありません!!

アイデアはありますか？アドバイスをいただければ幸いです。

4

2 に答える 2

1

最初の例では、 を使用する必要がありますnumpy.linalg.lstsq。いずれにせよ、行列の行と列を混同しているようです。のようなものを使用して修正しnumpy.linalg.solve(zip(dT, dTa, Constant), dV)ます。

于 2012-08-14T23:33:59.830 に答える
0

ということで、こんな写真を発見

定数を含むように調整しました。ここに私の解決したコードがあります:

    X = 0
    XX = 0
    XY = 0
    XZ = 0

    Y = 0
    YY = 0
    YZ = 0

    Z = 0

    for j in range(0, len(dTemp)):
        X = X + dTemp[j]
        XX = XX + (dTemp[j] * dTemp[j])
        XY = XY + (dTemp[j] * dTempA[j])
        XZ = XZ + (dTemp[j] * Applied[j])

        Y = Y + dTempA[j]
        YY = YY + (dTempA[j] * dTempA[j])
        YZ = YZ + (dTempA[j] * Applied[j])

        Z = Z + Applied[j]


    lhs = numpy.array([[XX, XY, X], [XY, YY, Y], [X, Y, 1]]) 
    rhs = numpy.array([XZ,  YZ, Z])

    coefficients = numpy.linalg (lhs, rhs)

    a = coefficients[0]
    b = coefficients[1]
    c = coefficients[2]

私はそれがうまくいったと思います！それでも、値は少しずれていますが、Matlab が別のアルゴリズムを使用していたためでしょうか?

于 2012-08-15T18:47:56.117 に答える