Дифференциация буфера с Delphi

Question

Периодически получаю буфер данных, который содержит ряд значений, которые являются фиксированным расстоянием во времени. Мне нужно их отличить. Это Су давно я сделал исчисление в школе ....

То, что я придумал это:

function DifferentiateBuffer(
    ABuffer: TDoubleDynArray; AVPS: integer): TDoubleDynArray; 
var 
    i: integer; 
    dt, dy: double; 
begin 
    if (AVPS = 0) then exit; 

    // calc the delta time 
    dt := 1/AVPS; 

    // run through the whole buffer 
    for i := 0 to high(ABuffer) do begin 
    if (i = 0) then 
     if (IsNan(FLastDiffValue) = false) then 
     dy := ABuffer[0] - FLastDiffValue 
     else 
     dy := ABuffer[0] 
    else 
     dy := Abuffer[i] - ABuffer[i -1]; 

    result[i] := dy/dt 
    end; 

    // remember the last value for next time 
    FLastDiffValue := ABuffer[high(ABuffer)]; 
end; 

AVPS это значения в секунду. Типичное значение для этого должно составлять от 10 до 100. Длина буферов обычно составляет от 500 до 1000 значений.

Я вызываю буфер время от времени с новыми данными, которые непрерывны с предыдущим блоком данных, следовательно, сохраняя последнее значение блока в следующий раз. Это последнее значение устанавливается в NAN в конструкторе.

Это то, что я сделал правильно? т. е. будет ли он правильно дифференцировать значения?

Мне также нужно интегрировать подобные данные ... что это может быть похоже?

Answer 1

В этом я не могу найти недостатки.

Первое значение, возвращенное, при первом вызове будет (ABuffer[0] - 0.0)/dt, которое основано на предположении, что сигнал начинается с нуля. Я полагаю, это то, что вы намереваетесь.

Теперь, вместо того чтобы попросить сообщество переполнения стека проверить ваш код, вы можете сделать гораздо лучше себя. Вы должны проверить код, чтобы доказать его точность. Протестируйте его с помощью модульной системы тестирования, например DUnitX. Подайте значения функций, для которых вы можете прогнозировать выход. Например, подайте значения из y = x , или y = sin (x).

Другим большим преимуществом написания тестов является то, что они могут быть выполнены снова и снова. По мере разработки кода вы рискуете ввести ошибки. Сегодня код может быть прав, но кто знает, будет ли он по-прежнему правильным после внесения изменений завтра. Если у вас есть сильный тест на месте, вы можете защитить от ошибок, возникших во время обслуживания.

Один комментарий к стилю заключается в том, что вы никогда не должны тестировать = false или = true. Оператор if работает с булевыми выражениями, поэтому сравнение с булевым значением всегда довольно бессмысленно. Я хотел бы написать тест, как это:

if not IsNan(FLastDiffValue) then 
    dy := ABuffer[0] - FLastDiffValue 
else 
    dy := ABuffer[0] 

или как это:

if IsNan(FLastDiffValue) then 
    dy := ABuffer[0] 
else 
    dy := ABuffer[0] - FLastDiffValue