qam64d.m

LDPC 64QAM仿真程序

function [f0,f1]=QAM64d(Z,Vn)
%Y=imag(Z);%AA=reshape(A,1,[])
%X=real(Z);
a10=[(-7:2:7)*1i-7,(-7:2:7)*1i-5,(-7:2:7)*1i-3,(-7:2:7)*1i-1];a11=[(-7:2:7)*1i+7,(-7:2:7)*1i+5,(-7:2:7)*1i+3,(-7:2:7)*1i+1];
a20=[(-7:2:7)*1i-7,(-7:2:7)*1i-5,(-7:2:7)*1i+5,(-7:2:7)*1i+7];a21=[(-7:2:7)*1i+3,(-7:2:7)*1i+1,(-7:2:7)*1i-1,(-7:2:7)*1i-3];
a30=[(-7:2:7)*1i-7,(-7:2:7)*1i-3,(-7:2:7)*1i+3,(-7:2:7)*1i+7];a31=[(-7:2:7)*1i-5,(-7:2:7)*1i-1,(-7:2:7)*1i+1,(-7:2:7)*1i+5];
a40=[(-7:2:7)-7i,(-7:2:7)-5i,(-7:2:7)-3i,(-7:2:7)-1i];a41=[(-7:2:7)+7i,(-7:2:7)+5i,(-7:2:7)+3i,(-7:2:7)+1i];
a50=[(-7:2:7)-7i,(-7:2:7)-5i,(-7:2:7)+5i,(-7:2:7)+7i];a51=[(-7:2:7)-3i,(-7:2:7)-1i,(-7:2:7)+1i,(-7:2:7)+3i];
a60=[(-7:2:7)-7i,(-7:2:7)-3i,(-7:2:7)+3i,(-7:2:7)+7i];a61=[(-7:2:7)-5i,(-7:2:7)-1i,(-7:2:7)+1i,(-7:2:7)+5i];  
f0=[];
for i=1:length(Z)
f10(i)=sum(exp(-abs(Z(i)-a10).^2/Vn));
f11(i)=sum(exp(-abs(Z(i)-a11).^2/Vn));
F10=f10(i)/(f10(i)+f11(i)+eps);

f20(i)=sum(exp(-abs(Z(i)-a20).^2/Vn));
f21(i)=sum(exp(-abs(Z(i)-a21).^2/Vn));
F20=f20(i)/(f20(i)+f21(i)+eps);

f30(i)=sum(exp(-abs(Z(i)-a30).^2/Vn));
f31(i)=sum(exp(-abs(Z(i)-a31).^2/Vn));
F30=f30(i)/(f30(i)+f31(i)+eps);

f40(i)=sum(exp(-abs(Z(i)-a40).^2/Vn));
f41(i)=sum(exp(-abs(Z(i)-a41).^2/Vn));
F40=f40(i)/(f40(i)+f41(i)+eps);

f50(i)=sum(exp(-abs(Z(i)-a50).^2/Vn));
f51(i)=sum(exp(-abs(Z(i)-a51).^2/Vn));
F50=f50(i)/(f50(i)+f51(i)+eps);

f60(i)=sum(exp(-abs(Z(i)-a60).^2/Vn));
f61(i)=sum(exp(-abs(Z(i)-a61).^2/Vn));
F60=f60(i)/(f60(i)+f61(i)+eps);

f0=[f0;F10,F20,F30,F40,F50,F60];
end
f0=reshape(f0,1,[]);
f1=1-f0;
%u=(sign(f1-0.5)+1)/2;