TA的每日心情 | 怒 2019-11-19 15:34 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
( b* j+ A- g1 x, N. N, P) |
msg=gf(randint(1,7),1); %information code1 h+ A5 m4 ~6 p" F5 S% E
c=bchenc(msg,15,7); %encode
% |, P- [7 n8 E# j7 w" s# i& j' ^e=gf([1 0 0 0 0 0 0 0 0 0 0 1 0 0 0],1);%2 bit errors8 X7 X$ F" m- T' ~( k% J6 F. G
r1=c+e;% receive code1 }( ]/ t9 K! F: N, H* n% F/ {
r=r1;
+ o# J7 e9 B* q* Z: {g=gf([1 1 1 0 1 0 0 0 1],1);%generator polynomial
6 I! c- y1 Y6 `0 L% i9 Q# mm=3;n=0;k=0;
$ A% r8 }1 U% M% t6 I8 vwhile m>2 & k<30
/ }8 U' {$ _/ Xk=k+1;* ? A+ j# U+ J) N/ `7 L* K& h
IF r(1)==1
" Z1 E% w; s- Wfor j=1:9
' B0 @1 u& I3 tr(j)=r(j)+g(j);
% e, X/ N# l% w; K( ]: Aend/ {2 b# j F0 H+ {3 Y! \
end# _; o% w5 N& r. ]: Q
m=sum(r==1);- j7 |9 O8 Z( |/ c
while r(1)==0 & m>2 & k<30 %left shift0 y% \7 w. K9 e9 p! P q t9 @" T
n=n+1; 6 ?, U" z, x5 M- S
for i=1:14
4 K+ i9 Q+ g6 {8 Er(i)=r(i+1);, u( ~6 W J, {. O$ y8 c& m
end( @) p6 C w. F) B8 ^* b2 Z
r(15)=0;
% g% V* x3 _+ f& j" ?+ Q: p2 bend
! t2 `& ]8 r* b+ o/ ?4 |4 z4 Jend Q" G0 }- O/ \+ l8 z: s! o
s=mod(n,15);2 P4 E: V3 i* A2 U- |
for k=1:s %cyclic right shift
7 u. C4 Y' w% J$ m1 qt=r(15);: u8 c4 I9 A! X- P& T
for i=0:13- }- k( y8 q W; x6 \' `
r(15-i)=r(14-i);
, P% n, B* ? } K0 P$ a7 Yend
' T4 B: V$ M2 j ^r(1)=t; d5 D- s+ V6 p
end( @, x- M1 E! P+ [( x, W
e1=r;
" {; h5 n* U$ R% t% N7 I2 ]m1=sum(e1==1);
* N+ I- j1 t) [$ I: q* sif m1<3
& o" Z# {! F7 m9 kdisp('OK e(x)=r(x) mod(g(x))');
! {, ~' @& R+ a. X* dend
: I6 I1 c* P5 v, x5 v) ? [1 [7 Z" N& Z: {3 P2 @
% Q. K0 c/ K/ i5 K$ X
想问一下这段代码是什么意思,最好能具体到每一个循环 是能够直接运行的嘛? 他的输出有哪些?最后disp那个输出的是什么意思? |
|