MATLAB 读取 Excel 文件并进行四次拟合 - 计算 R 方

MATLAB 读取 Excel 文件并进行四次拟合 - 计算 R 方

本代码使用 MATLAB 读取给定路径的 Excel 文件数据，并对指定列进行四次项函数拟合，最后计算调整后的 R 方。

% 指定文件路径和文件名
filePath = 'D:\桌面\暑假数模\第六次模拟赛题\B\111.xlsx';
sheetName = '题目1-1';

% 读取Excel文件中的数据
[num, txt, raw] = xlsread(filePath, sheetName);

figure;

for tar=1:21
    x1=[];y1=[];y2=[];
    for i=1:114
        if raw{i,1}==tar
            tmp=[x1,raw{i,2}];x1=tmp;
            tmp=[y1,raw{i,3}];y1=tmp;
            tmp=[y2,raw{i,4}];y2=tmp;
        end
    end
    % 拟合四次项函数
    p = polyfit(x1, y1, 4);

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

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

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

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

    disp(R_square)
end

使用说明:

1. 确保文件路径和文件名是正确的，并且 Excel 文件中的数据与代码中的索引匹配。
2. 将代码中的 filePath 和 sheetName 替换为您的实际文件路径和工作表名称。
3. 运行代码即可获得每个目标 tar 对应的调整后的 R 方值。

注意:

• 由于文件路径中的反斜杠需要转义，因此代码中使用了双反斜杠 \。
• 该代码假设 Excel 文件中第一列的数据对应于 tar 值，第二列对应于 x1 值，第三列对应于 y1 值，第四列对应于 y2 值。
• 如果您的 Excel 文件数据结构不同，请根据实际情况修改代码中的索引。
• 代码中使用了 xlsread 函数读取 Excel 文件，请确保您的 MATLAB 版本支持该函数。

该代码可以帮助您快速分析 Excel 文件中的数据，并使用四次项函数进行拟合，从而获得更准确的结果。

如果您有任何问题或需要进一步改进，请随时提出。