【MIMO无线编码】在瑞利衰落信道中比较2x2和4x4的BPSK调制系统（Matlab代码实现）

最新推荐文章于 2024-09-13 23:25:03 发布

然哥依旧

最新推荐文章于 2024-09-13 23:25:03 发布

阅读量422

点赞数 5

文章标签： matlab 算法网络

本文链接：https://blog.csdn.net/m0_64583023/article/details/139208677

版权

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

⛳️赠与读者

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

在瑞利衰落信道中比较2x2和4x4的BPSK调制系统，在MIMO无线编码系统中，对于更高的天线系统，BER方面编码的选择并不重要。在这里，我们对2x2和4x4的BPSK调制系统进行了比较研究，针对线性复域编码A和B（即LCFA和LCFB）在瑞利衰落信道中的比特误码概率。

我们发现，在2x2天线系统中，LCFA和LCFB在比特误码概率方面确实很重要。在ML（最大似然）解码的情况下，LCFB比LCFA结果要好得多。但当我们转向4x4天线系统时，我们发现这并不重要。因为随着系统接收器变得更高，对于决定原始信息，接收器会得到更多的编码信号，而不像低阶系统那样少。

📚2 运行结果

部分代码：

N = 10^4; % Block size
theta=randn(1);
snr = 1:2:15; % SNR values

bit_error=zeros(1,length(snr));
bit_error1=zeros(1,length(snr));
bit_error2=zeros(1,length(snr));
bit_error3=zeros(1,length(snr));

%%%lcfA

MT = 4; % Number of transmitters
MR = 4;

phi=(1/2)*[1 exp((-1i*pi)/8) exp((-1i*2*pi)/8) exp((-1i*3*pi)/8);1 exp((-1i*5*pi)/4) exp((-1i*10*pi)/8) exp((-1i*15*pi)/8);1 exp((-1i*9*pi)/8) exp((-1i*18*pi)/8) exp((-1i*27*pi)/8);1 exp((-1i*13*pi)/8) exp((-1i*26*pi)/8) exp((-1i*39*pi)/8)];
for i = 1:length(snr)
sum1=0;
for in=1:N
p = rand(4,1)>0.5;
s = 2*p-1 ;
x=phi*s;
lcf=[x(1,1) 0 0 0; 0 x(2,1) 0 0;0 0 x(3,1) 0; 0 0 0 x(4,1)];
r1=normrnd(1*cos(theta),1,[MT,MR]);
r2=normrnd(1*sin(theta),1,[MT,MR]);
h = 1/sqrt(2)*sqrt(r1^2+r2^2);

r1=normrnd(1*cos(theta),1,[MT,MR]);
r2=normrnd(1*sin(theta),1,[MT,MR]);
n = 1/sqrt(2)*sqrt(r1^2+r2^2);

y=h*lcf ;

y = y + 10^(-snr(i)/20)*n ;

%decoding
sHat0 = [-1 ; -1; -1 ;-1];
x=phi*sHat0;
lcf=[x(1,1) 0 0 0; 0 x(2,1) 0 0;0 0 x(3,1) 0; 0 0 0 x(4,1)];
y0000=h*lcf;
j0000=abs(y - y0000);
j0000=sum(sum(j0000,2),1);

sHat1=[-1 ;-1;-1 ; 1];
x=phi*sHat1;
lcf=[x(1,1) 0 0 0; 0 x(2,1) 0 0;0 0 x(3,1) 0; 0 0 0 x(4,1)];
y0001=h*lcf;
j0001=sum(sum(abs(y-y0001),2),1);