clc;clear?all;
bit_number=200;
P=32;?
D=68;
B=2;
d=P/B;
p=D/B;
a=1/2;
M=4;
symbol_number=bit_number/log2（M）;??
pi=3.14159;??????
id=sqrt（-1）;
T=1/10^（4）;
fs=1/T;
n_t=[-3?3];
v=100;

%-------------------------------------
%?add?noise?and?average?
snr=-8:1:8;??
times=2000;?????%?average?times
sigma=zeros（1length（snr））;
ve1=zeros（timeslength（snr））;

avve1=zeros（1length（snr））;

var1=zeros（1length（snr））;
for?t1=1:length（snr）????????
??sigma（t1）=sqrt（1/（2*10^（snr（t1）/10）））;?
?for?k=1:times
?????
???Data_input=randint（12*D）;?%生成BPSK信號。。
???
???[Data]=qpsk（Data_input）;????%??qpsk調制
???
???Pilot_input=zeros（12*P）;?
???[Pilot]=qpsk（Pilot_input）;?
???
???%ck=[Pilot（1:P/2）?Data?Pilot（P/2+1:P）];
???ck=[Pilot?Data?];

for?i=1:length（ck）
????r（i）=ck（i）*exp（id*（2*pi*v*i/fs））;
end
%?add?noise
y=r+sigma（t1）*（randn（1length（r））+id*randn（1length（r）））;
%?choose?point?from?y1
%?calculate?z
z=zeros（1length（y））;%一般賦值時先給出空間在賦值。。。
for?i=1:length（y）
????z（i）=y（i）*conj（ck（i））;
end
%Ip=[1:16?85:100];%前置-后置
Ip=[1:32];%前置
%?Fitz??mothed
for?m=1:P-1
????sum0=0;
????for?k=m:P
???????sum0=sum0+z（Ip（k））*conj（z（Ip（k）-m+1））;
????end
????RM（m）=sum0/（P-m）;
end
sum1=0;
for?i=1:P-1
????sum1=sum1+angle（RM（i））;
end
f_est=sum1/（pi*（P-1）*P*T）;?
var1（t1）=var1（t1）+（（f_est-v）*T）^2/times;

%??Kay?methed
%?gama=zeros（1（L0-1））;
%?for?i=1:L0-1
%??????gama（i）=3/2*L0/（L0^2-1）*（1-（（2*i-L0）/L0）^2）;
%?end
%?
%?for?i=2:L0?%不懂。。
%?????ve1（（t1-1）*times+k）=ve1（（t1-1）*times+k）+gama（i-1）*angle（z（i）*conj（z（i-1）））;
%??end
%?ve1（（t1-1）*times+k）=ve1（（t1-1）*times+k）/（2*pi*T）;
%?avve1（t1）=avve1（t1）+ve1（（t1-1）*times+k）/times;
%?var1（t1）=var1（t1）+（（ve1（（t1-1）*times+k）-v）*T）^2/times;


?end
end
%?MCRB
????MCRB=zeros（1length（snr））?;
????
????for?k=1:length（snr）
????????EsN0（k）=10^（snr（k）/10）;%snr和ESN之間的轉化。。
????????MCRB（k）=3/（2*pi^2*P^3*EsN0（k））??;?%T^2*
????end
????
figure?（1）

%?semilogy（snrMCRB‘-?k‘）;
%?hold?on
semilogy（snrvar1‘-^?k‘）;
hold?on

xlabel（‘???????Es/N0?dB‘）;
ylabel（‘??頻差歸一化均方差‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2123??2012-11-07?18:22??Fitz\Fitz.m

?????文件????????523??2012-07-12?08:37??Fitz\qpsk.m

?????目錄??????????0??2014-06-14?11:42??Fitz

-----------?---------??----------?-----??----

?????????????????2646????????????????????3