%% 編一個 2fsk 的調制方式的程序；
%% 基帶信號是一串隨機碼，每個基帶碼元用一個周期的載波發送；
%% 設置采樣頻率，每個碼元的采樣點數，載波頻率；
function [sig_M]=fsk_2(M)
     fs=1000; N=50; f=20; t=(1:N)'/fs; 
     % sig=sin(2*pi*f*t);
     % 設置發送的原始隨機碼序列
     M=M; Ji=rand(M,1); Ji=ceil(Ji-0.5); 
   %% 設計調制信號；共有2種頻率,每個碼元的發送時間50個點；
   %% sig_M表示調制后的信號;    
         sig_M=zeros(M*N,1);
         for i=1:M
             if Ji(i)==0
                 sig_M((i-1)*N+1:1:(i*N))=sin(2*pi*f*t);
             elseif Ji(i)==1
                 sig_M((i-1)*N+1:1:(i*N))=sin(2*pi*(f+f)*t);
             end;
         end;
         sig_M=sig_M;
         %%%% 開始畫圖；
   figure;
   subplot(321);
   stem(Ji);title('2進制基帶時域波形');grid;axis tight;

%    subplot(322);
%    stem(Ji_2);title('4進制基帶時域波形');grid;axis tight;

   subplot(323);
   plot(sig_M);title('2fsk時域波形');grid;axis tight;

   sig_fft1=abs(fftshift(fft(sig_M)));
   subplot(324);
   plot(fs*((0:M*N-1)-(M*N)/2)/(M*N),(sig_fft1));
   title('2fsk頻譜');grid;axis tight;

   sig_noise=awgn(sig_M,10,'measured');
   subplot(325);
   plot(sig_noise);
   title('2fsk時域噪聲波形');grid;axis tight;

   sig_fft2=abs(fftshift(fft(sig_noise)));
   subplot(326);
   plot(fs*((0:M*N-1)-(M*N)/2)/(M*N),(sig_fft2));
   title('2fsk含噪聲的頻譜');grid;axis tight;