bpsk_analysis_function.m

产生BPSK分析信号,可以设置基本S/N比等参数,并对信号画图显示.

%-----------二相相位编码信号频率估计子函数---------------%
%二相相位编码信号频率估计子函数
%f=BPSK_analysis_function(signal,N,deltt)
%输入参数：signal 为二相相位编码信号采样值
%          N 为采样点数
%          deltt 为采样间隔
%          figureon 为是否输出调频信号图，1为输出，0为不输出 
%输出参数：freq_BPSK为信号频率估计值

%参考资料： 夏阳松 刘渝 ，信号脉冲特征分析，南京航空航天大学硕士论文，1999年
%设计者：CJF
%设计时间：2005年1月19日

function freq_BPSK=BPSK_analysis_function(signal,N,deltt)

fs=1/deltt;
signalBPSK=zeros(1,N);
signalBPSK=signal;
signalN=zeros(1,N);
signalN=signalBPSK.^2;                  %仿真信号经平方运算，变成频率为2fc的复正弦波
                                             
freq_BPSK=SFMfunction(signalN,N,deltt)/2;       %调用子函数计算频率为2fc的复正弦波函数频率,除以2为二项编码信号载频估计值

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%采用I、Q相位差(PPP算法）的算法计算频率估计值
signal_angle=angle(signalBPSK)/pi;   %求信号的相位，修正到-1～+1之间
signal_angle_cha=signal_angle(2:length(signal_angle)) -  signal_angle(1:length(signal_angle)-1)   ;     %计算相位差

%将相位差修正到-1～+1之间
signal_angle_cha=signal_angle_cha-(signal_angle_cha>1)*2;
signal_angle_cha=signal_angle_cha+(signal_angle_cha<-1)*2;

f_ppp=mean(signal_angle_cha)*(fs/2) +fs;

%经试验，PPP算法的精度没有综合法高

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if abs(f_ppp-freq_BPSK)>fs/4
    freq_BPSK=freq_BPSK+fs/2;
end