%功能：生成psk調製訊號
%建立日期：2016/7/27
%建立人：Alice（[email protected]）
clear all; 
close all;
clc;
%%
max=100;
g=zeros(1,max);
g=randint(1,max);     %長度為max的隨機二進位制序列
cp=[];
mod1=[];
f=2*pi;                    
fc=10000;                  %載波頻率
fs=90000;                  %取樣率
Sp=100;                     %每個值100個取樣點

t=0:1/fs*2*pi/99:1/fs*2*pi;
 
%%
for n=1:length(g);
    if g(n)==0; 
        A=zeros(1,Sp);   %每個值100個點          
    else g(n)==1;
        A=ones(1,Sp);          
    end
    cp=[cp A];                   %碼元寬度100  
    c=cos(f*t*fc);                   %載波訊號  
    mod1=[mod1 c];         %與s(t)等長的載波訊號,變為矩陣形式
end
 
cm=[];mod=[];
%%
for n=1:length(g);
    if g(n)==0; 
        B=ones(1,Sp);      %每個值100個點 
        c=cos(f*t*fc);            %載波訊號
    else g(n)==1;
        B=ones(1,Sp); 
        c=cos(f*t*fc+pi);      %載波訊號
    end
    cm=[cm B];              %s(t)碼元寬度100   
    mod=[mod c];          %與s(t)等長的載波訊號
end
tiaoz=cm.*mod;          %e(t)調製
%%

ask=cp.*mod;
subplot(2,1,1);plot(cp,'LineWidth',1.5);grid on;
title('Binary Signal');
axis([0 Sp*length(g) -2.5 2.5]);

subplot(2,1,2);
plot(tiaoz,'LineWidth',1.5);grid on;
title('PSK modulation');
axis([0 Sp*length(g) -2.5 2.5]);