function bpskseq=BPSK(fbconnection);
fbconnection=[1 0 0 1 0 1 0 0 0 0 0 1];
n=numel(fbconnection);
x=m_sequence(fbconnection);
N=2^n-1;
for i=1:N
if x(i)==1
    bpskseq(i)=1;
else
    bpskseq(i)=-1;
end
end
bpskseq
BPSKSEQ=fft(bpskseq);
figure(2);
subplot(2,1,1);
stem([0:1:(N-1)],bpskseq);axis([0 100 -2 2]);
xlabel('n');ylabel('BPSK');title('BPSK信號');grid;
subplot(2,1,2);
n=0:length(bpskseq)-1;
plot(2.*n/length(bpskseq),abs(BPSKSEQ));
xlabel('frequency(*pi)');ylabel('abs(BPSKSEQ)');title('BPSK信號頻譜');grid;