TA的每日心情 | 怒 2019-11-19 15:34 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
: `3 W7 E# P2 E3 j
msg=gf(randint(1,7),1); %information code
+ ~+ y# o# A1 O B4 v+ p7 Yc=bchenc(msg,15,7); %encode; O2 ]5 @% L. k
e=gf([1 0 0 0 0 0 0 0 0 0 0 1 0 0 0],1);%2 bit errors0 L+ P2 w# n5 @5 N6 G( x
r1=c+e;% receive code
6 @7 c2 _* k5 ^' C5 j4 er=r1;& [4 o0 V: Q7 c! m* k. a/ }
g=gf([1 1 1 0 1 0 0 0 1],1);%generator polynomial
$ C" A) Y: e# V% V+ p6 @5 J# um=3;n=0;k=0;
$ G4 B- B3 t2 G! c& @8 c) Ywhile m>2 & k<303 T9 m: J1 I7 y# a! w. J u+ V
k=k+1;& ]3 R$ r! M9 O$ R$ F* P; J* m; h' x
IF r(1)==1
* k5 i# U" ~- @) W3 n! Lfor j=1:9
3 n# ]7 q/ h0 B9 d: Rr(j)=r(j)+g(j);6 h _/ x7 [" s
end. s5 s% P( D* g3 a9 {$ \
end
* j; U) Z0 X; N% b3 g" nm=sum(r==1);* S7 n8 S/ e0 R) j
while r(1)==0 & m>2 & k<30 %left shift
! v$ X. {4 \. f) B2 Bn=n+1;
8 ]0 j$ O4 [$ U1 g8 |. J1 |+ cfor i=1:141 X7 @( v+ y( }. Y3 ~+ M0 N5 X9 v
r(i)=r(i+1);9 K9 D/ W8 w& [ ^$ ]$ e- V( K
end' ?% V/ B. l5 x) H7 d
r(15)=0;* F) f3 J" T# h" ^4 J
end
4 ^- r: L5 C+ I8 P1 ~! M8 X Wend0 F9 R" I6 o- F8 e
s=mod(n,15);( ?) }1 V. D2 J
for k=1:s %cyclic right shift' P( W+ M& V2 F* A% c. ]9 `
t=r(15);& g8 f* Y s7 }- h& _% e7 \
for i=0:13
' [/ i& c* L; C+ I5 h& Nr(15-i)=r(14-i);
3 H, Q6 d! m: R9 Q0 rend" ~0 a) T% M4 S$ O3 W
r(1)=t; # O* o d' L2 x' A+ e; z; V! b% }
end
0 n' q: y9 _- ke1=r;- N' {9 [) l5 o! Q$ z
m1=sum(e1==1);( u: ]! y# e% N- C1 u' Z& S
if m1<3
2 S! u$ t9 k& Y% jdisp('OK e(x)=r(x) mod(g(x))');, _# a) a$ J3 t7 R0 G2 ~0 ^
end
c5 ~# t7 A; J8 Q9 [- k
$ V: O/ i5 Y& m: j6 G; Y; |. v, ?" \6 W6 c& @
想问一下这段代码是什么意思,最好能具体到每一个循环 是能够直接运行的嘛? 他的输出有哪些?最后disp那个输出的是什么意思? |
|