matdisp.m

function matdisp(R,R_bc,t,lambdap_H11,deg_lambda_H11,lambda_H,deg_lambda_H1,deg_lambda_H2);

% Initialization

lambdap_H22 = lambda_H(length(deg_lambda_H1)+1:length(deg_lambda_H1)+length(deg_lambda_H2))./(deg_lambda_H2');lambdap_H22 = lambdap_H22/sum(lambdap_H22);

lambdap_H11 = lambdap_H11/sum(lambdap_H11);

cs_lambdap_H11 = cumsum(lambdap_H11);
lambdap_H1 = (lambda_H(1:length(deg_lambda_H1))./deg_lambda_H1');lambdap_H1 = lambdap_H1/sum(lambdap_H1);
cs_lambdap_H1 = cumsum(lambdap_H1);
cs_lambdap_H21 = sort(unique([cs_lambdap_H1;cs_lambdap_H11]));
lambdap_H21 = [cs_lambdap_H21(1);cs_lambdap_H21(2:end)-cs_lambdap_H21(1:end-1)];
deg_lambda_H21 = zeros(size(lambdap_H21));for n = 1:length(deg_lambda_H21);deg_lambda_H21(n) = deg_lambda_H1(min(find(cs_lambdap_H1>=cs_lambdap_H21(n)-(1e-6)))) - deg_lambda_H11(min(find(cs_lambdap_H11>=cs_lambdap_H21(n)-(1e-6))));end;
deg_lambda_H2 = [reshape(deg_lambda_H21,1,length(deg_lambda_H21)),reshape(deg_lambda_H2,1,length(deg_lambda_H2))];
lambdap_H2 = [t*reshape(lambdap_H21,1,length(lambdap_H21)),(1-t)*reshape(lambdap_H22,1,length(lambdap_H22))];
deg_lambda_max = max([max(deg_lambda_H1),max(deg_lambda_H2)]);

% display

N = 500;M = round((1-R)*N);Ma = round(t*(1-R_bc)*N);mat_disp = zeros(M,N);
vNa = [1;round(cumsum(N*t*lambdap_H11)/sum(lambdap_H11))];
vNb = [1,round(cumsum(N*lambdap_H2)/sum(lambdap_H2))];
vM = [round((1-R_bc)*t*N);M];
for n = 1:length(vNa)-1;
	mat_disp(1:vM(1),vNa(n):vNa(n+1)) = deg_lambda_H11(n);
end;
for n = 1:length(vNb)-1;
	mat_disp(vM(1)+1:vM(2),vNb(n):vNb(n+1)) = deg_lambda_H2(n);
end;
%mat_disp = mat_disp/max(max(mat_disp));
figure;imagesc(linspace(0,1,N),linspace(0,1-R,M),mat_disp);
h=line([0,1],[t*(1-R_bc),t*(1-R_bc)]);set(h,'LineWidth',3,'Color',[0 0 0]);h=line([t,t],[0,t*(1-R_bc)]);set(h,'LineWidth',3,'Color',[0 0 0]);
%h=line([0,1],[0,1-R]);set(h,'LineWidth',2,'Color',[0 0 0]);
h=text(t/2,t*(1-R_bc)/2,'H_{11}');set(h,'FontSize',20);h=text(1/2,t*(1-R_bc)/2+(1-R)/2,'H_{2}');set(h,'FontSize',20);h=text((t+1)/2,t*(1-R_bc)/2,'0');set(h,'FontSize',30);
colormap('jet');colorbar;

% To compare with real matrix H
%figure;imagesc(H);colormap('gray');