^(-(kk*l)); end end end %q=zeros(L*Nt,L*Nt); %for ii=1:Nt % for jj=1:Nt % for l1=0:L-1 % for l2=0:L-1 % for kk=0:carrier_count-1 % q(l2+(ii-1)*L+1,l1+(jj-1)*L+1)= q(l2+(ii-1)*L+1,l1+(jj-1)*L+1)+tx_training_symbols(kk+1,ii)*conj(tx_t raining_symbols(kk+1,jj))*Wk^(-(kk*(-l1+l2))); % end % end % end % end %end %h=inv(q)*p; h=p/carrier_count; H_buf=zeros(carrier_count,Nt); for ii=1:Nt for kk=0:carrier_count-1 for l=0:L-1 H_buf(kk+1,ii)=H_buf(kk+1,ii)+h(l+(ii-1)*L+1,1)*Wk^(kk*l); end end end H_buf=conj(H_buf'); RRR=[]; for kk=1:carrier_count Y=Y_buf(:,kk); H=H_buf(:,kk); for co_ii=1:num_X for co_tt=1:size(eta,2) if eta(co_ii,co_tt)~=0 if coj_mt(eta(co_ii,co_tt),co_ii)==0 r_til(eta(co_ii,co_tt),:,co_ii)=Y(eta(co_ii,co_tt),:); a_til(eta(co_ii,co_tt),:,co_ii)=conj(H(epsilon(eta(co_ii,co_tt),co_ii)