北斗三号 B1C BOC信号仿真

本程序为Matlab程序，根据北斗三号卫星导航B1C特点，画了BOC（1，1）MBOC（6，1，4/33）以及北斗二号BPSK-R（2）频频分布。及信号自相关峰值图，可以看出，BOC调制的多峰值等特点，以及频频带宽特性。

%BOC（alfabeta）  PSD
%sin modulation .n is even
%create by bmj 2014.06.23
function boc
close all;
clear all;
clc;
alfa = 1;
beta = 1;
beta2 = 2;
f0 = 1.023e6;
pi = 3.1415926535898;
fs = alfa*f0;
n = 2*alfa/beta;
ts = 1/（2*fs）;
tc = 1./（beta.*f0）;
tc2 = 1./（beta2.*f0）;
f = -20./tc:0.001./tc:20./tc;
f2 = -20./tc2:0.001./tc2:20./tc2;
% f = -2./tc:0.1./tc:2./tc;
% f_x_boc = 1./（n.*ts）.*（sin（pi.*f.*ts）.*sin（pi.*f.*n.*ts）./（pi.*f.*cos（pi.*f.*ts）））.^2;
f_x_boc = n.*ts.*sinc（f.*ts）.^2.*sin（pi.*f.*ts）.^2;
f_x = tc2.*sinc（f2.*tc2）.^2;
figure
plot（f210.*log10（f_x）‘g-‘‘LineWidth‘2）;axis（[-8./tc8./tc-110-60]）;
hold on
plot（f10.*log10（f_x_boc）‘b-‘‘LineWidth‘2）;axis（[-8./tc8./tc-110-60]）;
xlabel（‘frequency  Hz‘）;
ylabel（‘PSD dBW/Hz‘）;
legend（‘BPSK-R（2）‘‘BOC（11）‘）;
grid on ;
grid minor;

alfa2 = 6;
beta2 = 1;
f0 = 1.023e6;
% pi = 3.1415926535898;1./（n.*ts）
fs = alfa2*f0;
n = 2*alfa2/beta2;
ts = 1/（2*fs）;
tc = 1./（beta2.*f0）;
f = -20./tc:0.001./tc:20./tc;
% f = -2./tc:0.1./tc:2./tc;
% f_x_boc2 = f0.*（sin（pi.*f.*ts）.*sin（pi.*f.*n.*ts）./（pi.*f.*cos（pi.*f.*ts）））.^2;
% f_x_boc2 = 16./（n.*ts.*pi.*pi.*f.*f）.*sin（pi.*f.*t