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