mimo_freqoffset.m

mimoofdm系统自适应功率分配算法

close all
tic
ee = 0.1;
mu = 0.006;
PRINTL=500;
figure (10)
colorsfoff = ['r' 'g' 'b'];
for ddd=1:3  %mu =0.01,0.005,0.0025
    mu = mu*0.5;
    foff = zeros(1,Mt );   %#
    Es = 1.0;
    akk = zeros(1,Mt );    %#
    delta = zeros(1,Mt );  %#
    mymse = zeros(1,Mt );  %#
    for ii = 1:Mt          %#
        if mod(ii,PRINTL)==0
            fprintf(2,'ii = %10d\n',ii);
        end
        Rx_spectrumCol = transpose(rec_symbol(ii,:)) ;
        if ii ~= 1
            for kk = 1:N
                Rx_spectrumCol(kk,1) = Rx_spectrumCol(kk,1)* exp(j*(-2.0)*pi*(kk-1)*foff(ii-1)/N);%%#
            end
        end
        cz = Rx_spectrumCol * conj(Rx_spectrumCol');
        for jj = 1:N
            %        cz(jj,jj)
            akk(ii) = akk(ii) + cz(jj,jj);
        end
        akk(ii) = akk(ii)/N;
        %    atan(1.0)
        %    atan(-1.0)
        delta(ii)=(1.0/pi)*atan( imag(N*akk(ii) - Es) / real(N*akk(ii) - Es) );
        if ii ~= 1
            foff(ii) = foff(ii-1)+mu*delta(ii-1);
        end
        mymse(ii) = (foff(ii)-ee)^2;
    end
    subplot(2,1,1)
    plot(1:Mt , foff,colorsfoff(ddd))%#
    %semilogy(1:symbols_per_carrier , foff,colorsfoff(ddd))%#
    grid
    ylabel('Frequency Offset')
    xlabel('Time index')
    hold on
    subplot(2,1,2)
    semilogy(1:Mt , mymse,colorsfoff(ddd))%#
    %plot(1:symbols_per_carrier , mymse,colorsfoff(ddd))%#
    grid
    ylabel('MSE')
    xlabel('Time index')
    hold on
    %title('MSE of the estimaotr')
end
hold off;
toc