-2

100 ポイントまで変化するデータ ポイントのセットを使用してa、b、c 値を決定する必要がありますy=ax2+bx+c。例 (-270,69) (-269,90) (-280,50) の場合。ポイントを使用して二次方程式を生成し、a、b、c値を決定するためにデータURLを補間しましたが、両方の方法でa、b、c値の違いを見つけました。注: Numpy を製品コードに使用することはできません。

def coefficent(x,y):
    x_1 = x[0]
    x_2 = x[1]
    x_3 = x[2]
    y_1 = y[0]
    y_2 = y[1]
    y_3 = y[2]

    a = y_1/((x_1-x_2)*(x_1-x_3)) + y_2/((x_2-x_1)*(x_2-x_3)) + y_3/((x_3-x_1)*(x_3-x_2))

    b = -y_1*(x_2+x_3)/((x_1-x_2)*(x_1-x_3))
    -y_2*(x_1+x_3)/((x_2-x_1)*(x_2-x_3))
    -y_3*(x_1+x_2)/((x_3-x_1)*(x_3-x_2))

    c = y_1*x_2*x_3/((x_1-x_2)*(x_1-x_3))
    + y_2*x_1*x_3/((x_2-x_1)*(x_2-x_3))
    + y_3*x_1*x_2/((x_3-x_1)*(x_3-x_2))

    return a,b,c

x = [1,2,3]
y = [4,7,12]

a,b,c = coefficent(x, y)
print a,b,c


> import numpy as np
>>> A, B, C = np.polyfit([1,2,3],[4,7,12],2)
>>> print A, B, C
1.0 -4.2727620148e-15 3.0
>>> print A, 'x^2 +', B, 'x +', C
1.0 x^2 + -4.2727620148e-15 x + 3.0
>>>
4

1 に答える 1

4

SOに投稿する前に計算bして行を分割しましたか？c質問に貼り付けられたコードはコンパイルされません。このバージョンは次のことを行います。

def coefficient(x,y):
    x_1 = x[0]
    x_2 = x[1]
    x_3 = x[2]
    y_1 = y[0]
    y_2 = y[1]
    y_3 = y[2]

    a = y_1/((x_1-x_2)*(x_1-x_3)) + y_2/((x_2-x_1)*(x_2-x_3)) + y_3/((x_3-x_1)*(x_3-x_2))

    b = (-y_1*(x_2+x_3)/((x_1-x_2)*(x_1-x_3))
         -y_2*(x_1+x_3)/((x_2-x_1)*(x_2-x_3))
         -y_3*(x_1+x_2)/((x_3-x_1)*(x_3-x_2)))

    c = (y_1*x_2*x_3/((x_1-x_2)*(x_1-x_3))
        +y_2*x_1*x_3/((x_2-x_1)*(x_2-x_3))
        +y_3*x_1*x_2/((x_3-x_1)*(x_3-x_2)))

    return a,b,c

x = [1,2,3]
y = [4,7,12]

a,b,c = coefficient(x, y)

print "a = ", a
print "b = ", b
print "c = ", c

出力は完璧です：

a =  1
b =  0
c =  3

これは、 からの回答よりも( の係数で4*10 -15ほど)正確です。3 つのデータ ポイントについても、数学的に正確です。bnumpy

あなたのコードがあなたに与える答えの問題は何ですか?

于 2013-10-04T07:24:59.760 に答える