mskk.m

来自「关于 8ask的 调制程序 用来做幅度调制的」· M 代码 · 共 50 行

clear;clc;
a=[1 -1 -1 1 1 1 -1 -1 -1];
N=length(a);
fa=zeros(1,N);
fa=0;
for k=2:N
    if a(k)==a(k-1)
        fa(k)=fa(k-1);
    else
        fa(k)=fa(k-1)+pi;
    end
end

lk=cos(fa);
Qk=-a.*cos(fa);

fs=10;     
tim=N/fs;

ff=1024;
t=1/ff:1/ff:tim;

ck=cos((pi*t*fs)/2);
sk=sin((pi*t*fs)/2);

ft=80;
ck0=cos(2*pi*(ft)*t);
sk0=sin(2*pi*(ft)*t);

%x=(lk).*(ck).*(ck0);
%y=(Qk).*(sk).*(sk0);
%smk=x+y;
%plot(t,a);