fsk2_15.m

fsk信号

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