clc
clear?all
close?all
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
L?=?64;
N?=?L/2;
n?=?2;
T?=?1/10^n;
fs?=?1/T;
f_bias?=?1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%添加AWGN%%%%%%%%%%%%%%%%%%%%%%%%%
SNR?=?-5:1:15;
times?=?400;?????%?設定循環次數
sigma?=?zeros（1length（SNR））;
var_Fitz?=?zeros（1length（SNR））;
var_Kay?=?zeros（1length（SNR））;
var_LR?=?zeros（1length（SNR））;
f=[];
for?t?=?1:length（SNR）
????sigma（t）?=?sqrt（1/（2*10^（SNR（t）/10）））;
????for?k?=?1:times
????????pilot_input?=?ones（12*L）;
????????[pilot]?=?qpsk（pilot_input）;
????????pilot_frequency?=?[pilot];
????????for?n?=?1:length（pilot_frequency）
????????????s（n）?=?pilot_frequency（n）*exp（j*2*pi*f_bias*n/fs）;%上變頻
????????end???%大概是因為沒有成形的緣故，處理的信號是離散時間的
????????r?=?s+sigma（t）*（randn（1length（s））+j*randn（1length（s）））;%疊加噪聲后的接收信號
????????%r_conj?=?zeros（1length（r））;%一般賦值時先給出空間在賦值也可有可無
????????Ip?=?1:L;%前置，IP是什么？??????
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Fitz%%%%%%%%%%%%%%%%%%%%%%%%%
????????sum_arg_cor?=?0;??
????????for?m?=?1:N
????????????sum_z?=?0;
????????????for?n?=?m:L
????????????????sum_z?=?sum_z+r（Ip（n））*conj（r（Ip（n）-m+1））;
????????????end;
????????????Correlation（m）?=?sum_z/（L-m）;
????????????sum_arg_cor?=?sum_arg_cor+angle（Correlation（m））;
????????end
????????f_est?=?sum_arg_cor/（pi*N*（N+1）*T）;
????????f_est?=?sum_arg_cor/（pi*N*（N-1）*T）;
????????var_Fitz（t）?=?var_Fitz（t）+（（f_est-f_bias）*T）^2/times;
?????????%——————————————————Kay算法—————————————————
????????f_est_Kay?=?0;
????????for?k?=?1:L-1
????????????R?=?3*L/2/（L^2-1）*（1-（（2*k-L）/L）^2）;
????????????f_est_Kay?=?f_est_Kay+R*angle（r（k+1）*conj（r（k）））/2/pi/T;
????????end
????????var_Kay（t）?=?var_Kay（t）+（（f_est_Kay-f_bias）*T）^2/times;
????%——————————————————L&R算法——————————————————
?????????sum_cor?=?0;??
????????for?m?=?1:N
????????????sum_z?=?0;
????????????for?n?=?m:L
????????????????sum_z?=?sum_z+r（Ip（n））*conj（r（Ip（n）-m+1））;
????????????end;
????????????Correlation（m）?=?sum_z/（L-m）;
????????????sum_cor?=?sum_cor+Correlation（m）;????????
????????end
?????????sum_arg?=angle（sum_cor）;
????????f_estimation?=?sum_arg/（pi*N*T）;
????????var_LR（t）?=?var_LR（t）+（（f_estimation-f_bias）*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*L^3*EsN0（k））;
????end
figure（1）
semilogy（SNRMCRB‘--‘‘LineWidth‘3）;hold?on
semilogy（SNRvar_Fitz‘-square‘‘LineWidth‘1.5）;hold?on
semilogy（SNRvar_Kay‘->r‘‘LineWidth‘1.5）;hold?on
semilogy（SNRvar_LR‘-dg‘‘LineWidth‘1.5）;hold?on
xlabel（‘Es/N0?dB‘）;
ylabel（‘歸一化載波頻率估計方差‘）;%最小均方差準則
title（‘三種時域頻率估計算法性能比較‘）;
legend（‘CRB‘‘Fitz算法‘‘Kay算法‘‘L&R算法‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-01-11?10:17??Fitz&Kay仿真性能比較\
?????文件????????2963??2016-10-25?19:08??Fitz&Kay仿真性能比較\main.m
?????文件?????????523??2012-07-12?08:37??Fitz&Kay仿真性能比較\qpsk.m