ZF-QAM-1.zip_64-QAM_MIMIO ZF_QAM ZF_ZF 64 _zf qam

clc clear all; %4X4 64-QAM ZF MIMIO DETECTION close all; N =4*10^4; %number of generated bits M=4; x = randi([0 63],1,N); %generating random integer values from 4 levels ,it can be changed for high modulation schemes xmod=qammod(x,64);%4qam modulation scheme for x values xmod=reshape(xmod,4,N/4); %inserting the values into (8 x n/8) matrix for 8x8 system xmod=kron(xmod,[1,1,1,1]);%repeting every modulated element 8 times for i=1:16 h(i,:)=1/sqrt(2)*(randn(1,N/4) + 1i*randn(1,N/4));% generating the cahnnel complex matrix(1x10000) values depending on gaussian sidtribution end H=reshape(h,4,N); y=reshape(sum(H.*xmod,1),4,N/4); %getting the values of received signals H=reshape(h,4,4,N/4);%the final H matrix 8X8 changed every time sample snr=linspace(0,18,19); %generatinig 0 to 18 values ser=zeros(1,length(snr));%generating 19 zeros vector for ii=1:length(snr) N1=1/sqrt(2)*(randn(1,N)+1i*randn(1,N)); N1=reshape(N1,4,N/4);%generating different gaussian matrix values of noise ynoisy=y+10^(-(snr(ii)-10*log10(16))/20)*N1;%adding additive white gaussian noise to the received signal %ynoisy=awgn(y,snr(ii),'measured'); ynoisy=reshape(ynoisy,4,1,N/4);%building a matrix of received symbols(8 by 8) B=[];recvd=[]; %in this part we generated the signals, channel matrix ,received signa for kk=1:N/4 Heq=transpose(H(:,:,kk)); %matrix transpose for each symbols or time slot H^T B=pinv(Heq); %(H^T)^-1 recvd=[recvd,B*ynoisy(:,:,kk)]; end finy=qamdemod(reshape(recvd,1,N),64); %demodulation process [num ty]=symerr(x,finy);%compare the error rate between x,finy ser(ii)=ty; end semilogy(snr,ser,'r-*'); grid on;hold on; title('Plot of Symbol error rate for Uncoded 4X4(64-QAM) System','FontSize',12); legend('sim (nTx=4, nRx=4, Uncoded(64-QAM)) Using ZF','location','southwest'); xlabel('SNR(dB) ---->','Color','k','FontSize',11);Ylabel('Symbol Error rate ---->','Color','k','FontSize',11);