#include?
#include?
#include?
#include?

using?namespace?itpp;
using?std::cout;
using?std::endl;

namespace?itpp?{
????template?Num_T*?begin（Vec?&v）?{
????????return?&v[0];
????}

????template?Num_T*?end（Vec?&v）?{
????????return?&v[v.size（）];
????}
};

int?main（int?argc?char?*argv[]）{
	int?iter=100;	??//?No.?of?iterations
	int?C=3;		??//?#?of?cells?in?a?network
	int?U=12;	??????//?#?of?active?users?in?each?cell?（worst?case）
	int?N=20;??????????//?#?of?maximum?transmitting?antennas?at?each?base?station
	int?P=19;?	??????//?total?length?of?past?samples?of?y/h_bcu?used?（Samples?for?operating）
	it_file?ff;
	Real_Timer?tt;?
	tt.tic（）;
	if?（argc?>?1）?{
	iter?=?atoi（argv[1]）;
	cout?<	}
	if?（argc?>?2）?{
	C?=?atoi（argv[2]）;
	cout?<	}
	if?（argc?>?3）?{
	U?=?atoi（argv[3]）;
	cout?<	}
	if?（argc?>?4）?{
	N?=?atoi（argv[4]）;
	cout?<	}
	if?（argc?>?5）?{
	P?=?atoi（argv[5]）;
	cout?<	}
	int?l=3;????????//?#?of?past?samples?considered
	vec?final_sum_rate=“0.0“;?//?stores?the?total?sum
	vec?avg_sum_rate=“0.0“;???//?average?of?data?rates?for?U?users
	ivec?Nt_vals?=?linspace_fixed_step（20?N?10）;
	int?Nt;	
	int?no_iterations_in_antenna;
	no_iterations_in_antenna=（（N-Nt_vals[0]）/10）+1;
	vec?Data_Rate[no_iterations_in_antenna];
	int?count=0;
	//Nt=N;?//?look?here?loop?terminating?guy

	QPSK?qpsk;?????????????????????//The?QPSK?modulator?class???
	AWGN_Channel?awgn_channel;?????//The?AWGN?channel?class

	/*?Generation?of?H_bc/h_bcu??which?is?coming?from?cell?to?base?station?*/
	RNG_randomize（）;
	double?Ec=1.0;
	double?Eb=Ec/2.0;
	vec?EbN0dBN0EbN0;

	EbN0dB=linspace（0.02010）;
	EbN0=inv_dB（EbN0dB）;
	N0=Eb/EbN0;???//This?will?also?work

	cout<<“No.?of?Cells?Selected?=“<	cout<<“No.?of?Users?Selected?=“<	

	for?（int?Nt?:?Nt_vals）?{
?????????final_sum_rate=“0.0“;
?????????cout<<“No.?of?antennas?at?each?base?station?=“<?????????for?（int?it?=?0;?it?			cout<<“The?iteration?number?is?=“<
	/*?QPSK?Symbol?generation?*/
	int?Number_of_bits;
	bvec?transmitted_bits;
	cvec?transmitted_symbols;
	Number_of_bits=2*U;??//?Focusing?on?one?cell?only?&?others?will?follow
	transmitted_bits?=?randb（Number_of_bits）;
	transmitted_symbols?=?qpsk.modulate_bits（transmitted_bits）;

????//?Corr?matrix?of?the?channel?coefficients
????cmat?R_bcu[C][C][U];
????cmat?Rsqr_bcu[C][C][U];??//?Rsqr_bcu?represents?the?square?root?of?the?corresponding?complex?matrix

????for（int?i=0;i????	for（int?j=0;j????	????for（int?k=0;k????	????????cmat?unsymmetricR_bcu?=?randn_c（NtNt）;??//dummy?cmat?to?construct?correlation?mat
????????????????R_bcu[i][j][k]=?unsymmetricR_bcu*hermitian_transpose（unsymmetricR_bcu）;
????????????????Rsqr_bcu[i][j][k]=?sqrtm（R_bcu[i][j][k]）;?
			}
		}
	}
	
????//?Declaring?channel?coeff.
???

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-04-05?02:42??Pilot_Contamination_Precoding-master\
?????文件?????????124??2018-04-05?02:42??Pilot_Contamination_Precoding-master\README.md
?????文件???????12658??2018-04-05?02:42??Pilot_Contamination_Precoding-master\downlink_precoding.cpp
?????文件???????12946??2018-04-05?02:42??Pilot_Contamination_Precoding-master\optimizing_lambda.ipynb
?????文件????????1382??2018-04-05?02:42??Pilot_Contamination_Precoding-master\pilot_contamination_data.txt
?????文件???????12056??2018-04-05?02:42??Pilot_Contamination_Precoding-master\prev_downlink_kalman.cpp