close?all;clear?all;?clc;
code7=[1110010];
code13=[11111-1-111-11-11];
tao=0.5e-6;
fc=10e6;
fs=40e6;
t_tao=0:1/fs:tao-1/fs;
N=length（code7）;%N=7
NN=length（code13）;%NN=13
pha=0;
t=0:1/fs:7*tao-1/fs;
tt=0:1/fs:13*tao-1/fs;
N2=length（t）;
NN2=length（tt）;
s=zeros（1length（t））;
????for?i=1:N
????????if?code7（i）==1
???????????pha=1;
????????else
???????????pha=0;
????????end???????
???????s（1（i-1）*length（t_tao）+1:i*length（t_tao））=sin（2*pi*fc*t_tao+pha*pi）;
????end

?ss=zeros（1length（tt））;
????for?i=1:NN
????????if?code13（i）==1
???????????pha=1;
????????else
???????????pha=0;
????????end???????
???????ss（1（i-1）*length（t_tao）+1:i*length（t_tao））=sin（2*pi*fc*t_tao+pha*pi）;
????end???
????
????
????

figure
subplot（211）
plot（ts）xlabel（‘t（單位：秒）‘）title（‘二相碼（7位巴克碼）‘）;

subplot（212）
plot（ttss）xlabel（‘t（單位：秒）‘）title（‘二相碼（13位巴克碼）‘）;


f=（-N2/2:N2/2-1）/N2*fs;
y_fft=fft（s）;
ff=（-NN2/2:NN2/2-1）/NN2*fs;
yy_fft=fft（ss）;
figure

subplot（211）
plot（fabs（y_fft））
xlabel（‘頻率f（單位:Hz）‘）title（‘7位巴克碼的頻譜‘）;

subplot（212）
plot（ffabs（yy_fft））
xlabel（‘頻率f（單位:Hz）‘）title（‘13位巴克碼的頻譜‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1212??2011-11-30?14:38??BARKER\Barker_13.m

?????文件????????592??2010-11-24?21:18??BARKER\Barker_7.m

?????目錄??????????0??2011-12-05?14:47??BARKER

-----------?---------??----------?-----??----

?????????????????1804????????????????????3