%16QAM和64QAM
M=[41664256];
for?i=1:4
????k=log2（M（i））;
????n=12000000;???%比特序列長度
????samp=1;?????%過采樣率
????x=randint（n1）;???%生成隨機二進制比特流
????%?????stem（x（1:50）‘filled‘）;??%畫出相應的二進制比特流
????%?????title（‘二進制隨機比特流‘）;
????%?????xlabel（‘比特序列‘）;ylabel（‘信號幅度‘）;
????x4=reshape（xklength（x）/k）;?????%將原始的二進制比特序列每四個一組分組，并排列成k行length（x）/k列的矩陣
????xsym=bi2de（x4.‘‘left-msb‘）;?????%將矩陣轉化為相應的16進制信號序列
????%?figure;
????%stem（xsym（1:50））;????????????????%畫出相應的16進制信號序列
????%title（‘64進制隨機信號‘）;
????%xlabel（‘信號序列‘）;ylabel（‘信號幅度‘）;
????y=modulate（modem.qammod（M（i））xsym）;?%用16QAM調制器對信號進行
????%?????scatterplot（y）;???????????%畫出16QAM信號的星座圖
????%?????text（real（y）imag（y）dec2bin（xsym））;
????%?????axis（[-5?5?-5?5]）;
????EbN0=-5:1:20;
????for?n=1:length（EbN0）;
????????snr（n）=EbN0（n）+10*log10（k）-10*log10（samp）;?%信噪比
????????yn=awgn（ysnr（n）‘me