function [out]=mapping_psdu(data,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];

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+12,127))*data((i-1)*50+k);
        sym((i-1)*200+100+101-k)=p(mod(i+12,127))*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+12,127))*imag(data((i-1)*100+k))+j*imag(data((i-1)*100+k));
        end
    end