Matlab通信仿真系列——幅度调制之抑制载波双边带调制DSBSC

微信公众号上线，搜索公众号小灰灰的FPGA,关注可获取相关源码，定期更新有关FPGA的项目以及开源项目源码，包括但不限于各类检测芯片驱动、低速接口驱动、高速接口驱动、数据信号处理、图像处理以及AXI总线等

本节目录

一、抑制载波双边带调制DSBSC
二、DSBSC信号的仿真
三、DSBSC信号的解调仿真
四、Matlab源码
1、DSBSC信号的仿真源码
2、DSBSC信号的解调仿真源码

本节内容
一、抑制载波双边带调制DSBSC
若在AM基础上，将载波抑制，同时将A0去掉，可输出抑制载波双边带信号，简称双边带信号DSB,时域和频域表达式：

波形和频谱：

DSBSC信号的带宽是基带信号带宽fH的两倍，即B(DSB)=2fH。
DSBSC信号功率为P(DSB)=1/2Ps,Ps为消息信号功率。
二、DSBSC信号的仿真

matlab代码示例：
消息信号是[-3 3]均匀分布的随机整数，
产生的时间间隔为1/10s,
用DSVBSC方法调制载波cos2πfct。
fc=100,0≤t≤10；求：
①消息信号和已调信号的频谱
②已调信号的功率和消息信号的功率

DSBSC信号的带宽是基带信号带宽fH的两倍，即B(DSB)=2fH。

DSBSC信号功率为P(DSB)=1/2Ps,Ps为消息信号功率。

三、DSBSC信号的解调仿真
DSB信号的包络不再与调制信号变化规律一致，需采用相干解调，即同步检波。
在DSBSC信号中，已调信号由m(t)和cos(2πfct)相乘(混频)给出

混频输出y(t)傅里叶变换为

DSBSC信号的解调模型

混频器中有一个低频分量m(t)/2和±2fc附加的高频分量。
当y(t)通过带宽为W的低通滤波器时，高频分量被滤除，低频分量被解调。

matlab代码示例：
消息信号是[-3 3]均匀分布的随机整数，
产生的时间间隔为1/2s,
用DSBSC方法调制载波cos2πfct。fc=100,0≤t≤5；求：
①用同步检波解调该信号，
低通滤波器的截止频率为100Hz,增益为2，
绘制原始信号和解调信号
②调制信号通过awgn信道，信噪比为20dB，
绘制解调后的信号与原始信号

四、Matlab源码
1、DSBSC信号的仿真源码

clear all;
clc;
%信号抽样时间间隔
ts=0.0025
%时间矢量
t=0:ts:10-ts;
%抽样频率
fs=1/ts;
%fft的频率分辨率
df=fs/length(t);
%生成随机数种子的消息序列
msg=randi([-3 3],100,1);
%扩展成取样信号形式
msg1=msg*ones(1,fs/10);
msg2=reshape(msg1.',1,length(t));
%求消息信号的频谱
Pm=fft(msg2)/fs;
f=-fs/2:df:fs/2-df;
subplot(2,1,1);
plot(f,fftshift(abs(Pm)));
title('消息信号频谱');

%载波频率
fc=100;
%已调信号
Sdsb=msg2.*cos(2*pi*fc*t);
%已调信号的频谱
Pdsb=fft(Sdsb)/fs;
subplot(2,1,2);
plot(f,fftshift(abs(Pdsb)));
title('DSBSC信号频谱');
axis([-200 200 0 2]);

%已调信号功率
Pc=sum(abs(Sdsb).^2)/length(Sdsb);
%消息信号功率
Ps=sum(abs(msg2).^2)/length(msg2);

2、DSBSC信号的解调仿真源码

clear all;
clc;
%信号抽样时间间隔
ts=0.0025
%时间矢量
t=0:ts:5-ts;
%抽样频率
fs=1/ts;
%fft的频率分辨率
df=fs/length(t);
f=-fs/2:df:fs/2-df;
%生成随机数种子的消息序列
msg=randi([-3 3],10,1);
%扩展成取样信号形式
msg1=msg*ones(1,fs/2);
msg2=reshape(msg1.',1,length(t));
subplot(3,1,1);
plot(t,msg2);
title('消息信号');


%载波频率
fc=100;
%已调信号
Sdsb=msg2.*cos(2*pi*fc*t);
%相干解调
y=Sdsb.*cos(2*pi*fc*t);
%解调后的频谱
Y=fft(y)./fs;
%低通滤波器的截止频率
f_stop=100;
n_stop=floor(f_stop/df);
%设计低通滤波器
Hlow=zeros(size(f));
Hlow(1:n_stop)=2;
Hlow(length(f)-n_stop+1:end)=2;
%解调信号通过低通滤波器
DEM=Y.*Hlow;
%最终得到的解调信号
dem=real(ifft(DEM))*fs;
subplot(3,1,2);
plot(t,dem);
title('无噪声的解调信号');

%调制信号经过awgn信道
y1=awgn(Sdsb,20,'measured');
%相干解调
y2=y1.*cos(2*pi*fc*t);
%解调后的频谱
Y2=fft(y2)./fs;
%解调信号通过低通滤波器
DEM2=Y2.*Hlow;
%最终得到的解调信号
dem2=real(ifft(DEM2))*fs;
subplot(3,1,3);
plot(t,dem2);
title('信噪比为20dB的解调信号');