f=1*300;???????					%采樣頻率?
g=（sign（rand（1300）-0.5+eps）+1）/2;		%產生二進制序列
sn=randn（1100*length（g））;			%產生高斯白噪聲
dt=2*pi/299;
t=0:dt:2*pi;
si=[];co=[];????					%si為正交分量，co為同相分量
sit=[];sqt=[];??					%sit為同相分量幅度，sqt為正交分量幅度
sb2=[];?????????					%輸入二進制序列
for?n=1:3:length（g）;					%一次取3個二進制數
if?g（n）==0?&&?g（n+1）==0?&&?g（n+2）==0	
%b1b2b3=000時正交分量和同相分量的幅值
????????it=-0.383*ones（1300）;
????????qt=-0.924*ones（1300）;??
????????b2=[zeros（1100）?zeros（1100）?zeros（1100）];
????elseif?g（n）==0?&&?g（n+1）==0?&&?g（n+2）==1		%b1b2b3=001時
????????it=-0.924*ones（1300）;
????????qt=-0.383*ones（1300）;
????????b2=[zeros（1100）?zeros（1100）?ones（1100）];
????elseif?g（n）==1?&&?g（n+1）==0?&&?g（n+2）==0
????????it=0.383*ones（1300）;
????????qt=-0.924*ones（1300）;?
????????b2=[ones（1100）?zeros（1100）?zeros（1100）?];
????elseif?g（n）==1?&&?g（n+1）==0??&&?g（n+2）==1
????????it=0.924*ones（1300）;
????????qt=-0.383*ones（1300）;??
????????b2=[ones（1100）?zeros（1100）?ones（1100）?];
????elseif?g（n）==0?&&?g（n+1）==1?&&?g（n+2）==0
????????it=-0.383*ones（1300）;
????????qt=0.924*ones（1300）;?
????????b2=[zeros（1100）?ones（1100）?zeros（1100）??];
????elseif?g（n）==0?&&?g（n+1）==1?&&?g（n+2）==1
????????it=-0.924*ones（1300）;
????????qt=0.383*ones（1300）;?
????????b2=[zeros（1100）?ones（1100）?ones（1100）??];
????elseif?g（n）==1?&&?g（n+1）==1?&&?g（n+2）==1
????????it=0.924*ones（1300）;
????????qt=0.383*ones（1300）;??
????????b2=[ones（1100）?ones（1100）?ones（1100）??];
????elseif?g（n）==1?&&?g（n+1）==1?&&?g（n+2）==0
????????it=0.383*ones（1300）;
????????qt=0.924*ones（1300）;
????????b2=[ones（1100）?ones（1100）?zeros（1100）??];
????end
????sb2=[sb2?b2];
????c=cos（f*t）;???s=sin（f*t）;
????sit=[sit?it];?sqt=[sqt?qt];???
????co=[co?c];????si=[si?s];??
end
figure（1）;
subplot（311）;
plot（sit）;%I路
title（‘I‘）;
subplot（312）;
plot（sqt）;%Q路
title（‘Q‘）;
psk=sit.*co+sqt.*si;							%調制后的8psk信號
%%%%%%%%%%%%%%%%%%%%%%%%%%
figure（2）;
subplot（311）;
?stem（sn）%輸出噪聲的時域圖
?title（‘高斯隨機噪聲‘）
zaosheng=fft（sn）;
abs（zaosheng）;
subplot（312）;%畫出噪聲的頻譜圖
abs（zaosheng）
plot（abs（zaosheng））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure（3）;
subplot（211）;plot（sb2‘LineWidth‘1.5）;
grid?on;title（‘二進制序列（信源）‘）;
axis（[0?3000?-1.5?1.5]）;
set（gca‘Xtick‘[3006009001200]）