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

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

   subplot(323);
   plot(sig_M);title('2ask時(shí)域波形');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('2ask頻譜');grid;axis tight;

   sig_noise=awgn(sig_M,10,'measured');
   subplot(325);
   plot(sig_noise);
   title('2ask時(shí)域噪聲波形');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('2ask含噪聲的頻譜');grid;axis tight;