How to set the lower bound of extrapolation?

Question

George McFadden 2012-7-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/43747-how-to-set-the-lower-bound-of-extrapolation

评论： charles alexander 2021-7-29

Hi!

Is it possible to set a lower bound for extrapolating using the interp2 function? I am dealing with nitrogen concentrations, and a negative concentration cannot exist. However, currently the extrapolation results have negative values. So I do not want the extrapolation to go less than 0. I do not desire to set an upperbound limit.

Thanks! George

3 个评论
显示 1更早的评论隐藏 1更早的评论

George McFadden 2012-10-30

The interpolation and extrapolation are both producing below 0 values. Thanks!

charles alexander 2021-7-29

yes, its is possible to set a lower bound.

you will need to find the value of c from the linear discriminant model.

extrapolating will be (c1+c2)/100*values less than 30,

find the percentage value of extrapolation, add the (percentage value to the value of extrapolation:c3).

(c1+c2)/c3*fill_value.

请先登录，再进行评论。

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

Answer 1

D G 2012-7-16

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/43747-how-to-set-the-lower-bound-of-extrapolation#answer_53687

在 MATLAB Online 中打开

You could set the boundary limit manually after the calculation like this.

concentration(concentration < 0) = 0;

I create a function that does just this.

function value = clip(value, clip, type)
  %CLIP Clip values
  % Adjust values in a vector such that no value extends beyond the
  % specified type.
  %
  %   Inputs:
  %       Clip min and max together:
  %           value (vector - 1 by x)
  %           clip  (vector - 1 by x > 1)
  %       Clip min and max seperately
  %           value (vector - 1 by x)
  %           clip  (scalar)
  %           type  (char - min/max)
  %
  %   Outputs:
  %       value - The clipped values
    if length(clip) > 1
        value(max(clip) < value) = max(clip);
        value(value < min(clip)) = min(clip);
    else
        switch lower(type)
            case 'min'
                value(value < clip) = clip;
            case 'max'
                value(clip < value) = clip;
        end %switch
    end
end %clip