limited_grassman_full.m

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

Precoder(:,:,185) = [-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(:,:,186) =  [-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(:,:,187) =   [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(:,:,188) =   [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(:,:,189) =   [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(:,:,190) =  [-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(:,:,191) =  [-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(:,:,192) =  [-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(:,:,193) =  [-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(:,:,194) =   [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(:,:,195) =   [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(:,:,196) =   [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(:,:,197) =  [-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(:,:,198) =  [-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(:,:,199) =  [-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(:,:,200) =   [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(:,:,201) =   [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(:,:,202) =   [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(:,:,203) =   [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(:,:,204) =  [-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(:,:,205) =  [-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(:,:,206) =  [-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(:,:,207) =   [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(:,:,208) =  [ 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(:,:,209) =   [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(:,:,210) =  [-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(:,:,211) =  [-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(:,:,212) =  [-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(:,:,213) =  [-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(:,:,214) =  [ 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(:,:,215) =  [ 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(:,:,216) =  [ 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(:,:,217) = [ -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(:,:,218) = [ -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(:,:,219) =  [-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(:,:,220) =   [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(:,:,221) =   [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(:,:,222) =  [ 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(:,:,223) =  [ 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(:,:,224) =  [-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(:,:,225) =  [-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(:,:,226) =  [-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(:,:,227) =   [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(:,:,228) =   [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(:,:,229) =   [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(:,:,230) = [  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(:,:,231) =  [-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(:,:,232) = [ -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(:,:,233) =  [-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(:,:,234) =  [ 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(:,:,235) =  [ 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(:,:,236) =  [ 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(:,:,237) =  [-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(:,:,238) =  [-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(:,:,239) =  [-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(:,:,240) =  [-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(:,:,241) =  [ 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(:,:,242) =   [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(:,:,243) =  [ 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(:,:,244) = [ -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(:,:,245) =  [-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(:,:,246) =  [-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(:,:,247) = [  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(:,:,248) =  [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(:,:,249) =  [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(:,:,250) =  [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(:,:,251) =  [-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(:,:,252) = [-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(:,:,253) =  [-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(:,:,254) =  [ 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(:,:,255) =  [ 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(:,:,256) =   [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];
%=============================================
%==========================================================================
% data initializing 
for k=1:1024
    index_feedback(1,k)=feeds(k+1024*4);
end
for k=1:1024
    G2_stbc(:,:,k)=[feeds(k),feeds(k+1024);feeds(k+1024*2),feeds(k+1024*3)];
    H(:,:,k)=[feeds(k+1024*5),feeds(k+1024*6),feeds(k+1024*7),feeds(k+1024*8)];
end
%==========================================================================
% precoding prossecedure 
%index_feedback=ones(1,1024);
for k=1:1024
    Grassman_precoder(:,:,k)=Precoder(:,:,index_feedback(1,k));
    G2_pre(:,:,k)=Grassman_precoder(:,:,k)*G2_stbc(:,:,k);
    H_pre(:,:,k)=H(:,:,k)*Grassman_precoder(:,:,k);
end
%==========================================================================
% output 
for k=1:1024
    G2_pre_11(k,1)=G2_pre(1,1,k);
    G2_pre_12(k,1)=G2_pre(1,2,k);
    G2_pre_21(k,1)=G2_pre(2,1,k);
    G2_pre_22(k,1)=G2_pre(2,2,k);
    G2_pre_31(k,1)=G2_pre(3,1,k);
    G2_pre_32(k,1)=G2_pre(3,2,k);
    G2_pre_41(k,1)=G2_pre(4,1,k);
    G2_pre_42(k,1)=G2_pre(4,2,k);
    H_pre_11(k,1)=H_pre(1,1,k);
    H_pre_12(k,1)=H_pre(1,2,k);
end
precoded=[G2_pre_11,G2_pre_21,G2_pre_31,G2_pre_41,G2_pre_12,G2_pre_22,G2_pre_32,G2_pre_42,H_pre_11,H_pre_12];