matlab怎么生成星座图,关于16QAM生成星座图的程序问题

%16QAM系统的调制、解调的仿真，信源为单极性不归零码(NRZ)

clear all  %清空工作空间的变量

close all  %关闭所有图形窗口

global dt df t f N T %定义全局变量

k=13;

N=2^k;    %总的取样点数

L=32;                %L为每个码元的取样点数

M=N/L;                %M码元总数

Rb=3;                %码元速率Rb=2Mb/s

Ts=1/Rb;                %码元宽度Ts

dt=Ts/L;        %时域取样间隔

df=1/(N*dt);%频域取样间隔

T=N*dt;     %时间截断长度

Bs=N*df/2;  %带宽

Na=4;      %示波器扫描宽度

alpha=0.5;   %升余弦滚降系数

t=linspace(-T/2,T/2,N); %频域横坐标

f=linspace(-Bs,Bs,N)+eps;%时域横坐标

%升余弦

hr1=sin(pi*t/Ts)./(pi*t/Ts);

hr2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);

hr=hr1.*hr2;

HR=abs(t2f(hr));

GT=sqrt(HR);

GR=GT;

figure(1)

set(1,'Position',[10,250,500,400])

%设定窗口位置及大小

figure(2)

set(2,'Position',[515,250,500,400])

%设定窗口位置及大小

figure(3)

set(3,'Position',[50,10,500,200])

%设定窗口位置及大小

figure(4)

set(4,'Position',[600,10,400,300])

%设定窗口位置及大小

A=1;

fc=1/Ts;

m=A*cos(2*pi*fc*t);     % 载波

n=-A*sin(2*pi*fc*t);

EP=zeros(size(f));     %输入信号功率谱

EPa=zeros(size(f));    %调制后信号功率谱

EPo=zeros(size(f));     %输出信号功率谱

for loop1=1:15    %误码率曲线

Eb_N0(loop1)=2*(loop1-1);  %Eb/N0 in dB

eb_n0(loop1)=10^(Eb_N0(loop1)/10);

Eb=0.