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

shyw559

本人菜鸟，第一次接触matlab程序画流程图，不知道该如何按照这个程序画流程图， 尤其是这个程序里用了很多for循环， 比较迷茫……麻烦大神指点一下，非常感谢！
%MRC_scheme.m
%接收分集-MRC
+ N- R/ v' O$ N! eclear,clf
3 @2 F' q" t5 T' pL_frame=130;
N_packet=4000;%130个编码后符号，仿真4000组
3 U1 B* d7 [" t! g# t! A: {b=4;%设置为1/2/3/4对应于BPSK/QPSK/8PSK/16-QAM
9 D( }# T# X; C0 B; ]2 A) SSNRdBs=[0:2:20];
sq2=sqrt(2);
for iter=1:3
' m4 q0 O7 U) ?% a7 G' q, Q    if iter==1
% p8 I: [" V. a1 P) M$ K        NT=1;0
        NR=1;%SISO
        gs='-kx';
- o( }- C" k/ h6 T+ U+ {6 _    elseif iter==2
        NT=1;
        NR=2;% 1发2收
6 d( N; n" v8 n8 O) l0 q        gs='-^';
    else
        NT=1;
        NR=4;% 1发4收
        gs='-ro';
    end
    sq_NT=sqrt(NT);
    for i_SNR=1:length(SNRdBs)
$ F8 b8 x/ A' f1 y( s        SNRdB=SNRdBs(i_SNR);
+ U) I6 u6 z! L        sigma=sqrt(0.5/(10^(SNRdB/10)));%瑞利分布参数，噪声标准差
        for i_packet=1:N_packet
            symbol_data=randint(L_frame*b,NT); %产生L_frame*b行NT列二进制信息
7 O* \: V/ x  q( t2 t1 m1 s            [temp,sym_tab,P]=modulator(symbol_data.',b);%调制
            X=temp.';%X为调制后序列（L_frame行1列）
- _8 j* D: h  |) \! {5 p            Hr=(randn(L_frame,NR)+j*randn(L_frame,NR))/sq2;%瑞利信道衰落系数
            H=reshape(Hr,L_frame,NR);%重置Hr行列（L_frame行NR列）
            Habs=sum(abs(H).^2,2);
+ p, [; i; b* L: H# O9 H! n" W8 X            Z=0;
            for i=1:NR
; S$ g; r$ Q" A& `) B9 @                R(:,i)=sum(H(:,i).*X,2)/sq_NT+...
                    sigma*(randn(L_frame,1)+j*randn(L_frame,1)); %经过瑞利信道并加复高斯白噪声
                Z=Z+R(:,i).*conj(H(:,i));%MRC合并，合并系数为conj(H(:,i))
            end
6 {/ J; M0 b/ a1 l1 v  Q            for m=1 %最大似然译码
! ^# q4 C' g1 o+ ?( a! @2 l                dl(:,m)=abs(sum(Z,2)-sym_tab(m)).^2+...
4 G7 D2 [1 j- ^* F& O6 E                    (-1+sum(Habs,2))*abs(sym_tab(m))^2;
            end
            [y1,i1]=min(dl,[],2);
            Xd=sym_tab(i1).';
            temp1=X>0;
            temp2=Xd>0;
            noeb_p(i_packet)=sum(sum(temp1~=temp2));%累加错误比特
        end
6 q: R$ [2 F8 h" {+ y  i        BER(iter,i_SNR)=sum(noeb_p)/(N_packet*L_frame*b);%计算误比特率
    end
- H! c/ D6 b/ t# q7 x- O% M* s: d    semilogy(SNRdBs,BER(iter,,gs),hold on
& k3 ~1 G  ~5 W( }    axis([SNRdBs([1 end]) 0.5*10^(-4) 1e0])
2 p  x) F1 ?  @* Oend
title('BER peRFoemancde of MRC Scheme')
xlabel('SNR[dB]'),ylabel('BER')
grid on,set(gca,'fontsize',9)
4 {" X- ?9 I% I! K* S9 X! Hlegend('SISO','MRC(Tx:1,Rx:2)','MRC(Tx:1,Rx:4)')
( T  ]/ Y+ y) n4 ^

wu68aq

这个其实很简单。