生成strassen矩阵乘法matlab代码

以下是Strassen矩阵乘法的Matlab代码实现：

function C = strassen(A, B)
% Strassen矩阵乘法算法
% 输入：A、B为n阶矩阵，n为2的幂次方
% 输出：C=A*B

% 确定矩阵A和B的大小
[n, ~] = size(A);

% 当矩阵规模小于等于64时，直接使用传统矩阵乘法
if n <= 64
    C = A * B;
    return;
end

% 划分A、B、C矩阵
mid = n / 2;
A11 = A(1:mid, 1:mid);
A12 = A(1:mid, mid+1:n);
A21 = A(mid+1:n, 1:mid);
A22 = A(mid+1:n, mid+1:n);
B11 = B(1:mid, 1:mid);
B12 = B(1:mid, mid+1:n);
B21 = B(mid+1:n, 1:mid);
B22 = B(mid+1:n, mid+1:n);
C11 = zeros(mid);
C12 = zeros(mid);
C21 = zeros(mid);
C22 = zeros(mid);

% 递归计算7个矩阵
M1 = strassen(A11+A22, B11+B22);
M2 = strassen(A21+A22, B11);
M3 = strassen(A11, B12-B22);
M4 = strassen(A22, B21-B11);
M5 = strassen(A11+A12, B22);
M6 = strassen(A21-A11, B11+B12);
M7 = strassen(A12-A22, B21+B22);

% 计算C矩阵
C11 = M1 + M4 - M5 + M7;
C12 = M3 + M5;
C21 = M2 + M4;
C22 = M1 + M3 - M2 + M6;
C = [C11, C12; C21, C22];
end

其中，我们将矩阵A和B划分成四个大小相等的子矩阵，然后递归计算7个矩阵，最后再根据这7个矩阵计算出矩阵C。当矩阵规模小于等于64时，直接使用传统矩阵乘法。