00001 function [a, b, sampling] = picoeff(N,R,C)
00002 % This function gives the coefficients "a" and "b" of the PI
00003 % controller at sampling frequency "sampling". N is the minimum
00004 % number of flows, R the maximum round trip delay and C the capacity
00005 % of the link in (average sized) packets/sec.
00006 %function [a, b, sampling] = picoeff(N,R,C)
00007 N = N*500;
00008 C = C*500;
00009 p_q = 1/R;
00010 w_g = 2*N/R/R/C;
00011 K = w_g*abs( (i*w_g/p_q+1)/((R*C)^3/(2*N)^2));
00012 num = K*[1/w_g 1];
00013 den = [1 0];
00014 [numn,denn] = bilinear(num,den,100*w_g);
00015 sampling = 100*w_g;
00016 format long
00017 a = numn(1);
00018 b = -numn(2);
