0

「解像度」が小さすぎるため、最近、私が持っているデータセットで問題に遭遇しました。データセットは、実際には最大値ではない最大値 (0) を持つように特徴付けられています。これは、すべての数値がその「解像度」内で見つかったより高い数値に正規化されているためです。これは、基本的に「X」の「1 x 1 単位」です。桁。Y 値は常に放物線に似ています。

   Y       X
-34,32   -93
-16,56   -92
-10,04   -91
-6,03    -90
-3,34    -89
-1,56    -88
-0,47    -87
 0       -86
-0,10    -85
-0,78    -84
-2,11    -83
-4,20    -82
-7,36    -81
-12,31   -80
-22,03   -79
-25,17   -78

これで、0 は 0 から -0.10 の間、または 0 から -0.47 の間のどこかにあるはずです。線形補間では誤差が大きすぎるため、3 次スプラインを使用すると思います。

私がする必要があるのは、Y REALLY = 0 である X パラメータを計算することです。

残念ながら、私が試したすべての関数は指定された X から Y を計算するため、インターネットで見つけたものを使用できませんでした (コードの例の下)。誰でも助けることができますか？

Function SpLine(PeriodCol As Range, RateCol As Range, x As Range)

Dim Period_Count As Integer
Dim Rate_Count As Integer
Dim c As Integer
Dim n As Integer
Dim i, k As Integer
Dim p, qn, sig, un As Single
Dim klo, khi As Integer
Dim h, b, a As Single

' shows the calculation of a cubic spline interpolated value given known values for fixed periods

Period_Count = PeriodCol.Rows.Count
Rate_Count = RateCol.Rows.Count

If Period_Count <> Rate_Count Then
    SpLine = "Error: Range count dos not match"
    GoTo endnow
End If

ReDim xin(Period_Count) As Single
ReDim yin(Period_Count) As Single

For c = 1 To Period_Count
    xin(c) = PeriodCol(c)
    yin(c) = RateCol(c)
Next c

ReDim u(Period_Count - 1) As Single
ReDim yt(Period_Count) As Single
n = Period_Count
yt(1) = 0
u(1) = 0

For i = 2 To n - 1
    sig = (xin(i) - xin(i - 1)) / (xin(i + 1) - xin(i - 1))
    p = sig * yt(i - 1) + 2
    yt(i) = (sig - 1) / p
    u(i) = (yin(i + 1) - yin(i)) / (xin(i + 1) - xin(i)) - (yin(i) - yin(i - 1)) / (xin(i) - xin(i - 1))
    u(i) = (6 * u(i) / (xin(i + 1) - xin(i - 1)) - sig * u(i - 1)) / p
Next i

qn = 0
un = 0
yt(n) = (un - qn * u(n - 1)) / (qn * yt(n - 1) + 1)

For k = n - 1 To 1 Step -1
    yt(k) = yt(k) * yt(k + 1) + u(k)
Next k

klo = 1
khi = n

Do
    k = khi - klo
    If xin(k) > x Then
        khi = k
    Else
        klo = k
    End If

    k = khi - klo
Loop While k > 1

h = xin(khi) - xin(klo)
a = (xin(khi) - x) / h
b = (x - xin(klo)) / h
SpLine = a * yin(klo) + b * yin(khi) + ((a ^ 3 - a) * yt(klo) + (b ^ 3 - b) * yt(khi)) * (h ^ 2) / 6

endnow:

End Function
4

1 に答える 1

0

誰も答えなかったので、他の方法を見つける必要があると思いました。

私の値は放物線の y に似ているので、この曲線を近似する 3 次多項式の係数を見つけました。

これは、セル 1x4 の配列に数式を入力することによって作成されます。

worksheet.MyRange1x4.formulaarray = "=LINEST(" & MyX & "," & MyY & "^{1,2,3})"

MyX と MyY は同じ次元 (1 xn) の配列です。

結果は、多項式の "a"、"b"、"c"、および "d" です。 ax^3+bx^2+cx+d

ここで、最大値の X を見つけたい場合は、デルタを解く必要がありました。これは、2 次方程式の 1 つに非常によく似ています。Delta = (-b (+/-) SQRT(b^2 - 3ac))/3a

デルタを強制的に 0 にすると、2 つの可能な X 値が得られました (+/-SQRT のため)。そのうちの 1 つは間違いなくすべてのほぼ最大値の範囲内にあり、もう 1 つはこの範囲では完全に間違っていました。

正しいものを選択すると、最終的にXが見つかりました。

于 2013-10-25T12:44:02.750 に答える