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

shyw559

本人菜鸟，第一次接触matlab程序画流程图，不知道该如何按照这个程序画流程图， 尤其是这个程序里用了很多for循环， 比较迷茫……麻烦大神指点一下，非常感谢！
+ J0 j- b. a1 M0 M2 c1 \%MRC_scheme.m
%接收分集-MRC
# L2 \) Q+ y$ {4 G% tclear,clf
L_frame=130;
N_packet=4000;%130个编码后符号，仿真4000组
b=4;%设置为1/2/3/4对应于BPSK/QPSK/8PSK/16-QAM
3 j) L5 [' c" e5 OSNRdBs=[0:2:20];
: J) @1 L: i/ z0 ksq2=sqrt(2);
for iter=1:3
# f% _) I; E" t) W    if iter==1
1 L9 {. }  o9 L1 |  f        NT=1;0
        NR=1;%SISO
* p6 V, P. n  F8 ?2 ~        gs='-kx';
    elseif iter==2
        NT=1;
) }' G3 I# m3 ]0 x        NR=2;% 1发2收
7 Q" L% x( m2 d7 `" }# ~        gs='-^';
    else
. l- e, x! ?  Y  W: {        NT=1;
        NR=4;% 1发4收
3 K" G- Y7 E& R" W6 R( W        gs='-ro';
    end
; c! y9 s7 J- L% J0 N/ R    sq_NT=sqrt(NT);
' o: J& g# l9 Z. Y0 m8 G+ O    for i_SNR=1:length(SNRdBs)
6 o! s! I8 {0 c* b6 w5 I        SNRdB=SNRdBs(i_SNR);
: @+ x* h# w! b2 ]2 a  ^        sigma=sqrt(0.5/(10^(SNRdB/10)));%瑞利分布参数，噪声标准差
& N. P. _% G* U        for i_packet=1:N_packet
9 D0 z6 Z9 O3 t& M$ ^( q9 U% g9 K: X5 d/ P            symbol_data=randint(L_frame*b,NT); %产生L_frame*b行NT列二进制信息
1 F( X6 z+ q6 |1 w            [temp,sym_tab,P]=modulator(symbol_data.',b);%调制
            X=temp.';%X为调制后序列（L_frame行1列）
            Hr=(randn(L_frame,NR)+j*randn(L_frame,NR))/sq2;%瑞利信道衰落系数
+ J- }4 c6 Z% X5 s9 ~% v            H=reshape(Hr,L_frame,NR);%重置Hr行列（L_frame行NR列）
5 v3 O& v0 r! Y( }            Habs=sum(abs(H).^2,2);
, q' r7 y4 S3 ^# a            Z=0;
! M0 R% Y/ U% h! v4 a            for i=1:NR
                R(:,i)=sum(H(:,i).*X,2)/sq_NT+...
                    sigma*(randn(L_frame,1)+j*randn(L_frame,1)); %经过瑞利信道并加复高斯白噪声
/ ^7 H# p) o+ U, [3 G6 c" W                Z=Z+R(:,i).*conj(H(:,i));%MRC合并，合并系数为conj(H(:,i))
1 l$ F" [# [2 q/ a! r7 H* O( U/ T) a$ d            end
: i/ M2 O) M3 t  u7 W            for m=1 %最大似然译码
) b5 @) t- {9 G7 x7 \/ d' |* T                dl(:,m)=abs(sum(Z,2)-sym_tab(m)).^2+...
& l2 K( l9 F- L2 ?! j6 I' Q' ~                    (-1+sum(Habs,2))*abs(sym_tab(m))^2;
3 g& N2 A: Q% `$ q! H            end
            [y1,i1]=min(dl,[],2);
1 N- j/ c' K/ V0 k0 [; \1 n            Xd=sym_tab(i1).';
            temp1=X>0;
            temp2=Xd>0;
' [! J  i1 {1 O  Q. D3 h+ n- L: b            noeb_p(i_packet)=sum(sum(temp1~=temp2));%累加错误比特
        end
$ [5 \8 G$ Z& q        BER(iter,i_SNR)=sum(noeb_p)/(N_packet*L_frame*b);%计算误比特率
; C: v8 ]  c1 B- o) u0 N    end
    semilogy(SNRdBs,BER(iter,,gs),hold on
    axis([SNRdBs([1 end]) 0.5*10^(-4) 1e0])
7 D  d' c4 j; {. s: |" a% y: K& Vend
8 i& {  u/ t' R% e$ t* xtitle('BER peRFoemancde of MRC Scheme')
xlabel('SNR[dB]'),ylabel('BER')
grid on,set(gca,'fontsize',9)
legend('SISO','MRC(Tx:1,Rx:2)','MRC(Tx:1,Rx:4)')
8 K4 M; T! u- m& }7 ?# V* ?

wu68aq

这个其实很简单。