test_cc.m

% demo script using plot_cc
%   to display singular vectors.
% Version 1.1 Sept.28,2009AD, L.P.Withers, Jr./E531
% Copyright 2009 The MITRE Corporation. ALL RIGHTS RESERVED. 
% 
%Redistribution and use in source and binary forms, with or without
%modification, are permitted provided that the following conditions are
%met:
%
%* Redistributions of source code must retain the above copyright
%notice, this list of conditions and the following disclaimer.
%
%* Redistributions in binary form must reproduce the above copyright
%notice, this list of conditions and the following disclaimer in the
%documentation and/or other materials provided with the distribution.
%
%* The US Government will not be charged any license fee and/or
%royalties related to this software.
%
%* Neither name of The MITRE Corporation; nor the names of its
%contributors may be used to endorse or promote products derived from
%this software without specific prior written permission.
%
%THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
%``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
%LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
%A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT
%HOLDERS OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
%INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
%BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS
%OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
%ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR
%TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE
%USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH
%DAMAGE.
%-------------------------------------------------------------------------

n=8;
R=(1/sqrt(2))*(randn(n)+i*randn(n)); %as for frequency-flat Rayleigh fading channels
dots=['r';'g';'b';'m';'c';'y'];
for k=1:2;
 [U,D,V]=svd(R);
 plot_cc(50*U(:,1),1,dots(mod(k-1,6)+1),1,1,1, 2.79,'time (us)',[],1,[]);
 drawnow;
 hold on;
 R= .9*R+.1*R';
end;

hold off;