clear all;
%function [y,Y] = fsk(a,f1,f2)
%2fsk調制，a為輸入碼流，T為碼元周期，f為采樣頻率，
%A為載波幅度，f1,f2為兩個載波的角頻率,f1對應0，f2對應1
f1=1e7;
f2=1.5e7;
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;