Как SolvStat Macro в Excel вычисляет частные производные? - PullRequest
Новая
       19

Как SolvStat Macro в Excel вычисляет частные производные?

0 голосов
/ 28 июня 2019

Я работаю с SolvStat Macro от E.J. Билло из книги «Excel для химиков», и я не могу понять, как вычисляются частные производные как что-либо, кроме 0, поскольку я не вижу явных назначений, которые бы изменили значение любого из массивов YCalc () значения или значения диапазона calc_ys. Следующий код извлечен из макроса: 

Dim YObsd(), YCalc(), ParmValu(), PartialDeriv(), Product(), ProdArray()

Set calc_ys = Application.InputBox("Podaj zakres wartości Y obliczonych z modelu." & msg1, _
            "SOLVER STATISTICS - Krok 2 z 4", , , , , , 8)

ReDim YObsd(N), YCalc(N)
x = 1
For Each F In calc_ys
   YCalc(x) = F.Value
   x = x + 1
Next

Set Parms = Application.InputBox("Zaznacz komórki zawierające  współczynniki obliczone metodą najmniejszych kwadratów przez Solvera." & msg2, _
"SOLVER STATISTICS – Krok  3 z 4", , , , , , 8)
N3 = Parms.Count
x = 1
For Each cell In Parms
   ParmValu(x) = cell.Value
   x = x + 1
Next
RMSD = Sqr(SSresiduals / (N - N3))

'Calculate table of partial differentials
increment = 0.000001
'1E-6 seems to be optimum value for increment for numerical differentiation.
'1E-3 is too large, 1E-12 is too small.  1E-9 gives results almost identical to 1E-6.
y = 1
For Each parm In Parms
   parm.Value = parm * (1 + increment)     'Increase regression coeffs by a small increment.
   If parm = 0 Then parm.Value = 1E-100    'If cell contains zero exactly, replace it with very small value
   CheckErrorSum = 0
   x = 1
   For Each cell In calc_ys
      PartialDeriv(x, y) = (cell - YCalc(x)) / (parm * increment)      'Partial deriv = delta(function)/delta(parameter)
      CheckErrorSum = CheckErrorSum + PartialDeriv(x, y)   'This sum used only for error checking.
      x = x + 1
   Next cell
   parm.Value = ParmValu(y)    'Restore original parameter value
   If CheckErrorSum = 0 Then
      MsgBox "Błąd w obliczeniach macierzy." & Chr(13) & Chr(13) & 
      "Najpowszechniejsze błędy: " & Chr(13) & Chr(13) & _
         "1.  Niepoprawny wybór komórek Y(obl)." & Chr(13) & Chr(13) & _
         "2.  Niepoprawny wybór komórek zawierających współczynniki regresji 
         " & Chr(13) & Chr(13) & "Zatrzymanie programu.", 16, "UNKNOWN ERROR"
      Exit Sub
   End If
   y = y + 1
Next parm

Значение YCalc (x) содержит значение, которое было назначено из диапазона в Excel; однако, ячейка (элемент в calc_ys) имеет немного другое значение, но из кода не похоже, что должно быть какое-либо различие. Чтобы вычислить частные производные, должна быть разница, но я не могу понять, как / почему этот код выдает разные значения для этих двух переменных.

Вот ссылка на полный код: http://ue.poznan.pl/data/upload/articles/20140213/330de4761384630218/solvstat.doc

1 Ответ

0 голосов
/ 01 июля 2019

Поскольку calc_ys - это диапазон ячеек, он не содержит отображаемые значения в электронной таблице, а формулы, которые вычисляют отображаемые значения. Эти формулы ссылаются на параметры, которые изменяются в разделе частных производных; следовательно, ячейки calc_ys в Range также модифицируются.

...