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