MATLAB自定义modifikuj函数异常求助：矩阵元素按规则替换失效

嘿，我帮你揪出这个MATLAB函数的问题啦！主要是素数判断的逻辑出了bug，咱们一步步拆解：

问题根源 #

你的内层循环判断素数的逻辑有两个关键错误：

1. 循环范围错误：你写的for k = 2:(A(i,j))会把素数本身算进去，比如数字2，k会取到2，这时候mod(2,2)==0会把prost设为0，导致程序把2错误判定为合数（赋值0），完全搞反了！
2. 冗余迭代+未及时终止：就算数字是合数，你也会循环到它本身，比如4会循环k=2、3、4，其实k=2的时候已经能判定是合数了，没必要继续跑后面的循环。

修正后的代码 #

我把你的代码调整了，修复了逻辑问题，还优化了效率：

function B = modifikuj(A)
[n,m] = size(A);
% 初始化输出矩阵，避免修改原输入矩阵（MATLAB函数的好习惯）
B = A;
for i = 1:n
    for j = 1:m
        val = B(i,j);
        if val == 1
            B(i,j) = 0.5;
        else
            prost = 1; % 默认标记为素数，除非找到因数
            % 素数判断只需要检查到sqrt(val)就够了，减少计算量
            for k = 2:floor(sqrt(val))
                if mod(val, k) == 0
                    prost = 0;
                    break; % 找到因数立刻跳出循环，不用白跑
                end
            end
            % 根据标记赋值
            B(i,j) = prost == 1 ? 1 : 0;
        end
    end
end

更简洁的MATLAB风格写法 #

其实MATLAB有内置的素数判断函数isprime，用向量化操作代替嵌套循环，代码更简洁效率也更高（尤其是处理大矩阵时）：

function B = modifikuj(A)
% 先初始化全0矩阵
B = zeros(size(A));
% 处理1的情况
B(A == 1) = 0.5;
% 处理素数的情况
B(isprime(A)) = 1;
% 剩下的元素就是合数，保持0即可
end

现在测试输入A = [1,2;3,4]，调用D = modifikuj(A)，就能得到你预期的[0.5, 1; 1, 0]啦！

内容的提问来源于stack exchange，提问作者yusano