function [out1,out2]=de_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];
sym=double(datain);
if(isequal(rate,[0 0 0 0 0])|isequal(rate,[0 0 0 0 1]))
    for i=1:(length/200)
        for k=1:50
        data((i-1)*50+k)=sym((i-1)*200+k);
        data((i-1)*50+k)=conj(sym((i-1)*200+101-k))+data((i-1)*50+k);
        data((i-1)*50+k)=sym((i-1)*200+100+k)*p(mod(i+11,127)+1)+data((i-1)*50+k);
        data((i-1)*50+k)=conj(sym((i-1)*200+100+101-k)*p(mod(i+11,127)+1))+data((i-1)*50+k);
        data((i-1)*50+k)=data((i-1)*50+k)/4;  
        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/200)
        for k=1:100
        data((i-1)*100+k)=sym((i-1)*200+k);
        data((i-1)*100+k)=j*real(p(mod(i+11,127)+1)*sym((i-1)*200+201-k))+...
        imag(p(mod(i+11,127)+1)*sym((i-1)*200+201-k))+data((i-1)*100+k);
        data((i-1)*100+k)=data((i-1)*100+k)/2;
        end
     end
elseif(isequal(rate,[0 0 1 0 1])|isequal(rate,[0 0 1 1 0])|isequal(rate,[0 0 1 1 1]))
     data=sym;
end
out1=data;