clear?all;
clc;
close?all;

M=16;
k=log2（M）;
n=100000;??????????????????????????%比特序列長度
samp=1;????????????????????????????%過采樣率
?x=randi（[01][1n]）;????????????????????%生成隨機(jī)二進(jìn)制比特流
?stem（x（1:50）‘filled‘）;????????????%畫出相應(yīng)的二進(jìn)制比特流信號(hào)
?title（‘二進(jìn)制隨機(jī)比特流‘）;
?xlabel（‘比特序列‘）;ylabel（‘信號(hào)幅度‘）;

x4=reshape（xklength（x）/k）;???????%將原始的二進(jìn)制比特序列每四個(gè)一組分組，并排列成k行l(wèi)ength（x）/k列的矩陣
?xsym=bi2de（x4.‘‘left-msb‘）;???????%將矩陣轉(zhuǎn)化為相應(yīng)的16進(jìn)制信號(hào)序列
figure;
?stem（xsym（1:50））;??????????????????%畫出相應(yīng)的16進(jìn)制信號(hào)序列
xlabel（‘信號(hào)序列‘）;ylabel（‘信號(hào)幅度‘）;
y=modulate（modem.pskmod（M）xsym）;??%用16PSK調(diào)制器對(duì)信號(hào)進(jìn)行調(diào)制
scatterplot（y）;????????????????????%畫出16PSK信號(hào)的星座圖
text