基于匯編語(yǔ)言的BCH解碼校驗(yàn)算法

在信號(hào)傳輸中，BCH碼以其獨(dú)特的優(yōu)點(diǎn)被廣泛應(yīng)用于微機(jī)級(jí)的通信中，但因其算法復(fù)雜，通常只用在動(dòng)態(tài)實(shí)時(shí)的無(wú)線通信中，而對(duì)更底層的單片機(jī)級(jí)的信號(hào)傳輸糾錯(cuò)，往往只采用奇偶校驗(yàn)等簡(jiǎn)單的校驗(yàn)方法。本文結(jié)合一些測(cè)控系統(tǒng)和監(jiān)控系統(tǒng)的開(kāi)發(fā)，摸索出了在實(shí)時(shí)動(dòng)態(tài)單片機(jī)中的BCH解碼檢糾方法，并通過(guò)匯編語(yǔ)言加以實(shí)現(xiàn)，取得了一定的效果。下面以BCH（15，7）碼為例進(jìn)行探討。

　　1 BCH碼在單片機(jī)中的放置結(jié)構(gòu)

　　BCH碼作為一種檢糾能力較強(qiáng)的循環(huán)碼，由信息多項(xiàng)式i(x)和監(jiān)督多項(xiàng)式j(luò)(x)組成。這里以c(x)表示整個(gè)BCH(15,7)碼的15位碼組多項(xiàng)式，則有：

　　在單片機(jī)中其放置的具體結(jié)構(gòu)如下：

　　其中，7位信息位放入寄存器R3中，8位監(jiān)督位放入寄存器R4中。

　　2 BCH解碼校驗(yàn)原理

　　二元BCH(15,7)碼的解碼校驗(yàn)原理是在時(shí)域上直接利用碼的代數(shù)結(jié)構(gòu)進(jìn)行解碼。首先，由于BCH(15,7)碼的糾錯(cuò)能力t=2，所以根據(jù)接收序列計(jì)算伴隨式sk=r(αk)，其在伽羅華域GF(24)上的規(guī)定連續(xù)根為α、α2、α3、α4。與其對(duì)應(yīng)的伴隨式分別為： s1=r(α)，s2=r(α2)，s3=r(α3)，s4=r(α4)。

　　然后，由伴隨式計(jì)算差錯(cuò)定位多項(xiàng)式[1]的系數(shù)。在二元BCH碼中，對(duì)于任何值都有s2k=s2k；同理可推，s4=s24=s41，s6=s23 等。所以在求差錯(cuò)定位多項(xiàng)式的系數(shù)時(shí)，僅須用到奇數(shù)下標(biāo)的伴隨式值。就BCH(15,7)碼而言，根據(jù)s1和s3這兩個(gè)伴隨式值便可計(jì)算出差錯(cuò)定位多項(xiàng)式的2個(gè)系數(shù)： σ1=s1和σ2=s3+s31 s1。

　　最后，依據(jù)Chien氏搜索算法對(duì)碼的每個(gè)位置逐位檢索，以確定其錯(cuò)誤位置。若s1=s3=0，則可判定無(wú)差錯(cuò)發(fā)生；若s31+s3=0，則有1個(gè)差錯(cuò)發(fā)生，錯(cuò)碼位置就是s1；若有2個(gè)或2個(gè)以上的差錯(cuò)發(fā)生，則可按σ1αi+σ2α2i=1（i=0,1,2,…，14）進(jìn)行搜索。若在搜索中找到的根少于2個(gè)，則說(shuō)明該多項(xiàng)式有的根在定位域之外，這表明發(fā)生的差錯(cuò)已超過(guò)2個(gè)；若找到的根恰好等于2個(gè)，則表示剛好有2個(gè)差錯(cuò)發(fā)生，可根據(jù)差錯(cuò)位置予以糾正。經(jīng)差錯(cuò)定位找到差錯(cuò)位置后，便可進(jìn)行糾錯(cuò)了。糾錯(cuò)的原理相對(duì)來(lái)說(shuō)比較簡(jiǎn)單，因?yàn)閱纹瑱C(jī)處理的是二進(jìn)制數(shù)，而二進(jìn)制數(shù)只有2個(gè)狀態(tài)，即不是“0”就是 “1”，因此糾錯(cuò)只須將對(duì)應(yīng)差錯(cuò)位取反。

　　3 BCH解碼校驗(yàn)算法的匯編語(yǔ)言實(shí)現(xiàn)

　　具體的解碼程序采用單片機(jī)的匯編語(yǔ)言實(shí)現(xiàn)，包括1個(gè)主程序和6個(gè)子程序。主程序的工作流程是整個(gè)程序的主線，決定著解碼的效率；而子程序則是為了提高主程序在伽羅華域上代數(shù)運(yùn)算的效率，優(yōu)化主程序的程序結(jié)構(gòu)。主程序的清單如下：

　　MOV03H，R3

　　MOV04H，R4

　　MOVR1，#60H；錯(cuò)誤位置初始地址

　　MOVR7，#00H；出錯(cuò)個(gè)數(shù)初始值

　　MOVR0，#00H；Chien氏搜索的初始值

　　LCALLS1；調(diào)用s1=r(α)子程序

　　MOVA，71H

　　CJNEA，#00H，L1；s1≠0

　　LCALLS3；調(diào)用s3=r(α3)子程序

　　MOVA，70H

　　CJNEA，#00H，L1；s3≠0

　　LJMPRIGHT；送至解碼輸出程序

　　L1：MOVA，71H

　　MOV78H，A；s1的矢量值

　　LCALLTAB2；s1的指數(shù)

　　MOVB，A

　　RLA

　　ADDA，B；得到s31的指數(shù)

　　LCALL DIV15；調(diào)用模15求余子程序

　　LCALL TAB1；s31的矢量值