matrix dimensions must agree

Question

Amanda McGovern 2019-11-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/489064-matrix-dimensions-must-agree

评论： Walter Roberson 2019-11-5

I am writing a lagrange interpolator function and when runnign the test code, my program fails for most tests and i get "matrix dimensions must agree" error. I am not sure where the error is. Can anybody tell? Also, is it appropriate for me to be using .* as opposed to just * operators? Thanks.

function p = lagrangeval(x,y,w)
%% lagrangeval evalutates the Lagrange interpolant for a set of knots (x,y)
%% Inputs: numbers x1,x2,..xn; values f(x1),f(x2)..f(xn)
%% Output: p - the value of the polynomial going through the n data points
 
function p = lagrangeval(x,y,w)
%% lagrangeval evalutates the Lagrange interpolant for a set of knots (x,y)
%% Inputs: numbers x1,x2,..xn; values f(x1),f(x2)..f(xn)
%% Output: p - the value of the polynomial going through the n data points
 
n=size(x,2);
k=size(w,2);
p=zeros(n,k);
for i=1:n
    p(i,1)=y(i);
end
for l=1:k
  for i=1:n-1
    for j=1:i
        p(i+1,j+1)=(((w(l)-x(i-j+1)).*p(i+1,j))-((w(l)-x(i+1)).*p(i,j)))./(x(i+1)-x(i-j+1));
    end
  end
end

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

Amanda McGovern 2019-11-4

x = 0:2*pi/5:2*pi; y = sin(x); w = [3 4]

x = 0:2*pi/8:2*pi; y = sin(x); w = [3 4 5];

Walter Roberson 2019-11-5

I do not encounter any error messages when I pass the above values to your function.

请先登录，再进行评论。

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

Answer 1

David Hill 2019-11-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/489064-matrix-dimensions-must-agree#answer_399816

在 MATLAB Online 中打开

Length of x and y must be equal since they are points. The following functions evaluations lagrange polynomial at each value of w (any length).

function p = lagrangeval(x,y,w)
a=length(x);
b=ones(1,a-1);
p=zeros(1,length(w));
for i=1:a
    p=p+arrayfun(@(z)y(i)*(prod(z*b-x(x~=x(i)))/prod(x(i)*b-x(x~=x(i)))),w);
end
end