У меня есть функция (сложный массив значений uul(:)), для которой я хотел бы вычислить ее производную. Я просто использую центральные конечные разности
uulder=(uul(3:end)-uul(1:end-2))./deltaeta(2:end-1)'/2;
uulderabs=(abs(uul(3:end))-abs(uul(1:end-2)))./deltaeta(2:end-1)'/2;
Я рисую это:
>> plot(eta,abs(uul))
>> hold on
>> plot(eta(1:end-2),abs(uulder))
>> hold on
>> plot(eta(1:end-2),uulderabs)
Функция и производные
Кажется, красная кривая неправильная, а желтая правильная. Однако мне понадобится вектор, абсолютное значение которого дает желтую кривую (и не принимает заранее абсолютные значения, как с uulderabs). Мне это нужно, потому что с желтой кривой я теряю реальную и мнимую части, и я хотел бы иметь возможность их использовать. Что я делаю не так?