polar-codes.rar

function [L,cc]=yima(M,n,e) %y=1:2^n; %Y=daoweixu(y,n); %LL=-ones(1,2^n); %for i=1:2^n % LL(i)=M(Y(i)); %end w0=1-e; w1=e; W0=e; W1=1-e; LL=daoweixu(M,n); for i=1:2^n if (LL(i)==0) LL(i)=w0/w1; else LL(i)=W0/W1; end end [~,P]=polarize(n,e); L=-ones(2^n,n+1)*0.1; L(:,n+1)=LL; u=-ones(1,2^n); S=-ones(2^n,n); %%%%%%%%% L S V 初始化为-1的矩阵 V=-ones(1,2^n-2); for j=1:n %指数项矩阵右上角初始赋值0 for i=1:2^(j-1) S(i,j)=0; end end l=1; while(l<=2^n) if (P(l)~=1) u(l)=0; else j=1; i=l; if (l>1) %%%%%%%% S(l,1)=u(l-1); end %%%%%%%% if (l>2) if (mod(l,2)==0) p=l-2; for t=1:p V(t)=u(t); end else p=l-1; for t=1:p V(t)=u(t); end end end while(j<n+1) if (rem(fix((i-1)/(2^(j-1))),2)==1) if (L(i-2^(j-1),j+1)==-0.1) %%%%%%%%%%%%%%%%%%%%%子程序2 这里 if (l>=3) p=fix((p+1)/2); for q=1:p V(q)=mod(V(2*q-1)+V(2*q),2); end i=i-2^(j-1); j=j+1; if(S(i,j)==-1) S(i,j)=V(p); %else % if (V(p)~=S(i,j)) % k=10000000 % end end i; j; V; else i=i-2^(j-1); j=j+1; end elseif (L(i-2^(j-1),j+1)~=-0.1&&L(i,j+1)==-0.1) %%%%%%%%%%%%%%%%%%%%%子程序4 if (l>=3) p=fix((p+1)/2); for q=1:p V(q)=V(2*q); end j=j+1; if(S(i,j)==-1) S(i,j)=V(p); % else % if (V(p)~=S(i,j)) % k=10000000 % end end i; j; V; else j=j+1; end elseif (L(i-2^(j-1),j+1)~=-0.1&&L(i,j+1)~=-0.1) L(i,j)=(L(i-2^(j-1),j+1)^(1-2*S(i,j)))*L(i,j+1); if L(i,j)==inf L(i,j)=1.000e+50; end if (j==1) %$$$$$$$$$$$$$$$$$$$$$$$ break; else j=1; i=l; if (l>2) if (mod(l,2)==0) p=l-2; for t=1:p V(t)=u(t); end else p=l-1; for t=1:p V(t)=u(t); end end end continue; end end else if (L(i,j+1)==-0.1) %%%%%%%%%%%%%%%%%%%%%子程序1 if (l>=3) p=fix((p+1)/2); for q=1:p V(q)=mod(V(2*q-1)+V(2*q),2); end j=j+1; if(S(i,j)==-1) S(i,j)=V(p); %else % if (V(p)~=S(i,j)) % k=10000000 % end end i; j; V; else j=j+1; end elseif (L(i,j+1)~=-0.1&&L(i+2^(j-1),j+1)==-0.1) %%%%%%%%%%%%%%%%%%%%%子程序3 if (l>=3) p=fix((p+1)/2); for q=1:p V(q)=V(2*q); end i=i+2^(j-1); j=j+1; if(S(i,j)==-1) S(i,j)=V(p); %else % if (V(p)~=S(i,j)) % k=10000000 % end end i; j; V; else i=i+2^(j-1); j=j+1; end elseif (L(i,j+1)~=-0.1&&L(i+2^(j-1),j+1)~=-0.1) L(i,j)=(L(i,j+1)*L(i+2^(j-1),j+1)+1)/(L(i,j+1)+L(i+2^(j-1),j+1)); if L(i,j)==inf L(i,j)=1.000e+50; end if (j==1) break; else j=1; i=l; if (l>2) if (mod(l,2)==0) p=l-2; for t=1:p V(t)=u(t); end else p=l-1; for t=1:p V(t)=u(t); end end end continue; end end end end if (L(l,1)>=1) u(l)=0; else u(l)=1; end end l=l+1; end L; S; cc=u;