modifiedwiener.m

PHS多天线基带调制与解调MATLAB程序

%function for wiener SAPSP

function [AirData, BitData] = ModifiedWiener (AirParam, BitParam)
global Z;
global alg;
global sys;
global env;
global n;

cisoid_S = exp(-j*2*pi*[1:alg.nRzzSym]*AirParam.FO/sys.rate).' ;

lower = 1 + alg.beginRzzSym*sys.over;
upper = lower + alg.nRzzSym*sys.over - 1;
uwstart = 8 + 1 + 6*sys.over;

if( alg.UWWeights == 1 )
   
    
    dUWZ = Z(uwstart+round(AirParam.TO):sys.over:(uwstart+10*sys.over)+round(AirParam.TO)-1,:) .* ...
             (cisoid_S(1:10) * ones(1,sys.nAnt)) ;
    RzzUW = dUWZ'*dUWZ;
    invRzzUW = inv(RzzUW);
%    invRzzUW = inv(RzzUW+(abs(trace(RzzUW))*10e-2)*eye(4));
    wStart = invRzzUW * (dUWZ'*BitParam.uwref) ;
end


wTrack = wStart ;

    cisoid_S = exp(-j*2*pi*[1:sys.nSym]*AirParam.FO/sys.rate).' * ones(1,sys.nAnt) ;
    dZ = Z(1+round(AirParam.TO):sys.over:1+round(AirParam.TO)+sys.over*(sys.nSym-1),:) ; 
    dZ = dZ .* cisoid_S ;

    zTrack = dZ * wTrack;


    BitData.bits = phs_demod_mex(zTrack.',1,sys.tracking) ;

    [BitData.S_uwe, BitData.S_crce, BitData.S_type, BitData.S_info] = ...
    PayloadExtract_mex(BitData.bits,0,-1,-1,sys.keyS) ;

    [BitData.I_uwe, BitData.I_crce, BitData.I_type, BitData.I_info] = ...
    PayloadExtract_mex(BitData.bits,0,-1,-1,sys.keyI) ;


    switch alg.FreqAlignEstMode

    case 'normal'


        cisoid_S = exp(-j*2*pi*[1:alg.nFreqEstSym*sys.over]*AirParam.FO / ...
                      (sys.over*sys.rate)).' ;
        lower = 1 + alg.beginFreqEstSym*sys.over;
        upper = lower + alg.nFreqEstSym*sys.over - 1;
        dZ = Z(lower:upper,:) .* (cisoid_S * ones(1,sys.nAnt));
  


        z = dZ * wTrack ;
      

        [AirData.TO, AirData.FO] = Freq_Time_Recovery(z, sys.over, sys.rate);
 
    case 'perfect'
        AirData.TO = AirParam.TO;
        AirData.FO = AirParam.FO;
    end


    AirData.w = wTrack;