%*******************modulation****************
D=[1 0 1 1 0 0 1 0];%產生隨機數
a=[];
N=length(D);%基帶信號符號數目
t=0:0.01:7.99;
y=cos(2*pi*t);
d0=D;d1=2*D-1;%單極性信號變雙極性1 -1 1 1 -1 -1 1 -1
D0=[];
x=ones(1,100);%100個1
for k=1:N,
        D0=[D0,d1(k)*x];
end%開關信號,為便于觀察
z=y.*D0;%bpsk的表達式
%*****************demodulation***************
recos=z.*y;
meg=find(recos>0);
ch(meg)=1;
meg=find(recos<0);
ch(meg)=-1;
demo=(ch+1)/2;
%***************drawing**********************
subplot(2,1,1);plot(z);grid ;title('BPSK modulation');axis([0,800,-2,2]);
subplot(2,1,2);plot(demo);grid ;title('BPSK demodulation');axis([0,800,-2,2]);

xlabel('1 0 1 1 0 0 1 0');