尝试在MATLAB中将值存储在向量中时重塑错误

Question

我使用LaGrange的插值方法在Matlab中编写了以下插值代码，但是当尝试将值存储在向量中（L（i）= ....）时，我会不断遇到整形错误（元素不得更改）。

我不确定该怎么做，如果在此之前有人问过这个问题，我深表歉意，但是我已经尝试过搜索，但是遇到此问题的人实际上正在尝试使用重塑功能。

我要做的就是将J循环的结果存储在向量中。

syms X
x = [-1, 0, 2]
y = [4, 1, -1]

i = 1
j = 1
P = 0

while i < 4
  L(i) = 1
  while j < 4
    if (i ~= j)
      L(i) = L(i) * (X - x(j)) / (x(i) - x(j))
    end
    j = j + 1
  end
  P = P + y(i) * L(i)
  i = i + 1
end

Answer 1

您将L初始化为具有L(i) = 1的双精度数组，但随后尝试在其中存储符号内容。 您需要将L初始化为符号数组，因为要在其中存储相同的内容。 因此，在循环之前添加以下行：

L = sym('L', [3 1]);

还要避免将i和j用作变量名，因为它们是为虚数保留的。

Answer 2

我在使用此代码时遇到了与您类似的错误，并且能够通过提前初始化L来解决它。 这是经过修改的原始代码（ clear all代码只是为了确保没有发生任何有趣的事情）。

clear all;
syms X L;
L(1:3) = [1 1 1];
x = [-1, 0, 2]
y = [4, 1, -1]

i = 1
j = 1
P = 0

while i < 4
  while j < 4
    if (i ~= j)
      L(i) = L(i) * (X - x(j)) / (x(i) - x(j))
    end
    j = j + 1
  end
  P = P + y(i) * L(i)
  i = i + 1
end