(DOC)-matlab与通信仿真实验指导书(下)

matlab与通信仿真实验指导书(下)

《Matlab与通信仿真》实验指导书(下)

刘毓杨辉徐健和煦黄庆东吉利萍编著

上课时间：学年第学期

系部：

班级：姓名：班内序号：指导教师：

实验课程成绩：

通信与信息工程学院

2009-12

目录

实验一MATLAB基础实验.........................................................1

实验二实验三实验四实验五实验六实验七实验八

实验一成绩

绘图和确知信号分析实验...............................................8

实验二成绩

随机信号与数字基带实验.............................................15

实验三成绩

模拟调制实验.................................................................24

实验四成绩

模拟信号数字传输实验（一）.....................................32

实验五成绩

模拟信号数字传输实验（二）.....................................39

实验六成绩

数字频带传输系统实验.................................................45

实验七成绩

通信系统仿真综合实验.................................................57

实验八成绩

Matlab与通信仿真

实验一MATLAB基础实验

一、实验目的

了解MATLAB程序设计语言的基本特点，熟悉MATLAB软件运行环境掌握创建、保存、打开m文件及函数的方法

掌握变量等有关概念，具备初步的将一般数学问题转化为对应的计算机模型并进行处理的能力

二、实验内容及步骤

1.在CommandWindow里面计算

①(358)510；

②sin(3)1

③A4

73

④A7.5

5.4

12i

⑤Z

56i

258

376，B4

911.26.63.4

852

9

6，计算：CAB,DAB,A\C,C/B；3

4

1

3.1，求A,A,A；6.1

34i

，输入复数矩阵；78i

12

2.建立.m文件，用for循环语句生成50×50的矩阵A：

50

2351

5051

,将A矩阵进行水99

平和垂直翻转得到矩阵B和C。将A矩阵的前10行，10列变成0并赋值给D。

3.建立.m文件，随机产生一个50×50的矩阵，元素值为从0到255，要求用0和255对该矩阵进行标记，元素值大于128的标记为255，元素值小于128的标记为0。4.产生一个均值为2.4方差为0.2大小为3×4的随机矩阵。

5.（选做）编写函数使用0.618搜索法（近似黄金分割法）求给定函数的极值：搜索法求解

min(t)的基本过程：给出[a,b]，使得t在[a,b]中。[a,b]称为搜索区间。迭代缩短[a,b]的长

t0

度。当[a,b]的长度小于某个预设的值，或者导数的绝对值小于某个预设的正数，则迭代终止。

1

Matlab与通信仿真

以函数min(t)t2t1t03，作为处理对象，其中搜索区间定为[0,3],精度定为0.5。

三、预习报告：（原理，流程图，编程思想，基本程序）

2

Matlab与通信仿真

3

Matlab与通信仿真

4

Matlab与通信仿真

四、实验报告：（调试好的程序，实验结果及分析）5

Matlab与通信仿真

6

Matlab与通信仿真

7

Matlab与通信仿真

实验二绘图和确知信号分析实验

一、实验目的

掌握二维平面图形的绘制方法，能够使用这些方法进行常用的数据可视化处理理解周期信号的傅里叶级数展开的物理意义掌握信号的傅里叶变换及其反变换

二、实验原理

1.周期信号的傅里叶级数

若一周期信号ftftkT，其中k为整数，T成为信号的周期。若周期信号在一个周期内可积，则可通过傅里叶级数对该信号进行展开。其傅里叶展开式如下：

ft

n

Fne

j2nfst

，Fn

1T

T/2T/2

fte

j2nfst

dt

其中，T为信号最小周期；fs1/T为信号的基波；Fn为傅里叶展开系数，其物理意义为频率分量nfs的幅度和相位。2.信号的傅里叶变换及其反变换

对于非周期信号st，满足绝对可积的条件下，可利用傅里叶变换对其进行频域分析。

Sf

ste

j2ft

dt，st

Sfe

j2ft

df

其中，Sf称为信号st傅里叶变换，表示了该信号的频谱特性。

三、实验内容

1．假设N=12。对于M=4，5，7，10，在0≤n≤2N-1区间上画出xM[n]sin(添上适当标注。用plot和stem分别绘制该信号，并比较。

2MnN

)，并

1,0tT/2

2．设周期信号一个周期[0，T]的波形为st，其中T=1。求该信号

0，T/2tT

傅里叶级数展开式，并用MATLAB画出傅里叶级数展开后的波形，并通过展开式项数的变化考察其对st的逼近程度，考察其物理意义。

1,0t1/23．设非周期信号st，求该信号的傅里叶变换，MATLAB画出傅里

0，1/2t1

叶变换后的频谱，并对频谱进行反变换，画出st的波形。

8

Matlab与通信仿真

9

Matlab与通信仿真

10

Matlab与通信仿真

11

Matlab与通信仿真

五、实验报告：（调试好的程序，实验结果及分析）

12

Matlab与通信仿真

13

Matlab与通信仿真

14

Matlab与通信仿真

实验三随机信号与数字基带实验

一、实验目的

掌握库函数产生随机数方法

基带信号波形生成和其功率谱密度

理解采用蒙特卡罗算法仿真的思想

二、实验原理

（一）.库函数产生随机数

①均匀分布的随机数

利用MATLAB库函数rand产生。rand函数产生（0，1）内均匀分布的随机数，使用方法如下：

1）x=rand(m)；产生一个m×m的矩阵，所含元素取值均为在（0，1）内均匀分布的随机数。

2）x=rand(m，n)；产生一个m×n的矩阵，所含元素取值均为在（0，1）内均匀分布的随机数。

3）x=rand；产生一个随机数。

4）调用rand(1,N)给出均值为0.5，功率为1/12=0.083的白噪声。

5）调用rand(1,N)给出均值为0.5，功率为n0.083的白噪声，现在要均值为0，功率为0.01。

②高斯分布的随机数

randn函数产生均值为0，方差为1的高斯分布的随机数，使用方法如下：

1）x=randn(m)；产生一个m×m的矩阵，所含元素都是均值为0，方差为1的高斯分布的随机数。

2）x=randn(m，n)；产生一个m×n的矩阵，所含元素都是均值为0，方差为1的高斯分布的随机数。

3）x=randn；产生一个均值为0，方差为1的高斯分布的随机数。

4）调用randn，生成均值为0，方差为1，服从高斯分布的白噪声信号u(n)。2

（二）.基带信号波形生成和信号的功率谱密度

1).要画出完整的基带信号波形，每一个码元要采n个样值。2）信号f(t)的功率谱密度为：P()limFT()

T2T。

（三）.蒙特卡罗算法

蒙特卡罗估计是指通过随机实验估计系统参数值的过程。蒙特卡罗算法的基本思想：由概率论可知，随机实验中实验的结果是无法预测的，只能用统计的方法来描述。故需进行大量的随机实验，如果实验次数为N，以NA表示事件A发生的次数。若将A发生的概率近似

15

Matlab与通信仿真

为相对频率，定义为NAN。这样，在相对频率的意义下，事件A发生的概率可以通过重复无限多次随机实验来求得，即：PA

lim

N

NAN

在二进制数字通信系统中，若N是发送端发送的总码元数，NA是差错发生的次数，则

总误码率可通过蒙特卡罗算法计算。

三、实验内容

1．产生随机数。

利用rand函数产生5×4的（0，2）内均匀分布的随机数；

利用randn函数产生5×4的均值为0，方差为2的高斯分布的随机数。

2．利用随机数产生单极性基带信号，每一个码元要采8个样值，并画出其波形和功率谱密。3．利用蒙特卡罗算法仿真二进制双极性基带通信系统的误码率。

假定通信系统满足以下条件：

①信源输出的数据符号是相互独立和等概的双极性基带信号②发送端没有发送滤波器，接收端没有接收滤波器③信道是加性高斯白噪声信道

数字基带信号传输系统模型如图1所示：

图1数字基带信号传输系统模型

Vd0，误码率：

当P1

P012

P

eerfc

2

1

1

erfc2时

利用Q(x)，

x

e

t/2

2

dt，erfc(x)

x

e

t

2

dt，

可得Q(x)

12

erfc，故可用Q函数表示误码率PeQ。抽样判决器输入信噪比：

rA

2

2

n

r1An1，发1

抽样判决器输入信号为：

r0An2，发0

2

A为判决器输入有用信号电压，n1，n2为信道输入的均值为0，方差为n高斯噪声。

通信系统的蒙特卡罗仿真模型如图2所示。编程实现二进制基带通信系统的误码率的蒙特卡罗仿真，并和理论误码率比较。

16

Matlab与通信仿真

图2通信系统的蒙特卡罗仿真模型

四、预习报告：（原理，流程图，编程思想，基本程序）

17

Matlab与通信仿真

18

Matlab与通信仿真

19

Matlab与通信仿真

五、实验报告：（调试好的程序，实验结果及分析）

20

Matlab与通信仿真

21

Matlab与通信仿真

22

Matlab与通信仿真

23

Matlab与通信仿真

实验四模拟调制实验

一、实验目的

掌握线性模拟调制信号的波形及产生方法掌握线性模拟调制信号的频谱特点

掌握线性模拟调制信号的解调方法

掌握线性模拟调制系统的MATLAB仿真实现

二、实验原理

1.AM调制

AM信号的时域表示式：

频谱：

sAM(t)[A0m(t)]cosctA0cosctm(t)cosctSAM()A0[(c)(c)]12[M(c)

M(c)]调制器模型如图1所示：

ct

图1AM调制器模型

AM的时域波形和频谱如图2所示：

时域频域

图2AM调制时、频域波形

2.DSB-SC调制

DSB信号的时域表示式

24

Matlab与通信仿真

频谱：

sDSB(t)m(t)cosct

SDSB()

12

[M(c)M(c)]

DSB的时域波形和频谱如图3所示：

sDSB

t

H

H

c

c时域频域

图3DSB调制时、频域波形

DSB的相干解调模型如图4所示：：

图4DSB调制器模型

3.SSB调制

SSB信号的时域表示式1

sSSB(t)

2

Amcosmtcosct

12

Amsinmtsinct

频谱：分为上边带和下边带，均为双边带的一半。

三、实验内容（任选一）

1.用matlab产生一个频率为1Hz，振幅为1的余弦信源，设载波频率为10Hz，A=2。①AM系统中，当相干解调器输入信噪比为15dB时，系统各点的波形及其频谱，并对调制信号和解调信号进行比较分析。

②分析比较不同信噪比下的解调信号。

2.用matlab产生一个频率为1Hz，振幅为1的余弦信源，设载波频率为10Hz。

①DSB-SC系统中，当相干解调器输入信噪比为15dB时，系统各点的波形及其频谱，并对调制信号和解调信号进行比较分析。

25

Matlab与通信仿真

②分析比较不同信噪比下的解调信号。

3.用matlab产生一个频率为1Hz，振幅为1的余弦信源，设载波频率为10Hz。

①SSB系统中，当相干解调器输入信噪比为15dB时，系统各点的波形及其频谱，并对调制信号和解调信号进行比较分析。

②分析比较不同信噪比下的解调信号。

四、预习报告：（原理，流程图，编程思想，基本程序）

26

Matlab与通信仿真

27

Matlab与通信仿真

28

Matlab与通信仿真

五、实验报告：（调试好的程序，实验结果及分析）

29

Matlab与通信仿真

30

Matlab与通信仿真

31

Matlab与通信仿真

实验五模拟信号数字传输实验（一）

一、实验目的

掌握低通信号抽样定理

理解13折线A率逐次比较型PCM编码仿真的思想

二、实验原理

1．低通信号的抽样定理

一频带限制在0,fH赫内的时间连续信号mt，若以fs2fH速率对mt等间隔

Tsfs12fH抽样，则mt将被所得抽样函数mst完全确定。

2．验证低通信号的抽样定理

抽样的过程是将输入的模拟信号与抽样信号相乘而得，通常抽样信号是一个周期为Ts的周期脉冲信号，抽样后得到的信号称为抽样序列。理想抽样信号定义如下：

T(t)

其中，P(t)

10

t0t0

,fs

1Ts

P(tnT

n

s

)

称为抽样速率。

因此抽样后的信号为

xs(t)x(t)T(t)

经带宽为fH低通滤波器后可恢复抽样信号。3．PCM的基本原理如图1所示：

k

x(kTs)P(tkTs)

图1PCM的基本原理

抽样是对模拟信号进行周期性的扫描，把时间上连续的信号变成时间上离散的信号。我们要求经过抽样的信号应包含原信号的所有信息，即能无失真地恢复出原模拟信号，抽样速率的下限由抽样定理确定。

量化是把经抽样得到的瞬时值进行幅度离散，即指定Q规定的电平，把抽样值用最接近的电平表示。

编码是用二进制码组表示有固定电平的量化值。实际上量化是在编码过程中同时完成的。4．逐次比较型PCM编码的产生原因：

①均匀量化时其量化信噪比随信号电平的减小而下降。产生这一现象的原因就是均匀量

32

Matlab与通信仿真

化时的量化级间隔Δ为固定值，而量化误差不管输入信号的大小均在(-Δ/2,Δ/2)内变化。故大信号时量化信噪比大，小信号时量化信噪比小。对于语音信号来说，小信号出现的概率要大于大信号出现的概率，这就使平均信噪比下降。同时，为了满足一定的信噪比输出要求，输入信号应有一定范围(即动态范围)，由于小信号信噪比明显下降，也使输入信号范围减小。要改善小信号量化信噪比，可以采用量化间隔非均匀的方法，即非均匀量化。

②压扩特性曲线用A律压扩特性曲线，工程上使用13折线来近似A律压扩特性曲线，目的是为了使用逐次比较算法来提高运算速度。

5．逐次比较型编码的编码方法

码位安排：

极性码段落码段内码

M1M2M3M4M5M6M7M8

其中：

①第一位M1表示量化值的极性正负。M1=1代表信号极性为正M1=0代表信号极性为负；后面7位分为段落码和段内码两部分,用于表示量化值的绝对值。

②M2M3M4为段落码，分别对应表1段落编码

③M5M6M7M8为段内码，分别对应表1PCM编码的后四位.

三、实验内容

1．利用MATLAB软件验证低通抽样定理，若低通信号为

x(t)0.1cos(0.15t)0.5cos(4t)

①画出该低通信号的波形；

②画出抽样速率为fs4Hz的抽样序列，画出抽样后的信号波形；

②画出经低通滤波器恢复的波形。

2．设输入一个样值x∈[-2048，+2048]，对x进行A律PCM编码。要求编写成函数，该函数输入变量为样值，输出变量为A律13折线逐次比较的八位PCM编码，设码元宽度为1，画出其波形。

四、预习报告：（原理，流程图，编程思想，基本程序）

33

Matlab与通信仿真

34

Matlab与通信仿真

35

Matlab与通信仿真

五、实验报告：（调试好的程序，实验结果及分析）

36

Matlab与通信仿真

37

Matlab与通信仿真

38

Matlab与通信仿真

实验六模拟信号数字传输实验（二）

一、实验目的

掌握13折线A律逐次比较型PCM编，译码原理

二、实验原理

参见实验五模拟信号的数字传输仿真（一）实验原理部分

三、实验内容

输入信号xtsin4t，对该信号进行采样，要求符合低通采样定理，然后将采样信号进行非均匀量化，对每个样值的量化结果利用模拟信号数字传输的仿真（一）所编写的编码函数进行PCMA律13折线编，译码。

（1）画出原始模拟信号的波形

（2）画出编码结果的波形

（3）试编写程序将编码结果进行译码并通过低通滤波器来恢复模拟信号，对恢复的信

号与原始模拟信号进行比较分析。

四、预习报告：（原理，流程图，编程思想，基本程序）

39

Matlab与通信仿真

40

Matlab与通信仿真

41

Matlab与通信仿真

五、实验报告：（调试好的程序，实验结果及分析）

42

Matlab与通信仿真

43

Matlab与通信仿真

44

Matlab与通信仿真

实验七数字频带传输系统实验

一、实验目的

掌握数字频带传输系统调制解调的仿真过程掌握数字频带传输系统误码率仿真分析方法

二、实验原理

数字频带信号通常也称为数字调制信号，其信号频谱通常是带通型的，适合于在带通型信道中传输。数字调制是将基带数字信号变换成适合带通型信道传输的一种信号处理方式，正如模拟通信一样，可以通过对基带信号的频谱搬移来适应信道特性，也可以采用频率调制、相位调制的方式来达到同样的目的。1.调制过程1）2ASK

如果将二进制码元“0”对应信号0，“1”对应信号Acos2fct，则2ASK信号可以写成如下表达式：

sTtangtnTsAcos2fct

n

10tTs

an0,1，gt。

0其他

可以看到，上式是数字基带信号mt框图如图1所示：

agtnT经过DSB调制后形成的信号。其调制

n

s

n

图12ASK信号调制框图

2ASK信号的功率谱密度为：Psf2）2FSK

A

2

4

[PmffcPmffc]

将二进制码元“0”对应载波Acos2f1t，“1”对应载波Acos2f2t，则形成2FSK信号，可以写成如下表达式：

sTt

gtnTAcos2

n

s

n

f1tn

agtnTAcos2

n

s

n

f2tn

45

Matlab与通信仿真

当an0时，对应的传输信号频率为f1；当an1时，对应的传输信号频率为f2。上式中，n、n是两个频率波的初相。2FSK也可以写成另外的形式如下：

sTtAcos2fct2hangtnTs

n

10tTs

其中，an1,1，fcf1f2/2，gt，hf2f1为频偏。

0其他

2FSK信号也可以看成是两个不同载波的2ASK信号的叠加：

sTts1tcos1t1s2tcos2t2

当这两项不相关时（如载波之间频率差足够大），它的功率谱密度为：

Psf

14

14

[Ps1ff1Ps1ff1][Ps2ff1Ps2ff2]

3）2PSK

将二进制码元“0”对应相位为的载波Acos2fct，“1”对应载波相位为0的载波

Acos2fct，则2PSK信号可以写成如下表达式：

sTtangtnTsAcos2fct

n

10tTs

其中an1,1，gt。其实现框图与2ASK相同，只是输入是双极性

0其他

的信号。

2PSK信号的功率谱密度为：

Psf

A

2

4

[PmffcPmffc]

2．解调过程

数字频带通信系统相干解调基本框图如图2所示。

图2数字频带通信系统相干解调基本框图

3.抗噪性能分析

二进制数字频带传输系统普通接收机误码率公式如表1所示。其中r为解调器输入端信噪比。

46

Matlab与通信仿真

表1二进制数字频带传输系统普通接收机误码率公式

三、实验内容

1.用Matlab产生独立等概的二进制信源。

①画出2ASK调制解调中信号波形及其功率谱变化过程；②画出2FSK调制解调中信号波形及其功率谱变化过程；（设f1f2③画出2DPSK调制解调中信号波形及其功率谱。

三个题目任选其一，并在不同信噪比下对其进行抗噪声性能分析。1