rate.m

这是一个超宽带的源码

%
%	r = rate(H,x)
%
%	Calculates rates of links. 
%
%	H(i,j)  - fading between source of link <i> and dst of link <j>
%  x       - vector [a1;..;a(N+1);p11;..;pN1;..;P1(N+1)1;..;PN(N+1)];
%
%  r(i)    - rate at link <i> calculated as
%            r = sum aj * fr(snr(pj))


function xx=rate(H,x)

N = size(H,1);

xx = zeros(N,N+1);
for i=1:N+1
    xx(:,i) = fr(snr(x(N+1+(i-1)*N+1:N+1+i*N),H));
end

xx = xx*x(1:N+1);