m06_ssbdom.m

MATLAB仿真技术与实例应用教程

t0=0.15;    %信号保持时间
ts=0.001;   %采样时间间隔
fc=250;     %载波频率
fs=1/ts;    %采样频率
df=0.3;     %频率分辨力
t=[0:ts:t0];%时间矢量
m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];  
    %定义信号序列
c=cos(2*pi*fc.*t);                       %载波信号
b=sin(2*pi*fc.*t);
ussb=m.*c-imag(hilbert(m)).*b;          %上边带调制信号
lssb=m.*c+imag(hilbert(m)).*b;          %下边带调制信号
[M,m,dfl]=fft_seq(m,ts,df);              %傅立叶变换
M=M/fs;
[U,ussb,dfl]=fft_seq(ussb,ts,df);          
U=U/fs;
[L,lssb,dfl]=fft_seq(lssb,ts,df);
f=[0:dfl:dfl*(length(m)-1)]-fs/2;        %频率矢量
clf
subplot(2,2,1);plot(t,ussb(1:length(t)));   %未调制信号
title('上边带信号');
subplot(2,2,2);plot(t,lssb(1:length(t)));   %已调制信号
title('下边带信号');
subplot(2,2,3);plot(f,abs(fftshift(U))); %未调制信号频谱
title('上边带调制信号频谱');
subplot(2,2,4);plot(f,abs(fftshift(L))); %已调制信号频谱
title('下边带调制信号频谱');