線性分組碼的編碼與譯碼-資料下載頁(yè)
2025-08-05 10:50本頁(yè)面
【正文】 for i=1:1:k %要求為輸入字符，否則默認(rèn)輸入的為數(shù)值 Q_charA(i,:)=input(sprintf(39。輸入一個(gè)校正子，以檢測(cè)C%d位的錯(cuò)碼(%d位)39。,k+ri,r),39。s39。)。 for j=1:1:r %將字符轉(zhuǎn)換為數(shù)值后寫(xiě)入矩陣 Q(i,j)=str2double(Q_charA(i,j))。 end endend% % 生成完全的校正子與錯(cuò)碼位置之間的對(duì)應(yīng)關(guān)系 %生成單位矩陣（字符型），以補(bǔ)齊校正子for i=r1:1:0 Qr(ri)=2^i。endQ_charB=dec2bin(Qr)。%將十進(jìn)數(shù)格式轉(zhuǎn)換為二進(jìn)制表示方式%生成校正子與錯(cuò)碼位置的對(duì)應(yīng)關(guān)系Q_charC(1:k,:)=Q_charA。Q_charC(k+1:k+r,:)=Q_charB。Q_charD=reshape(sprintf(39。C%02d39。,k+r1:1:0),8,k+r)39。disp(39。校正子與錯(cuò)碼位置的對(duì)應(yīng)關(guān)系39。)。Q_char=[Q_charC Q_charD]% % 線性分組碼內(nèi)部的函數(shù)關(guān)系及計(jì)算 Ik=eye(k)。%生成k*k的單位矩陣，與Q合成為生成矩陣GG=[Ik Q]。%生成矩陣P=Q39。Ir=eye(r)。%生成r*r的單位矩陣，與P合成為監(jiān)督矩陣HH=[P Ir]。%監(jiān)督矩陣% % 計(jì)算所有的可用碼組 valid_yesOrNo=input(39。是否要顯示所有可用碼組?Y/N:39。,39。s39。)。if(valid_yesOrNo==39。Y39。) %初始化2^k行k+r列，與1行k列的矩陣 valid_codes=zeros(2^k,k+r)。 valid_buffer=zeros(1,k)。 %初始化每行碼重?cái)?shù)組 weight_array=zeros(1,2^k)。 %遍歷可用信息碼并輸出所有可用碼組 for i=1:1:2^k %將第i行的十進(jìn)制數(shù)轉(zhuǎn)制為4位二進(jìn)制數(shù)（字符矩陣） valid_binary=dec2bin(i1,k)。 for j=1:1:k %將字符轉(zhuǎn)換為數(shù)值放入矩陣（行向量）中 valid_buffer(j)=str2double(valid_binary(j))。 end %將得出的第i行的行向量valid_codes,并得出當(dāng)前信息碼對(duì)應(yīng)全長(zhǎng)碼字 valid_codes(i,:)=mod(valid_buffer*G,2)。 %求每行可用碼字的碼重 weight_array(i)=length(nonzeros(valid_codes(i,:)))。 end disp(39。最終可用碼組：39。)。 valid_codes disp(39。最小碼重：39。)。 min_weight=min(weight_array(2:2^k)) %count=1。 for i=1:1:2^k1 for j=i+1:1:2^k %兩兩的碼距存入數(shù)組中 %distance_array(count)=length(nonzeros(valid_codes(i,:)valid_codes(j,:)))。 distance_array(i,j)=length(nonzeros(valid_codes(i,:)valid_codes(j,:)))。 %count=count+1。 end end %調(diào)整碼距矩陣 disp(39。最小碼距：39。)。 %min_distance=min(min(distance_array)) distance_array(find(distance_array==0))=NaN。 min(min(distance_array))end% % 輸入或自動(dòng)生成碼字及校驗(yàn)碼字的正確性 check_yesOrNo=input(39。是否自動(dòng)生成要校驗(yàn)的碼字？Y/N:39。,39。s39。)。if(check_yesOrNo==39。Y39。) disp(39。校驗(yàn)碼字39。)。 check_code=randint(1,k+r,[0 1])else %直接輸入矩陣更好，但此處就當(dāng)輸入的字符 check_char=input(sprintf(39。請(qǐng)輸入要校驗(yàn)的碼字(%d位): 39。,k+r),39。s39。)。 %初始化要校驗(yàn)的碼字 check_code=zeros(1,k+r)。 for i=1:1:k+r %將字符陣轉(zhuǎn)換為數(shù)值陣 check_code(i)=str2double(check_char(i))。 endenddisp(39。得出校驗(yàn)子為:39。)。S=mod(H*check_code39。,2)39。%新建最終校驗(yàn)矩陣，檢驗(yàn)錯(cuò)誤碼位置Q_final(1:k,:)=Q。Q_final(k+1:k+r,:)=Ir。flag=1。%校驗(yàn)標(biāo)記，默認(rèn)為1，檢測(cè)不出，為2，無(wú)錯(cuò)，為0，有錯(cuò)碼for i=1:1:k+r %判斷兩矩陣是否相同，是則返回1，否則返回0 if(isequal(S,Q_final(i,:))) errors=sprintf(39。錯(cuò)碼位置：C%d，原碼應(yīng)為：39。,k+ri) flag=0。 %糾錯(cuò) modify_code=check_code。 modify_code(i)=mod(check_code(i)+1,2)。 modify_code elseif(isequal(S,zeros(1,r))) flag=2。 endend%有兩位及以上出錯(cuò)的碼字，檢測(cè)不出，或者無(wú)錯(cuò)if(flag==1) errors=sprintf(39。檢測(cè)不出，冏...39。)elseif(flag==2) errors=sprintf(39。無(wú)錯(cuò)39。)end
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1