matlab实现DSB信号调制解调：相干解调器解调，画出原始信号和解调信号；通过AWGN信号，信噪比20db，画出原始信号和解调信号进行对比

以下是 Matlab 实现 DSB 信号调制解调的代码：

%% 生成原始信号
fs = 100; % 采样率
t = 0:1/fs:1; % 时间范围
f1 = 10; % 原始信号频率
f2 = 40;
s1 = sin(2*pi*f1*t); % 原始信号1
s2 = sin(2*pi*f2*t); % 原始信号2
s = s1 + s2; % 合成原始信号

%% DSB 调制
fc = 20; % 载波频率
c = cos(2*pi*fc*t); % 载波信号
m = s .* c; % DSB 调制信号

%% DSB 相干解调
x = m .* c; % 相干解调器输入信号
[b, a] = butter(6, 2*fc/fs); % 设计低通滤波器
y = filter(b, a, x); % 低通滤波
s_rec = y - mean(y); % 解调信号

%% 绘图
subplot(2,1,1)
plot(t,s,'b','LineWidth',1.5);
title('原始信号');
xlabel('时间 (s)');
ylabel('幅度');

subplot(2,1,2)
plot(t,s_rec,'r','LineWidth',1.5);
title('解调信号');
xlabel('时间 (s)');
ylabel('幅度');

以下是加入 AWGN 信道后的代码：

%% 生成原始信号
fs = 100; % 采样率
t = 0:1/fs:1; % 时间范围
f1 = 10; % 原始信号频率
f2 = 40;
s1 = sin(2*pi*f1*t); % 原始信号1
s2 = sin(2*pi*f2*t); % 原始信号2
s = s1 + s2; % 合成原始信号

%% DSB 调制
fc = 20; % 载波频率
c = cos(2*pi*fc*t); % 载波信号
m = s .* c; % DSB 调制信号

%% 加入 AWGN 信道
SNR = 20; % 信噪比
Pm = mean(m.^2); % 信号功率
Pn = Pm / 10^(SNR/10); % 噪声功率
n = sqrt(Pn) * randn(size(m)); % 加性高斯白噪声
m_noise = m + n; % 加入噪声后的信号

%% DSB 相干解调
x = m_noise .* c; % 相干解调器输入信号
[b, a] = butter(6, 2*fc/fs); % 设计低通滤波器
y = filter(b, a, x); % 低通滤波
s_rec = y - mean(y); % 解调信号

%% 绘图
subplot(2,1,1)
plot(t,m,'b','LineWidth',1.5);
hold on;
plot(t,m_noise,'r','LineWidth',1.5);
title('原始信号和加噪声后的信号');
xlabel('时间 (s)');
ylabel('幅度');
legend('原始信号','加噪声后的信号');

subplot(2,1,2)
plot(t,s,'b','LineWidth',1.5);
hold on;
plot(t,s_rec,'r','LineWidth',1.5);
title('原始信号和解调信号');
xlabel('时间 (s)');
ylabel('幅度');
legend('原始信号','解调信号');

需要注意的是，由于 AWGN 信号是随机的，每次运行代码得到的结果可能会有所不同。