hamming_simulation.m

哈弗曼编码

%Hamming code simulation 
%Simulation for encoding and decoding of a [7,4] Hamming code. The decoder
%can correct one error as shown and as theory states. The table at the end
%of the file shows the various outputs with different error positions and
%message bits. One error can be placed at any of the 7 bit locations and
%corrections made.
clear
n = 7%# of codeword bits per block
k = 4%# of message bits per block
A = [ 1 1 1;1 1 0;1 0 1;0 1 1 ];%Parity submatrix-Need binary(decimal combination of 7,6,5,3)            
G = [ eye(k) A ]%Generator matrix
H = [ A' eye(n-k) ]%Parity-check matrix
% ENCODER%
msg = [ 1 1 1 1 ] %Message block vector-change to any 4 bit sequence
code = mod(msg*G,2)%Encode message
% CHANNEL ERROR(add one error to code)%
%code(1)= ~code(1);
%code(2)= ~code(2);
code(3)= ~code(3);
%code(4)= ~code(4);%Pick one,comment out others
%code(5)= ~code(5);
%code(6)= ~code(6);
%code(7)= ~code(7);
recd = code %Received codeword with error
% DECODER%
syndrome = mod(recd * H',2)
%Find position of the error in codeword (index)
find = 0;
for ii = 1:n
    if ~find
        errvect = zeros(1,n);
        errvect(ii) = 1;
        search = mod(errvect * H',2);
        if search == syndrome
            find = 1;
            index = ii;
        end
    end
end
disp(['Position of error in codeword=',num2str(index)]);
correctedcode = recd;
correctedcode(index) = mod(recd(index)+1,2)%Corrected codeword
%Strip off parity bits
msg_decoded=correctedcode;
msg_decoded=msg_decoded(1:4)



%Error position     Syndrome       Decimal   4 bit word   codeword  dmin
%    1               111             7        0000       0000000
%    2               110             6        0001       0001011     3
%    3               101             5        0010       0010101     4
%    4               011             3        0011       0011110     3
%    5               100             4        0100       0100110     3
%    6               010             2        0101       0101101     3
%    7               001             1        0110       0110011     4
%No error will give syndrome of 000           0111       0111000     3
%                                             1000       1000111     4
%                                             1001       1001100     3
%                                             1010       1010010     4
%                                             1011       1011001     3
%                                             1100       1100001     3
%                                             1101       1101010     3
%                                             1110       1110100     4
%                                             1111       1111111     3
%Any exclusive or additions of any two codewords should give another
%codeword.