【无标题】RFID实验

实验三

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');

运行结果为：

可以看到，随着立体角度范围的缩小，立体方向图也逐渐变窄。