interp1 monotonic increasing error

2 次查看（过去 30 天）
Brasco , D.
Brasco , D. 2014-7-16
编辑： Roger Parkyn 2014-10-8
hi guys i got two arrays to interpolate (interp1). but i got an error :
>>Ct_related=interp1(J_00_B,J_00_Ct,10)
Error using griddedInterpolant
The grid vectors are not strictly
monotonic increasing.
Error in interp1 (line 191)
F =
griddedInterpolant(X,V,method);
how can i solve this ? any help?

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

采纳的回答

Star Strider
Star Strider 2014-7-16
If you run this code:
difB = find([NaN; diff(J_00_B)] <= 0);
Q = [[difB-2:difB+2]' J_00_B(difB-2:difB+2)]
you will see that the values of J_00_B for indices 138 and 139 are the same. The independent variable has to be strictly monotonically increasing, so consecutive duplicate values are not allowed. See ‘Strictly Monotonic’ in the ‘More About’ section of the interp1 documentation. (The NaN in the diff call is to make the index references correct.)
Correcting that by adding a small value to the second duplicate:
J_00_B(difB) = J_00_B(difB)+1E-10;
Ct_related=interp1(J_00_B,J_00_Ct,10)
produces:
Ct_related =
96.1895e-003
  5 个评论
显示 3更早的评论
Star Strider
Star Strider 2014-8-14
My pleasure!
I apologise for not including the flipud lines the first time. I was interrupted while writing my Answer, and forgot them.
Roger Parkyn
Roger Parkyn 2014-10-8
编辑：Roger Parkyn 2014-10-8
I have used that same trick (i.e. using diff to find duplicated x-data and then adding a very small increment to the duplicates) but yesterday I found a weakness in this approach: if you have more than two identical x-values then adding the increment to all the duplicates (i.e. two or more duplicates) still leaves some of them the same. With my data it was OK just to eliminate duplicates so that is what I did (may not suit all purposes of course), e.g.:
idx = diff([-9.99; x_data]) == 0.0; % Find duplicate x-values
zz = interp1(x_data(~idx), y_data(~idx), x_i); % Only use non-duplicated data

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Resizing and Reshaping Matrices 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by