BOC调制的matlab仿真

matlab编写的BOC调制的仿真，包括信号的发生，频谱问题，BPSK的对比，信号的捕获

%*******************************************************算法性能评估********
%********************
%% Clean up the environment first =========================================
clear; close all; clc;

format （‘compact‘）;
format （‘long‘ ‘g‘）;              %数字显示格式设置

%--- Include folders with functions ---------------------------------------
addpath include             % The software receiver functions
addpath geoFunctions        % Position calculation related functions  %此处addpath为引用文件夹功能

settings = initSettings（）;
Pfa=0.001;

CN0			=	[35:55];  
times       =   [1:100];
dir_miss=zeros（size（CN0））;
B_l_miss=zeros（size（CN0））;
ASP_miss=zeros（size（CN0））;
CCRW_miss=zeros（size（CN0））;
AAC_miss=zeros（size（CN0））;
SCPC_miss=zeros（size（CN0））;

acqResults.AAC=zeros（size（CN0））;
acqResults.SCPC=zeros（size（CN0））;
acqResults_dir = zeros（size（CN0）size（times））;
acqResults_B_l = zeros（size（CN0）size（times））;
acqResults_ASP = zeros（size（CN0）size（times））;
acqResults_CCRW = zeros（size（CN0）size（times））;
acqResults_AAC = zeros（size（CN0）size（times））;
acqResults_SCPC = zeros（size（CN0）size（times））;

%% Initialization =========================================================
                      %推导应为10*lg（Ps/N0）
for CN0index=   1:21
    for times=   1:100

    
% Find number of samples per spreading boc code period
samplesPerCode = round（settings.samplingFreq * （settings.BOCcodeLength/ ...
                        settings.BOCcodeFreq））;  %每个伪码周期采样点的数目，round为取整函数
% samplesPerCode = round（settings.samplingFreq * （settings.CAcodeLength/ ...
%                         settings.codeFreq））;
                    
% Create three vectors of data
signalI = zeros（1settings.msToProcess*1e-3*settings.samplingFreq）;
signalQ = zeros（1settings.msToProcess*1e-3*settings.samplingFreq）;
              %产生初始化的同步与正交信号零向量序列，括号右侧为500ms内采样点数目
% signal  = zeros（1settings.msToProcess*settings.samplingFreq*1e-3）;

% Find sampling period
ts = 1 / settings.samplingFreq;  %采样周期

% Find phase points of the carrier wave 
phasePoints = （1 : （settings.msToProcess*1e-3*settings.samplingFreq）） * 2 * pi * ts;
% phasePoints = （0 : （settings.msToProcess*samplesPerCode*settings.BOCn-1）） * 2 * pi * ts;
                                 %500ms内所有采样点对应的相位
% Set the number of visible satellites
% Visible_sv = [1 3 6 15 20 26];   %可见星
Visible_sv = 1;
frqshift = zeros（size（Visible_sv））;% used to save different Doppler frqshifts carrier frqs code phases of different visible satellites
frqBins = zeros（size（Visible_sv））;
codePhase = zeros（size（Visible_sv））;%每个可见星分别对应的多普勒频移，载波频率，码相位偏移

% navdatLength = settings.msToProcess/20; %每比特数据20个码周期，所以导航数据长度为25



% Generate all PRN numbers of visible satellites...
for ii = 1:length（Visible_sv）
    
    PRN = Visible_sv（ii）;
    fprintf（‘%02d ‘ PRN）;   %输出即将赋值的可见星号码
    
%% Correlate signals ======================================================   
%--- Generate carrier wave frequency grid -----------
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-04-28 21:01  researchonBOCsignalcapturing\
     文件       11782  2013-03-13 18:10  researchonBOCsignalcapturing\ASPacquisition.asv
     文件       11834  2013-03-26 08:59  researchonBOCsignalcapturing\ASPacquisition.m
     文件       13091  2013-03-21 14:02  researchonBOCsignalcapturing\B_lacquisition.asv
     文件       13109  2013-04-18 19:11  researchonBOCsignalcapturing\B_lacquisition.m
     文件       14548  2013-05-11 15:46  researchonBOCsignalcapturing\acqdetpro.asv
     文件       14550  2013-05-11 15:51  researchonBOCsignalcapturing\acqdetpro.m
     文件       30029  2013-06-02 13:17  researchonBOCsignalcapturing\acqperformencce.m
     文件       12254  2013-01-22 11:50  researchonBOCsignalcapturing\acquisition.m
     文件       10865  2013-06-07 22:56  researchonBOCsignalcapturing\acquisition2.m
     文件       13601  2013-03-21 08:54  researchonBOCsignalcapturing\acuisitionbocmn.asv
     文件       13603  2013-03-21 08:55  researchonBOCsignalcapturing\acuisitionbocmn.m
     文件     1745749  2013-05-16 11:59  researchonBOCsignalcapturing\aspect.fig
     文件        3370  2006-08-10 01:20  researchonBOCsignalcapturing\calculatePseudoranges.m
     文件        6724  2006-08-14 19:38  researchonBOCsignalcapturing\findPreambles.m
     文件       16625  2013-05-30 08:24  researchonBOCsignalcapturing\genenavsignal.asv
     文件       16848  2013-06-07 22:06  researchonBOCsignalcapturing\genenavsignal.m
     文件         170  2012-04-13 09:04  researchonBOCsignalcapturing\generateBOCCode.m
     目录           0  2015-04-28 20:31  researchonBOCsignalcapturing\geoFunctions\
     文件        1720  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\cart2geo.m
     文件        5055  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\cart2utm.m
     文件         679  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\check_t.m
     文件         945  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\clksin.m
     文件         619  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\clsin.m
     文件        1135  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\e_r_corr.m
     文件        2650  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\findUtmZone.m
     文件        1543  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\geo2cart.m
     文件        4461  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\leastSquarePos.m
     文件        5201  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\satpos.m
     文件        2892  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\togeod.m
     文件        1418  2006-08-22 21:46  researchonBOCsignalcapturing\geoFunctions\topocent.m
............此处省略38个文件信息