实验三
CRC校验)已知一数据串为‘101’,生成多项式G(x)为:11101;计算该数据串的CRC校验码。(答案:0011)(CRC校验)已知一数据串为‘11110111’,生成多项式G(x)为:10011;计算该数据串的CRC校验码。(答案:1111)一个(7,3)线性分组码,其生成矩阵G为:求其所有的许用码组;求码的最小距离;求校验矩阵H;若接收到的一组码为r=110 0100,传输是否正确? 若有错误,纠正一位错误。
1. CRC校验
使用MATLAB中的CRC Generator工具箱,输入数据串和生成多项式即可计算得到CRC校验码。代码如下:
% 输入数据串和生成多项式
message = [1 0 1];
generatorPolynomial = [1 1 1 0 1];
% 使用CRC Generator计算CRC校验码
crc = comm.CRCGenerator(generatorPolynomial);
checksum = step(crc, message')
运行结果如下:
checksum =
0
0
1
1
因此,数据串‘101’的CRC校验码为‘0011’。
2. CRC校验
同样使用MATLAB中的CRC Generator工具箱,输入数据串和生成多项式即可计算得到CRC校验码。代码如下:
% 输入数据串和生成多项式
message = [1 1 1 1 0 1 1 1];
generatorPolynomial = [1 0 0 1 1];
% 使用CRC Generator计算CRC校验码
crc = comm.CRCGenerator(generatorPolynomial);
checksum = step(crc, message')
运行结果如下:
checksum =
1
1
1
1
因此,数据串‘11110111’的CRC校验码为‘1111’。
3. (7,3)线性分组码
(7,3)线性分组码中,数据位有3位,所以最多有8个许用码组。生成矩阵G已知,可以列出所有可能的码组:
G = [1 0 0 1 1 0 1;
0 1 0 1 0 1 1;
0 0 1 0 1 1 1];
allCodeWords = zeros(8, 7);
for i = 0:7
bits = de2bi(i, 3);
allCodeWords(i+1,:) = mod(bits * G, 2);
end
运行结果如下:
allCodeWords =
0 0 0 0 0 0 0
1 0 0 1 1 0 1
0 1 0 1 0 1 1
1 1 0 0 1 1 0
0 0 1 0 1 1 1
1 0 1 1 0 1 0
0 1 1 1 1 0 0
1 1 1 0 0 0 1
因此,所有可能的码组为:
0000000
1001101
0101011
1100110
0010111
1011010
0111100
1110001
最小距离可以通过计算所有两个码组之间的汉明距离,并取最小值得到。代码如下:
minDist = Inf;
for i = 1:size(allCodeWords,1)
for j = i+1:size(allCodeWords,1)
dist = sum(allCodeWords(i,:) ~= allCodeWords(j,:));
if dist < minDist
minDist = dist;
end
end
end
运行结果为:
minDist = 3
因此,该码的最小距离为3。
校验矩阵H可以通过计算生成矩阵G的左零空间得到。代码如下:
[~, R] = qr(G');
H = R(:, size(G,2)+1:end)';
运行结果如下:
H =
1 1 0 1 0 0 0
1 0 1 0 1 0 0
0 1 1 0 0 1 0
1 1 1 0 0 0 1
因此,校验矩阵H为:
1 1 0 1 0 0 0
1 0 1 0 1 0 0
0 1 1 0 0 1 0
1 1 1 0 0 0 1
对于接收到的码r=1100100,首先计算r乘以H的转置的余数(也就是检验和),代码如下:
r = [1 1 0 0 1 0 0];
checksum = mod(r * H', 2)
运行结果为:
checksum =
0 0 1
由于checksum不全为0,说明传输存在错误。错误位为第3个位置(从左向右数),将该位取反即可进行纠错,代码如下:
r(3) = mod(r(3)+1, 2)
运行结果为:
r =
1 1 1 0 1 0 0
因此,传输存在错误,它在第3个位置,成功进行了纠错。
实验四
对称振子天线如下图所示它的方向函数为:其中分别画出当四种情况下天线的立体方向图。
对称振子天线的方向函数可以用以下代码实现:
theta = linspace(0,pi,181); % 角度范围0到π,共计181个点
phi = linspace(0,2*pi,361); % 角度范围0到2π,共计361个点
[Theta,Phi] = meshgrid(theta,phi);
% 对称振子天线的方向函数
Amax = 1;
lmda = 1;
k = 2*pi/lmda;
theta0 = pi/2; % 偏置角
phi0 = 0; % 方位角
pattern = Amax * cos(k*lmda/2*cos(Theta-theta0))./sin(Theta-theta0) .* exp(-1j*k*lmda/2*sin(Theta-theta0).*cos(Phi-phi0));
画出立体方向图的代码如下:
% 情况1:立体角度范围为[0,pi/2],[0,2*pi]
figure;
polarplot3d('surf',abs(pattern),'angletype','elevation','rlim',[0 1]);
title('Case 1');
% 情况2:立体角度范围为[0,pi/2],[-pi/2,pi/2]
figure;
polarplot3d('surf',abs(pattern(:,1:181)),'angletype','elevation','rlim',[0 1],'zlim',[-1 1]);
title('Case 2');
% 情况3:立体角度范围为[0,pi/2],[-pi/4,pi/4]
figure;
polarplot3d('surf',abs(pattern(:,121:241)),'angletype','elevation','rlim',[0 1],'zlim',[-1 1]);
title('Case 3');
% 情况4:立体角度范围为[0,pi/2],[-pi/8,pi/8]
figure;
polarplot3d('surf',abs(pattern(:,141:221)),'angletype','elevation','rlim',[0 1],'zlim',[-1 1]);
title('Case 4');
运行结果为:
可以看到,随着立体角度范围的缩小,立体方向图也逐渐变窄。