x = [-5:-1:5]; %% 画图用点 z = @(x) 1 + x.^2; m = @(x) 1 ./ z(x); X1 = [-5:10/5:5]; Y1 = m(X1); %% 5次多项式的插值点 y1 = Lagrange(X1,Y1,5,x); X2 = [-5:10/10:5]; Y2 = m(X2); %% 10次多项式的插值点 y2 = Lagrange(X2,Y2,10,x); X3 = [-5:10/20:5]; Y3 = m(X3); %% 20次多项式的插值点 y3 = Lagrange(X3,Y3,20,x); plot(x,m(x),'r',x,y1,'b',x,y2,'m',x,y3,'c'); legend('1/(1+X^2)','L_5(x)','L_{10}(x)','L_{20}(x)'); %% Larange Interpolation % 输入： X 插值点 % Y 插值点函数值 % n 插值次数 % x 逼近点 % y 逼近值 function y = Lagrange(X,Y,n,x) if n >= length(X) fprintf('错误：插值点不够\n'); return end m = length(x); y = zeros(m,1); for k = 0 : n for i = 1:m y(i) = y(i) + Y(k+1)*prod(x(i)-X([1:k,k+2:end]))/prod(X(k+1)-X([1:k,k+2:end])); end end end legend 函数的输入参数那里不正确

根据警告信息，可以看出是 legend 函数的输入参数不正确。在 legend 函数的第一个输入参数中，应该是一个包含每个图形的标签的 cell 数组，而不是将所有标签放在一个字符串中。

因此，可以将 legend 函数的输入参数修改如下：

legend({'1/(1+X^2)','L\_5(x)','L\_{10}(x)','L\_{20}(x)'});

其中，用花括号 {} 将所有标签放在一个 cell 数组中，并且用反斜杠转义特殊字符，如下划线。