graph.m

数字通信调制解调

% Program 6-5
% graph.m
%
% The function of drawing the graph of simulation result.
%
% Input argument
%   filename : name of the file which simulation result was stored.
%
% Output argument
%   nothing
%
% Programmed by M.Okita
% Checked by H.Harada
%

function graph(filename)

mtitle1 = {'Throughput of Pure ALOHA system'        ...     % definition of the title
           'Throughput of Slotted ALOHA system'     ...
           'Throughput of np CSMA system'           ...
           'Throughput of Slotted np ISMA system'   ...
          };
mtitle2 = {'Average Delay time of Pure ALOHA system'    ... % definition of the title
           'Average Delay time of Slotted ALOHA system' ...
           'Average Delay time of np CSMA system'       ...
           'Average Delay time of Slotted ISMA system'  ...
          };

fid      = fopen(filename,'r');
nul      = fscanf(fid,'%s',2);
protocol = fscanf(fid,'%d',1);
nul      = fscanf(fid,'%s',2);
caputre  = fscanf(fid,'%d',1);
nul      = fscanf(fid,'%s',2);
dtime    = fscanf(fid,'%f',1);
nul      = fscanf(fid,'%s',2);
brate    = fscanf(fid,'%d',1);
nul      = fscanf(fid,'%s',2);
srate    = fscanf(fid,'%d',1);
nul      = fscanf(fid,'%s',2);
plen     = fscanf(fid,'%d',1);
nul      = fscanf(fid,'%s',1);
mnum     = fscanf(fid,'%d',1);
nul      = fscanf(fid,'%s',2);
mcn      = fscanf(fid,'%f',1);
nul      = fscanf(fid,'%s',2);
tcn      = fscanf(fid,'%f',1);
nul      = fscanf(fid,'%s',2);
spend    = fscanf(fid,'%d',1);

idx = 0;
while 1
    data0 = fscanf(fid,'%f',1);
    data1 = fscanf(fid,'%f',1);
    data2 = fscanf(fid,'%f',1);
    data3 = fscanf(fid,'%f',1);
    data4 = fscanf(fid,'%f',1);
    data5 = fscanf(fid,'%f',1);
    data6 = fscanf(fid,'%f',1);
    
    if feof(fid) == 1                                           % eof
        break
    end
    
    idx = idx + 1;
    
    tg(idx) = data0;                                            % offered traffic
    g(idx)  = data3 / data1;                                    % actual offered traffic
    s(idx)  = data2 / data1;                                    % actual throughput
    w(idx)  = data4 / spend * srate / plen;                     % average delay
end

fclose(fid);

ts = theorys(protocol,g,dtime);                                 % calculation of the throughput

if protocol < 3                                                 % Pure ALOHA & Slotted ALOHA
    plot(g,s,'bo:',g,ts,'r-');                                  % normal graph
    legend('result','theory',0);                                % legend
else                                                            % np-CSMA & Slotted np-ISMA
    semilogx(g,s,'bo:',g,ts,'r-');                              % semi log graph
    if protocol == 3                                            % legend
        legend(strcat('result (a=',num2str(dtime),')'),'theory',0);
    else
        legend(strcat('result (d=',num2str(dtime),')'),'theory',0);
    end
end
title(char(mtitle1(protocol)),'FontSize',16);                   % title
xlabel('Traffic(Simulation result)','FontSize',14);             % x axis label
ylabel('Throughput','FontSize',14);                             % y axis label

figure(2);                                                      % preparation of the new graph windos
if protocol < 3                                                 % Pure ALOHA & Slotted ALOHA
    plot(g,w,'bo:');                                            % normal graph
    legend('result');                                           % legend
else                                                            % np-CSMA & Slotted np-ISMA
    semilogx(g,w,'bo:');                                        % semi log graph
    if protocol == 3                                            % legend
        legend(strcat('result (a=',num2str(dtime),')'),0);
    else
        legend(strcat('result (d=',num2str(dtime),')'),0);
    end
end
title(char(mtitle2(protocol)),'FontSize',16);                   % title
xlabel('Traffic(Simulation result)','FontSize',14);             % x axis label
ylabel('Average Delay time(packet)','FontSize',14);             % y axis label
 
%%%%%%%%%%%%%%%%%%%%%% end of file %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%