10用户CDMA matlab实现

clear;
clc;
close all;


%%     基带信号产生 扩频码产生

N=20;                 %码元个数                         
xn1=round（rand（1N））; %单极性码，码元周期为100ms，频率为10Hz 
xn2=round（rand（1N））;
xn3=round（rand（1N））;
xn4=round（rand（1N））;
xn5=round（rand（1N））;
xn6=round（rand（1N））;
xn7=round（rand（1N））;
xn8=round（rand（1N））;
xn9=round（rand（1N））;
xn10=round（rand（1N））;
xn1=1-2*xn1;           %双极性码
xn2=1-2*xn2;
xn3=1-2*xn3;
xn4=1-2*xn4;
xn5=1-2*xn5;
xn6=1-2*xn6;
xn7=1-2*xn7;
xn8=1-2*xn8;
xn9=1-2*xn9;
xn10=1-2*xn10;
Xn1=zeros（1567000）;
Xn2=zeros（1567000）;
Xn3=zeros（1567000）;
Xn4=zeros（1567000）;
Xn5=zeros（1567000）;
Xn6=zeros（1567000）;
Xn7=zeros（1567000）;
Xn8=zeros（1567000）;
Xn9=zeros（1567000）;
Xn10=zeros（1567000）;
Cn1=zeros（1567000）;
Cn2=zeros（1567000）;
Cn3=zeros（1567000）;
Cn4=zeros（1567000）;
Cn5=zeros（1567000）;
Cn6=zeros（1567000）;
Cn7=zeros（1567000）;
Cn8=zeros（1567000）;
Cn9=zeros（1567000）;
Cn10=zeros（1567000）;
for i=1:20
    for j=1:28350
        Xn1（j+28350*（i-1））=xn1（i）;             %采样，采样周期为1.59/450ms，一个基带信号码元采样次28350，频率为283500Hz
        Xn2（j+28350*（i-1））=xn2（i）;
        Xn3（j+28350*（i-1））=xn3（i）;
        Xn4（j+28350*（i-1））=xn4（i）;
        Xn5（j+28350*（i-1））=xn5（i）;
        Xn6（j+28350*（i-1））=xn6（i）;
        Xn7（j+28350*（i-1））=xn7（i）;
        Xn8（j+28350*（i-1））=xn8（i）;
        Xn9（j+28350*（i-1））=xn9（i）;
        Xn10（j+28350*（i-1））=xn10（i）;
    end
end
cn1=m_creat（[1 0 0 0 0 1 1][0 0 0 0 0 1]）;   %63扩频码，双极性码，码元周期为1.59ms频率为630Hz一个扩频码采样450次
cn2=[cn1（6:63）cn1（1:5）];cn3=[cn1（11:63）cn1（1:10）];cn4=[cn1（16:63）cn1（1:15）];cn5=[cn1（21:63）cn1（1:20）];cn6=[cn1（26:63）cn1（1:25）];
cn7=[cn1（31:63）cn1（1:30）];cn8=[cn1（36:63）cn1（1:35）];cn9=[cn1（41:63）cn1（1:40）];cn10=[cn1（46:63）cn1（1:45）];
cn11=[];cn22=[];cn33=[];cn44=[];cn55=[];cn66=[];cn77=[];cn88=[];cn99=[];cn1010=[];
for i=1:N;
    cn11=[cn11cn1];
    cn22=[cn22cn2];
    cn33=[cn33cn3];
    cn44=[cn44cn4];
    cn55=[cn55cn5];
    cn66=[cn66cn6];
    cn77=[cn77cn7];
    cn88=[cn88cn8];
    cn99=[cn99cn9];
    cn1010=[cn1010cn10];
end
for i=1:1260
    for j=1:450
        Cn1（j+450*（i-1））=cn11（i）;
        Cn2（j+450*（i-1））=cn22（i）;
        Cn3（j+450*（i-1））=cn33（i）;
        Cn4（j+450*（i-1））=cn44（i）;
        Cn5（j+450*（i-1））=cn55（i）;
        Cn6（j+450*（i-1））=cn66（i）;
        Cn7（j+450*（i-1））=cn77（i）;
        Cn8（j+450*（i-1））=cn88（i）;
        Cn9（j+450*（i-1））=cn99（i）;
        Cn10（j+450*（i-1））=cn1010（i）;
    end
end
%Xn1
figure （1）;
subplot（211）;
plot（0:2000/567000:2000-2000/567000Xn1）;
axis（[0 2000 -1.5 1.5]）;
xlabel（‘ms‘）
title（‘Xn1基带信号‘）;
f_Xn1=fft（Xn12835000）*2/567000;  %频谱分辨率为0.1Hz
subplot（212）;
plot（0:0.1:283500-0.1abs（f_Xn1））;
axis（[0 150 0 1]）;
xlabel（‘Hz‘）;
title（‘Xn1基带信号频谱‘）;
figure （2）;
subplot（211）;
plot（0:2000/567000:2000-2000/567000Cn1）;
axis（[0 500 -1.5 1.5]）;
xlabel（‘ms‘）;
title（‘Xn1-m序列‘）;
f_Cn1=fft（Cn12835000）*2/567000;
subplot（212）;
plot（0:0.1:283500-0.1abs（f