clc
clear all;
close all;
N = 10^5;
M=4;
x = randi([0 3],1,N);
xmod=qammod(x,4);
xmod=reshape(xmod,8,N/8);
xmod=kron(xmod,[1,1,1,1,1,1,1,1]);
for i=1:64
h(i,:)=1/sqrt(2)*(randn(1,N/8) + 1i*randn(1,N/8));
end
H=reshape(h,8,N);
y=reshape(sum(H.*xmod,1),8,N/8);
H=reshape(h,8,8,N/8);
snr=linspace(0,18,19);
ser=zeros(1,length(snr));
for ii=1:length(snr)
N1=1/sqrt(2)*(randn(1,N)+1i*randn(1,N));
N1=reshape(N1,8,N/8);
ynoisy=y+10^(-(snr(ii)-10*log10(16))/20)*N1;
%ynoisy=awgn(y,snr(ii),'measured');
ynoisy=reshape(ynoisy,8,1,N/8);
finy=[];
for kk=1:N/8
c=zeros(8,1);
r=ynoisy(:,:,kk);
Heq=transpose(H(:,:,kk));
for jj=1:8
B=pinv(Heq);
b=abs(B.*B);
b=b(:,1)+b(:,2)+b(:,3)+b(:,4);
b=b+c;
[ss, dd] = min(b,[],1);
c(dd)=100;
temp=qamdemod(B(dd,:)*r,4);
ydemod(dd)=temp;
r=r-qammod(temp,4).*Heq(:,dd);
Heq(:,dd)=0;
end
finy=[finy,ydemod];
end
[num ty]=symerr(x,finy);
ser(ii)=ty;
ty
end
semilogy(snr,ser,'g-*');
grid on;hold on;
title('Plot of Symbol error rate for Uncoded 8X8(4-QAM) System','FontSize',12);
legend('sim (nTx=8, nRx=8, Uncoded(4-QAM)) Using V-Blast/ZF','location','southwest');
xlabel('SNR(dB) ---->','Color','k','FontSize',11);
ylabel('Symbol Error rate ---->','Color','k','FontSize',11);
vblastZF.rar_blast_capacity
版权申诉
43 浏览量
2022-07-15
18:20:24
上传
评论
收藏 804B RAR 举报
小贝德罗
- 粉丝: 70
- 资源: 1万+
最新资源
- keil2 + proteus + 8051.exe
- 1961ee27df03bd4595d28e24b00dde4e_744c805f7e4fb4d40fa3f695bfbab035_8(1).c
- mediapipe-0.9.0.1-cp37-cp37m-win-amd64.whl.zip
- windows注册表编辑工具
- mediapipe-0.9.0.1-cp37-cp37m-win-amd64.whl.zip
- 校园通行码预约管理系统20240522075502
- 车类型数据集6250张VOC+YOLO格式.zip
- The PyTorch implementation of STGCN.STGCN-main.zip
- 092300108.cpp
- 车类型数据集6000张VOC+YOLO格式.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈