% 实现2ASK调制以及相干解调，并画出调制解调后的信号波形。 文件名为：twoask.m
% 并且含有T2F、F2T、LPF这三个子函数
clear all
global dt t f df N T % 全局变量
close all
N = 2^20;
dt = 0.001; % ms
df = 1/(N*dt); % kHz
T = N*dt; % 截短时间
Bs = N*df/2; % 系统带宽
t = linspace(-T/2,T/2,N); % 时域横坐标
f = linspace(-Bs,Bs,N)+eps; % 频域横坐标
fm = 1; % 模拟基带信号频率为1kHz
fc = 6; % 正弦载波信号频率为6kHz
mt = cos(2*pi*fm*t); % 模拟基带信号
m(mt>0)=1; % 发送比特
m(mt<0)=0;
c = cos(2*pi*fc*t); % 正弦载波信号
s = m.*c; % DSB-SC AM已调信号
S = T2F(s);
n = awgn(s,20); % AWGN信道噪声
r = s+n; % 接收信号
y = r.*c;
Lpfil = LPF(fc); % 低通滤波器
Y = T2F(y);
yrt = real(F2T(Lpfil.*Y));
yr = yrt-mean(yrt);
yr(yr>0)=1; % 发送比特
yr(yr<0)=0;
subplot(3,1,1)
plot(t,m,'LineWidth',1.2)
plot(t,m,'LineWidth',1.2);
subplot(3,1,1)
plot(t,m,'LineWidth',1.2);
subplot(3,1,1)
plot(t,m,'LineWidth',1.2);
title('输入信号');
axis([-1,+1,-0.2,1.2*max(m)])
xlabel('t (ms)')
ylabel('s(t) (V)')
subplot(3,1,2)
plot(t,s,'LineWidth',1.2)
title('2ASK 调制后波形');
axis([-1,+1,1.2*min(c),1.2*max(c)])
xlabel('t (ms)')
ylabel('s(t) (V)')
subplot(3,1,3)
plot(t,yr,'LineWidth',1.2)
axis([-1,+1,-0.2,1.2*max(yr)])
title('2ASK 解调后波形');
xlabel('t (ms)')
ylabel('s(t) (V)')
