模拟幅度调制系统抗干扰性能仿真分析

1.引言

一个客观存在的系统多多少少都会受到噪声的干扰，这些噪声来自方方面面，其中最常见的就是存在于信道的高斯白噪声。系统抗干扰性的好坏决定着系统的可靠性，这里我们对一个通信系统进行仿真，模拟信号从产生到接收方将其还原的整个过程，通过计算讨论和分析模拟幅度调制系统抗干扰性。

2.系统模型和仿真

仿真过程中使用的信源为：

　　　　　　　　　　　　　　　　图 1 m(t)
其中fm=10Hz；
频谱M(f):

　　　　　　　　　　　　　　　　图 2 M(f)
载波：

　　　　　　　　　　　　　　　　图 3 c(t)
其中fc=40Hz;
C(f)频谱：

　　　　　　　　　　　　　　　　图 4 C(f)
仿真过程中使用的滤波器:
带通：

　　　　　　　　　　　　　　图 5 带通滤波器
低通：

　　　　　　　　　　　　　　　　　　　　图 6 低通滤波器

仿真过程中高斯白噪声产生：

n=wgn(1,(T/dt),-15);

1.DSB-SC调制和解调

1.调制过程

　　　　　　　　　　　　　　　　　　图 7 DSB-SC 调制系统

上图就是双边带抑制载波调幅信号的产生过程，调制过程中信号的变化如下图所示：
1.m(t)和M(f)

　　　　　　　　　　　　　　图 8 原始信号和频谱

m(t)与载波相乘调制：

　　　　　　　　　　　　　　图 9 调制信号

可见，时域下s(t)高频信号的边缘即是原始信号的体现。而频域下，原始信号的频谱为冲激函数，其冲激频率点即是fm=10Hz，载波对应的频域特性也是一个冲激函数，冲激点在fc=40Hz处，由于时域相乘对应频域卷积，而信号与冲激函数在频域的卷积对应频谱搬移，则最终经过调制可以看到原始信号的频谱被搬移到了f=+fc和f=-fc处。至此，信号的调制已经完成，s(t)下一步应该进入信道进行传输。

2.解调过程

在上述调制过程的结尾提到信号经过信道进行传输，接收方接收到信号后进行解调后即可获得原始信号，但是接收方收到的信号并非经过调制器输出的s(t)，而是在其基础上叠加了噪声，解调过程的模型如下：

　　　　　　　　　　　　　　　　　　　　　　图 10 解调系统框图

其中，nw（t）即噪声，一般为零均值高斯白噪声，白噪声功率谱密度为均匀分布的值，带通滤波器的作用是将无限长的功率谱密度截断为一个截断信号，经过解调后再通过一个低通滤波器即可获得原始信号的近似。
解调过程信号变化如下：

　　　　　　　　　　　　　　　　　　　　图 11 接收端信号

从上图可以很明显的看出，在时域和频域下噪声作用在信号上的体现。
通过带通滤波器处理后频谱图如下：

　　　　　　　　　　　　　　　　　　　　　图 12 滤波后的频谱 　　　　　　　　　　　　　　　　

可见，m(t)频谱两边的噪声已经被滤除，但是m(t)的频谱也发生了失真，这是由于滤波器不理想导致的，在理想滤波器下滤波后的频谱应是上面信号在通带的截断，目的在于限制噪声功率。
相干解调：

　　　　　　　　　　　　　　图 13 相干解调
通过低通滤波器：

　　　　　　　　　　　　　　图 14 还原信号
至此DSB-SC调制解调完成。

3.信噪比分析

　　　　　　　　　　　　　　　图 15 信噪比分析
输入基带信号平均功率为0.5，输出基带信号平均功率为0.2566，
通过计算可得到 :
输入信噪比 (SNR)i= 1.7519e+03，
输出信噪比 (SNR)o=3.4920e+03。
输入信噪比 /输出信噪比=1.9936。

4.matlab代码

T_start=-0.5;%开始时间
T_stop=0.5;%截止时间
T=T_stop-T_start;%仿真持续时间
T_sample=1/1000;%采样间隔
f_sample=1/T_sample; % 采样速率
N=T/T_sample;% 采样点数
f_res=f_sample/N_sample;%频率分辨率
fm=10;fc=40;
am=sqrt(2);
Fs=1000; 
wc=2*pi*fc;
wm=fm*2*pi;
n=-N/2+1:N/2;
t=n*T_sample;             %时间序列
f=n*Fs/N;
figure(1);
%%%%%%%%
mt=am*cos(2*pi*t*fm);  %信源
Mf=fftshift(abs(fft(mt)));
mtt=modulate(mt,fc,Fs,'amdsb-sc');
MTT=fftshift(abs(fft(mtt)));
yt=demod(mtt,fc,Fs,'amdsb-sc');
Yf=fftshift(abs(fft(yt)));
subplot(211);
plot(t,mt);
title('m(t)');
subplot(212);
plot(f,Mf);axis([-2*fc 2*fc 0 max(Mf)]);   
title('M(f)');
 figure(2);
%载波
sam=mt.*cos(wc*t);
SAM=fftshift(abs(fft(sam)));
subplot(211)
plot(t,sam);
hold on;    %画出AM信号波形
plot(t,mt,'r--');
title('s(t)');
subplot(212);
plot(f,SAM);axis([-2*fc 2*fc 0 max(SAM)]);title('S(f)');
%d叠加噪声
n=wgn(1,(T/dt),-15);
si=n+sam;
SI=fftshift(abs(fft(si)));
figure(3);
subplot(211);
plot(t,si);
title('si(t)');
subplot(212);
plot(f,SI);axis([-2*fc 2*fc 0 max(SI)]);   
title('SI(f)');
%输入信号通过带通滤波
s=filter(band_dsb,si);
S=fftshift(abs(fft(s)));
figure(4);
subplot(2,1,1);
plot(f,SI);axis([-2*fc 2*fc 0 max(SI)]);title('SI(f)');
subplot(2,1,2);
plot(f,S);axis([-