a=randint(1,5500);%a為1行5500列的矩陣
T=bianma(a);%編碼
t=reshape(T',1,7500);%t為編碼后的1行7500列的矩陣
T1=zaosheng(t,50);%加噪聲
t=jiucuo(T1);%對漢明編碼進行一位糾錯解碼
T=reshape(t,15,500);
T=T';%行列變換
T=T(:,1:11);%解碼
a1=reshape(T,5500,1);
a1=a1';%行列變換，得到a1為糾錯解碼后的輸出碼
n=xor(a1,a)*ones(5500,1);%n為誤碼個數
f=n/5500;%f為誤碼率