ask_2.m

matlab的调制解调

%% 编一个 2ask 的调制方式的程序；
%% 基带信号是一串随机码，每个基带码元用一个周期的载波发送；
%% 设置采样频率，每个码元的采样点数，载波频率；
function [sig_M]=ask_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))=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进制基带时域波形');grid;axis tight;

%    subplot(322);
%    stem(Ji_2);title('4进制基带时域波形');grid;axis tight;

   subplot(323);
   plot(sig_M);title('2ask时域波形');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时域噪声波形');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;