【通信原理】数字基带传输的线路码型

数字基带传输的线路码型

简单介绍数字基带传输的线路码型的信号波形的特点，以及生成方法。注意观察频谱。文末附Matlab代码。

以下包括双极性NRZ、单极型NRZ、双极型RZ、单极型RZ、差分码，曼切斯特码/数字双相码、密勒码、CMI码、AMI码、HDB3码。

参数：Rb = 20 bps；fs = 200 Hz；仿真时间2s。

PART1

1. 双极性NRZ

   • “1”对应“+1V”，“0”对应“-1V”
   • 双极性的优点是抗噪能力强一些
   • 缺点是生成电路需要正负双电源供电。
   • 不归零信号抗噪声能力较强

2. 单极性NRZ

   • “1”对应“+1V”，“0”对应“0V”
   • 单极性的优点是可以采用单电源供电
   • 缺点是具有直流分量，只能在直流耦合的电路中使用。如果在交流耦合的电路中，比如通过一个电容器，则会过滤掉直流分量。

3. 双极性RZ

   • “1”对应“+1V”，“0”对应“-1V”。维持其电平半个bit时间（$T_b/2$），将电平回到零电平。
   • 归零信号跳变边缘丰富，易于定时

4. 单极性RZ

   • “1”对应“+1V”，“0”对应“0V”。“1”维持正电平半个bit时间（$T_b/2$），将电平回到零电平。

5. 差分码

   • “1”对应“上跳电平”或者“下跳电平”，“0”保持上一bit对应的电平。
   • 第一个bit如果为1，上跳或者下跳则由电路复位后保持的电平确定。本仿真中初相（或者说第0个bit）为0。
   • 生成方法是进行$m\_ne{w}_i=xor(m_i,m_{i-1})$

PART2

6. 数字双相码/曼切斯特码
   • “1”对应“下跳电平”，“0”对应“上跳电平”。
   • 数字双相码没有直流分量，包含丰富的定时信息。
   • 但占用的带宽加大。用于10Mbps的以太网中。
7. 密勒码/延迟调制码
   • “1”对应“上跳电平”或者“下跳电平”，“0”对应负或正电平，连0时要交替。密勒码是数字双相码使用下降沿触发一级触发器后的结果，克服了数字双相码的定时相位不确定性。
8. 传号反转码(CMI)
   • “1”对应交替的正或负电平，“0”对应“上跳脉冲”。
   • 没有直流分量，跳变丰富，容易接收、易于检错。
9. 交替传号反转码（AMI）
   • “1”——交替的(+1)或(-1)，“0”对应“0V”。
   • 长连“0”串时，无法提取时钟。
   • AMI码3电平信号，是北美电话系统中的接口标准之一。
10. HDB3码
    • 规则：先将码改为AMI码，然后将四连0调整为“0 0 0 V”或者“B 0 0 V”，保证V和前面的B同极性（不同极性的时候就调整为“B 0 0 V”），同时V、B各自正负交替出现。
    • HDB3具有AMI码的优点，译码简单，且容易提取定时时钟。

PART3 Matlab代码

% ---------------------------------------------
% 作者：wsg
% 日期：2022年12月4日
% 邮箱：shouguo19@163.com
% 主题：基带传输的线路码型的仿真
% 参考：李晓峰《通信原理（第二版）》，第四章第九节。
% ---------------------------------------------

% 基本参数
Rb = 20;              % 信息bit速率
t_span = 2;           % 仿真时间
fs = 200;            % 采样频率

unitBitSa = fs/Rb;          % 单位原始信息采样点数
bit_length = Rb*t_span;     % 信息序列长度

t = transpose(0:1/fs:(t_span-1/fs));    % 仿真时间轴 
N_points = length(t);
f = (1/N_points:1/N_points:1/2)*fs;     % 频率轴

m = randi([0 1], bit_length, 1);        % 随机产生的bit序列

name_list = ["双极性NRZ","单极性NRZ","双极性 RZ", "单极性 RZ", "差分码",...
             "曼切斯特码", "密勒码", "传号反转码CMI", "交替传号反转码AMI","HDB3码"];

%% part1
rz_gt = zeros(unitBitSa,1);
rz_gt(1:unitBitSa/2) = ones(unitBitSa/2, 1);    % 归零码生成滤波器
nrz_gt = ones(unitBitSa,1);                     % 不归零码生成滤波器

% 1. 双极性NRZ
s1 = filter(nrz_gt, 1, upsample(m*2-1, unitBitSa));

% 2. 单极性NRZ
s2 = filter(nrz_gt, 1, upsample(m, unitBitSa));

% 3. 双极性 RZ
s3 = filter(rz_gt, 1, upsample(m*2-1, unitBitSa));

% 4. 单极性 RZ
s4 = filter(rz_gt, 1, upsample(m, unitBitSa));

% 5. 差分码
m5 = zeros(bit_length,1);
tmp = 0;    % 初相
for i = 1:bit_length
    tmp = xor(tmp,m(i));
    m5(i) = tmp;
end
s5 = filter(nrz_gt, 1, upsample(m5, unitBitSa));


%% part2
down_pulse = zeros(unitBitSa,1)-1;
down_pulse(1:unitBitSa/2) = ones(unitBitSa/2, 1);
up_pulse = flipud(down_pulse);

% 6. 曼切斯特码/数字双相码
m_tmp = [m,~m]';
m6 = m_tmp(:);
s6 = filter(ones(unitBitSa/2,1),1,upsample(m6*2-1,unitBitSa/2));


% 7. 密勒码/延迟调制码
bph = [m,~m]';
bph = bph(:);

m7 = zeros(bit_length*2,1);
tmp = 0;
for i = 1:bit_length*2
    if isequal([tmp,bph(i)], [1,0])
        m7(i) = ~m7(i-1);
    else
        if i==1  
            m7(i) = tmp;
        else
            m7(i) = m7(i-1);
        end
    end
    tmp = bph(i);
end
s7 = filter(ones(unitBitSa/2,1),1,upsample(m7*2-1,unitBitSa/2));


% 8. CMI
s8 = zeros(N_points,1);
flag = 1;
for i = 1:bit_length
    if m(i)==1
        if flag == 1 
            s8((1:unitBitSa)+(i-1)*unitBitSa) = 1;
        else 
            s8((1:unitBitSa)+(i-1)*unitBitSa) = -1;
        end
        flag = ~flag;
    else
        s8((1:unitBitSa)+(i-1)*unitBitSa) = up_pulse;
    end
end


% 9.AMI
s9 = zeros(N_points,1);
flag = 1;
for i = 1:bit_length
    if m(i)==1
        if flag == 1 
            s9((1:unitBitSa)+(i-1)*unitBitSa) = 1;
        else 
            s9((1:unitBitSa)+(i-1)*unitBitSa) = -1;
        end
        flag = ~flag;
    else
        s9((1:unitBitSa)+(i-1)*unitBitSa) = 0;
    end
end

% 10. HDB3
m10 = zeros(bit_length,1);
flagb = 0;
flagv = 1;
for i = 1:bit_length
    if m(i)==1
        flagb = ~flagb;
        if flagb == 1 
            m10(i) = 1;
        else 
            m10(i) = -1;
        end
    else
        m10(i) = 0;
        if i>3 && isequal(m10(i-3:i),zeros(4,1))    
            if flagv == flagb
                if flagv == 1
                    m10(i) = 1;
                else
                    m10(i) = -1;
                end
            else
                if flagv == 1
                    m10(i-3:i) = [1;0;0;1];
                else
                    m10(i-3:i) = [-1;0;0;-1];
                end
                flagb = ~flagb;
            end
            flagv = ~flagv;
        end
    end
end
s10 = filter(ones(unitBitSa,1),1,upsample(m10,unitBitSa));

%% 绘图
figure(3)
for i = 1:5
    s = eval("s"+num2str(i));
    S = abs(fft(s));

    subplot(5,2,i*2-1)
    plot(t,s);xlabel("t/s");ylim([-1.2 1.2]);title(name_list(i))
    subplot(5,2,i*2)
    plot(f,S(1:N_points/2));xlabel("f/Hz");
end

figure(2)
for i = 1:5
    s = eval("s"+num2str(i+5));
    S = abs(fft(s));

    subplot(5,2,i*2-1)
    plot(t,s);xlabel("t/s");ylim([-1.2 1.2]);title(name_list(i+5))
    subplot(5,2,i*2)
    plot(f,S(1:N_points/2));xlabel("f/Hz");
end
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174