0 个投票

Why does:
diffpressure = gradient(pressure)./gradient(density);
and
diffpressure = diff(pressure)./diff(density);
give drastically different results? The derivative using each method are completely different. See attached figure

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

 采纳的回答

Matt J
Matt J 2018-11-5
编辑：Matt J 2018-11-5

2 个投票

gradient uses central differences while diff(x) uses right-hand differences x(i+1)-x(i).

12 个评论
显示 10更早的评论 隐藏 10更早的评论

Benjamin
Benjamin 2018-11-5
Look at my attached figure, should they really be that different?
Matt J
Matt J 2018-11-5
I don't see why not, but on the other hand the picture isn't all that informative. You should provide the input data in a .mat file.
Benjamin
Benjamin 2018-11-5
Check file now. It shows the original function, and the derivative using diff and derivative using grad. Surely, the derivative is not that sensitive to the method. I have hundreds of points on the curve.
Matt J
Matt J 2018-11-5
I see no change in the attachment.
Benjamin
Benjamin 2018-11-5
Sorry, check now
Matt J
Matt J 2018-11-5
编辑：Matt J 2018-11-5
You have only plotted one of the original functions. We cannot tell if this is density or pressure and without both we also cannot gain any insight into the ratio of their gradients.
Once again, the most helpful thing would be to attach your numeric data in a .mat file.
Matt J
Matt J 2018-11-5
编辑：Matt J 2018-11-5
My preliminary guess would be that the gradient of density is tending toward zero, where the ratio calculation gradient(pressure)./gradient(density) becomes numerically unstable.
Benjamin
Benjamin 2018-11-5
The y-axis is pressure, x-axis is density. Ok, I'll try and attach, but quick question: If I have x values and yvalues, how would I differentiate y wrt to x? do I do diff(y)/diff(x) or just do diff(y) and plot it as plot(x,diff(y))
Benjamin
Benjamin 2018-11-5
Perhaps the better question is, if I have numerical X and Y data, what way should I differentiate y with respect to x?
Benjamin
Benjamin 2018-11-5
Do you think diff(y) / diff(x) is the way to go here?
Matt J
Matt J 2018-11-5
编辑：Matt J 2018-11-5
If the X-data are not equi-spaced, use diff().

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by