NR 物理层编码: CRC循环冗余码

最新推荐文章于 2022-11-14 18:21:13 发布

明朝百晓生

最新推荐文章于 2022-11-14 18:21:13 发布

阅读量836

点赞数

分类专栏：芯片开发通讯技术文章标签： java servlet 数据库

4AM_明朝百晓生

本文链接：https://blog.csdn.net/chengxf2/article/details/126638220

版权

线性分组码循环冗余码编码理论模二除法错误检测

关键词由CSDN通过智能技术生成

芯片开发同时被 2 个专栏收录

94 篇文章 35 订阅

订阅专栏

通讯技术

54 篇文章 9 订阅

订阅专栏

前言：

循环码是线性分组码的一个重要分支

1：线性分组码编解码思考

2： CRC循环冗余码讲解

通信原理概论8.6线性分组码和8.7循环码_哔哩哔哩_bilibili

信道编码/差错控制码/编码理论/数字通信/线性分组码、循环码、卷积码_哔哩哔哩_bilibili

一线性分组码编解码思考

发送方编码

发送信息为u ：1 行k列向量

$c=uG$ : 1行 n列向量

系统编码矩阵为

$G=[I_k,Q_{k,m}]$ [k,n]矩阵

m=n-k

Q主要是编码规则，I 取信息位

例：

$p_1=u_1+u_2,p_2=u_2+u_3,p_3=u_3+u_4$

$c=[u_1,u_2,u_3,u_4.p_1,p_2,p_3,p_4]$

接收方译码

接收方收到的信息位

$R=C+E$ $=>E=R+C$

E: 错误图样，理想值为全0向量

伴随式,只跟E有关系

$S=RH^T=(C+E)H^T=CH^T+EH^T=EH^T$

当S不全为0,我们解码出E的时候，再根据

$C=R+E$ 可以解码出发送码字（上面+都是模二加法）

其中H为校验矩阵

$H=[Q^T,I_m]_{m,n}$

原理是

$u_1+u_2+p_1=0,u_2+u_3+p_2=0,u_3+u_4+p_3=0$

我们要解码出E,本质上就是奇偶校验矩阵H中的

行重除以2种可能，也就是后面LDPC码为什么要强调稀疏

以上面H 为例， $E=[e_1,e_2,e_3,e_4,e_5,e_6,e_7]$

$e_1+e_2+e_5=0,e_2+e_3+e_6=0,e_3+e_4+e_7=0$

就是解上面的线性方程组，结果可能为0，或1，对应的e组合为

$2^7=128$ ，可以先做一个码表，比如要使得e_1+e_2+e_5=1,这种组合

相对于整体是一个固定的。

# -*- coding: utf-8 -*-
"""
Created on Thu Sep  1 11:09:22 2022

@author: chengxf2
"""
import torch
import numpy as np
import copy 


Bits =[]
'''
生成N个二进制位的组合
'''
def getBit(k,ret):
    
    retLeft =  copy.deepcopy(ret)
    retRight = copy.deepcopy(ret)
    retLeft.append(0)
    retRight.append(1)
    
    N = len(retRight)
    if N<k:
     
        getBit(k,retLeft)
        getBit(k,retRight)
    else:
       
        Bits.append(copy.deepcopy(retLeft))
        Bits.append(copy.deepcopy(retRight))
        
        
def Check(p1,p2,p3):
    
    N = len(Bits)
    num = 0
    for index in range(N):
        item = Bits[index]
        #print(item)
        
        c1  =  item[0]+item[1]+item[4]
        c2 =   item[1]+item[2]+item[5]
        c3 =   item[1]+item[3]+item[6]
        
        #print(item,"\t  index%d 奇偶校验: "%index,c1%2,c2%2,c3%2,((c3%2)and p3))
        if ((c1%2) == p1)  and ((c2%2) == p2)  and ((c3%2)== p3):
            print("\n index ",item)
            num= num+1
    return num
        
     
            
'''
p1: e1+e2+e5 =0
p2: e2+e3+e6 =0
p3  e3+e4+e7 = 0
'''
def Concomitant( P1,P2, P3):

       num = Check(p1,p2,p3)
     
       print("\n num ",num)

if __name__ =="__main__":    
    ret =[]
    k=7
    getBit(k,ret)
    
    p1 = 1
    p2 = 0
    p3 = 1
    Concomitant(p1,p2,p3)

    #print("\n 所有的bit流组合为 ",Bits)

二 CRC循环冗余码讲解

2.1 循环除法

码字转换为生成多项式 ,生成多项式作模二除法,得到的结果转换为码字

例： 10111 对110的模2除法

$c=[10111],c(x)=x^4+x^2+x+1$