a=randint（1360）;%a為1行360列的矩陣
T1=bianma（a）;?????%編碼（120*7矩陣）
t0=mod（T12）%把T化成二進制矩陣
t1=reshape（t01840）;%t1為交織后生成的1*840的矩陣
for?p=1:20
????T2=awgn（t1p4）;
T3=T2*2-1;?%把加完噪聲的序列變成1?-1序列以便判決
for?i=1:840;
????if?T3（i）>0
????????receive（i）=1;
????else?
????????receive（i）=0;
????end
end
T4=reshape（receive1207）%去交織
error=0;
for?i=1:840
????if??T4（i）==t1（i）;
????else?error=error+1;
????end
end
ERROR（p）=error/840;