%%%%最大似然ML??最大似然方法聯(lián)合實(shí)現(xiàn)符號定時(shí)同步和載波同步仿真
clear?all;
%**********************?preparation?part?***************************
para=128;???%?Number?of?parallel?channel?to?transmit?（points）?%并行信道傳輸個(gè)數(shù)
fftlen=128;?%?FFT?length?%fft?長度
noc=1024;????%?Number?of?carrier?%載波的個(gè)數(shù)
nd=6;???????%?Number?of?information?OFDM?symbol?for?one?loop%OFDM?symbol?循環(huán)個(gè)數(shù)的說明？
ml=2;???????%?Modulation?level?:?QPSK?%調(diào)制方式：QPSK
sr=250000;??%?Symbol?rate??%符號的比特率
br=sr.*ml;??%?Bit?rate?per?carrier?%每個(gè)載波的比特率?
gilen=128;???%?Length?of?guard?interval?（points）%%保護(hù)間隔的長度
ebn0=3;?????%?Eb/N0??%Eb為單位比特的平均信號的能量；n0為噪聲的單邊功率普密度
SNR=15;?????%?ENR??%信噪比
T=10^（-6）;???%?Sampling?time?:?1us%%采樣時(shí)間1us?抽樣頻率1024kHz
%**************************?transmitter發(fā)射機(jī)?*********************************
%**************************?Data?generation?數(shù)據(jù)生成程序****************************?
seldata=rand（1para*nd*ml）>0.5;??%??rand?:?built?in?function??%產(chǎn)生1x（128*6*2）=1X1536的數(shù)據(jù)
M=length（seldata）/noc;
%******************?Serial?to?parallel?conversion串并轉(zhuǎn)換?***********************
paradata=reshape（seldataparand*ml）;?%??reshape?:?built?in?function
%**************************?QPSK?modulation?%QPSK?調(diào)制*****************************?
[ichqch]=qpskmod（paradataparandml）;
kmod=1/sqrt（2）;?%??sqrt?:?built?in?function
ich1=ich.*kmod;%實(shí)部的數(shù)據(jù)
qch1=qch.*kmod;%虛部的數(shù)據(jù)
%*******************?IFFT?************************
x=ich1+qch1.*i;
y=ifft（x）;??????%??ifft?:?built?in?function
ich2=real（y）;???%??real?:?buil