psk_telek.m

Hierarchical－QAM的程序

clear;
theta=input('Enter Theta=');
d=2*(pi/360)*theta;
m=length(d);
n=2^m; 
p=zeros(1,n);
for i=0:n-1
   b=dec2bin(i,m);
   for j=1:m
      p(i+1)=p(i+1)+((b(j)-'0')*2-1)*d(j);
   end
   g(i+1,:)=char(graycode(i,m)+'0');
end
polar(p,ones(1,n),'o');
for i=1:n
   text(cos(p(i)),sin(p(i)),g(i,:));
end