%16apsk調(diào)制
fc=0.002;%載波頻率即中心頻率
Tc=1/fc;%載波周期
fs=1;%采樣頻率
dt=1/fs;%采樣間隔
fo=0.001;%碼元速率
To=1/fo;%碼元寬度
N=6;%碼元個數(shù)
t=0:dt:(N*To-dt);%時間長度
a(1)=0;
a(2)=5;
a(3)=0;
a(4)=9;
a(5)=6;
a(6)=7;
for k=1:N;
    for m=1:To;
       c(To*(k-1)+m)=a(k);
    end;
end;
figure(1);
subplot(2,1,1); 
plot(c);
title('原學(xué)號數(shù)據(jù)');
b=[1,10,4,7,0,11,5,6,2,9,3,8,12,15,13,14];
for n=1:N;
    s=a(n)+1;
    d(n)=b(s);
end;
for k=1:N;
    for m=1:To;
        if (d(k)<12)
        B(To*(k-1)+m)=3*sin(2*pi*fc*t(To*(k-1)+m)+d(k)*pi/6+pi/12);
        else
        B(To*(k-1)+m)=3/3*2.^0.5*sin(2*pi*fc*t(To*(k-1)+m)+(d(k)-12)*pi/2+pi/4);
       end;
    end;
end;
subplot(2,1,2); 
plot(B(1:N*To));
title('16APSK調(diào)制后數(shù)據(jù)');