xiangweiguji.m

该程序是同步中的相位估计方法

%本程序对频偏估计算法进行了仿真，给出了不同信噪比
%条件下频偏估计均方差的曲线图
Clear;
C1C;
M=4; % Q P SK 调 制
Es =1; %每符号能量
Eb =1/2; %每比特能量
T =1/16000; %符号宽度
delta f=0.05*(1/T); %载波频偏(以符号速率的百分比表示)
theta=pi/8; %载波初始相位的初始化
H = 65; %仿真中用于频偏和相位估计的符号数
EbNOdB = [2:10]; %信噪比
loop-nu m一10; %统计平均的循环次数
delta-f m se=zeros(l,length(EbNOdB)); %估计方差初始化
delta-f m ean= zeros(l,leng(h(EbNOdB)); %估计均值初始化
%首先对本地载波频偏进行估计
for卜l:length(EbNOdB)
for n = 1 :loopnun
EbN O = 1 " Eb NWB(1y10); %信噪比的真值
Std _de v= sg rt( 1/EbNO); %噪声的方差
eta =2 *S td d e v*randn(1,H);
% e ta 是 均 值 为 0,方差为4*(NO/Eb)的高斯白噪声
for m = - (H -1)/ 2:(H-1)/2
phi (m 城H + 1冲 )=
8*p i*d elt af *T *(m 十( H+ 1Y 2)+4*thetaO+eta(m+(H十1)/2);
0/64 倍 the ta( k)表 达 式 ，己 经去 除调制数据
开环结构突发PSK信号载波同步算法研究
end;
for k=-(H-1)/2:(H-3)/2
bet a(k 十( H+1)/2)= ( 3*H/2)*(1一《(k十1/2)/(H/2))^2Y(H^2-1);
%窗 函 数
del ta_ phi (k+(H+1Y2 )= phi(k十1+(H+1)/2}phi(k+(B十1Y2);
%求 相 邻 码元相差
end;
alphaO estim = 0;
for k二一(H-1羚:(H-3户
alph aO e s tim= a lphaOe stim+ d elta_phi(k+(H+1Y2)*beta(k+(H+1Y2 );
%相 邻 码 元相差估计
end;
delta fe stimate= alphaOeses tim/(8*pi*T);%估计deltaf
delta_f m ean(1)=delta
%先累加delta-f
fr口ean(1)+delta f estimate;
delta_f mse(1)=delta f mse(1)+delta f estimate^2;
%先累加delta f平方
end;
delta fmean(lyloop一 num;%恢复delta-f 的均值
de lta f m set heoretical=
sgrt ((1/(T^2))*3./(pi^2*a*H*(R^2-1)*10.^(EbNOdB/10))/10);
% 频偏估计误差的均方差的理论值
delt af m se(1)=sgrt(deltaf m se(1)lloop_num-deltaf m ean(1)^2);
%恢 复delta-f 标准差
end;
%绘图命令
figure(1);
plot(EbNOdB,delta f mean,，一b o');
hold on;

plot(EbNOdB,delta f*ones(l,length(EbNOdB)),，一r "');
hold o氏
axis([2,10,797.5,802.5]);
xlabel('信噪比(dB)');ylabel('频偏估计均a);
title(.频偏估计均值与输入信号信噪比的关系曲线图，);
set(gca,'XTickMode','manual','XTick',[2:1:10]);
set(gca,'YTickMode','manual','nick',[797.5:0.5:1802.5]);
text(7,801.25,['符号速率:',num2str(1/T/1000),'KH z']);
text(7,800.75,[.实际频偏:',num2str(deltaes0, 'Hz ']);
legend('仿真值，’理论值，);
grid;
figure(2);
plot(EbNOM,deltaf rose,'-b o ');
hold on;
plot(EbNOdB,deltaf m se_theoretical,，一r") ;
hold o氏
axis([2,10,0,10]);
xlabel('信噪比(dB)');ylabeIC频偏估计标准差(Hz)');
title('频偏估计标准差与输入信号信噪比的关系曲线图)
sct(gca,'XTickMode','manual',XT ick',[2:1:10]);
text(7,7.5,['符号速率:',num2str(1/T/1000);K Hzj);
text(7,6.5,['实际频偏:',num2str(delta f),'Hzj);
legend('仿真值.，f理论旬;
grid;