paf_preprocess.m

来自「Mtlab toolbox containing useful m files 」· M 代码 · 共 47 行

% paf_preprocess.m
% File to load LPIT signal and plot PAF, PACF, ACF using Levanon's PAF code
% (available on his Web site).
% See Appendix A and B.
% Distributed with LPIT

fd=input('Signal filename (.mat) to be analyzed   ','s');
load(fd);
sw=input('(1) Phase Modulation, (2) Frequency Modulation, (3) Polytime, (4) Frequency hopping');
periods=input('number of periods to include N  ');
if sw==1
ffss=input('sampling frequency used fs (Hz) ');
ffcc=input('carrier frequency fc (Hz) ');
cpp=input('cycles per phase cpp ');
    nncc=input('number of phase codes Nc  ');
    numsamps=ffss*cpp/ffcc*nncc;
    u_amp1=sqrt(I(1:periods*numsamps,1).^2+Q(1:periods*numsamps,1).^2);
    u_phase1=atan2(Q(1:periods*numsamps,1),I(1:periods*numsamps,1))./pi;
    u_phase=u_phase1';
    u_amp=u_amp1';
elseif sw==2
ffss=input('sampling frequency used fs (Hz) ');
mp=input('modulation period (in seconds)  ');
    numsamps=ffss*mp*2;
    u_amp1=sqrt(I(1:periods*numsamps,1).^2+Q(1:periods*numsamps,1).^2);
    u_phase1=atan2(Q(1:periods*numsamps,1),I(1:periods*numsamps,1))./pi;
    u_phase=u_phase1';
    u_amp=u_amp1';
elseif sw==3
ffss=input('sampling frequency used fs (Hz) ');
mp=input('code period ');
    numsamps=ffss*mp;
    u_amp1=sqrt(I(1:periods*numsamps,1).^2+Q(1:periods*numsamps,1).^2);
    u_phase1=atan2(Q(1:periods*numsamps,1),I(1:periods*numsamps,1))./pi;
    u_phase=u_phase1';
    u_amp=u_amp1';
elseif sw==4
ffss=input('sampling frequency used fs (Hz) ');
mp=input('Number of samples in a single FH sequence ');
    numsamps=mp;
    u_amp1=sqrt(I(1:periods*numsamps,1).^2+Q(1:periods*numsamps,1).^2);
    u_phase1=atan2(Q(1:periods*numsamps,1),I(1:periods*numsamps,1))./pi;
    u_phase=u_phase1';
    u_amp=u_amp1';

end