%-----雙極性BPSK(相乘法)-----%
clc;
clear;
M=10;%發送符號數
scode=randint(1,M) %發送bit
flag=find(scode==0);
scode(flag)=-1  %轉換成雙極性碼
Tb=1/5;%碼元寬度
f0=1600;%載波頻率
fs=10*f0;%采樣頻率
t=0:1/fs:Tb;%碼元采樣時間
N=length(t);%碼元采樣點數
x=cos(2*pi*f0*t);%產生載波信號
r=[]; %調制信號
for i=1:M   
      dd=(repmat(scode(i),1,N));
      source=dd(:)';
      r=[r,x.*source] ;
end
plot(r)

%-----相位選擇法BPSK-----%
M=10; %發送符號數
scode=randint(1,M) ;%發送bit
Tb=1/5;%碼元寬度
f0=1600;%載波頻率
fs=10*f0;%采樣頻率
phi1=0;%載波相位
phi2=pi;%載波相位
t=1/fs:1/fs:Tb;%碼元采樣時間
N=length(t);%碼元采樣點數
r=[];%調制信號
for i=1:M   
     if scode(i)==1           
        r=[r,cos(2*pi*f0*t+phi1)];   
     else       
        r=[r,cos(2*pi*f0*t+phi2)];  
     end
end
plot(r)