function fskdigital(s,f0,f1)
%本程序?qū)崿F(xiàn)FSK調(diào)制
%s--輸入二進(jìn)制序列，f0,f1-兩個載波信號的頻率
%調(diào)用舉例：（f0和f1必須是整數(shù)）fskdigital（[1 0 1 1 0],1,2)
t=0:2*pi/99:2*pi;                 %初始定義
s=[1 0 1 1 0 1 0 1 0];
f0=1;
f1=2;
cp=[];mod=[];bit=[];                
for n=1:length(s);                %調(diào)制過程
    if s(n)==0;
        cp1=ones(1,100);
        c=sin(f0*t);
        bit1=zeros(1,100);
    else s(n)==1
        cp1=ones(1,100);
        c=sin(f1*t);
        bit1=ones(1,100);
    end
    cp=[cp cp1];
    mod=[mod c];
    bit=[bit bit1];
end 
fsk=cp.*mod;
subplot(2,1,1);                 %分別畫出原信號、已調(diào)信號示意
plot(bit,'LineWidth',1.5);grid on;
ylabel('Binary Signal');
axis([0 100*length(s) -2.5 2.5]);
subplot(2,1,2);
plot(fsk,'LineWidth',1.5);grid on;
ylabel('FSK modulation');
axis([0 100*length(s) -2.5 2.5]);