function[T]=zaosheng(t,m)
%函數用于添加噪聲
%T為輸出一個500行15列的加噪聲后的矩陣
%t為輸入一個1行7500列的矩陣
%m為長度7500的噪聲序列中1的個數
r1=randerr(1,7500,m);%r1為模擬信道矩陣，誤碼率為m/7500
t1=r1+t;%t1為第一次加過噪聲后的編碼
T1=reshape(t1,15,500);
T1=T1';%T1為變換后的矩陣
T=mod(T1,ones(500,15)*2);%把T1化成二進制矩陣