通信原理课程设计
设计题目：AM 及SSB 调制与解调及抗噪声性能分析
班级：学生姓名：学生学号：指导老师：
目录
一、引言 (3)
1.1概述 (3)
1.2课程设计的目的 (3)
1.3课程设计的要求 (3)
二、A M调制与解调及抗噪声性能分析 (4)
2.1AM 调制与解调 (4)
2.1.1AM 调制与解调原理 (4)
2.1.2调试过程 (6)
2.2相干解调的抗噪声性能分析 (9)
2.2.1 抗噪声性能分析原理 (9)
2.2.2调试过程 (10)
三、S SB调制与解调及抗噪声性能分析 (12)
3.1 SSB 调制与解调原理 (12)
3.2SSB 调制解调系统抗噪声性能分析 (13)
3.3调试过程 (15)
四、心得体会 (19)
五、参考文献 (19)
一、引言
1.1概述
《通信原理》是通信工程专业的一门极为重要的专业基础课，但内容抽象，基本概念较多，是一门难度较大的课程，通过MATLAB仿真能让我们更清晰地理解它的原理，因此信号的调制与解调在通信系统中具有重要的作用。本课程设计是AM及SSB 调制解调系统的设计与仿真，用于实现AM及
SSB 信号的调制解调过程，并显示仿真结果，根据仿真显示结果分析所设计的系统性能。在课程设计中，幅度调制是用调制信号去控制高频载波的振幅，使其按调制信号的规律变化，其他参数不变。同时也是使高频载波的振幅载有传输信息的调制方式。
1.2课程设计的目的
在此次课程设计中，我需要通过多方搜集资料与分析：
(1)掌握模拟系统AM和SSB调制与解调的原理；
(2)来理解并掌握AM和SSB调制解调的具体过程和它在MATLAB中的实现方法；
(3)掌握应用MATLAB分析系统时域、频域特性的方法，进一步锻炼应用MATLAB进行编
程
仿真的能力。
通过这个课程设计，我将更清晰地了解AM和SSB的调制解调原理，同时加深对MATLAB
这
款《通信原理》辅助教学操作的熟练度。
1.3课程设计的要求
(1)熟悉MATLAB的使用方法，掌握AM信号的调制解调原理，以此为基础用MATLAB编程
实现信号的调制解调；
(2)设计实现AM调制与解调的模拟系统，给出系统的原理框图，对系统的主要参数
进行设计说明；
(3)采用MATLAB语言设计相关程序，实现系统的功能，要求采用一种方式进行仿真，即
直接采用MATLAB语言编程的静态方式。要求采用两种以上调制信号源进行仿真，并记录各个输出点的波形和频谱图；
(4)对系统功能进行综合测试，整理数据，撰写课程设计论文。
二、AM调制与解调及抗噪声性能分析
2.1A M 调制与解调
2.1.1AM 调制与解调原理
幅度调制是由调制信号去控制高频载波的幅度，使正弦载波的幅度随着调制信号而改变的调制方案，属于线性调制。
AM信号的时域表示式：
频谱：
调制器模型如图所示：
A0 cos c t
图1-1 调制器模型
AM的时域波形和频谱如图所示：
图 1-2 调制时、频域波形
AM 信号的频谱由载频分量、上边带、下边带三部分组成。它的带宽是基带信号带宽
的
在波形上，调幅信号的幅度随基带信号的规律而呈正比地变化，在频谱结构上，它的频谱完全是基 带信号频谱在频域内的简单搬移。
所谓相干解调是为了从接受的已调信号中，不失真地恢复原调制信号，要求本地载波和接收信 号的载波保证同频同相。相干载波的一般模型如下：
c t cos c t
将已调信号乘上一个与调制器同频同相的载波，得
2
S AM (t) cosw c t [A 0 m(t )] cos 2
w c t
11
2
[A 0 m(t)] 2[A 0 m(t)]cos2w c t 由上式可知，只要用一个低通滤波器，就可
以将第
1 项与第 2项分离，无失真的恢
复出原始的调
制信号
1
M 0(T) [A 0 M(T)]
2 相干解调的关键是必须产生一个与调制器同频同相位的载波。如
果同频同相位的条件得不到满 足，则会破坏原始信号的恢复。
频
s
p t
LPF
s
d t
2.1.2调试过程
t=-1:0.00001:1; % 定义时长
A1=6; %调制信号振幅
A2=10; %外加直流分量f=3000; %载波频率w0=2*f*pi; %角频率
Uc=cos(w0*t); %载波信号subplot(5,2,1);
plot(t,Uc); %画载波信号title(' 载波信号');
axis([0,0.01,-1,1]); %坐标区间
T1=fft(Uc); %傅里叶变换
subplot(5,2,2);
plot(abs(T1));%画出载波信号频谱
title(' 载波信号频谱'); axis([5800,6200,0,200000]); % 坐标区间mes=A1*cos(0.002*w0*t); %调制信号
subplot(5,2,3);
plot(t,mes);% 画出调制信号title(' 调制信号');
T2=fft(mes); % 傅里叶变换subplot(5,2,4);
plot(abs(T2)); %画出调制信号频谱title(' 调制信号频谱');
axis([198000,202000,0,1000000]); % 坐标区间
Uam1=A2*(1+mes/A2).*cos((w0).*t); %AM 已调信号subplot(5,2,5);
plot(t,Uam1);% 画出已调信号title(' 已调信号');
T3=fft(Uam1); %已调信号傅里叶变换subplot(5,2,6);
plot(abs(T3)); ;%画出已调信号频谱title(' 已调信号频谱');
axis([5950,6050,0,900000]); % 坐标区间
sn1=20; %信噪比db1=A1^2/(2*(10^(sn1/10))); % 计算对应噪声方差
n1=sqrt(db1)*randn(size(t)); % 生成高斯白噪声
Uam=n1+Uam1; %叠加噪声后的已调信号
Dam=Uam.*cos(w0*t); % 对AM 已调信号进行解调subplot(5,2,7);
plot(t,Dam);% 滤波前的AM 解调信号title(' 滤波前的AM 解调信号波形');
T4=fft(Dam); % 求AM信号的频谱subplot(5,2,8);
plot(abs(T4));% 滤波前的AM解调信号频谱title(' 滤波前的AM 解调信号频谱');
axis([187960,188040,0,600000]);
Ft=2000; %采样频率
fpts=[100 120]; % 通带边界频率fp=100Hz 阻带截止频率fs=120Hz mag=[1 0];
dev=[0.01 0.05]; %通带波动1%，阻带波动5%
[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev,Ft);
%kaiserord 估计采用凯塞窗设计的FIR滤波器的参数b21=fir1(n21,wn21,Kaiser(n21+1,beta));
%由fir1 设计滤波器
z21=fftfilt(b21,Dam); %FIR低通滤波
subplot(5,2,9);
plot(t,z21,'r');% 滤波后的AM 解调信号title(' 滤波后的AM 解调信号波形');
axis([0,1,-1,10]);
T5=fft(z21); %求AM 信号的频谱subplot(5,2,10);
plot(abs(T5),'r');% 画出滤波后的AM 解调信号频谱title(' 滤波后的AM 解调信号频谱'); axis([198000,202000,0,500000]);
运行结果:
2.2 相干解调的抗噪声性能分析 2.2.1 抗噪声性能分析原理
AM 线性调制系统的相干解调模型如下图所示。
图 3.5.1 线性调制系统的相干解调模型
图中 S
m (t)
可以是 AM 调幅信号，带通滤波器的带宽等于已调信号带宽 [10]
。下面讨论 AM 调制系
统的抗噪声性能 [11]
。
AM 信号的时域表达式为
S AM (t) [A 0 m(t )] cosw c t
通过分析可得 AM 信号的平均功率为
2
2 (S )
A 0 m 2
(t)
(S
i )
AM
i AM
2 2
又已知输入功率
Ni n0B
, 其中 B 表示已调信号的带宽。 由此可得 AM 信号在解调器的输
入信噪比为
因此解调器的输出噪声功率为
(S i N i ) AM
A 02
m 2
(t) A 02
m 2
(t)
2n 0B AM
4n 0f H
AM 信号经相干解调器的输出信号为
1
m 0(t) 1
2m(t)
因此解调后输出信号功率为
(S 0) AM
2
1 2
m 20 (t) m 2
(t)
4 变成窄带噪声 n
i (t)
，经乘法器相乘后的输出噪声为
在上图中输入噪声通过带通滤波器之后， n p (t) n i (t)cosw c t [n c (t)cosw c t-n (t)sinw t]cosw t
11 n c (t) [n c (t)cos2w c t-经 LPF 后，
n 0(t) 1
2n c (t)
2
N0 n20 (t) 1n c2(t) 1 N i 44
可得AM信号经过解调器后的输出信噪比为
(S0 N0 ) AM m2(t)
n0B
m2(t)
2n0f H
由上面分析的解调器的输入、输出信噪比可得AM信号的信噪比增益为
G AM
2
S0 N0 2m2(t) S i N i A
2 m2(t)
2.2.2 调试过程
clf;
t=0:0.01:2;
fc=50;
A=10;
fa=5;
mt=A*cos(2*pi*fa.*t); xzb=5;
snr=10.^(xzb/10);
db=A^2./(2*snr);
nit=sqrt(db).*randn(size (mt));
psmt=(A+mt).*cos(2*pi*fc .*t); psnt=psmt+nit; xzb1=30;
snr1=10.^(xzb1/10);
%清除窗口中的图形
%定义变量区间%给出相干载
波的频率%定义输入信号幅
度%定义调制信号频率
%输入调制信号表达式
%输入小信躁比(dB)
%由信躁比求方差
%产生小信噪比高斯白躁
声
%输出调制信号表达式%输出叠加小信噪比已调信号
波形%输入大信躁比(dB)
db1=A^2./(2*snr1);
nit1=sqrt(db1).*randn(size(mt) ) ; psnt1=psmt+nit1;
subplot(2,2,1);
plot(t,nit,'g');
title(' 小信噪比高斯白躁声'); xlabel(' t');
ylabel(' nit');
subplot(2,2,2);
plot(t,psnt,'b');
title(' 叠加小信噪比已调信号波形'); %由信躁比求方差%产生大信噪比高斯白躁声%输
出已调信号波形%划分
画图区间%画出输入信
号波形
xlabel('时间');
ylabel(' 输出调制信号');
subplot(2,2,3);
plot(t,nit1,'r');
title(' 大信噪比高斯白躁声
');
xlabel(' t');
ylabel('nit');
subplot(2,2,4);
plot(t,psnt1,'k');
title(' 叠加大信噪比已调信号
波形'); xlabel('时间');
ylabel(' 输出调制信号');
运行结果：
由上图可见，当输入信号一定时，随着噪声的加强，接收端输入信号被干
扰得越严重。而变
非线性元件如滤波器等的存在。非线性失真也会随噪声加大而变大。
%length用于长度匹配
%画出输入信号与噪声叠加
%画出输出信号波形
、 SSB 调制与解调及抗噪声性能分析
3.1 SSB 调制与解调原理
单边带调制信号是将双边带信号中的一个边带滤掉而形成的。根据方法的不同，产生 SSB 信号 的方法有：滤波法和相移法。
由于滤波法在技术上比较难实现 , 所以在此我们将用相移法对 SSB 调制与解调系统进行讨论与 设计。
相移法和 SSB 信号的时域表示
设单频调制信号为
m(t) A m cos m t 载波为
c(t) cos c t
则其双边带信号 DSB 信号的时域表示式为
s DSB (t) A m cos m tcos c t 11
2
A m cos( c m )t 2 A m cos( c m )t
若保留上边带，则有
1
s USB (t ) A m cos( C
m )t 2
若保留下边带，则有
将上两式合并得：
由希尔伯特变换
A m co?s m t A m sin m t
故单边带信号经过希尔伯特变换后得 :
把上式推广到一般情况，则得到
m?(t)的傅里叶变换 M ?
( )为
式中
m?(t)是m(t)的希尔伯特变换
若 M ( )是 m ( t ) 的傅里叶变换，则
11
2 A m cos m cos c t 2 A m sin m sin c t
1
1
s LSB (t) 21
A m cos( C m )t 12
1 A m cos m tcos c t A m sin m tsin c t
M?( ) M ( ) j sgn
上式中的[-jsgn ] 可以看作是希尔伯特滤波器传递函数，即H h( ) M?( )/M ( ) jsgn
移相法SSB调制器方框图
相移法是利用相移网络，对载波和调制信号进行适当的相移，以便在合成过程中将其中的一个边带抵消而获得SSB信号。相移法不需要滤波器具有陡峭的截止特性，不论载频有多高，均可一次实现SSB调制。
SSB信号的解调
SSB 信号的解调不能采用简单的包络检波，因为SSB 信号是抑制载波的已调信号，
它
的包络不能直接反映调制信号的变化，所以仍需采用相干解调。
SSB信号的性能
SSB信号的实现比AM、DSB要复杂，但SSB调制方式在传输信息时，不仅可节省发射功率，而且它所占用的频带宽度比AM、DSB减少了一半。它目前已成为短波通信中一种重要的调制方式。
3.2S SB 调制解调系统抗噪声性能分析
噪声功率
11
N o N i n0B
44
这里， B = fH 为SSB 信号的带通滤波器的带宽。信号功率
SSB 信号
与相干载波相乘后，再经低通滤波可得解调器输出信号
1
m o (t) m(t)
4
因此，输出信号平均功率
2
1 2 S o m o 2(t) 16m 2
(t)
输入信号平均功率为
因m?(t)与m(t)的幅度相同，所以具有相同的平均功率，故上式
12
S i m 2
(t)
4
单边带解调器的输入信噪比为
N i
n 0B 4n 0B
单边带解调器的输出信噪比为
制度增益
因为在 SSB 系统中，信号和噪声有相同表示形式，所以相干解调过程中，信号和噪声中的正交分量 均被抑制掉，故信噪比没有改善。
S i 1
m 2
(t) 4 m 2
(t)
S
o
N o
1 16
m 2(t) m 2
(t) 1
n 0B
4
4n 0B G
SSB
S o /N o S i /N i
3.3调试过程
先建立 3 个M 文件
1.afd_butt
function [b,a] = afd_butt(Wp,Ws,Rp,As);
if Wp <= 0 error( 'Passband edge must be larger than 0' )
end
if Ws <= Wp
error( 'Stopband edge must be larger than Passband edge' ) end
if (Rp <= 0) | (As < 0)
error( 'PB ripple and/or SB attenuation ust be larger than 0'
end
N = ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(Wp/Ws)));
fprintf( '\n*** Butterworth Filter Order = %2.0f \n' ,N)
OmegaC = Wp/((10^(Rp/10)-1)^(1/(2*N))); [b,a]=u_buttap(N,OmegaC);
2.u_buttap
function [b,a] = u_buttap(N,Omegac); [z,p,k] = buttap(N);
p = p*Omegac;
k = k*Omegac^N;
B = real(poly(z));
b0 = k;
b = k*B;
a = real(poly(p));
3.imp_invr
%脉冲响应不变法子程序function [b,a]=imp_invr(c,d,T) [R,p,k]=residue(c,d);
p=exp(p*T);
[b,a]=residuez(R,p,k); b=real(b).*T;
a=real(a);
程序：
>> t0=0.1;fs=12000; %t0 采样区间，fs 采样频率
>> fc=1000;Vm0=2.5;ma=0.25; %fc 载波频率，Vm0 输出载波电压振幅，ma 调幅度>> n=-t0/2:1/fs:t0/2; %定义变量区间
>> N=length(n);
>> A=4; %定义调制信号幅度
>> x1=A*cos(150*pi*n); %调制信号
>> x2=hilbert(x1,N); %对x1 做希尔伯特变换
>> y=(Vm0*x1.*cos(2*pi*fc*n)-Vm0*x2.*sin(2*pi*fc*n))/2; %保留上边带的已调波信号>> xzb=2; %输入小信噪比(dB)
>> snr=10.^(xzb/10);
>> [h,l]=size(x1); %求调制信号的维度
>> fangcha=A*A./(2*snr); %由信噪比求方差
>> nit=sqrt(fangcha).*randn(h,l); %产生高斯白噪声
>> yn=y+nit; %叠加小信噪比噪声的已调波信号
>> xzb=10; %输入小信噪比(dB)
>> snr1=10.^(xzb/10);
>> [h,l]=size(x1); %求调制信号的维度
>> fangcha=A*A./(2*snr1); %由信噪比求方差
>> nit1=sqrt(fangcha).*randn(h,l); %产生高斯白噪声
>> yn1=y+nit1; % 叠加小信噪比噪声的已调波信号
>> figure(1)
>> subplot(2,2,1) %划分画图区间
>> plot(n,x1) %画出调制信号的波形
>> title('调制信号');
>> subplot(2,2,2) %划分画图区间
>> plot(n,y) %画出已调波信号波形
>> title(' 已调波信号');
>> subplot(2,2,3) %划分画图区间
>> plot(n,yn) %画出叠加噪声的已调波信号波形
>> title(' 叠加小信噪比噪声的已调波信号');
>>subplot(2,2,4) %划分画图区间
>> plot(n,yn1) %画出叠加噪声的已调波信号波形
>> title(' 叠加大信噪比噪声的已调波信号');
>> X=fft(x1); %调制信号x1 的傅里叶变换
>> Y=fft(y); %已调信号y 的傅里叶变换
>> Yn=fft(yn); %叠加小信噪比噪声的已调信号yn 的傅里叶变换
>> Yn1=fft(yn1); %叠加大信噪比噪声的已调信号yn 的傅里叶变换
>> w=0:2*pi/(N-1):2*pi; %定义变量区间
>> figure(2)
>> subplot(2,2,1) %划分画图区间
>> plot(w,abs(X)) %画出调制信号频谱波形
>> axis([0,pi/4,0,3000]); %给出横纵坐标的范围
>> title(' 调制信号频谱');
>> subplot(2,2,2) %划分画图区间
>> plot(w,abs(Y)) %画出已调波信号频谱
>> axis([pi/6,pi/4,0,2500]); %给出横纵坐标的范围
>> title(' 已调波信号频谱');
>> subplot(2,2,3) %划分画图区间
>> plot(w,abs(Yn)) %画出叠加小信噪比噪声的已调波信号频谱
>> axis([pi/6,pi/4,0,2500]); %给出横纵坐标的范围
>> title(' 叠加小信噪比噪声的已调波信号频谱');
>> subplot(2,2,4) %划分画图区间
>> plot(w,abs(Yn1)) %画出叠加大信噪比噪声的已调波信号频谱
>> axis([pi/6,pi/4,0,2500]); %给出横纵坐标的范围
>> title(' 叠加大信噪比噪声的已调波信号频谱');
>> y1=y-2*cos(1500*pi*n);
>> y2=Vm0*y1.*cos(2*pi*fc*n); %将已调幅波信号的频谱搬移到原调制信号的频谱处>> wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1; %滤波器参数设计
>> OmegaP=wp/T;OmegaS=ws/T;
>> [cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);
>> [b,a]=imp_invr(cs,ds,T);
>> y=filter(b,a,y2);
>> yn=y+nit;
>> figure(3)
>> subplot(2,1,1) %划分画图区间
>> plot(n,y) %画出解调波波形
>> title('解调波');
>> Y=fft(y); %解调波y 的傅里叶变换
>> subplot(2,1,2) %划分画图区间
>> plot(w,abs(Y)) %画出解调信号频谱
>> axis([0,pi/6,0,2500]); %给出横纵坐标的范围>> title(' 解调信号频谱');
运行结果:
*** Butterworth Filter Order = 6
10
5
-5
-10
调制信号
叠加小信噪比噪声的已调波信号
已调波信号
-0.05 0 0.05
调制信号频谱
3000 2000 1000
0 0.2 0.4 0.6
叠 加 小 信 噪 比 噪 声 的 已 调 波 信 号 频 谱叠 加 大 信 噪 比 噪 声 的 已 调 波 信 号 频 谱
2500 2000 1500 1000 500
10 -0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05
解调信号频谱
2500
2500 0
0.55 0.6 0.65 0.7 0.75
5 0 -5
-10
已调波信号频谱
2000 1500 1000 500 0.55 0.6 0.65 0.7 0.75
解调波
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
2000
1500
1000
500
可以清晰地看出，加大噪声后，信号的波形杂乱无章，起伏远大于加小噪声时的波形。
造成此现象的原因是当信噪比较小时，噪声的功率在解调信号中所占比重较大，所以会造成杂波较多的情况；而信噪比很大时，噪声的功率在解调信号中所占比重就很小了，噪声部分造成的杂乱波形相对就不是很明显，甚至可以忽略。
综上所述，叠加噪声会造成解调信号的失真，信噪比越小，失真程度越大。所以当信噪比低于一定大小时，会给解调信号带来严重的失真，导致接收端无法正确地接收有用信号。所以在解调的实际应用中，应该尽量减少噪声的产生。
四、心得体会
这次的课程设计时间虽短但收获很多。我们用MATLAB进行了AM及SSB调制与解调的研究。不
但又加深了课本的知识，而且也对matlab 的基本知识有了一定掌握。
本次课程设计中实现了通信基本知识与MATLAB的结合，并在实际中设计并仿真AM及SSB调制
与解调的过程。
这次课程设计中我们不得不对AM原理以及SSB原理进行更深一层次的理解，对书中原来学到的只知其果不懂其因的理论，在设计中也有了更深刻的认识。
这次程序需要自己写而我们MATLAB的基础不是很好。这次课程设计虽然很简单，由于没有基础，查找了很多相关的资料，而且我们在编译和调试过程中除了很多次错误，这是过程，学习就是在过程中进行的，经过自己几天的劳动，再加上同学们和老师的帮助，不仅对读程序有了很大提高，而且自己的编译水平也上了一个新台阶，更加熟系了MATLAB的应用，也对其中的函数有了大概的了解
总之这次课程设计收获很大，是一场难得的锻炼机会。
五、参考文献
[1]樊昌信，曹丽娜。通信原理（第六版）。国防工业出版社。