clc;
close?all;
clear?all;
%?20181120
%?----------------
%?by?呂冬治
%?-----------------
%?此程序用于現(xiàn)代通信原理中的
%?模擬線性調(diào)制部分雙邊帶調(diào)幅中的
%?常規(guī)雙邊帶調(diào)幅（AM）的仿真實(shí)現(xiàn)
%?仿真給出了調(diào)制與解調(diào)過程
%?并分析了在不同的信噪比的條件下
%?利用demod（）函數(shù)實(shí)現(xiàn)解調(diào)的解調(diào)性能分析
%%?信號(hào)參數(shù)
fm?=?1e3;??????%?調(diào)制信號(hào)頻率
fc?=?10e3;?????%?載波頻率
fs?=?200e3;????%?采樣頻率
Am?=?1;????????%?調(diào)制信號(hào)幅度
A?=?2;?????????%?直流分量
N?=?1e3;???????%?采樣點(diǎn)數(shù)
K?=?N-1;
n?=?0:K;
f?=?n*fs/N;????%?頻率
t?=?0:1/fs:K/fs;?????????%?信號(hào)時(shí)長(zhǎng)
wm?=?2*pi*fm;??%?調(diào)制信號(hào)角頻率
wc?=?2*pi*fc;??%?載波信號(hào)角頻率
thm?=?0;???????%?起始相位
thc?=?0;
%%?信號(hào)調(diào)制
yt?=?Am*cos（wm*t+thm）;??%?調(diào)制信號(hào)
fft_yt?=?abs（fft（yt））;
figure（1）
subplot（211）;
plot（tyt）;
xlabel（‘時(shí)間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘頻率為100Hz的調(diào)制信號(hào)LV‘）;
subplot（212）;
plot（f（1:N/2）fft_yt（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘調(diào)制信號(hào)頻譜‘）;
y0?=?A+yt;???????????????%?帶有直流分量的調(diào)制信號(hào)
yc?=?cos（2*pi*fc*n/fs）;??%?載波信號(hào)
fft_yc?=?abs（fft（yc））;
figure（2）
subplot（211）;
plot（tyc）;
xlabel（‘時(shí)間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘頻率為10kHz的載波信號(hào)LV‘）;
subplot（212）;
plot（f（1:N/2）fft_yc（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘載波信號(hào)頻譜‘）;
sig_AM?=?y0.*yc;????????????????%?已調(diào)制信號(hào)
fft_sig_AM?=?abs（fft（sig_AM））;
figure（3）
subplot（211）;
plot（tsig_AM）;
xlabel（‘時(shí)間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘已調(diào)制信號(hào)LV‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_AM（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘已調(diào)制信號(hào)頻譜‘）;
%%?解調(diào)
sig?=?demod（sig_AMfcfs‘a(chǎn)m‘）;??%?對(duì)AM調(diào)制信號(hào)進(jìn)行解調(diào)?
fft_sig?=?abs（fft（sig））;?????????%?對(duì)AM解調(diào)信號(hào)進(jìn)行傅里葉變換??
figure（4）
subplot（211）;
plot（tsig）;
xlabel（‘時(shí)間/s‘）;
ylabel（‘幅值‘）;
axis（[0?5e-3?0?2]）;
grid?on;?
title（‘AM解調(diào)信號(hào)波形LV‘）;
subplot（212）;
plot（f（1:N/2）fft_sig（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
axis（[0?1e5?0?500]）;
grid;
title（‘AM解調(diào)信號(hào)頻譜‘）;
%%?加噪聲
sig_noise_AM?=?awgn（sig_AM3）;?????%?添加信噪比為6dB的高斯白噪聲
fft_sig_noise_AM?=?abs（fft（sig_noise_AM））;?????????%?對(duì)AM解調(diào)信號(hào)進(jìn)行傅里葉變換??
sig_noise_AM2?=?awgn（sig_AM18）;?????%?添加信噪比為15dB的高斯白噪聲
fft_sig_noise_AM2?=?abs（fft（sig_noise_AM））;?????????%?對(duì)AM解調(diào)信號(hào)進(jìn)行傅里葉變換??
figure（5）
subplot（211）;
plot（tsig_noise_AM）;
xlabel（‘時(shí)間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘添加信噪比為3dB噪聲后的AM信號(hào)LV‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_noise_AM（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid;
title（‘添加信噪比為3dB噪聲后的AM信號(hào)頻譜‘）;
figure（7）
subplot（211）;
plot（tsig_noise_AM2）;
xlabel（‘時(shí)間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘添加信噪比為18dB噪聲后的AM信號(hào)LV‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_noise_AM2（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid;
title（‘添加信噪比為18dB噪聲后的AM信號(hào)頻譜‘）;
%%?解調(diào)
sig_noise?=?demod（sig_noise_AMfcfs‘a(chǎn)m‘）;??%?對(duì)AM調(diào)制信號(hào)進(jìn)行解調(diào)?
fft_sig_noise?=?abs（fft（sig_noise））;?????????%?對(duì)AM解調(diào)信號(hào)進(jìn)行傅里葉變換??
sig_noise2?=?demod（sig_noise_AM2fcfs‘a(chǎn)m‘）;??%?對(duì)AM調(diào)制信號(hào)進(jìn)行解調(diào)?
fft_sig_noise2?=?abs（fft（sig_noise2））;?????????%?對(duì)AM解調(diào)信號(hào)進(jìn)行傅里葉變換??
figure（6）
subplot（211）;
plot（tsig_noise）;
xlabel（‘時(shí)間/s‘）;
ylabel（‘幅值‘）;
axis（[0?5e-3?0?2]）;
grid?on;?
title（‘添加3dB噪聲后AM解調(diào)信號(hào)LV‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_noise（1:N

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????3883??2018-11-20?20:03??AM_LV4.m

-----------?---------??----------?-----??----

?????????????????3883????????????????????1