LDPC BP算法

function decision=SPA（LLRrow_indexcol_indexrwcwMaxiterationH）;
LLRout=zeros（size（LLR））;%分配空间
ErrorFlag = 1;
times=0;
node_to_check=zeros（size（col_index））;%分配空间
check_to_node=zeros（size（row_index））;%分配空间
%% 初始化
for i=1:size（node_to_check1） %矩阵第一维的size，也就是显示行123...
    for j=1:cw（i）
        node_to_check（ij）=LLR（i）;
    end
end
%% 
while and（times < MaxiterationerrorFlag == 1）

		times=times+1;
			
		%pass information from checks to Variable
        for i=1:size（row_index1）
        for j=1:rw（i）
            temp=1;
            for k=1:rw（i）
                if row_index（ik）~=row_index（ij）
                    for t=1:cw（row_index（ik））
                    if col_index（row_index（ik）t）==i
                        temp=temp*tanh（0.5*node_to_check（row_index（ik）t））;
                    end
                    end
                end
            end
            check_to_node（ij）=2*atanh（temp）;
        end
        end
       %LLR total 
       for i=1:size