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

shyw559

本人菜鸟，第一次接触matlab程序画流程图，不知道该如何按照这个程序画流程图， 尤其是这个程序里用了很多for循环， 比较迷茫……麻烦大神指点一下，非常感谢！
+ ^# t" i% F- s: U! Z%MRC_scheme.m
%接收分集-MRC
clear,clf
. F# P7 P" s+ \) ^& DL_frame=130;
N_packet=4000;%130个编码后符号，仿真4000组
9 d! i+ }1 p  ]% M* F8 X$ L. sb=4;%设置为1/2/3/4对应于BPSK/QPSK/8PSK/16-QAM
SNRdBs=[0:2:20];
+ ^$ A3 O+ I3 m5 r7 Jsq2=sqrt(2);
5 B. Y* G/ }" p% afor iter=1:3
    if iter==1
        NT=1;0
        NR=1;%SISO
) l% d! x8 L1 x4 `$ x% J        gs='-kx';
4 q" S, f) T+ o) ]" G5 Y  y" ~    elseif iter==2
! s+ D+ n3 |* `3 o! u$ C2 k! p        NT=1;
        NR=2;% 1发2收
3 m. s; T, C; J. Y- t3 F' J8 `        gs='-^';
    else
        NT=1;
        NR=4;% 1发4收
, Z7 B0 f( _9 @2 C3 n. l/ B        gs='-ro';
" Y- O/ z2 R# C' o) g) \    end
8 i, l5 z; P* _! h    sq_NT=sqrt(NT);
" ]/ _) N! ]: E. {    for i_SNR=1:length(SNRdBs)
1 B% n7 R* }6 V        SNRdB=SNRdBs(i_SNR);
8 u0 b! @9 y/ w" U0 Y        sigma=sqrt(0.5/(10^(SNRdB/10)));%瑞利分布参数，噪声标准差
: ~8 m4 l: b* m. e* e: U7 p0 H        for i_packet=1:N_packet
* v$ `$ O9 j5 y            symbol_data=randint(L_frame*b,NT); %产生L_frame*b行NT列二进制信息
            [temp,sym_tab,P]=modulator(symbol_data.',b);%调制
; P/ S( M% J$ F( H6 t            X=temp.';%X为调制后序列（L_frame行1列）
            Hr=(randn(L_frame,NR)+j*randn(L_frame,NR))/sq2;%瑞利信道衰落系数
            H=reshape(Hr,L_frame,NR);%重置Hr行列（L_frame行NR列）
* G7 j2 ~$ I- J3 m            Habs=sum(abs(H).^2,2);
4 J7 [5 W0 P2 J) m+ C            Z=0;
            for i=1:NR
# T+ ^  n. _5 a: {                R(:,i)=sum(H(:,i).*X,2)/sq_NT+...
2 T9 R" o, o7 O) X, }" t9 R                    sigma*(randn(L_frame,1)+j*randn(L_frame,1)); %经过瑞利信道并加复高斯白噪声
2 q" h% h1 n  F% z. s2 u2 [                Z=Z+R(:,i).*conj(H(:,i));%MRC合并，合并系数为conj(H(:,i))
) P- N/ C% s, M3 n8 v            end
            for m=1 %最大似然译码
# E4 g0 [- b+ k0 s0 Q2 {. Y                dl(:,m)=abs(sum(Z,2)-sym_tab(m)).^2+...
. ^1 C$ T. o& E3 R& b* x                    (-1+sum(Habs,2))*abs(sym_tab(m))^2;
$ u+ X& @3 @$ k- h6 @: Z5 g; e% e            end
            [y1,i1]=min(dl,[],2);
( \0 `' Z" t& i1 L; r/ z            Xd=sym_tab(i1).';
            temp1=X>0;
            temp2=Xd>0;
            noeb_p(i_packet)=sum(sum(temp1~=temp2));%累加错误比特
, s9 r* i. i4 @! N+ T        end
        BER(iter,i_SNR)=sum(noeb_p)/(N_packet*L_frame*b);%计算误比特率
    end
    semilogy(SNRdBs,BER(iter,,gs),hold on
2 |' b# S! }$ n0 {    axis([SNRdBs([1 end]) 0.5*10^(-4) 1e0])
/ _' h- n* z$ B) ]6 R' Bend
& k1 C- _5 n  D. h! |+ T/ ytitle('BER peRFoemancde of MRC Scheme')
6 U2 h2 g6 d; r6 |5 Wxlabel('SNR[dB]'),ylabel('BER')
grid on,set(gca,'fontsize',9)
# m: {( E+ d1 Dlegend('SISO','MRC(Tx:1,Rx:2)','MRC(Tx:1,Rx:4)')
8 a3 ]4 [* w- M6 J! I

wu68aq

这个其实很简单。