系数矩阵的列压缩存储CCS的代码matlab

下面是一个用 MATLAB现的系数矩阵列压缩存储（CCS）的示例代码：

function [values, row_ind, col_ptr] = compress_matrix(matrix)
    [m, n] = size(matrix);
    nnz = nnz(matrix);
    
    % 初始化存储压缩矩阵所需的数组
    values = zeros(nnz, 1);
    row_ind = zeros(nnz, 1);
    col_ptr = zeros(n+1, 1);
    
    % 压缩矩阵
    k = 1;
    for j = 1:n
        col_ptr(j) = k;
        for i = 1:m
            if matrix(i,j) ~= 0
                values(k) = matrix(i,j);
                row_ind(k) = i;
                k = k + 1;
            end
        end
    end
    col_ptr(n+1) = nnz + 1;
end

这个函数接受一个输入矩阵 `matrix`，并返回压缩后的系数矩阵的值数组 `values`、行索引数组 `row_ind` 和列指针数组 `col_ptr`。其中，`values` 存储非零元素的值，`row_ind` 存储非零元素的行索引，`col_ptr` 存储每一列第一个非零元素在 `values` 数组中的索引。

你可以将你的系数矩阵作为输入传递给这个函数，然后使用返回的压缩后的数组进行存储和处理。