课程作业记录10：63位PN码序列扩频通信Matlab仿真

该过程中的扩频与解扩均采用63位PN码序列。
给出源码如下：

clear all;
close all;

SNR = -30:1:-20;
Eb_No = SNR-10*log10(4);

T=100;
Fs=12600;
bit=20;

%生成随机码元，设前两个码元为1、0
%并进行BPSK调制
rs=randint(1,bit*T);
rs(1)=1;
rs(2)=0;
rs=rs-0.5;
rs=-2*rs;

%扩频
%生成用于扩频的伪随机序列
d=round(rand(1,Fs*T/6));
pn_seq=[];
carrier=[];
%生成的级数为6
t=[0:2*pi/5:2*pi];
for k=1:Fs*T/6
    if d(1,k)==0
        sig=-ones(1,6);
    else
        sig=ones(1,6);
    end
    c=cos(t);   
    carrier=[carrier c];
    pn_seq=[pn_seq sig];
end
r_sp=pn_seq.*r;

%采样
for i=1:1:Fs*T;
    r(i)=rs(ceil(i/630));
end

%开始流程图
%设定抽样时间
t=1/Fs:1/Fs:T;

for aa=1:1:length(SNR)
    
        %BPSK载波
        gi=r.*sin(2*pi*1000.*t)*sqrt(2);
        
        %加性高斯白噪声    
        N0 = 1/10^(SNR(aa)/10);
        gii=gi+sqrt(2/N0)*(randn(1,length(gi)));
        
        %生成基本干扰信号
        J1=20*randn(1,Fs*T);