constbuilder.m

离散多载波传输(DMT)的MATLAB程序及注释

function y = constbuilder(N)
b=ceil(log2(N));
p=zeros(N,b);
for k=0:(N-1)
        u=dec2bin(k);
        w=zeros(1,length(u)); 
        for i=1:length(u)
            w(i)=str2num(u(i));
        end
        v=[zeros((b-length(w)),1)',w];
        p(k+1,:)=v;
end
if N==16
    %a=[.354+.354*j;sqrt(2)/2;sqrt(2)/2*j;-.354+.354*j;-sqrt(2)/2*j;.354-.354*j;-.354-.354*j;-sqrt(2)/2];
    %b=exp(j*2*pi/8.*(0:7)');
    %M=[p,[a;b]];  
    %a=[-.25-.25*j;-.75-.25*j;.-.25-.75*j;-.75-.75*j;.25-.25*j;.25-.75*j;.75-.25*j;.75-.75*j;-.25+.25*j;-.25+.75*j;-.75+.25*j;-.75+.75*j;
       %.25+.25*j;.75+.25*j;.25+.75*j;.75+.75*j];
     a=0.5*[-0.5-0.5*j;-1.5-0.5*j;-0.5-1.5*j;-1.5-1.5*j;0.5-0.5*j;0.5-1.5*j;1.5-0.5*j;1.5-1.5*j;-0.5+0.5*j;-0.5+1.5*j;-1.5+0.5*j;-1.5+1.5*j;
         0.5+0.5*j;1.5+0.5*j;0.5+1.5*j;1.5+1.5*j]; 
   M=[p,a];
else
end    
y=matlaberr(M);