limited_grassman_full.m

线形内插在MIMO OFDM上的实现,包括格拉斯曼域的仿真

function [ precoded ] = Limited_grassman_full( feeds )
%==========================================================================
% data type definition
index_feedback=zeros(1,1024);
H=zeros(1,4,1024);
Grassman_precoder=zeros(4,2,1024);
G2_stbc=zeros(2,2,1024);
G2_pre=zeros(4,2,1024);
G2_pre_11=zeros(1024,1);
G2_pre_12=zeros(1024,1);
G2_pre_21=zeros(1024,1);
G2_pre_22=zeros(1024,1);
G2_pre_31=zeros(1024,1);
G2_pre_32=zeros(1024,1);
G2_pre_41=zeros(1024,1);
G2_pre_42=zeros(1024,1);
H_pre_11=zeros(1024,1);
H_pre_12=zeros(1024,1);
H_pre=zeros(1,2,1024);
%=================================================================
Precoder(:,:,1) =   [0.0000 + 0.5000i  -0.5000 + 0.0000i;
  -0.5000 + 0.0000i   0.5000 - 0.0000i;
  -0.0000 - 0.5000i  -0.5000 + 0.0000i;
   0.5000 - 0.0000i   0.5000 - 0.0000i];
Precoder(:,:,2) =  [-0.4016 + 0.2978i  -0.2978 - 0.4016i;
   0.3536 - 0.3536i  -0.3536 + 0.3536i;
   0.4466 + 0.2248i   0.2248 - 0.4466i;
  -0.4998 - 0.0123i  -0.4998 - 0.0123i];
Precoder(:,:,3) =  [-0.4785 - 0.1451i   0.1451 - 0.4785i;
   0.0000 + 0.5000i  -0.0000 - 0.5000i;
  -0.4016 + 0.2978i   0.2978 + 0.4016i;
   0.4994 + 0.0245i   0.4994 + 0.0245i];
Precoder(:,:,4) =  [-0.1684 - 0.4708i   0.4708 - 0.1684i;
  -0.3536 - 0.3536i   0.3536 + 0.3536i;
  -0.0855 - 0.4926i  -0.4926 + 0.0855i;
  -0.4986 - 0.0368i  -0.4986 - 0.0368i];
Precoder(:,:,5) =   [0.2778 - 0.4157i   0.4157 + 0.2778i;
   0.5000 - 0.0000i  -0.5000 + 0.0000i;
   0.4785 + 0.1451i   0.1451 - 0.4785i;
   0.4976 + 0.0490i   0.4976 + 0.0490i];
Precoder(:,:,6) =   [0.4994 - 0.0245i   0.0245 + 0.4994i;
  -0.3536 + 0.3536i   0.3536 - 0.3536i;
  -0.3448 + 0.3621i   0.3621 + 0.3448i;
  -0.4962 - 0.0612i  -0.4962 - 0.0612i];
Precoder(:,:,7) =   [0.3172 + 0.3865i  -0.3865 + 0.3172i;
  -0.0000 - 0.5000i   0.0000 + 0.5000i;
  -0.1684 - 0.4708i  -0.4708 + 0.1684i;
   0.4946 + 0.0734i   0.4946 + 0.0734i];
Precoder(:,:,8) =  [-0.1215 + 0.4850i  -0.4850 - 0.1215i;
   0.3536 + 0.3536i  -0.3536 - 0.3536i;
   0.4962 + 0.0612i   0.0612 - 0.4962i;
  -0.4926 - 0.0855i  -0.4926 - 0.0855i];
Precoder(:,:,9) =  [-0.4619 + 0.1913i  -0.1913 - 0.4619i;
  -0.5000 + 0.0000i   0.5000 - 0.0000i;
  -0.2778 + 0.4157i   0.4157 + 0.2778i;
   0.4904 + 0.0975i   0.4904 + 0.0975i];
Precoder(:,:,10) =  [-0.4289 - 0.2571i   0.2571 - 0.4289i;
   0.3536 - 0.3536i  -0.3536 + 0.3536i;
  -0.2464 - 0.4350i  -0.4350 + 0.2464i;
  -0.4879 - 0.1096i  -0.4879 - 0.1096i];
Precoder(:,:,11) =  [-0.0490 - 0.4976i   0.4976 - 0.0490i;
   0.0000 + 0.5000i  -0.0000 - 0.5000i;
   0.4994 - 0.0245i  -0.0245 - 0.4994i;
   0.4850 + 0.1215i   0.4850 + 0.1215i];
Precoder(:,:,12) =  [ 0.3705 - 0.3358i   0.3358 + 0.3705i;
  -0.3536 - 0.3536i   0.3536 + 0.3536i;
  -0.2026 + 0.4571i   0.4571 + 0.2026i;
  -0.4819 - 0.1334i  -0.4819 - 0.1334i];
Precoder(:,:,13) =   [0.4904 + 0.0975i  -0.0975 + 0.4904i;
   0.5000 - 0.0000i  -0.5000 + 0.0000i;
  -0.3172 - 0.3865i  -0.3865 + 0.3172i;
   0.4785 + 0.1451i   0.4785 + 0.1451i];
Precoder(:,:,14) =   [0.2138 + 0.4520i  -0.4520 + 0.2138i;
  -0.3536 + 0.3536i   0.3536 - 0.3536i;
   0.4879 - 0.1096i  -0.1096 - 0.4879i;
  -0.4748 - 0.1568i  -0.4748 - 0.1568i];
Precoder(:,:,15) =  [-0.2357 + 0.4410i  -0.4410 - 0.2357i;
  -0.0000 - 0.5000i   0.0000 + 0.5000i;
  -0.1215 + 0.4850i   0.4850 + 0.1215i;
   0.4708 + 0.1684i   0.4708 + 0.1684i];
Precoder(:,:,16) =  [-0.4946 + 0.0734i  -0.0734 - 0.4946i;
   0.3536 + 0.3536i  -0.3536 - 0.3536i;
  -0.3786 - 0.3266i  -0.3266 + 0.3786i;
  -0.4665 - 0.1799i  -0.4665 - 0.1799i];
Precoder(:,:,17) =  [-0.3536 - 0.3536i   0.3536 - 0.3536i;
  -0.5000 + 0.0000i   0.5000 - 0.0000i;
   0.4619 - 0.1913i  -0.1913 - 0.4619i;
   0.4619 + 0.1913i   0.4619 + 0.1913i];
Precoder(:,:,18) =  [ 0.0734 - 0.4946i   0.4946 + 0.0734i;
   0.3536 - 0.3536i  -0.3536 + 0.3536i;
  -0.0368 + 0.4986i   0.4986 + 0.0368i;
  -0.4571 - 0.2026i  -0.4571 - 0.2026i];
Precoder(:,:,19) =  [ 0.4410 - 0.2357i   0.2357 + 0.4410i;
   0.0000 + 0.5000i  -0.0000 - 0.5000i;
  -0.4289 - 0.2571i  -0.2571 + 0.4289i;
   0.4520 + 0.2138i   0.4520 + 0.2138i];
Precoder(:,:,20) =  [ 0.4520 + 0.2138i  -0.2138 + 0.4520i;
  -0.3536 - 0.3536i   0.3536 + 0.3536i;
   0.4224 - 0.2675i  -0.2675 - 0.4224i;
  -0.4466 - 0.2248i  -0.4466 - 0.2248i];
Precoder(:,:,21) =   [0.0975 + 0.4904i  -0.4904 + 0.0975i;
   0.5000 - 0.0000i  -0.5000 + 0.0000i;
   0.0490 + 0.4976i   0.4976 - 0.0490i;
   0.4410 + 0.2357i   0.4410 + 0.2357i];
Precoder(:,:,22) =  [-0.3358 + 0.3705i  -0.3705 - 0.3358i;
  -0.3536 + 0.3536i   0.3536 - 0.3536i;
  -0.4665 - 0.1799i  -0.1799 + 0.4665i;
  -0.4350 - 0.2464i  -0.4350 - 0.2464i];
Precoder(:,:,23) =  [-0.4976 - 0.0490i   0.0490 - 0.4976i;
  -0.0000 - 0.5000i   0.0000 + 0.5000i;
   0.3705 - 0.3358i  -0.3358 - 0.3705i;
   0.4289 + 0.2571i   0.4289 + 0.2571i];
Precoder(:,:,24) =  [-0.2571 - 0.4289i   0.4289 - 0.2571i;
   0.3536 + 0.3536i  -0.3536 - 0.3536i;
   0.1334 + 0.4819i   0.4819 - 0.1334i;
  -0.4224 - 0.2675i  -0.4224 - 0.2675i];
Precoder(:,:,25) =  [ 0.1913 - 0.4619i   0.4619 + 0.1913i;
  -0.5000 + 0.0000i   0.5000 - 0.0000i;
  -0.4904 - 0.0975i  -0.0975 + 0.4904i;
   0.4157 + 0.2778i   0.4157 + 0.2778i];
Precoder(:,:,26) =  [ 0.4850 - 0.1215i   0.1215 + 0.4850i;
   0.3536 - 0.3536i  -0.3536 + 0.3536i;
   0.3076 - 0.3942i  -0.3942 - 0.3076i;
  -0.4088 - 0.2879i  -0.4088 - 0.2879i];
Precoder(:,:,27) =   [0.3865 + 0.3172i  -0.3172 + 0.3865i;
   0.0000 + 0.5000i  -0.0000 - 0.5000i;
   0.2138 + 0.4520i   0.4520 - 0.2138i;
   0.4016 + 0.2978i   0.4016 + 0.2978i];
Precoder(:,:,28) =  [-0.0245 + 0.4994i  -0.4994 - 0.0245i;
  -0.3536 - 0.3536i   0.3536 + 0.3536i;
  -0.4998 - 0.0123i  -0.0123 + 0.4998i;
  -0.3942 - 0.3076i  -0.3942 - 0.3076i];
Precoder(:,:,29) =  [-0.4157 + 0.2778i  -0.2778 - 0.4157i;
   0.5000 - 0.0000i  -0.5000 + 0.0000i;
   0.2357 - 0.4410i  -0.4410 - 0.2357i;
   0.3865 + 0.3172i   0.3865 + 0.3172i];
Precoder(:,:,30) =  [-0.4708 - 0.1684i   0.1684 - 0.4708i;
  -0.3536 + 0.3536i   0.3536 - 0.3536i;
   0.2879 + 0.4088i   0.4088 - 0.2879i;
  -0.3786 - 0.3266i  -0.3786 - 0.3266i];
Precoder(:,:,31) =  [-0.1451 - 0.4785i   0.4785 - 0.1451i;
  -0.0000 - 0.5000i   0.0000 + 0.5000i;
  -0.4946 + 0.0734i   0.0734 + 0.4946i;
   0.3705 + 0.3358i   0.3705 + 0.3358i];
Precoder(:,:,32) =   [0.2978 - 0.4016i   0.4016 + 0.2978i;
   0.3536 + 0.3536i  -0.3536 - 0.3536i;
   0.1568 - 0.4748i  -0.4748 - 0.1568i;
  -0.3621 - 0.3448i  -0.3621 - 0.3448i];
Precoder(:,:,33) =   [0.5000 - 0.0000i   0.0000 + 0.5000i;
  -0.5000 + 0.0000i   0.5000 - 0.0000i;
   0.3536 + 0.3536i   0.3536 - 0.3536i;
   0.3536 + 0.3536i   0.3536 + 0.3536i];
Precoder(:,:,34) =   [0.2978 + 0.4016i  -0.4016 + 0.2978i;
   0.3536 - 0.3536i  -0.3536 + 0.3536i;
  -0.4748 + 0.1568i   0.1568 + 0.4748i;
  -0.3448 - 0.3621i  -0.3448 - 0.3621i];
Precoder(:,:,35) =  [-0.1451 + 0.4785i  -0.4785 - 0.1451i;
   0.0000 + 0.5000i  -0.0000 - 0.5000i;
   0.0734 - 0.4946i  -0.4946 - 0.0734i;
   0.3358 + 0.3705i   0.3358 + 0.3705i];
Precoder(:,:,36) =  [-0.4708 + 0.1684i  -0.1684 - 0.4708i;
  -0.3536 - 0.3536i   0.3536 + 0.3536i;
   0.4088 + 0.2879i   0.2879 - 0.4088i;
  -0.3266 - 0.3786i  -0.3266 - 0.3786i];
Precoder(:,:,37) =  [-0.4157 - 0.2778i   0.2778 - 0.4157i;
   0.5000 - 0.0000i  -0.5000 + 0.0000i;
  -0.4410 + 0.2357i   0.2357 + 0.4410i;
   0.3172 + 0.3865i   0.3172 + 0.3865i];
Precoder(:,:,38) =  [-0.0245 - 0.4994i   0.4994 - 0.0245i;
  -0.3536 + 0.3536i   0.3536 - 0.3536i;
  -0.0123 - 0.4998i  -0.4998 + 0.0123i;
  -0.3076 - 0.3942i  -0.3076 - 0.3942i];
Precoder(:,:,39) =  [ 0.3865 - 0.3172i   0.3172 + 0.3865i;
  -0.0000 - 0.5000i   0.0000 + 0.5000i;
   0.4520 + 0.2138i   0.2138 - 0.4520i;
   0.2978 + 0.4016i   0.2978 + 0.4016i];
Precoder(:,:,40) =  [ 0.4850 + 0.1215i  -0.1215 + 0.4850i;
   0.3536 + 0.3536i  -0.3536 - 0.3536i;
  -0.3942 + 0.3076i   0.3076 + 0.3942i;
  -0.2879 - 0.4088i  -0.2879 - 0.4088i];
Precoder(:,:,41) =  [ 0.1913 + 0.4619i  -0.4619 + 0.1913i;
  -0.5000 + 0.0000i   0.5000 - 0.0000i;
  -0.0975 - 0.4904i  -0.4904 + 0.0975i;
   0.2778 + 0.4157i   0.2778 + 0.4157i];
Precoder(:,:,42) = [ -0.2571 + 0.4289i  -0.4289 - 0.2571i;
   0.3536 - 0.3536i  -0.3536 + 0.3536i;
   0.4819 + 0.1334i   0.1334 - 0.4819i;
  -0.2675 - 0.4224i  -0.2675 - 0.4224i];
Precoder(:,:,43) =  [-0.4976 + 0.0490i  -0.0490 - 0.4976i;
   0.0000 + 0.5000i  -0.0000 - 0.5000i;
  -0.3358 + 0.3705i   0.3705 + 0.3358i;
   0.2571 + 0.4289i   0.2571 + 0.4289i];
Precoder(:,:,44) =  [-0.3358 - 0.3705i   0.3705 - 0.3358i;
  -0.3536 - 0.3536i   0.3536 + 0.3536i;
  -0.1799 - 0.4665i  -0.4665 + 0.1799i;
  -0.2464 - 0.4350i  -0.2464 - 0.4350i];
Precoder(:,:,45) =   [0.0975 - 0.4904i   0.4904 + 0.0975i;
   0.5000 - 0.0000i  -0.5000 + 0.0000i;
   0.4976 + 0.0490i   0.0490 - 0.4976i;
   0.2357 + 0.4410i   0.2357 + 0.4410i];


Precoder(:,:,46) =   [0.4520 - 0.2138i   0.2138 + 0.4520i;
  -0.3536 + 0.3536i   0.3536 - 0.3536i;
  -0.2675 + 0.4224i   0.4224 + 0.2675i;
  -0.2248 - 0.4466i  -0.2248 - 0.4466i];


Precoder(:,:,47) =  [ 0.4410 + 0.2357i  -0.2357 + 0.4410i;
  -0.0000 - 0.5000i   0.0000 + 0.5000i;
  -0.2571 - 0.4289i  -0.4289 + 0.2571i;
   0.2138 + 0.4520i   0.2138 + 0.4520i];


Precoder(:,:,48) =   [0.0734 + 0.4946i  -0.4946 + 0.0734i;
   0.3536 + 0.3536i  -0.3536 - 0.3536i;
   0.4986 - 0.0368i  -0.0368 - 0.4986i;
  -0.2026 - 0.4571i  -0.2026 - 0.4571i];


Precoder(:,:,49) =  [-0.3536 + 0.3536i  -0.3536 - 0.3536i;
  -0.5000 + 0.0000i   0.5000 - 0.0000i;
  -0.1913 + 0.4619i   0.4619 + 0.1913i;
   0.1913 + 0.4619i   0.1913 + 0.4619i];


Precoder(:,:,50) =  [-0.4946 - 0.0734i   0.0734 - 0.4946i;
   0.3536 - 0.3536i  -0.3536 + 0.3536i;
  -0.3266 - 0.3786i  -0.3786 + 0.3266i;
  -0.1799 - 0.4665i  -0.1799 - 0.4665i];


Precoder(:,:,51) =  [-0.2357 - 0.4410i   0.4410 - 0.2357i;
   0.0000 + 0.5000i  -0.0000 - 0.5000i;
   0.4850 - 0.1215i  -0.1215 - 0.4850i;
   0.1684 + 0.4708i   0.1684 + 0.4708i];


Precoder(:,:,52) =   [0.2138 - 0.4520i   0.4520 + 0.2138i;
  -0.3536 - 0.3536i   0.3536 + 0.3536i;
  -0.1096 + 0.4879i   0.4879 + 0.1096i;
  -0.1568 - 0.4748i  -0.1568 - 0.4748i];


Precoder(:,:,53) =   [0.4904 - 0.0975i   0.0975 + 0.4904i;
   0.5000 - 0.0000i  -0.5000 + 0.0000i;
  -0.3865 - 0.3172i  -0.3172 + 0.3865i;
   0.1451 + 0.4785i   0.1451 + 0.4785i];


Precoder(:,:,54) =   [0.3705 + 0.3358i  -0.3358 + 0.3705i;
  -0.3536 + 0.3536i   0.3536 - 0.3536i;
   0.4571 - 0.2026i  -0.2026 - 0.4571i;
  -0.1334 - 0.4819i  -0.1334 - 0.4819i];


Precoder(:,:,55) =  [-0.0490 + 0.4976i  -0.4976 - 0.0490i;
  -0.0000 - 0.5000i   0.0000 + 0.5000i;
  -0.0245 + 0.4994i   0.4994 + 0.0245i;
   0.1215 + 0.4850i   0.1215 + 0.4850i];


Precoder(:,:,56) =  [-0.4289 + 0.2571i  -0.2571 - 0.4289i;
   0.3536 + 0.3536i  -0.3536 - 0.3536i;
  -0.4350 - 0.2464i  -0.2464 + 0.4350i;
  -0.1096 - 0.4879i  -0.1096 - 0.4879i];


Precoder(:,:,57) =  [-0.4619 - 0.1913i   0.1913 - 0.4619i;
  -0.5000 + 0.0000i   0.5000 - 0.0000i;
   0.4157 - 0.2778i  -0.2778 - 0.4157i;
   0.0975 + 0.4904i   0.0975 + 0.4904i];


Precoder(:,:,58) =  [-0.1215 - 0.4850i   0.4850 - 0.1215i;
   0.3536 - 0.3536i  -0.3536 + 0.3536i;
   0.0612 + 0.4962i   0.4962 - 0.0612i;
  -0.0855 - 0.4926i  -0.0855 - 0.4926i];


Precoder(:,:,59) =   [0.3172 - 0.3865i   0.3865 + 0.3172i;
   0.0000 + 0.5000i  -0.0000 - 0.5000i;
  -0.4708 - 0.1684i  -0.1684 + 0.4708i;
   0.0734 + 0.4946i   0.0734 + 0.4946i];


Precoder(:,:,60) =   [0.4994 + 0.0245i  -0.0245 + 0.4994i;
  -0.3536 - 0.3536i   0.3536 + 0.3536i;
   0.3621 - 0.3448i  -0.3448 - 0.3621i;
  -0.0612 - 0.4962i  -0.0612 - 0.4962i];


Precoder(:,:,61) =   [0.2778 + 0.4157i  -0.4157 + 0.2778i;
   0.5000 - 0.0000i  -0.5000 + 0.0000i;
   0.1451 + 0.4785i   0.4785 - 0.1451i;
   0.0490 + 0.4976i   0.0490 + 0.4976i];


Precoder(:,:,62) =  [-0.1684 + 0.4708i  -0.4708 - 0.1684i;
  -0.3536 + 0.3536i   0.3536 - 0.3536i;
  -0.4926 - 0.0855i  -0.0855 + 0.4926i;
  -0.0368 - 0.4986i  -0.0368 - 0.4986i];


Precoder(:,:,63) =  [-0.4785 + 0.1451i  -0.1451 - 0.4785i;
  -0.0000 - 0.5000i   0.0000 + 0.5000i;
   0.2978 - 0.4016i  -0.4016 - 0.2978i;
   0.0245 + 0.4994i   0.0245 + 0.4994i];


Precoder(:,:,64) =  [-0.4016 - 0.2978i   0.2978 - 0.4016i;
   0.3536 + 0.3536i  -0.3536 - 0.3536i;
   0.2248 + 0.4466i   0.4466 - 0.2248i;
  -0.0123 - 0.4998i  -0.0123 - 0.4998i];


Precoder(:,:,65) =  [-0.0000 - 0.5000i   0.5000 - 0.0000i;
  -0.5000 + 0.0000i   0.5000 - 0.0000i;
  -0.5000 - 0.0000i  -0.0000 + 0.5000i;
   0.0000 + 0.5000i   0.0000 + 0.5000i];


Precoder(:,:,66) =   [0.4016 - 0.2978i   0.2978 + 0.4016i;
   0.3536 - 0.3536i  -0.3536 + 0.3536i;
   0.2248 - 0.4466i  -0.4466 - 0.2248i;
   0.0123 - 0.4998i   0.0123 - 0.4998i];


Precoder(:,:,67) =   [0.4785 + 0.1451i  -0.1451 + 0.4785i;
   0.0000 + 0.5000i  -0.0000 - 0.5000i;
   0.2978 + 0.4016i   0.4016 - 0.2978i;
  -0.0245 + 0.4994i  -0.0245 + 0.4994i];


Precoder(:,:,68) =   [0.1684 + 0.4708i  -0.4708 + 0.1684i;
  -0.3536 - 0.3536i   0.3536 + 0.3536i;
  -0.4926 + 0.0855i   0.0855 + 0.4926i;
   0.0368 - 0.4986i   0.0368 - 0.4986i];


Precoder(:,:,69) =  [-0.2778 + 0.4157i  -0.4157 - 0.2778i;
   0.5000 - 0.0000i  -0.5000 + 0.0000i;
   0.1451 - 0.4785i  -0.4785 - 0.1451i;