根据H在有限域GF(2^m)上求解生成矩阵G

原理

有时间再补充。

注1:使用高斯消去法。如果Py不为单位阵,则说明进行了列置换,此时G不是系统形式。

注2:校验矩阵H必须是行满秩才存在对应的生成矩阵G,且生成矩阵G通常不唯一。

matlab实现:只做列置换,不做行置换

Matlab 复制代码
function [G, Px, Py] = Gaussian_Elimination_in_GFq(H, q)
    
    % initial
    [m, n] = size(H);
    H_stair = mod(H,q);
    Px = eye(m); 
    Py = eye(n);
    
    for i=1:m
        for j=i:n
            if gcd(H_stair(i,j), q) == 1
                break;
            else
                if j==n
                    error('Gaussian_Elimination_in_GFq: The H is not full rank in GF(%d).',q);
                else
                    continue;
                end
            end
        end

        if j ~= i
            H_stair(:, [i, j]) = H_stair(:, [j, i]);
            Py(:, [i, j]) = Py(:, [j, i]);
        end
    end
    
    for i=1:m
        [~, x, ~] = gcd(H_stair(i,i), q);
        inv_i = mod(x, q);

        H_stair(i, :) = mod(inv_i * H_stair(i, :), q);
        Px(i, :) = mod(inv_i * Px(i, :), q);

        for j=1:m
            if j~=i && H_stair(j,i)~=0
                factor=H_stair(j, i);
                H_stair(j, :) = mod(H_stair(j, :) - factor * H_stair(i, :), q);
                Px(j, :) = mod(Px(j, :) - factor * Px(i, :), q);
            end
        end
    end

    parity_matrix = mod((q-1).*H_stair(:,m+1:n),q);
    G_fixed=[parity_matrix',eye(n-m)];
    G=G_fixed*Py';
end

matlab实现:先做行置换,若无逆元,再做列置换

Matlab 复制代码
function [G, Px, Py] = Gaussian_Elimination_in_GFq(H, q)
    
    [m, n] = size(H);
    H_stair = mod(H,q);
    Px = eye(m); 
    Py = eye(n);
    
    %% column find
    for j=1:m
        for i=j:m
            if gcd(H_stair(i,j), q) == 1
                break;
            else
                if i==m
                    %% row find
                    for k=j+1:n
                        if gcd(H_stair(j,k), q) == 1
                            break;
                        else
                            if k==n
                                error('Gaussian_Elimination_in_GFq: The H is not full rank in GF(%d).',q);
                            else
                                continue;
                            end
                        end
                    end
                    %% column exchange
                    H_stair(:,[j,k]) = H_stair(:,[k,j]);
                    Py(:,[j,k]) = Py(:,[k,j]);
                    i=j;
                    break;
                else
                    continue;
                end
            end
        end
        %% row exchange
        if i ~= j
            H_stair([j, i],:) = H_stair([i, j],:);
            Px([j, i],:) = Px([i, j],:);
        end
    end

    for i=1:m
        [~, x, ~] = gcd(H_stair(i,i), q);
        inv_i = mod(x, q);

        H_stair(i, :) = mod(inv_i * H_stair(i, :), q);
        Px(i, :) = mod(inv_i * Px(i, :), q);

        for j=1:m
            if j~=i && H_stair(j,i)~=0
                factor=H_stair(j, i);
                H_stair(j, :) = mod(H_stair(j, :) - factor * H_stair(i, :), q);
                Px(j, :) = mod(Px(j, :) - factor * Px(i, :), q);
            end
        end
    end

    parity_matrix = mod((q-1).*H_stair(:,m+1:n),q);
    G_fixed=[parity_matrix',eye(n-m)];
    G=G_fixed*Py';
    
end
相关推荐
海蓝可知天湛6 分钟前
利用Genspark自定义智能体:算法竞赛测试数据反推与生成工具
算法·aigc
BothSavage10 分钟前
Qwen3-VL-8B-Instruct推理测试transformer+sglang双版本
算法
尤超宇27 分钟前
YOLOv3 目标检测算法核心技术
算法·yolo·目标检测
cyclel1 小时前
散列表的小想法
算法
Code小翊1 小时前
堆的基础操作,C语言示例
java·数据结构·算法
余俊晖1 小时前
如何让多模态大模型学会“自动思考”-R-4B训练框架核心设计与训练方法
人工智能·算法·机器学习
Emilia486.1 小时前
【Leetcode&nowcode&数据结构】顺序表的应用
数据结构·算法·leetcode
一水鉴天1 小时前
整体设计 逻辑系统程序 之27 拼语言整体设计 9 套程序架构优化与核心组件(CNN 改造框架 / Slave/Supervisor/ 数学工具)协同设计
人工智能·算法
小年糕是糕手1 小时前
【数据结构】双向链表“0”基础知识讲解 + 实战演练
c语言·开发语言·数据结构·c++·学习·算法·链表
PyHaVolask2 小时前
数据结构与算法分析
数据结构·算法·图论