function [U, S, V] = svd_decompose_channel(Mt, Mr, h_f, N);% [U S V] = svd_decompose_channel(Mt, Mr, h_f, N);%% Function decomposes the channel at each subcarrier into its SVD components%%   Mt - # Tx antennas%   Mr - # Rx antennas%   h_f - MIMO impulse response - Mr rows, Mt*L columns, where L is the number of%         channel taps%   N - # subcarriersU = [];S = [];V = [];for i = 1:N    [Utmp Stmp Vtmp] = svd(h_f(:,(i-1)*Mt+1:i*Mt));    U=[U Utmp];    V=[V Vtmp];    S=[S Stmp];endS = sum(S,1);