第一次接触matlab程序画流程图，不知道该如何按照这个程序画流程图。

shyw559

本人菜鸟，第一次接触matlab程序画流程图，不知道该如何按照这个程序画流程图， 尤其是这个程序里用了很多for循环， 比较迷茫……麻烦大神指点一下，非常感谢！
+ \, a8 \; G& j8 C! ]) a) E$ q%MRC_scheme.m
, G) A; Z# K- Y. F7 z%接收分集-MRC
8 ~9 X! v9 ]3 }5 r# Qclear,clf
L_frame=130;
N_packet=4000;%130个编码后符号，仿真4000组
6 J% n# u7 T8 u3 kb=4;%设置为1/2/3/4对应于BPSK/QPSK/8PSK/16-QAM
SNRdBs=[0:2:20];
sq2=sqrt(2);
for iter=1:3
    if iter==1
6 [3 j# h' [3 B6 \        NT=1;0
$ f9 e. p! ]7 v, P( d) R$ e1 ?        NR=1;%SISO
        gs='-kx';
# \' T  S2 T2 j' v8 a    elseif iter==2
7 u2 f9 W1 `& `. w3 F        NT=1;
, u" T1 \+ H6 }+ K5 C        NR=2;% 1发2收
        gs='-^';
    else
        NT=1;
        NR=4;% 1发4收
        gs='-ro';
    end
# t0 R9 r4 C2 e, `5 A9 s# @; I    sq_NT=sqrt(NT);
    for i_SNR=1:length(SNRdBs)
. x3 |/ t2 y) d        SNRdB=SNRdBs(i_SNR);
        sigma=sqrt(0.5/(10^(SNRdB/10)));%瑞利分布参数，噪声标准差
5 m1 }5 {& a+ F* h8 _% N  X7 x        for i_packet=1:N_packet
            symbol_data=randint(L_frame*b,NT); %产生L_frame*b行NT列二进制信息
. S$ ^5 }4 x" Q- R) A, f6 y            [temp,sym_tab,P]=modulator(symbol_data.',b);%调制
& s; x. ^, z3 h' W' I$ l            X=temp.';%X为调制后序列（L_frame行1列）
2 v' e% P7 y5 `( F7 K( j- ^- [            Hr=(randn(L_frame,NR)+j*randn(L_frame,NR))/sq2;%瑞利信道衰落系数
- m5 F) @9 X0 L$ s% T0 O            H=reshape(Hr,L_frame,NR);%重置Hr行列（L_frame行NR列）
            Habs=sum(abs(H).^2,2);
, ^) d' U6 S) v0 J3 d, A: {! N            Z=0;
            for i=1:NR
                R(:,i)=sum(H(:,i).*X,2)/sq_NT+...
- {- m+ M0 }  ^                    sigma*(randn(L_frame,1)+j*randn(L_frame,1)); %经过瑞利信道并加复高斯白噪声
0 ~1 v2 Y4 s8 ~3 h% B, f                Z=Z+R(:,i).*conj(H(:,i));%MRC合并，合并系数为conj(H(:,i))
            end
' u% W1 q8 s" _' \0 a            for m=1 %最大似然译码
1 W" j4 }& e3 b; S1 c  r- q                dl(:,m)=abs(sum(Z,2)-sym_tab(m)).^2+...
                    (-1+sum(Habs,2))*abs(sym_tab(m))^2;
            end
            [y1,i1]=min(dl,[],2);
            Xd=sym_tab(i1).';
- Z: f1 z8 J8 q0 P6 e: e1 @            temp1=X>0;
            temp2=Xd>0;
            noeb_p(i_packet)=sum(sum(temp1~=temp2));%累加错误比特
        end
        BER(iter,i_SNR)=sum(noeb_p)/(N_packet*L_frame*b);%计算误比特率
    end
    semilogy(SNRdBs,BER(iter,,gs),hold on
, x( `" `' c% t  Y5 u2 s    axis([SNRdBs([1 end]) 0.5*10^(-4) 1e0])
end
' l5 V. @! \" m* ftitle('BER peRFoemancde of MRC Scheme')
0 y* Q: c2 @* U+ Q( K  }0 q5 Fxlabel('SNR[dB]'),ylabel('BER')
. Y  t  `' H5 c# n. ~2 G) B, Egrid on,set(gca,'fontsize',9)
legend('SISO','MRC(Tx:1,Rx:2)','MRC(Tx:1,Rx:4)')
* |0 ^2 M- W: r5 _+ ?+ B

wu68aq

这个其实很简单。