写出matlab代码求一组数据的四次项拟合函数的调整后的R方

假设已有一组数据 x 和 y，要求四次项拟合函数的调整后的 R 方。

% 输入数据 x 和 y
x = [1, 2, 3, 4, 5];
y = [2, 4, 8, 16, 32];

% 拟合四次项函数
p = polyfit(x, y, 4);

% 计算拟合值
y_fit = polyval(p, x);

% 计算总平方和
SS_total = sum((y - mean(y)).^2);

% 计算残差平方和
SS_residual = sum((y - y_fit).^2);

% 计算调整后的 R 方
R_square = 1 - SS_residual / (length(y) - 1) / (SS_total / (length(y) - length(p) - 1));

这段代码首先使用 polyfit 函数拟合了一个四次项多项式，然后使用 polyval 函数计算拟合值。接下来，计算总平方和 SS_total 和残差平方和 SS_residual，最后使用这两个值计算调整后的 R 方。