T=100;????????%原始信息序列長度
samp=100;?????%采樣點數
N=samp*T;
m=randint（1T）;?????%生成序列
data1=2*m-1;??????????????????%單/雙極性變換

for?i=1:T
????data1_orig（samp*（i-1）+1:samp*i）=data1（i）;??%信息序列擴展
end

%?信源序列圖及其頻譜圖
figure（1）;
subplot（212）;
pwelch（data1_orig[][][]10）;%data1_orig功率譜密度
title（‘信源序列功率譜密度‘）;
subplot（211）;
x1=0+0.01:0.01:100;
plot（x1data1_orig）;axis（[0?100?-1.2?1.2]）;
title（‘信源序列‘）;

%?將信源分成兩路
data11=zeros（1N/2）;????????????
for?i1=1:T/2??
????data11（samp*（i1-1）+1:samp*i1）=data1（2*i1-1）;
end
data12=zeros（1N/2）;??????????????
for?i2=1:T/2??
????data12（samp*（i2-1）+1:samp*i2）=data1（2*i2）;
end

%??對載波抽樣
for?j1=1:N/2??????????????????????
????a1（j1）=cos（（j1-1）*pi/20）;
end
for?j2=1:N/2??????????????????????
????a2（j2）=sin（（j2-1）*pi/20）;
end

%?調制
data21=data11.*a1;??%用余弦調制?????????????????????????
data22=data12.*a2;??%用正弦調制

data3_modul=data21+data22;?????????????????%經過理想信道的QPSK信號
data3_gauss=awgn（data3_modul10‘measured‘）;??%經過AWGN信道后的QPSK信號


%調制后的信號波形和頻譜圖
figure（2）;??????
subplot（211）;
x2=0+0.01:0.01:50;

plot（x2data3_modul）;???title（‘調制后無衰落QPSK信號波形圖‘）;

subplot（212）;
pwelch（data3_modul[][][]10）;??title（‘調制后無衰落QPSK信號頻譜圖‘）;????????????????????????????????????


%只經過高斯信道后的信號波形和頻譜圖

figure（3）;?

subplot（211）;
x2=0+0.01:0.01:50;

plot（x2data3_gauss）;title（‘經AWGN信道后的信號波形圖（信噪比10dB）‘）;?
???????????????????????????????
subplot（212）;
pwelch（data3_gauss[]