My interpolation works for one set, but does not work for the other

Question

Ismail Qeshta 2017-11-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/366566-my-interpolation-works-for-one-set-but-does-not-work-for-the-other

评论： Ismail Qeshta 2017-11-12

Hi,

I am using the following code for interpolation from data in files. The problem is that the code runs very nicely in the first set, X1, Y1, while it shows that the X and V values don't have the same length in the second set: X2, Y2. I can't figure out the problem. Both codes present the same type of results and have same column length. The files of both sets are attached herewith.

clear; clc;
y=load ('Y1.txt');
x=load ('X1.txt');
plot(x(:,2),-sum(y(:,2:11),2))
y1=-sum(y(:,2:11),2);
x1=x(:,2);
y3=unique(y1);
x3=unique(x1);
y2 = [1308.9, 2514.9, 4797.9]*1000;
x2 = interp1(y3, x3, y2, 'linear');
figure(1)
plot(x1, y1, '-g')
hold on
plot(x2, y2, 'bp')
hold off
grid
legend('Data', 'Interpolated Points', 'Location', 'NW')

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

Ismail Qeshta 2017-11-12

编辑：Ismail Qeshta 2017-11-12

在 MATLAB Online 中打开

Error message is the following:

Error using interp1>reshapeAndSortXandV (line 423)
X and V must be of the same length.
Error in interp1 (line 92)
[X,V,orig_size_v] = reshapeAndSortXandV(varargin{1},varargin{2});
Error in Interpolate (line 15)
x2 = interp1(y3, x3, y2, 'linear');

请先登录，再进行评论。

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

Answer 1

Jan 2017-11-12

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/366566-my-interpolation-works-for-one-set-but-does-not-work-for-the-other#answer_290665

编辑：Jan 2017-11-12

在 MATLAB Online 中打开

The message is clear:

X and V must be of the same length.

Set a breakpoint in the failing line and check the dimensions:

size(y3)
size(x3)

Or:

y = load('Y2.txt');
x = load('X2.txt');
y1 = -sum(y(:,2:11),2);
x1 = x(:,2);
y3 = unique(y1);
x3 = unique(x1);
size(y3)
size(x3)
 ans =
        5504           1
 ans =
        5501           1

Then an interpolation is not possible. A solution might be:

[x3, ix] = unique(x1);
y3       = y1(ix);

2 个评论
显示无隐藏无

Ismail Qeshta 2017-11-12

在 MATLAB Online 中打开

Hi Jan,

Yes. You are right. I have just checked it. Do you have any idea on how to solve this issue? Thanks.

a= size(y3)
b= size(x3)
  a =
          1220           1
  b =
          1215           1

Ismail Qeshta 2017-11-12

Thank you very much Jan. Your solution worked.

请先登录，再进行评论。

My interpolation works for one set, but does not work for the other

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

采纳的回答

2 个评论
显示无隐藏无

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

My interpolation works for one set, but does not work for the other

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

采纳的回答

2 个评论 显示 无隐藏 无

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

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

2 个评论
显示无隐藏无