How can I differentiate without decreasing the length of a vector?

63 次查看（过去 30 天）
Javad
Javad 2014-7-18
评论： John D'Errico 2020-3-10
I have some vectors and want to differentiate them up to second order. I don't want to use "diff" because it reduces the length of vector in higher orders! Is there any other function or method that I differentiate and keep the length of vector constant?

请先登录，再回答此问题。

采纳的回答

Jan
Jan 2014-7-18
gradient is smarter for calculating derivatives:
x = rand(1, 100);
d2 = gradient(gradient(x));
The Savtizky Golay smoothing filter can be applied to calculate a smoothed derivative by fitting polynmials to local parts of the signal. Look in the FileExchange for many different submissions:
http://www.mathworks.com/matlabcentral/fileexchange/index?utf8=%E2%9C%93&term=savitzky
  1 个评论
John D'Errico
John D'Errico 2020-3-10
Jan is correct, of course. I might only add one idea, to fit the data using a smoothing spline, then differentiate the spline, and evaluate the derivative spline at the original data points.
spl = csaps(x,y);
spld = fnder(spl);
yprimepred = fnval(spld,x);
As I've done it here, this uses tools from the curve fitting toolbox, though there are alternative ways to implement it too.

请先登录，再进行评论。

更多回答（1 个）

Daniel kiracofe
Daniel kiracofe 2014-7-18
My standard approach is to use 2nd order centered difference for the main part of the vector, and use first order forward and backward difference at the boundaries:
function d = cdiff(x, dt)
if (nargin<2)
dt =1 ;
end
d(1) = (x(2) - x(1)) / dt;
d(length(x)) = ( x(end) - x(end-1) ) / dt;
ndx = 2:(length(x)-1);
d(ndx) = (x( ndx+1) - x(ndx-1)) / (2 * dt);

类别

在 Help Center 和 File Exchange 中查找有关 Smoothing 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by