% pltSetB.m % 11 Sept 2004 function pltSetB(vrt, TB, ttl) for k = 1:4 Trb = TB(:,k) ; Tr = reshape((Trb+1)/2, 29, 3) ; xx = pow2(6:-1:-1)*Tr(1:8,:) ; yy = pow2(6:-1:-1)*Tr((1:8)+8,:) ; zt = xx+j*yy ; fi = pow2(6:-1:0)*Tr((1:7)+2*8,:) ; rot = exp(j*pi*(90-fi)/180) ; r = [2 1]*Tr((1:2)+23,:) ; g = [2 1]*Tr((1:2)+25,:) ; b = [2 1]*Tr((1:2)+27,:) ; cc = [r ; g ; b]/3 ; % rotation and translation of three tiangles vrt z = vrt.*rot([1 1 1],:) + zt([1 1 1],:) ; X = real(z) ; Y = imag(z) ; subplot(2,2,k) fill(X(:,1), Y(:,1), cc(1,:), ... X(:,2), Y(:,2), cc(2,:), ... X(:,3), Y(:,3), cc(3,:)) axis([0 80 0 80]), axis equal, grid on if k ==1, title(ttl), end end