ssb.m

SSB单边带调制信号 设定原始信号并对其进行单边带调制

%显示模拟调制的波形及解调方法SSB
    %信源
    close all;
    
    clear all;
    
    dt=0.001;       %采样时间间隔
    
    fm=1;           %信源最高频率
    
    fc=10;          %载波中心频率
    
    T=5;            %信号时长
 
    am=sqrt(2);     %频谱

    sn=4;         %信噪比
    
    t=0:dt:T;
    
%%%%%%%%% 画画调制信号波形 %%%%%%%%%%%%%%
    
    mt=am*cos(2*pi*fm*t);     %信源
    
    subplot(5,2,1);
    
    plot(t,mt);         %画调制信号波形
    
    axis([0,5,-2,2]);
    
    xlabel('t');
    
    ylabel('幅度 V');
    
    title('调制信号');
    
%%%%%%%%% 画画调制信号频谱 %%%%%%%%%%%%%%

    [f,sf]=T2F(t,mt);
    
    subplot(5,2,2);
    
    plot(f,sf);         %画调制信号频谱
    
    axis([-2*fm,2*fm,0,4]);
    
    xlabel('f');
    
    axis([-3,3,0,4]);
    title('调制信号频谱');
    

    
%%%%%%%%% 画载波信号波形 %%%%%%%%%%%%%%
    
    mtt=am*cos(2*pi*fc*t);     %信源
    
    subplot(5,2,3);
    
    plot(t,mtt);                    %画载波信号波形
    
    axis([0,1,-2,2]);
    
    xlabel('t');
    
    title('载波信号');

%%%%%%%%% 画画载波信号频谱 %%%%%%%%%%%%%%

    [f,sff]=T2F(t,mtt);
    
    subplot(5,2,4);
    
    plot(f,sff);         %画调制信号频谱
    
    axis([-2*fm,2*fm,0,4]);
    
    xlabel('f');
    
    axis([-12,12,0,2]);
    title('调制信号频谱');
    
%%%%%%%%% 画出已调信号波形 %%%%%%%%%%%%%%    
    
    s_ssb=0.5*am*cos(2*pi*(fc-fm)*t);
    
    subplot(5,2,5);
    
    plot(t,s_ssb);                  %画出SSB信号波形
    
    title('已调调制信号');
    
    xlabel('t');
    
    axis([0,1,-1,1]);
    
    
%%%%%%%%% 画画已调信号频谱 %%%%%%%%%%%%%%

    [f,sff]=T2F(t,s_ssb);
    
    subplot(5,2,6);
    
    plot(f,sff);         %画SSB信号频谱
    
    axis([-12,12,0,2]);
    
    xlabel('f');
    
    
    title('已调信号频谱');
%%%%%%%%% 画出SSB信号波形 %%%%%%%%%%%%%%    
     
   
   s_ssb=0.5*am*cos(2*pi*(fc-fm)*t); 
   db=2*sn/am^2;      
   noise=sqrt(db)*randn(size(t));
   s_ssb=s_ssb+noise;
   subplot(5,2,7);
   plot(t,s_ssb);                  %画出SSB信号波形
   title('SSB调制信号');
   xlabel('t');
   axis([0,1,-1,1]);

%%%%%%%%% 画画SSB信号频谱 %%%%%%%%%%%%%%

    [f,sfff]=T2F(t,s_ssb);
    
    subplot(5,2,8);
    
    plot(f,sfff);         %画SSB信号频谱
    title('SSB调制信号频谱');
    axis([-12,12,0,2]);
    xlabel('f');
    
 %%%%%%%%% 画出解调波形 %%%%%%%%%%%%%%     
    %SSB  demodulation
    rt=4*s_ssb.*cos(2*pi*fc*t);
    
    rt=rt+mean(rt);
        
    [f,rf]=T2F(t,rt);
    
    [t,rt]=lpf(f,rf,fm);
    
    subplot(5,2,9);
    
    plot(t,rt);
    
    title('相干解调后的波形'); 
     axis([0,5,-2,2]);
    xlabel('t');
%%%%%%%%% 画画解调信号频谱 %%%%%%%%%%%%%%

    [f,sffff]=T2F(t,rt);
    
    subplot(5,2,10);
    
    plot(f,sffff);         %画SSB信号频谱
    title('解调信号频谱');
    xlabel('f');
    axis([-3,3,0,4]);