function [out]=mapping_psdu(datain,length,rate);
p=[1,1,1,1,-1,-1,-1,1,-1,-1,-1,-1,1,1,-1,1,-1,-1,1,1,-1,1,1,-1,1,1,1,1,1,1,-1,1,...
1,1,-1,1,1,-1,-1,1,1,1,-1,1,-1,-1,-1,1,-1,1,-1,-1,1,-1,-1,1,1,1,1,1,-1,-1,1,1,...
-1,-1,1,-1,1,-1,1,1,-1,-1,-1,1,1,-1,-1,-1,-1,1,-1,-1,1,-1,1,1,1,1,-1,1,-1,1,-1,1,...
-1,-1,-1,-1,-1,1,-1,1,1,-1,1,-1,1,1,1,-1,-1,1,-1,-1,-1,1,1,1,-1,-1,-1,-1,-1,-1,-1];
data=double(datain);

if(isequal(rate,[0 0 0 0 0])|isequal(rate,[0 0 0 0 1]))
    for i=1:(length/50)
        for k=1:50
        sym((i-1)*200+k)=data((i-1)*50+k);
        sym((i-1)*200+101-k)=conj(data((i-1)*50+k));
        sym((i-1)*200+100+k)=p(mod(i+11,127)+1)*data((i-1)*50+k);
        sym((i-1)*200+100+101-k)=p(mod(i+11,127)+1)*conj(data((i-1)*50+k));
        end
    end
elseif(isequal(rate,[0 0 0 1 0])|isequal(rate,[0 0 0 1 1])|isequal(rate,[0 0 1 0 0]))
     for i=1:(length/100)
        for k=1:100
        sym((i-1)*200+k)=data((i-1)*100+k);
        sym((i-1)*200+201-k)=p(mod(i+11,127)+1)*(imag(data((i-1)*100+k))+j*real(data((i-1)*100+k)));
        end
     end
elseif(isequal(rate,[0 0 1 0 1])|isequal(rate,[0 0 1 1 0])|isequal(rate,[0 0 1 1 1]))
     sym=data;
end
out=sym;