www.gusucode.com > MC-CDMA系统的仿真matlab源码程序 > mc-cdma/rayleigh(移动信道)/dopplercomp.m

    %%*********************************************************************
%%%%     M.1225信道下，ML估计，不同多普勒频移的比较
%%%********************************************************************
clear

fc=3.5e9;                % operation frequency
%V=75;                    % moving speed, V km/h
%fdmax=V*fc/3e8/3.6;      % maximum frequency shift: fmax=V*fc/C

numusers=4;
wordsize=2;   % represent the mode of modulation, wordsize=2:QPSK, wordsize=4:16QAM, wordsize=6:64QAM;  
Numsymb=100;  % the number of symbols
[ber1,berreal1,mse1]=ML_20(numusers,wordsize,Numsymb,3*fc/3e8/3.6,sqrt(2));  % ML estimation using 16 pilots
[ber3,berreal3,mse3]=ML_20(numusers,wordsize,Numsymb,75*fc/3e8/3.6,sqrt(2));  % ML estimation using 20 pilots
[ber2,berreal2,mse2]=ML_20(numusers,wordsize,Numsymb,125*fc/3e8/3.6,sqrt(2));  % ML estimation using 24 pilots
figure
subplot(1,2,1)
k=0:5:30;
semilogy(k,ber1,'-rx')
hold on
semilogy(k,ber2,'-r*')
semilogy(k,ber3,'-ro')
%semilogy(k,berreal1,'-bx')
%semilogy(k,berreal2,'-b*')
%semilogy(k,berreal3,'-bo')
grid on
axis([0 30 10^(-4) 1])
xlabel('SNR/dB'),ylabel('BER')
legend('V=3','V=75','V=125')
set(gcf,'color',[1,1,1])

subplot(1,2,2)
semilogy(k,mse1,'-rx')
hold on
semilogy(k,mse2,'-r*')
semilogy(k,mse3,'-ro')
grid on
axis([0 30 10^(-4) 1])
xlabel('SNR/dB'),ylabel('MSE')
legend('V=3','V=75','V=125')