svd_decompose_channel.m

本程序利用改进的CHOW算法

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  L为信道有效时间
%   N - # subcarriers子载波数

U = [];
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];
end

S = sum(S,1);