C(f,:)=diag(S);
0042         VGxxf=Gxxf*U1;                  <span class="comment">% Virtual cross spectrum</span>
0043         Gxxf=real(diag(Gxxf));          <span class="comment">% Reduce to autospectra</span>
0044         <span class="keyword">for</span> x_sig=1:R
0045             <span class="keyword">for</span> pc_sig=1:R
0046                 VCxx(f,x_sig,pc_sig)=abs(VGxxf(x_sig,pc_sig))^2/Gxxf(x_sig)/PC(f,pc_sig);
0047             <span class="keyword">end</span>
0048         <span class="keyword">en