encode_伍坚.m

来自「关于LDPC编码的MATLAB实现仿镇,求关于LDPC码的FPGA实现」· M 代码 · 共 20 行

function [x]=encode(H,s)
% H为校验矩阵，s为信息序列
[n,m]=size(H);
G=H(1:n,[m]);
g=6;
A=H(1:(n-g),1:(m-n));
B=H(1:(n-g),(m-n+1):(m-n+g));
C=H((n-g+1):n,1:(m-n));
D=H((n-g+1):n,(m-n+1):(m-n+g));
E=H((n-g+1):n,(m-n+g+1):m);
T=H(1:(n-g),(m-n+g+1):m);
p=inv(-E*inv(T)*B+D)*(-E*inv(T)*A*s'+C*s');
p1=round(p);            %将p转换为离其最近的整数
p2=mod(p1,2);
p2=p2';
q=inv(T)*(A*s'+B*p2');
q1=round(q);
q2=mod(q1,2);
q2=q2';
x=[s,p2,q2];                                          %编码结果;