请用MATLAB编写全局多项式拉格朗日插值函数

好的，以下是MATLAB代码：

function [P] = global_poly_lagrange(x,y)
% 全局多项式拉格朗日插值函数
% x: 插值节点
% y: 插值节点对应的函数值

n = length(x); % 节点个数
syms t; % 符号变量

% 计算基函数
L = ones(n,1);
for i = 1:n
    for j = 1:n
        if i ~= j
            L(i) = L(i) * (t - x(j)) / (x(i) - x(j));
        end
    end
end

% 构造插值多项式
P = 0;
for i = 1:n
    P = P + y(i) * L(i);
end

% 化简多项式
P = simplify(P);
end

使用方法：输入要插值的节点和对应函数值，即可得到插值多项式P。