OFDM同步ML算法MATLAB仿真.doc

m1=３；c1＝zeros(m１,2０4８);
fｏr zhen1=1：m1
M1=1992;
ａ1＝raｎd（1，７968)；
ａ1=sign(a1－0。5)；
b１=reｓhape(a１，199２,4)；
　 xxx1＝2＊b１(：,１)＋ｂ1(：，2)＋j＊(2*ｂ1(:，3)+b1（:,4)）;
xx１=iffｔ(xxx1’）＊19９2；
　　 geshu1=length（xx１）；
　 s11＝1；s2１=56；ｓ３1=１;
　　w11＝s1１：s３1:s21;
　 addcp１=ones(1，ｇeshu１+s21);
　fｏr　i1=w11
　 aｄdｃｐ1(1,i1）=xx１(1992-s２１＋i１)；
　 ｅnd
　 　 w21=1:１:gｅｓhu1；
ｆor　k１＝w21
　 addｃp1（1，s２1＋k1)=xx１（k1)；
　　 end
　
　　ｃ1(zｈen1，1:20４８)=addcp１;
　　
　 c１1＝c1＇；
　 adｄcp１1＝reshapｅ（c1１，1，6144）；
end

foｒ　n=1：3
　 for m=１:2048 　 　 　　
　 　 　addcp1１（(n－1）＊2048+m)＝addcp11((ｎ—1)＊204８+m）*ｅxp（i＊2*pi＊0。５＊ｍ/204８);　 　　 ％加入频偏V以后的各样点值
　 end
end


　％信道模型：带多普勒频移的瑞利衰落信道
　 ｆade=［0 0 ０．2 0．3]； ％多径的幅度
　 　 Ｌen=ｌength(aｄｄcp11）；
　 f＿Ｌen=leｎgtｈ(fade）；
　 ｓch=ａddcp11；
　　 　　　fｏｒ m=１：f_Len
　　　　 　 sch(１+m：Leｎ)=sch（1+m：Ｌen）+ｆade(m)*ａddcp11（１:Len-m)；
　 　 　 end
　　　 　
　 　 fｄ＝；％多普勒频移

　 foｒ m=1:lengtｈ（sch） 　　 　　 　 　 　 　 　 　　
　 　　 ａdｄcp１１（ｍ）=scｈ(ｍ）＊ｅxｐ（i＊2＊pi＊fd＊ｍ/2048）; 　　
％加入频偏V以后的各样点值
　 　 　　eｎd
　 　 　
ＳNR=15；
noise=awgｎ（ａddcｐ1１,SＮR）;　
　
aｄｄcp11＝addcｐ1１＋noｉse；

n1＝lｅnｇth(aｄdcp1１）；
ｆoｒ h1=1:n1—1９92
　 as1=0；
　　　 if h１<４098
　 　　 fｏｒ ｌ1=h１：h1+ｓ21—1
　 　 　 as１=as1+abｓ（aｄｄｃｐ１１（l1）＊coｎj（aｄdcp11（ｌ1+M1)）)—