搜索
%C/A码BPSK调制%C/A码时钟为1.023MHz%BOC(n,m),取m=1,可以对比n=2,4调制结果;ca_length=10;%C/A码长度m_boc=1;%BOC调制系数mn_boc=2;%BOC调制系数n,可以取n_boc=4对比结果Nc=32;%一个周期载波信号采样点数Rs=1023000;%C/A码采样时钟N=ca_length*n_boc*Nc;%采样数n=1:N;%%%%%%%%%%%%%%%%%%%%%%%%ca_code_generator产生卫星1的C/A码，ca_code_generator2产生卫星2的C/A码%根据需求，选择不同卫星输出C/A码%产生C/A码，对应卫星1ca_code=ca_code_generator(ca_length)%%产生C/A码，对应卫星2%ca_code=ca_code_generator2(ca_length)%%%%%%%%%%%%%%%%%%%%%%%%将C/A码转变成波形ca_wave=zeros(1,N);fori=0:ca_length-1forj=1:n_boc*Ncca_wave(n_boc*Nc*i+j)=ca_code(i+1);endend%画C/A码波形n_cn=0:ca_length/N:ca_length-ca_length/N;figure(1);subplot(2,1,1);%画C/A码时域波形plot(n_cn,ca_wave);axis([0,ca_length,-0.2,1.2]);xlabel('Ts=1/Rs');title('C/A码');%Frequencyanalysisfs2=Rs;B2=0.5*fs2;df2=fs2/N;f2=-B2:df2:B2-df2;CA_wave=fft(ca_wave);subplot(2,1,2);%画C/A码频域波形plot(f2,fftshift(abs(CA_wave)));xlabel('Hz');ylabel('Aplitude');title('C/A码频谱');%%%%%%%%%%%%%%%%%%%%%%%%BPSK调制%产生正弦载波x_sin=sin(2*pi*n/Nc);n_sin=0:1/Nc:n_boc*ca_length-1/Nc;figure(2);subplot(3,1,1);%画C/A码时域波形plot(n_cn,ca_wave);axis([0,ca_length,-0.1,1.1]);xlabel('Ts=1/Ra');title('C/A码');subplot(3,1,2);%画正弦波载波时域波形plot(n_sin,x_sin);xlabel('Ts/n_boc=1/(n_boc*Rs)');title('正弦载波');%modulationmod_x_sin=zeros(1,N);forn_mod=1:Nmod_x_sin(n_mod)=(2*ca_wave(n_mod)-1)*x_sin(n_mod);endsubplot(3,1,3);%画BPSK调制信号时域波形plot(n_sin,mod_x_sin);xlabel('Ts/n_boc=1/(n_boc*Rs)');title('C/A码BPSK调制结果');%频域分析fs1=n*Rs;B1=0.5*fs1;df1=fs1/N;f1=-B1:df1:B1-df1;X_sin=fft(x_sin);MOD_X_sin=fft(mod_x_sin);figure(3);subplot(3,1,1);%画C/A码频域波形plot(f2,fftshift(abs(CA_wave)));xlabel('Hz');ylabel('Aplitude');title('C/A码波形');subplot(3,1,2);%画正弦波载波信号频域波形plot(f1,fftshift(abs(X_sin)));xlabel('Hz');ylabel('|X_sin|');title('正弦载波频域分析');subplot(3,1,3);%画BPSK调制信号频域波形plot(f1,fftshift(abs(MOD_X_sin)));xlabel('Hz');ylabel('|MOD_X_sin|');title('C/A码BPSK调制频域分析');