考虑载波和采样频率的2PSK调制 MATLAB仿真

%功能：生成psk调制信号
%创建日期：2016/7/27
%创建人：Alice（764499604@qq.com）
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]);