BPSK.rar_For Each Other_Gram Schmidt_Gram-Schmidt

clear; clc; %---------Input Fields------------------------ N=10000000; %Number of input bits EbN0dB = -6:2:10; % Eb/N0 range in dB for simulation %--------------------------------------------- data=randn(1,N)>=0; %Generating a uniformly distributed random 1s and 0s bpskModulated = 2*data-1; %Mapping 0->-1 and 1->1 M=2; %Number of Constellation points M=2^k for BPSK k=1 Rm=log2(M); %Rm=log2(M) for BPSK M=2 Rc=1; %Rc = code rate for a coded system. Since no coding is used Rc=1 BER = zeros(1,length(EbN0dB)); %Place holder for BER values for each Eb/N0 index=1; for k=EbN0dB, %------------------------------------------- %Channel Noise for various Eb/N0 %------------------------------------------- %Adding noise with variance according to the required Eb/N0 EbN0 = 10.^(k/10); %Converting Eb/N0 dB value to linear scale noiseSigma = sqrt(1./(2*Rm*Rc*EbN0)); %Standard deviation for AWGN Noise noise = noiseSigma*randn(1,length(bpskModulated)); received = bpskModulated + noise; %------------------------------------------- %Threshold Detector estimatedBits=(received>=0); %------------------------------------------ %Bit Error rate Calculation BER(index) = sum(xor(data,estimatedBits))/length(data); index=index+1; end %Plot commands follows plotHandle=plot(EbN0dB,log10(BER),'r--'); set(plotHandle,'LineWidth',1.5); title('SNR per bit (Eb/N0) Vs BER Curve for BPSK Modulation Scheme'); xlabel('SNR per bit (Eb/N0) in dB'); ylabel('Bit Error Rate (BER) in dB'); grid on; hold on; theoreticalBER = 0.5*erfc(sqrt(10.^(EbN0dB/10))); plotHandle=plot(EbN0dB,log10(theoreticalBER),'k*'); set(plotHandle,'LineWidth',1.5); legend('Simulated','Theoretical'); grid on;