?
T=1;??????????????????????????%?基帶信號寬度也就是頻率?
fc=10/T;??????????????????????%?載波頻率?
ml=2;?????????????????????????%?調制信號類型的一個標志位?
nb=100;???????????????????????%?傳輸的比特數?
delta_T=T/200;????????????????%?采樣間隔?
fs=1/delta_T;?????????????????%?采樣頻率?
SNR=0;????????????????????????%?信噪比?
t=0:delta_T:nb*T-delta_T;?????%?限定t的取值范圍?
N=length（t）;??????????????????%?采樣數???
%-------------------------??調制部分?
%?基帶信號的產生?
data=randn（1nb）>0.5;?????????
datanrz=data.*2-1;?????????????%變成極性碼
data1=zeros（1nb/delta_T）;?????%創建一個1*nb/delta_T的零矩陣
for?q=1:nb?
????data1（（q-1）/delta_T+1:q/delta_T）=datanrz（q）;??%將極性碼變成對應的波形信號
end??
%將基帶信號變換成對應波形信號
?data0=zeros（1nb/delta_T）;
for?q=1:nb
???????data0（（q-1）/delta_T+1:q/delta_T）=data（q）;
end
figure（1）;
plot（data0）title（‘原始信號‘）;
axis（[0?20000?-2?2]）;
%串并轉換，將奇偶位數據分開
idata=datanrz（1:ml:（nb-1））;???%將奇偶位分開，因此間隔為2
qdata=datanrz（2:ml:nb）;
%將I相和Q相信號畫出
figure（2）;
subplot（211）;
stairs（idata）title（‘I路信號‘）;
axis（[0?50?-2?2]）;
subplot（212）;
stairs（qdata）;title（‘Q路信號‘）;
axis（[0?50?-2?2]）;
%QPSK信號的調制
ich=zeros（1nb/delta_T/2）;????%創建一個1*nb/delta_T/2的零矩陣，以便后面存放奇偶位數據
for?i=1:nb/2?
????ich（（i-1）/delta_T+1:i/delta_T）=idata（i）;?
end??
?
for?ii=1:N/2?
????a（ii）=sqrt（2/T）*cos（2*pi*fc*t（ii））;??
end?
idata1=ich.*a;???????%將奇數位數據與余弦函數相乘，得到一路的調制信號?????????
qch=zeros（1nb/2/delta_T）;