svd_decompose_channel.m

来自「mimoofdm系统自适应功率分配算法」· M 代码 · 共 24 行

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
% [U,S,V] = SVD(X);;;;;X = U*S*V'.( U,V is V'=inv(V),S is diag)
%
%   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 - # 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);