Matlab：非单调插值

Question

我在使用非单调函数进行插值时遇到问题。

x=[34,35,36,37,38,39,41,42,43,44,45,46]
y=[0.33 0.166 0.25 0.33 0.5 0.58 0.66 0.66 0.5 0.58 0.8 0.66]

我希望能够对对应于 0.25 的 Y 值的 X 值进行插值。

绘制(x,y)时，我可以看到有两个点对应于0.25的Y值，因此正常插值（使用interp1）将不起作用。

Error: "The grid vectors are not strictly monotonic increasing."

有没有办法对两个点进行插值然后计算它们的平均值，我找不到适合我的问题的方法。

Answer 1

您可以select独特点数为：

% x and y are your data available
[xN, index] = unique(x);
yN = y(index);

% Now you have only unique x and y values stored in xN and yN

请注意，这将丢弃您的重复数据点。

Answer 2

您可以先计算重复的平均值，然后进行插值：

%test data
x=[34,35,36,37,38,39,41,42,43,44,45,46];
y=[0.33 0.166 0.25 0.33 0.5 0.58 0.66 0.66 0.5 0.58 0.8 0.66];

% sort values
[~,idx] = sort(y);
x = x(idx);
y = y(idx);

% use all values only once. Average the values
y_unique = unique(y);
x_unique = zeros(size(y_unique));
for i = 1:length(y_unique)
   x_unique(i) =  mean(x(y==y_unique(i)));
end

% test
interp1(y_unique,x_unique,0.33)

Answer 3

遍历数组并检测是否传递了查询值。然后取平均值。

%test data
x=[34,35,36,37,38,39,41,42,43,44,45,46];
y=[0.33 0.166 0.25 0.33 0.5 0.58 0.66 0.66 0.5 0.58 0.8 0.66];

%value to query
yq = 0.25;
result = [];
for i = 1:length(y)-1
    if (y(i) <= yq && yq < y(i+1)) || (y(i) > yq && yq >= y(i+1))
        result(end+1) = interp1([y(i),y(i+1)],[x(i),x(i+1)],yq);
    end
end

xq = mean(result);

Matlab：非单调插值

Matlab: non-monotonic interpolation

matlab

interpolation