fsk2_05.m

fsk信号

clear all;
%function [y,Y] = fsk(a,f1,f2)
%2fsk调制，a为输入码流，T为码元周期，f为采样频率，
%A为载波幅度，f1,f2为两个载波的角频率,f1对应0，f2对应1
f1=1e7;
f2=1.25e7;
data=randint(400,1,2);
f=80*f1;%f为采用频率
k=1;
for i=1:length(data)
    if data(i)==1
        for j=1:f/f1
            y(k)=cos(16*pi*f2*j/f);
            k=k+1;
        end
    elseif data(i)==0
        for j=1:f/f1
            y(k)=cos(16*pi*f1*j/f);
            k=k+1;
        end
    end
end
signal=awgn(y,5,'measured',[],'db');  
long=8;