apsk.m

来自「用16apsk方式调制对学号数据进行调制」· M 代码 · 共 41 行

%16apsk调制
fc=0.002;%载波频率即中心频率
Tc=1/fc;%载波周期
fs=1;%采样频率
dt=1/fs;%采样间隔
fo=0.001;%码元速率
To=1/fo;%码元宽度
N=6;%码元个数
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('原学号数据');
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调制后数据');