QPSK信号生成代码

clc;
clear;

t=0:0.0001:1-0.0001;%时间步长
f=30;%载波
M=500;
m=4;%4相调制
fai=0:2*pi/m:2*pi;%相位分配
Carrier1=sin(2*pi*f*t+fai(1));%4种相位分别代表不同的值
Carrier2=sin(2*pi*f*t+fai(2));
Carrier3=sin(2*pi*f*t+fai(3));
Carrier4=sin(2*pi*f*t+fai(4));
P=length(t)/M;
randNum=rand(1,P);
Signal=zeros(1,length(t));
for i=1:P
    if randNum(i)>=0.5
        randNum(i)=1;
    else
        randNum(i)=0;
    end
end
for j=0:P-1
    for p=1:M
        Signal(j*M+p)=randNum(j+1);
    end
end
chuanToBing1=zeros(1,length(t));
chuanToBing2=zeros(1,length(t));
chuanToBing3=zeros(1,length(t));