ISDN网络模型

clear all;
clc;
s1=0;
for u=1:20
    for t=1:20
        ton([u],[t])=0.35*6*randint(1,1,[0,1]);
        
    end;
end;
for u=1:20
    for t=1:20
        s1=ton([u],[t])+s1;
        totalonvoice(t)=s1;
    end;
    s1=0;
end;
[d]=1000*[0 128 512 1024];
s2=0;
for u=1:20
    for t=1:20
        ton([u],[t])=0.35*randint(1,1,[0,1]);
    end;
end;
for u=1:20
    for t=1:20
        if((ton([u],[t])>=0)&(ton([u],[t])<=2))
            tond([u],[t])=ton([u],[t])*d(1);
        end;
        if((ton([u],[t])>=3)&(ton([u],[t])<=5))
            tond([u],[t])=ton([u],[t])*d(2);
        end;
        if((ton([u],[t])>=6)&(ton([u],[t])<=8))
            tond([u],[t])=ton([u],[t])*d(3);
        end;
        if((ton([u],[t])>=9)&(ton([u],[t])<=10))
            tond([u],[t])=ton([u],[t])*d(4);
        end;
    end;
end;
for t=1:20
    for u=1:20
        s2=(tond([u],[t]))+s2;
        totalondata(t)=s2;
    end;
    s2=0;
end;
[x]=randint(1,1,[1,5]);
s3=0;
a=20;
[p]=[0 0 0.5 0.5; 0.117 0 0.33 0.5; 0.4 0.2 0 0.4; 2 0.3 0.5 0];
ts1=exppdf(x,7);
ts2=exppdf(x,1);
ts3=exppdf(x,3);
ts4=exppdf(x,4);
[d]=[3.75 2.9 3.4 6]*power(10,6);
for u=1:20
    for t=1:5
        p1([u],[t])=ts1*p(1,3)*d(1)*a;
    end;
    for t=6:10
        p1([u],[t])=ts2*p(2,4)*d(2)*a;
    end;
    for t=11:15
        p1([u],[t])=ts3*p(3,2)*d(3)*a;
    end;
    for t=16:20
        p1([u],[t])=ts4*p(4,2)*d(4)*a;
    end;
end;
s5=0;
for u=1:20
    for t=1:20
        s3=p1([u],[t])+s3;
        totalonvideo(t)=s3;
    end;
    s3=0;
end;
buffer(1)=0;
for t=1:20
    totaltraffic(t)=0.05*totalondata(t)+.15*totalonvoice(t)+0.8*totalonvideo(t);
    buffer(t+1)=totalondata(t)-0.05*totalondata(t)+totalonvoice(t)-0.15*totalonvoice(t)+totalonvideo(t)-0.8*totalonvideo(t);
    totaltraffic(t+1)=totaltraffic(t)+buffer(t+1);
    if totalonvideo(t)<23.25*10^6
        totaltraffic(t)=totalondata(t)+totalonvoice(t)+0.8+totalonvideo(t)+buffer(t+1);
    end;
end;
bar(totaltraffic);
totaltraffic
for i=1:20
    if(totaltraffic(i)<=155*power(10,6))
        rem(i)=0;
    else
        rem(i)=totaltraffic(i)-(155*power(10,6))
    end;
end;
rem
xlabel('time');
ylabel('traffic');
title('ISDN traffic model');