[data1, text, raw] = xlsread('振型','Sheet1'); ZX = data1; m=size(ZX,1); c=5 ; d=[1:m,]; for i=1:m-c E=ZX*(ZX'*ZX)^(-1)*ZX'; EA=diag(E); [a,b]=min(EA); ZX(b,:)=[]; d(:,b)=[]; end M = ZX ; [m, n] = size(M); MAC = zeros(n, n); for i = 1:n for j = 1:n a = M(:, i); b = M(:, j); MAC(i, j) = (a' * b)^2 / ((a' * a) * (b' * b)); end end disp(MAC); diag_MAC = diag(MAC); % MAC矩阵的对角元素 non_diag_MAC = MAC - diag(diag_MAC); % MAC矩阵的非对角元素 max_non_diag_MAC = max(non_diag_MAC(:)); % 最大非对角元素值 disp(max_non_diag_MAC);

如果你的目标是将上述代码中的`c`变为`c=1:15`，你需要使用一个外部循环来遍历每个`c`的值，并在每个循环迭代中执行相同的操作。下面是修改后的代码示例：

[data1, text, raw] = xlsread('振型','Sheet1');
ZX = data1; 
m = size(ZX, 1);
c_range = 1:15;

for c = c_range
    d = 1:m;
    
    for i = 1:m-c
        E = ZX * (ZX' * ZX)^(-1) * ZX';
        EA = diag(E);
        [a, b] = min(EA);
        ZX(b, :) = [];
        d(:, b) = [];
    end
    
    M = ZX; 
    [m, n] = size(M);
    MAC = zeros(n, n);
    
    for i = 1:n
        for j = 1:n
            a = M(:, i);
            b = M(:, j);
            MAC(i, j) = (a' * b)^2 / ((a' * a) * (b' * b));
        end
    end
    
    disp(['当c=', num2str(c), '时，计算得到的MAC矩阵：']);
    disp(MAC);
    
    diag_MAC = diag(MAC); % MAC矩阵的对角元素
    non_diag_MAC = MAC - diag(diag_MAC); % MAC矩阵的非对角元素
    max_non_diag_MAC = max(non_diag_MAC(:)); % 最大非对角元素值
    
    disp(['当c=', num2str(c), '时，最大非对角元素值：']);
    disp(max_non_diag_MAC);
end

在修改后的代码中，我们使用外部循环`for c = c_range`来遍历每个`c`的值。在每个循环迭代中，我们执行与原始代码相同的操作，并在每次迭代之后输出MAC矩阵和最大非对角元素值。

请注意，以上代码仅提供了一个示例，你可能需要根据具体情况进行适当的修改和调整。如果有任何进一步的问题，请随时提问。