Long back in time we discussed the BER (bit error rate) for BPSK modulation in a simple AWGN channel (time stamps states August 2007). Almost an year back! It high time we discuss the **BER for BPSK** in a **Rayleigh multipath channel**.

In a brief discussion on Rayleigh channel, wherein we stated that a circularly symmetric complex Gaussian random variable is of the form,

,

where real and imaginary parts are zero mean independent and identically distributed (iid) Gaussian random variables with mean 0 and variance .

The magnitude which has a probability density,

is called a** Rayleigh random variable. **This model, called **Rayleigh fading channel model**, is reasonable for an environment where there are large number of reflectors.

## System model

The received signal in **Rayleigh fading channel** is of the form,

, where

is the received symbol,

is complex scaling factor corresponding to **Rayleigh** multipath channel

is the transmitted symbol (taking values +1′s and -1′s) and

is the Additive White Gaussian Noise (AWGN)

**Assumptions**

1. The channel is flat fading – In simple terms, it means that the multipath channel has only one tap. So, the convolution operation reduces to a simple multiplication. For a more rigorous discussion on flat fading and frequency selective fading, may I urge you to review Chapter 15.3 **Signal Time-Spreading** from [DIGITAL COMMUNICATIONS: SKLAR]

2. The channel is randomly varying in time – meaning each transmitted symbol gets multiplied by a randomly varying complex number . Since is modeling a Rayleigh channel, the real and imaginary parts are Gaussian distributed having **mean 0** and **variance 1/2**.

3. The noise has the Gaussian probability density function with

with and .

4. The channel is known at the receiver. Equalization is performed at the receiver by dividing the received symbol by the apriori known i.e.

where

is the additive noise scaled by the channel coefficient.

## Bit Error Rate

The equations listed below refers Chapter 14.3 in [DIGITAL COMMUNICATION: PROAKIS]

If you recall, in the post on** BER computation in AWGN**, the probability of error for transmission of either +1 or -1 is computed by integrating the tail of the Gaussian probability density function for a given value of bit energy to noise ratio . The bit error rate is,

.

However in the presence of channel , the effective bit energy to noise ratio is . So the bit error probability for a given value of is,

,

where .

To find the error probability over all random values of , one must evaluate the conditional probability density function over the probability density function of .

**Probability density function of **

From our discussion on chi-square random variable, we know that if is a Rayleigh distributed random variable, then is chi-square distributed with two degrees of freedom. since is chi square distributed, is also chi square distributed. The probability density function of is,

.

**Error probability**

So the error probability is,

.

Somehow, this equation reduces to

.

Note:

1. I have not yet figured out the math to reduce the above integral to the answer. If some one knows, kindly drop in a comment. 2. Another way for finding the bit error rate might be to find the pdf of . However, I do not know how to find pdf following the division of two random variables.## Simulation Model

It will be useful to provide a simple Matlab/Octave example simulating a BPSK transmission and reception in Rayleigh channel. The script performs the following

(a) Generate random binary sequence of +1′s and -1′s.

(b) Multiply the symbols with the channel and then add white Gaussian noise.

(c) At the receiver, equalize (divide) the received symbols with the known channel

(d) Perform hard decision decoding and count the bit errors

(e) Repeat for multiple values of and plot the simulation and theoretical results.

Click here to download Matlab/Octave script for BER computation of BPSK in Rayleigh fading channel

**Figure: BER plot of BPSK in Rayleigh fading channel**

When compared to the AWGN case, around **25dB degradation due to the multipath** channel (at the point). This is both good and bad: bad because we need to spend so much energy to get a reliable wireless link up (in this era of global warming), and good because we signal processing engineers are trying to figure out ways for improving the performance.

## Reference

[DIGITAL COMMUNICATION: PROAKIS] Digital Communications by John Proakis

D id you like this article? Make sure that you do not miss a new article
by subscribing to RSS feed OR subscribing to e-mail newsletter.
* Note: Subscribing via e-mail entitles you to download the free e-Book on BER of BPSK/QPSK/16QAM/16PSK in AWGN.*

{ 279 comments… read them below or add one }

This website truly has all oof the info I needed about this

subject and didn’t know wwho to ask.

sir,

i am generating a channel with the modelling (outdated channel model)

g=a*h+sqrt(1-a^(2))*n

where ‘a’ is the correlation coefficient with 0<a<1, h is the Rayleigh channel, n has the same distribution as h i.e Rayleigh.

for a=1, g=h, perfect, and i can model 'g' as 'h'. h=sqrt(1/2)*(randn(1,N)+j*randn(1,N)).

but for the values a<1, how to model the channel 'g'.

hello sir,

In ur code, you have written

EsN0dB = EbN0dB + 10*log10(nDSC/nFFT) + 10*log10(64/80);

Can you pls explain how the third term “10*log10(64/80)” comes??

@Simran: The 64/80 is to handle the cyclic prefix in the OFDM symbol. In this case we are assuming that the ofdm symbol is of duration 3.2us and has an extra 0.8us of cyclic prefix.

So the ratio is 3.2us/4us which with 20Msps sampling translates to 64/80.

Hi, why is the rayleigh channel same as the noise channel? i thought the rayleigh channel should be simulated using the rayleigh random variable posted in the other article ?

@Bai: Rayleigh channel is formed using complex gaussian random variable.

http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/

Hi Krishna,

You can find the proof of bit error rate formula in David Tse’s book problem 3.1

http://snag.gy/2yRO0.jpg

@hmilan: thanks

relay network with feedback delay over rayleigh fading channel

Hi,Khishma

you said the PDF of r was :p(r) = 1/(Eb/No)*exp(-r/(Eb/No)) ,

should it be :p(r) = 1/(r’)*exp(-r/r’) ?

where , r’ = E(h^2)*Eb/No ;

.

.

@Fan: When trying to look back at the derivation of the BER @

http://www.dsplog.com/2009/01/22/derivation-ber-rayleigh-channel/

could see the following comment – “The resulting BER in a communications system in the presence of a channel h, for any random values of |h|^2, must be calculated evaluating the conditional probability density function Pb|hover the probability density function of gamma.”

Am unable to full recall the details, will ponder on this a bit, and update.

Hi Krishna,

I derived PDF of SNR in Rayleigh fading channel.

How can I share here? It is set of mathematics and if I paste it is not appearing neat.

I wish to present equations as you presented in your write up.

It is just 8 lines of equations.

Please give hints to write equations as you have written.

Thank you.

@Abhijith: You can mail me the pdf/doc and I can share it.

OR upload to an online site like dropbox and paste the public url

Your posts help me a lot,thks.

And I saw that someone had asked you some questions in Chinese.

I am curious that you could read in Chinese.Is that really ? haha

@Fan: I know as much Chinese as I know Spanish Xie xie

krishna plz give the code for ber of bpsk over rician channel

@payel saha: I have not tried modeling rician channel

Hi sir,

have recently gone through the matlab code for ” BER computation of BPSK in Rayleigh fading channel”. There, after equalization , we perform hard decision decoding by

% receiver – hard decision decoding

ipHat = real(yHat)>0;

which effectively converts all incoming complex signals into 1 having real part greater than zero. That means at this stage we are not utilizing the imaginary part.

So my question is, why should we add complex ” noise ” OR complex ” h ” ? as finally we are not utilizing imaginary part of complex signal atall ?

@Skm: In the BPSK case, as you rightly said we are not using the information on the imaginary dimension and that component is ignored. I added the imaginary part in the simulations to keep it as a general case.

n = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)]; % white gaussian noise, 0dB variance

h = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)]; % Rayleigh channel

y = h.*Signal + 10.^(-SNR_db(i)/20)*n;

hello Sir! what does those above equations mean? i will very happy if you could explain them as i am to write a code which calculate the BER of OFDM using BPSK.

Thanks

@Charly Tchouadou-Ndalleu: randn generates zero mean unity variance Gaussian distributed random variables. To add on the real and imaginary dimensions, we use two instances of randn and scale by 1/sqrt(2) to keep the variance at unity. The term n is an additive noise (modeling the thermal noise) and h is a multiplicative noise (modeling the phase and amplitude changes introduced by channel)

@Krishna: Thank you very much for your informative articles.

I have few questions slightly related to this article. I need your opinion on how you would have approached this scenario. This scenario is hypothetical.

Specification:

Assume I have to design a communication system for a specific area and specific purpose. My customer specified me minimum BER (QoS) and transmit power. I have liberty to choose rest (for eg: delay, modulation, complexity).

Approach:

Now I have to choose modulation technique, error correction algorithm. This depends on, what channel and bandwidth is.

As it is specific place, I don’t know how channel looks like . I would have approached this way,

Step 1: Measure and model channel:

Assume I am able to measure channel. Find PDF of the channel somehow.

Using this PDF, I can find outage probability or Bit error rate.

Step 2: To meet this BER, choose modulation technique and error correction algorithm.

Step3: Decide how much Sync you need, based on Channel.

Step 4: Start design and test back.

My question is Q1: How to measure real time channel if I don’t have built communication system?

Q2. Is my approach is right?

Q3. If I am right, have I overlooked something.

Thank you very much.

I am waiting for your inputs.

With regards,

Abhijith Gopalakrishna.

@Abhijith:

Nice approach. One would also need to look at the range requirement of the communication link. To start off, assuming a free-space path loss model, transmit power and receiver noise bandwidth one can get a good estimate for the range.

Now depending on whether the system is going to be indoor or outdoor, the free-space path loss model needs to be tweaked. Doing a channel measurement in the target deployment environment is definitely a good idea, however it may be expensive and time consuming. Alternately, one can look at typical channel models available in the literature and pick one which can closely match your desired use case. Hope this helps.

hello krishna,

I want to design a transceiver for wireless sensor networks for the application of home automation or for the indoor environment. So, I want to know about the energy efficient modulation schemes which is suitable for my appliction and what are the characteristics i have to consider for choosing these modulation schmes, so under these conditions which schemes is more suitable for my application..

@DHARMESH: The key aspects to look out are – data rate requirement and the range (since it is indoor, am assuming that it is of the order of 50meters). Typically, simpler the modulation scheme (i.e BPSK ish), the RF and associated circuits can be cheaper and power efficient.

The Zigbee standard seems to enable a cheap and power efficient transceivers

http://en.wikipedia.org/wiki/ZigBee

How can I do this using “Gamma-Gamma Channel”. Please help me out.

@Rupom: what is gamma-gamma channel – any pointers?

Hi Krishna, thanks for the blog and the MATLAB demo..

I have some questions about the commands

1. This line below, h is normalized by sqrt(2) to make the average power gain=1, am I right?

h = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)]; % Rayleigh channel

2. In the line below Why is the SNR divided by 20, not 10?

y = h.*qamSig + 10^(-Eb_N0_dB(ii)/20)*n;

3. what about the case of M-QAM? Let’s say it’s then Es_No_dB instead of Eb_N0_dB, do I need to change anything else for y in the equation above???

thank you

@Aiyu: replies

1. yes

2. the scaling is done on a voltage signal. hence using 1/20

3. Es/N0 should be straight forward. For Eb/N0 you might need to scale by M.

Es/N0 = k*Eb/N0, where k= log2(M)

Please see http://www.dsplog.com/2008/06/05/16qam-bit-error-gray-mapping/

Hey!

Based on what can you assume that channel is known in the receiver? When does the rx obtain this information? I mean channel changes a lot in a very short period, so when is the CSI supposed to be obtained?

What I noticed is that the curve obtained here matches the one obtained by the command berfading, so my q is basically about justifying how and why Rx knows CSI.

Thanks

@filip: In most cases, there will be training sequence in the transmission which will help the receiver estimate the channel. This will be a noisy estimate of the channel (which results in a performance loss), but useful nevertheless…

hello

thank u alot for your notes it was really usefull for me!

I am student of ms in system communications engineering would you mind introducing me some fresh, exciting and excellent topic in wireless communication (esp mimo and stc) for my ceminar?

thank u again and excuse me for my bad english

@sayyed: Please have a look at http://www.complextoreal.com/tutorial.htm

Very nice list of tutorials

hai,

can you please provide me a source for understanding jakes model…

@pavan: will add to the TO-DO list

Hi Krishna

DO you know the codes for the plot of BER for a BPSK modulation in a RICIAN channel? Plz reply asap

thanks

@Ritesh: Have not tried to simulate Rician channel

Mr. Sankar probably you will be the right person to solve a problem i have encountered in my MTech. project. Let us consider a Rayleigh frequency selective fading channel,while using BPSK. Now at the receiver i want to put a digital phase lock loop.

In this case how i will calculalte the BER performance. Will the DPLL degrade or enhance the perfomance of the BER

@Pranjal: Why do you want the DPLL in place? Am sure it’s not because of the channel, rather it might be necessitated by any phase errors.

Hi Krishna

you are a very helpful man .. i just want to ask you if you deal with matlab function that creat rayleigh channel >> name (rayleighchan)

its a good function .. and have a many optional parameter

thank you >> i will check your blog daily

@ammar: Thanks. I do not have the function rayleighchan() handy

Hello Mr. Krishna,

Why the theory and the simulation of Rayleigh channel in “ber-bpsk-rayleigh” is same but in ” ber-16qam-rayleigh” is so different?

@Nam: Well, the 16QAM is more error prone when compared to BPSK.

Perhaps you found this out a while back, as this is an older post, but a description of the evaluation of that integral is found in section 5.1 (pg 99-101) of the book “Digital Communication Over Fading Channels” by Marvin Simon.

@Colin: Thanks. Yes, indeed found the derivation.

I do not have the book handy, but will try to get it

Hi sir ,

I saw ur Ber_bpsk_reyleigh program i gone through that link but dont know how to send the data in the qpsk can u please help me and tell some hints about the rician channel.please sir

@seshu: My replies:

1/ For QPSK, send data on I channel as well as Q channel

2/ I have not discussed Rician channel model

hi krishna sir

I saw ur example of ber -bpsk-rayleigh channel but i need the the qpsk for the rician channel I dont know coding sir so please help me sir give some hints where to change thats enough for me sir .

@shu: My replies:

1/ For QPSK, send data on I channel as well as Q channel

2/ I have not discussed Rician channel model

thank you Mr krishna

Hello Mr. Krishna,

I’m trying to simulate frequency non-selective rayleigh fading channel like you did, but including the path loss (PL in dB) between two nodes.

So instead of having complex number h (channel) with the real and imaginary parts are Gaussian distributed random variable having mean 0 and variance 1/2,

the variance became

var = 1/(2*10^(PL/20))

Is this correct? Thank you.

@nico: The variance should use factor of 1/10. var = 1/(2*10^(PL/10)). Agree?

in the expression for y, I think you should use y=|h|x+n instead of using y=hx+n. Right?

@Hassan Moradi: No. The channel can distort the phase and amplitude of the transmit signal

what is meant by a power limited in communication system.

@Tom: A system where we cannot transmit very high power (due to battery/size limitations etc)

I simulate 16QAM modulation (no coding) over Rayleigh channel, I want to change profile of channel, so I write:

chan = stdchan(1/Rb,150,’cost207RAx4′);

chan.PathDelays = [0, 0.2*1e-6, 0.4*1e-6, 0.5*1e-6];

chan.AvgPathGaindB = [0, -2, -10, -16];

chan.NormalizePathGains = 1;

chan.ResetBeforeFiltering = 0;

I need your help, could you show me I am wrong or right?

@Tuyen Tran: I do not have matlab, hence unable to comment on the stdchan() function. But in general, your code seems to be correct. Make sure that you convolve your time domain samples with the channel taps.

Thank you for your respone!

Hi Krishna,

Do u have any topic related to BER for BPSK in Rician channel? I’d like to know how to write (matlab code) for Rician… Tks.

@Rebecca: Sorry, I have not discussed Rician channel

hii sir

how to get rayleigh channel using BPSK modulation & adding AWGN noise by usingg taps.actually what is tap? &

why it should be used in rayleigh channel?

@dolly: A signal transmitted over air can reach the receiver through different paths. The gain and phase information corresponding to each path is notionally stored as tap. In this post, I have assumed a single tap Rayleigh channel

I have discussed briefly on

a) Rayleigh channel @ http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/

b) 3 Tap fixed ISI channel with ZF equalization http://www.dsplog.com/2009/11/29/ber-bpsk-isi-channel-zero-forcing-equalization/

hey i just had a look on the matlab code that you posted..

http://www.dsplog.com/db-install/wp-content/uploads/2008/08/script_ber_bpsk_rayleigh_channel.m

% Transmitter

ip = rand(1,N)>0.5; % generating 0,1 with equal probability

s = 2*ip-1; % BPSK modulation 0 -> -1; 1 -> 0

correct me if i am wrong..cos i am new to matlab and to this topic i have very less knowledge..

i checked the output for n=10

ip =1 0 0 1 0 1 1 1 0 0

s =1 -1 -1 1 -1 1 1 1 -1 -1

please help me out..

cheers bijoy

@bijoy babu: Yes

Hi friend, please, I do not understand why the BER of BPSK in Rayleigh channel is equal to the BPSK OFDM in Rayleigh channel. Thank you very much.

@FERNANDO: In this simulation model, even though we have a 10-tap multipath channel, as the channel duration is less than the cyclic prefix (of 16 samples), we do not have inter symbol interference. Each subcarrier experiences flat fading as in the BER of BPSK in Rayleigh channel case. Hence the BER is identical.

hi,

I have similar kind of project to do but in my project prof is asking to plot BER of BPSK i.e. BER vs γ. Does it look different from when we plot it against Eb/No?

@Saurabh: What is y?

Can you distinguish clearly correlated and uncorrelated Raykeigh Channel and how are they formed from awgn? How to implement BER simulation for correlated and uncorrelated Rayleigh Channel?

@kadir: AWGN is an additive noise whereas as the channel is a multiplicative noise. For uncorrelated channel, the channel varies for each symbol. For correlated channel, the channel varies slowly across symbols

hello mr. krishna,

i would ask you a question regarding the rayleigh channel. how can i generate Rayleigh channel impulse response as ready numbers to be used in a program without using the rand()+j*randn() commands. in other words, i don’t want to have a loop or large number of symbols to get the plot of the prob. of error. I want to have a table of values for channel impulse response to be used for simulating a MIMO system. some papers do list down in a table the channel as complex values between each transmit and receive antenna to get a MIMO channel matrix, how did they get such values? this is exactly my question

Hi Krishna,

I have a little query..In your MATLAB program of calculating BER for BPSK in raleigh fading Channel,while calculating you have used

yhat=y./h;

Is it OK to divide y by h as above, as it may lead to severe noise amplification as h is a random variable and may take any value?

Thanks.

@Sanket: You are right, it can indeed result in noise amplification…

Hi Sir,,

in book Fundamental of Wireless Communication (by Tse and Viswanath)

Rician fading modeled as:

h = sqrt(K/(K+1))*line_of_sight_signal + sqrt(1/(K+1))*reflected_and_scattered_signal

K = K-factor

From the above equation, i have modified your code (rayleigh channel) in order to simulate ricean channel.

% generating 0,1 with equal probability

ip = rand(1,N)>0.5;

% BPSK modulation 0 -> -1; 1 -> 0

s = 2*ip-1;

% Rician channel

h = 1/sqrt(2)*(sqrt(K/(K+1))*s + sqrt(1/(K+1))*(randn(1,N) + j*randn(1,N)))

is this code correct?

please help me..

thanks

@sinto: Well, from a quick look it does not seem to be right.

hi kishna..

I am working in turbo decoding.my recieved signal is effected by AWGN noise and reyleigh fading. which is in thae same formate which u have taken like hS+W.

My problem is that i want the exact expession of probability of error for QPSK.

Whatever u have derieved is for BPSK case.Can u derive it or please give me some hint.

Anyways as u asked the full mathematical concept of ur expression of BER for BPSK,u can go through my thesis supervisors book,named by DIGITAL COMMUNICATION AND SIGNAL PROCESSING BY K.VASUDEVAN.

hope this will help u and waiting for ur responce.

Hi Krishna,

I hope you are doing well. I am simulating 16 QAM in a Rayleigh channel and am trying to compute the SER(& BER). I need to verify my results using the theoretical formula. Can you kindly tell me the formula of SER for 16 QAM in Rayleigh channel?

I shall be grateful

Thanking in advance,

Zaffar

hello mr. krishna,

i would ask you a question regarding the rayleigh channel. how can i generate Rayleigh channel impulse response as ready numbers to be used in a program without using the rand()+j*randn() commands. in other words, i don’t want to have a loop or large number of symbols to get the plot of the prob. of error. I want to have a table of values for channel impulse response to be used for simulating a MIMO system. some papers do list down in a table the channel as complex values between each transmit and receive antenna to get a MIMO channel matrix, how did they get such values? this is exactly my question.

Hi

How r u ? hope fine. I m very satisfy to this site, because it is very helpful to me to make my project. if u have QPSK matlab code for BER in Rayleigh channel, pls send me.

thanks.

Hi

I would like to know that how IFFT makes signls to be orthogonally modulated in OFDM..also when we give 1s and 0s as input to IFFT block the o/p is complex value…and some of them are very close to zero..for eg for 1 it is 0.012+j0.032…suppose if this is transmitted thro’ severely faded channel still there is a possibility that these signal might get into deep fade..so at the receiver there might be a problem in decoding the data as detector would decide upon the quantized value…

@prabag: When we are doing IFFT, we are loading information on subcarriers exp(jwt), exp(j2wt) and so on… you may find a detailed discussion at http://www.dsplog.com/2008/02/03/understanding-an-ofdm-transmission/

At the receiver, we are not going to look at the time domain samples and attempt demodulation. We take the fft and look at the fft output prior to attempting demodulation

http://www.dsplog.com/2008/08/26/ofdm-rayleigh-channel-ber-bpsk/

Hello

I am working on a paper in which i am trying to plot ber for bpsk using rayleigh channel. but, the thing is, i need to use estimated SNR including the pilot symbols and noise samples as parameters. I tried with the corresponding SNR estimate equations in the paper, but, i am experiencing constant BER.

I have followed the same steps posted by you in “BER for BPSK in Rayleigh channel” to produce the rayleigh output and later used estimated SNR to calculate BER, so that i can vary pilot symbols(N) and Noise samples(L). But, i am stuck. can you please refer my paper and guide me, whether i am in right way..??

http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=01437347&tag=1

@student: Sorry, due to time constraints I wont be able to help with the coding part. Good luck.

10^(-Eb_N0_dB(ii)/20): help me please

wy you used teh signe negatif – in( -Eb_N0_dB(ii)/20), and d’ont Eb_N0_dB(ii)/20

thanks

10^(-Eb_N0_dB(ii)/20): help me please

wy you used teh signe negatif – in( -Eb_N0_dB(ii)/20), and d’ont Eb_N0_dB(ii)/20

thanks

@aissou: The negative sign came as I am scaling the noise voltage. I am keeping the signal swing the same and reducing the swing of noise voltage to simulate various Eb/N0 values.

Krishna Pillai,

Hellos sir,

I am told to assume Rayleigh fading channel with BPSK modulation. Using MATLAB

plot bit error probability (BEP) under coherent and non-coherent detection

when receiver is equipped with three antennas to exploit diversity.

Your figures will include plots from simulation. Use average SNR (complex)

from -5 to 20 dB.

Please help me resolve this, thank you so much

Krishna Pillai,

Hello Sir, I was asked to assume Rayleigh fading channel with BPSK modulation. Using MAT-LAB plot bit error probability (BEP) under coherent and non-coherent de-tection. Your ¯gures should include plots from both analysis and simulation.Use average SNR (complex) from -5 to 20 dB.

I am having difficulty writing the code for non coherent, I have been able to write the code for coherent detection.

Please assist me in this.

@Obinna O: Hope you have finished the project by now.

Hello sir,

The posts here are very useful. great job.

Presently,I am working on a IEEE paper “On the Impact of SNR Estimation Error on Adaptive Modulation” . Here is th link

http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=01437347

But, I got stuck with the silmulation of the imperfect SNR Estimation.

I need to plot the fig.1 i.e; Estimated SNR vs pdf basing on the theoritical equations stated under imperfect SNR Estimation. can u please help me with the simulation ??

Thank You.

@usha: Sorry, due to time constraints, may I decline to help with the simulation. Good luck.

hi

i was trying to simulate BER for BPSK using Rayleigh fading channel.

But the difference is, I need to follow the steps,

1. Generate BPSK-modulated random data for the desired signal and

interfering users.

2. Generate the independent complexed Gaussian random vectors and multiply them with the corresponding covariance matrix represented by J0(x) (i.e., Jo(2*pi*d) where d is the ratio of the spacing between the two adjacent antennas to the wavelength) to get the correlated Rayleigh fading channel.

3. Generate complexed AWGN

4. Make hard-decision on the desired signal to get the BER.

And the algorithm you have followed says:

Simulation Model

(a) Generate random binary sequence of +1’s and -1’s.

(b) Multiply the symbols with the channel and then add white Gaussian noise.

(c) At the receiver, equalize (divide) the received symbols with the known channel

(d) Perform hard decision decoding and count the bit errors

(e) Repeat for multiple values of and plot the simulation and theoretical results.

My question is,

in the 2nd step, i am multipying the symbols with covariance matrix to obtain Rayleigh fading channel, and then in step 3, i am generating AWGN.

whereas, you are multiplying the symbols with the channel and then adding AWGN to it.

Do both the simulations yield similar results?

Can you please make me understand the difference.

@Student: In my simulation model, the channel is independent flat fading Rayleigh channel which is applied to a single users transmission and then we add awgn to it.

For your case, its bit more complex. We need to have

a) Generate output of channel for desired user by using flat fading independent Rayleigh channel

b) Generate output of channel for undesired user by using correlated flat fading Rayeligh channel

c) Add (a) and (b), and then add awgn.

So, your results will be poorer due to the additional interference from the undesired user.

Good luck.

10^(-Eb_N0_dB(ii)/20) ; pourquoi utiliser le moin ‘-’

@bopuhafs: To make the noise power to be lower than the signal power

Hi Krishna Sankar

My name is Puripong, I’ a student from Thailand.

I’m looking for the introduction to simulate BER performance in Rayleigh fading channel and this is a very nice post {BER for BPSK in Rayleigh channel}.

I found 25dB degradation when compared to BPSK-AWGN

Thank you very much.

I have some question about h,

Why does the real and imaginary parts must have variance equal to 1/2 ?

Could you give me an explaination or some reference textbooks for this question ?

Hi Krishna

I have a problem in my m.file.I produce a symbol strean by x=randint(N,1,M) where M is the order of modulation, the fading channel by H=sqrt(0.5)*(randn(N,1)+j*randn(N,1)).In this state the channel is changing every one bit(if M=2), is this true?I think it is not possible in real systems.What should i do if I want to keep the channel fixed in every coherence time of channel?How can i write the m-file of this issue?

help me please.

thank u.

Sara

@Sara: Well, the most simplest way would be to keep the same channel for a group of symbols, then use an independent realization for the next group of symbols. If we want t make the channel evolve over time, then it becomes slightly more sophisticated…. I have not yet discussed that.

Hi Krishna

I have a problem in my m.file that i sent to your email.I want you to help me about this problem. thank you.

@Sara: I will take a look. Sorry for the delay.

i m very satisfy to this site, bcause it is very hepful to me to make mmy project. if u have qpsk matlab code for BER then plz send me .

@amit: Thanks. For QPSK, plz refer

http://www.dsplog.com/2007/11/06/symbol-error-rate-for-4-qam/

HI Krishna,

I sent you an email regarding my project.

Please reply to that mail.

thanks very much for replying

i doo agree with ur comment

however i dint get my answer yet

can u plz tell me the code for channel model for bpsk

tht is frequency selective channel (u have done it with flat faing) ? with the help of this i can compare the ber of bpsk with rayleigh as compared to ber of bpsk rayleigh ofdm .so i can show how effectiev ofdm is in frequency slective channel.

can i use the same 10tap model used in ur code for bpskofdm in the case of bpsk only..

or i shuld use some other code..?

plz help me

thanks

@mak_m:

1/ Yes, with OFDM, the 10 tap channel is in effect gives same performance as flat fading channel

2/ You can use 10 tap channel with BPSK modulation. However, you need to do pulse shaping etc at the transmitter.

thanks , hi krishna can u plz tel me…i m really confused .can u plz tell me ray leigh channel used in BER for BPSK(single tap) is different from ber for bpsk ofdm(10tap)… but their graphs are same… they both are f lat fadding?right

but if the ber of bpsk is same as ber of bpsk in ofdm in , then y did we use ofdm ..can u plz tell me which channel model shuld i use so i can show the performance improvement due to ofdm which i cant see in both of ur programs…

tel me which channel model shuld i use in case of bpsk without ofdm…ie to degrade the ber of ur simulation..

ihave gone through all the comments but didnt find the answers..

thanks in advance

waiting for ur reply

@mak_m:

1/ Thanks to OFDM, even though we have a 10-tap frequency selective channel, each subcarrier experienced a flat fading channel. Hence the performance of flat fading BPSK is comparable to OFDM performance.

2/ If we do not use OFDM, the performance for BPSK in a 10 tap channel will be much poorer due to inter symbol interference.

Hope this helps.

can u plz tel me how did u drive the formula…for ber of bpsk in rayleigh channel…can i get further explanation

can u plz telme the source of it

thanks in advance

gud work keep it up

@mak_m: I am referring Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt

Dear Krishna Sankar,

your posts helps me a lot. I want to how can i find the BER equation for 16-QAM and QPSK over the RAYLEIGH channel ? please help me….

@M_abs: You may modify the modulated variable s to correspond to QPSK/16QAM respectively. Make sure that you have the corresponding demodulators too. You may refer to

http://www.dsplog.com/2007/11/06/symbol-error-rate-for-4-qam/

http://www.dsplog.com/2007/12/09/symbol-error-rate-for-16-qam/

http://www.dsplog.com/2008/06/05/16qam-bit-error-gray-mapping/

hi…

I have a doubt….within coherence time channel coffecients are constant…

if I am correct than why we are multiplying each data bit by rayleigh variable…..??

h = 1/sqrt(2)*[randn(nTx,N) + j*randn(nTx,N)]; % Rayleigh channel

sr = (1/sqrt(nTx))*kron(ones(nTx,1),s);

% Channel and noise Noise addition

hEff = h.*exp(-j*angle(h));

y1 = sum(h.*sr,1) + 10^(-Eb_N0_dB(ii)/20)*n;

please reply

@vikas: My replies:

1/ “within coherence time channel coffecients are constant…” [krishna] yes, channel coefficients are almost the same

2/ “why we are multiplying each data bit by rayleigh variable…..??” [krishna] rayleigh variable corresponds to channel coefficients. these coefficients remain the same for the coherence time duration

ok. What you said is one of solution for it. After I read sklar’s book again for fading, I know the anwser , thanks!

Hi Krishna,

“When compared to the AWGN case, around 25dB degradation due to the multipath channel (at the point). This is both good and bad: bad because we need to spend so much energy to get a reliable wireless link up (in this era of global warming), and good because we signal processing engineers are trying to figure out ways for improving the performance.”

1. I think ” because we signal processing engineers are trying to figure out ways for improving the performance” can be undestandbe that if channel estimation is perfect, we can try to combat effect of rayleigh channel and get Ber vs EbN0 much closer to that of AWGN.

2. in your code for “BER for BPSK modulation in a Rayleigh fading channel” , there is only one rayleigh path, so the channel is flat.

% Channel and noise Noise addition

y = h.*s + 10^(-Eb_N0_dB(ii)/20)*n;

% equalization

yHat = y./h;

you try to simulate as having a perfect channel estimation and then equalizes it. But the performance is still too far away from that of AWGN. What is the better channel estimation and equalization method for it to get close to performance in AWGN?

@lumingui: You are right in saying that, in this current simulations where we know the ‘perfect channel’ and equalization, we are having optimal performance. However, think about the case where we have

(a) Multiple receive antennas.

http://www.dsplog.com/2008/09/28/maximal-ratio-combining/

Maximal Ratio Combining and bring down the BER by quite a bit. And theoretically, if we have infinite receive antennas, I think we reach close to AWGN performance

(b) Multiple transmit antennas with Beamforming

http://www.dsplog.com/2009/04/13/transmit-beamforming/

(c) Multiple transmit antennas with Space Time Block Coding (STBC)

http://www.dsplog.com/2008/10/16/alamouti-stbc/

Does this answer your concern?

Hi Krishna,

When you derived an expression for the average probability of symbol error in a rayleigh channel, you obtained a linear relationship(approximation for higher Eb/N0′s) between the logBER and the average Eb/N0(as you have shown in the plots). This derivation gives the average BER for a rayleigh channel with no mention of an EQ anywhere in the derivation. But when you simulate, you are getting the theoretical result of a linear dependence after having used an EQ. Not using an EQ would give very poor BER since the phase of the received points is now uniformly distributed and when you use a BPSK, you would have a only a 50% chance that the transmitted point would be received correctly. So although i understand the importance of using an EQ, what i dont follow is the gap between the theoretical derivation for average BER with no mention of the effects of an EQ, and the simulation which uses it. Please tell me what you think about what im saying. Thanks.

@Hemanth: You are not correct when you mentioned that the theoretical derivation of BER in Rayleigh channel did not include the effect of equalizer. Recall that, we started of the derivation by assuming that the instantaneous SNR is |h|^2Eb/N0. Its onlu after equalization the instantaneous SNR is |h|^2Eb/N0.

Note: y = hx + n

After equalization, xhat = y/h = (h*/|h^2|)y = h*(hx+n)/(|h|^2) = (|h|^2x + n)/(|h|^2)

Do you agree?

Hi Krishna,

I agree with your explanation for where the EQ comes in. I was looking at it in a different way which lead to my confusion. It stems from the fact that the instantaneous SNR is (|h|^2)Eb/No, even when the EQ is not used(dont agree with you on this point )

Since y=hx+n,

Received signal power = E[(hx)(hx)^h](^h is hermittian),

= E((hx)(x^hh^h)],

=h(Px)(h^h)

=Px(|h|^2)

Noise power remains the same as before.

So, you can see that the signal energy is scaled by a factor of |h|^2, because of which the received SNR is going to be (|h|^2)Eb/No. So whether or not you use this kind of an EQ(ZFE), your SNR remains the same.

I think where we have used the EQ implicitly, is when we say the average BER is obtained by taking an expectation over the ‘old’ expression for a BPSK ber, i.e. 0.5*erfc(sqrt(gamma))(the instantaneous ber in this case). This has been derived assuming noise adds only along a single direction, while the effect of a random channel(with uniform phase distribution) would be to throw the constellation point anywhere in the signal space. So if you say the instantaneous BER is 0.5*erfc(sqrt(gamma)), then the channel model is x+n’ or utmost kx+n’ where k is real and n’ is Gaussian and using a zfe restores our situation to this model. Tell me if this sounds right to you. Thanks.

@Hemanth: Some comments:

“It stems from the fact that the instantaneous SNR is (|h|^2)Eb/No, even when the EQ is not used”

[krishna] We can have different types of receivers. The most simplest type wont even consider that there was a channel and just do hard decision decoding on the received symbols. We might loose around 50% of the symbols, but still thats one type of receiver.

Another receiver type, which does division by the channel h, results in instantaneous SNR of (|h|^2)Eb/No.

So your statement that the instantaneous SNR is independent of equalizaton is not correct.

“This has been derived assuming noise adds only along a single direction”

[krishna] Its a kinda loose sentence, which I am finding difficult to comprehend. But, yes the noise is additive and not multiplicative. Note that additive noise can result in the transmitted signal to be anywhere in the signal space.

Dear Krishna,

I have a question related to Rayleigh Multipath Channel. Lets assume that there is no white noise caused by the channel. Additionaly the transmitted signal will be up-mixed and then down-mixed, exactly as you described in the IQ-Imbalance article, but without the phase imbalance. So the structure looks like this :

Guard Insertion -> Up-Convertion -> Multipath Channel -> Down-Convertion -> Guard Removal -> Equalizer

The problem is that the received signal is disturbed. Since we have no AWGN, the signal is ideally equalized, right ? So, what can be resposible for the signal distortion , down-converter ? How could I repair the signal then ? Can you see any problems with the transmission/reception path I described ?

Regards

Joel

@joel: The multipath channel causes a phase and amplitude distortion to the channel. That has to be removed (aka equalization) prior to demodulation. I do no see issues with the chain which you have described.

The thing is, that the equalizer removes the multipath distortion by dividing the output by the channel impule response. Since there is no noise, the qualization is ideal , i.e. Y = X*H/H = X. I actually suspected that the down-convertion introduces additional phase rotation which spoils the equalization. I hoped that you negate or confirm my assumptions.

Regards

@joel: What you said is probable. How are you estimating the channel. If you are estimating the channel by defining a known sequence (in frequency domain) at the transmitter, makes it pass through the model which you have described above, then estimate the channel @ the receiver, then the estimated channel will include the effects of multipath + phase rotations which may have occurred.

Hi Krishna,

I used your script for BER with FSK modulation in AWGN channel and modified to get BER with FSK modulation in Rayleigh channel. But there is some problem inside, i dont know that, following code which I changed is

% Matlab Script for computing the BER for Binary FSK modulation in a Rayleigh fading channel

clear

N = 10^5 % number of bits or symbols

T = 10; % symbol duration

t = [0:1/T:0.99]; % sampling instants

tR = kron(ones(1,N),t); % repeating the sampling instants

Eb_N0_dB = [0:11]; % multiple Eb/N0 values

for ii = 1:length(Eb_N0_dB)

% generating the bits

ip = rand(1,N)>0.5; % generating 0,1 with equal probability

freqM = ip+1; % converting the bits into frequency, bit0 -> frequency of 1, bit1 -> frequency of 2

freqR = kron(freqM,ones(1,T)); % repeating

x = (sqrt(2)/sqrt(T))*cos(2*pi*freqR.*tR); %generating the FSK modulated signal

% noise

n = 1/sqrt(2)*[randn(1,N*T) + j*randn(1,N*T)]; % white gaussian noise, 0dB variance

h = 1/sqrt(2)*[randn(1,N*T) + j*randn(1,N*T)]; % Rayleigh channel

% coherent receiver

y = h.*x + 10^(-Eb_N0_dB(ii)/20)*n; % additive white gaussian noise

equalization

yHat = y./h;

op1 = conv(yHat, sqrt(2/T)*cos(2*pi*1*t)); % correlating with frequency 1

op2 = conv(yHat, sqrt(2/T)*cos(2*pi*2*t)); % correlating with frequency 2

% demodulation

ipHat = [real(op1(T+1:T:end)) < real(op2(T+1:T:end))]; %

nErr(ii) = size(find([ip - ipHat]),2); % counting the number of errors

end

simBer = nErr/N;

theoryBerAWGN = 0.5*erfc(sqrt(10.^(Eb_N0_dB/10))); % theoretical ber

EbN0Lin = 10.^(Eb_N0_dB/10);

theoryBer = 0.5.*(1-sqrt(EbN0Lin./(EbN0Lin+2))); %theoretical BER

% plot

close all

figure

semilogy(Eb_N0_dB,theoryBerAWGN,'r-','LineWidth',2);

hold on

semilogy(Eb_N0_dB,theoryBer,'b-');

hold on

semilogy(Eb_N0_dB,simBer,'mx-');

axis([0 11 10^-4 0.5])

grid on

legend('theory:AWGN','theory:fsk-coh', 'sim:fsk-coh');

xlabel('Eb/No, dB')

ylabel('Bit Error Rate')

title('Bit error probability curve')

But why the theory and simulation curves are separated each other, they should cling tightly. Can you do me a favour to check and find a problem to get them closely?

Thanks,

In the simulation model above, every fading coefficient is uncorrelated and affects only 1 bit. So, is it consider as uncorrelated fast fading? What are some of the ways to estimate the channel since we need it for equalization?

@bluray: Yes, this is uncorrelated fast fading

Typically, there will be a preamble sequence which is known by both transmitter and receiver, using which receiver can form an estimate of the channel.

Dear Krishna,

can you tell me how to plot the power angular profile for an incoming signal.

@Ideal: Sorry, am not familiar with modeling angular profile.

hi……i saw youe code for alamouti STBC…..it’s a great one…..but i need a simpler one eventhough it’s too long…

i don’t know how to use functions equivelent to those you used with your code…. i need a very simple code using very simole function…..

so can anyone help me please??

@fof: Good luck.

hi krishna. First i would like to thankyou for answering the questions of all the people like me..

I have written a code for modelling OFDM using BPSK where the channel taps are given by h1=1 h2=0.5 h3=0.3 and i have used awgn channel. I am unable to figure out the mistake i have done in the code.Can u please help me??here is the code

nsubc=4;

nsym=1;

N=nsubc*nsym;

SNR=4;

a=sign(rand(1,N)-0.5);

a_vec=reshape(a,nsubc,nsym)

a_fft=sqrt(nsubc)*ifft(a_vec);

h=[1 0.5 0.3];

hnorm=h./norm(h);

chanipac=[a_fft(nsubc-1,:); a_fft(nsubc,:); a_fft];

for j=1:nsym

p=chanipac(:,j);

chanresp(:,j)=conv(p’,hnorm);

end;

n=1/sqrt(2)*(randn(nsubc+4,nsym)+i*randn(nsubc+4,nsym));

y=chanresp+10^(-SNR/20)*n;

for i=1:nsubc

cprem(i,:)=y(i+2,:);

end;

hfreq=fft(h,nsubc);

hfreqinv=1./hfreq;

hdiag=diag(hfreqinv);

o_fft=(1/sqrt(nsubc))*fft(cprem);

fde=hdiag*o_fft;

o=fde>0;

op_vect=2*o-1;

op=reshape(op_vect,1,N);

ber=size(find(a-op),2)

berrat=ber/N

@nikitha: I did not probe deeply. However, from a quick look it seems that there are issues in subcarrier assignment. You may refer to the post on BER for BPSK in OFDM with Rayleigh channel @

http://www.dsplog.com/2008/08/26/ofdm-rayleigh-channel-ber-bpsk/

Hope this helps.

hi all,

i’m modeling a SUI channel,i’ve already found a model that gives a matrix of the channel coefficients, but i don’t know how to add the channel effect to my OFDM modulated signal

thanks for anyone gives me a hand…

@wafaa: You can just convolve the ofdm time domain signal with the channel taps. Agree? You may use the code in the post

http://www.dsplog.com/2008/08/26/ofdm-rayleigh-channel-ber-bpsk/ for reference.

Thank you Krishna for all your efforts. The information you have posted about matched filter are very intersting for me

hello sir

if we were to extende ofdm to ofdma how do i divide the sub carriers to say two users

thanks

@hayle: Simplistically, I can think of two approaches

- interleaved : where all even subcarriers are assigned to one user and odd subcarriers to another user

- block : where -ve subcarriers are asisgned to one user and +ve subcarriers to another user.

And between these two approaches, there can be lots of permutations possible.

Hello Sir,

I am working on multi-relay cooperative communication system. I consider two hop communication and rayleigh fading channel with decode and forward relaying. Could you please send me the matlab program to determine the end-to-end (i.e., from source to relay and from relay to destination)outage probability under different relay selection strategy where all the relay nodes are distributed according to homogeneous Poisson point process in a circle of infinity radius.

@Alam: Sorry, I have not worked on co-operative communications.

Dear S,Alam,

I am also worikng on multi-hop relay network. I am new in this area and need some help from you. if you can give meur email address please. my email is assad.abbasi@yahoo.com.

i’ll wait for ur reply,plz contzct me on my email.

thanks

Dear S.Alam

Did you find material for this problem ?

I am also worikng on multi-hop relay network. I am new in this area and need some

help from you. if you can give meur email address please. my email nguyenkimhieuha@gmail.com

i’ll wait for ur reply,plz contzct me on my email.

your simulation on the top was for BPSK in multipath and has 35dB at 10-4 which is the same value also for OFDM for multipath.if we were to say OFDM is better we expect to get a value less than this

thanks

@mohammed: As I said in my previous comment – in both the cases the signal experiences flat fading. Hence the BER performance is comparable.

thanks kirshina

modifying BPSK do you mean using channels other than flat fadding for bpsk (give me an example)and OFDM i think is assumed flat fading .

thanks once againe

@mohammed : Yes, channels other than flat fading for BPSK.

hi kirshina

for both bpsk and ofdm using bpsk on multipath your BER curve is alomst the same,and ofdm is said to superior in multipath can you explaine this to me

thanks

@mohammed: The BPSK simulation was using flat fading channel and the OFDM simulations, thanks to cylcic prefix is also having a flat fading channel. Hence both the simulations are having identical performance.

However, if you modified the BPSK simulation to have a multipath channel, the performance will be much poorer than with OFDM.

i got the answer: for which i have posted previously:

% For Slow Fading, Channel coefficients are the same over block transmission;

h_10 = sqrt(1 / 2) * abs(randn + j * randn) * ones(1, N);

h_20 = sqrt(1 / 2) * abs(randn + j * randn) * ones(1, N);

% For Fast Fading, Channel coefficients changes from one symbol to the other

h_10 = sqrt(1 / 2) * abs(randn(1, N) + j * randn(1, N));

h_20 = sqrt(1 / 2) * abs(randn(1, N) + j * randn(1, N));

But my question is:

As we all know small scale fading (based on Doppler spread) is divided in to

1.Slow Fading

2.fast fading

for fast fading, there is high doppler spread , channel variations are ofcourse faster than signal variations, coherance time symbol period , channel variations are slower than signal variations

How can one write code for simulating doppler spread ????

Can U help me by giving hint..

thanks a lot in advance

@R.Ramya: I have not tried simulating doppler, hence unable to help. Sorry.

Hi,

I want to change it from slow fading to a fast fading.

% For Slow Fading, Channel coefficients are the same over block transmission;

h_10 = sqrt(1 / 2) * abs(randn + j * randn) * ones(1, N);

h_20 = sqrt(1 / 2) * abs(randn + j * randn) * ones(1, N);

what should i do for fast fading?

(For Fast Fading, Channel coefficients changes from one symbol to the other )

thanks a lot in advance.

Greetings Krishna. Thank you for all these efforts. I want to ask you if there is a matlab code about using matched filter as an optimum binary detection. How can this simulated?Also equivalently the correlator reciever. How can I simulate the integrator?

Thank you very much.

Ahmed

@Ahmed: Thanks.

In most of the simulation models, I have assumed flat fading channels and did not use filtering. However, there are two post on transmit pulse shaping filters:

http://www.dsplog.com/2008/04/14/transmit-pulse-shape-nyquist-sinc-rectangular/

http://www.dsplog.com/2008/04/22/raised-cosine-filter-for-transmit-pulse-shaping/

You may modify one of these to adapt to BER computation using matched filtering at receiver. You may simulate the integrator as a convolution function.

% example matlab code snippet to show matched filtering with a rectangular tx pulse shaping filter

clear all

N = 7

am = 2*(rand(1,N)>0.5)-1 + j*(2*(rand(1,N)>0.5)-1); % generating random binary sequence

fs = 10; % sampling frequency in Hz

% recatangular filter

gt1 = ones(1,fs);

amUpSampled = [am;zeros(fs-1,length(am))];

amU = amUpSampled(:).’;

% received seqeunce

rt = conv(amU,gt1);

% matched filter

yt = conv(rt,gt1/10);

% extracting the samples

at = yt(10:10:end-10)

err = at-am

Hope this helps.

hello every body,

I am doing search on multiuser detection and I need some help to simulate CDMA system over multipath rayleigh fading channel,

@rim: We have not yet written posts on CDMA systems. However, the multiuser detection problem can in general be viewed as a MIMO V-BLAST transmission scenario, where each spatial dimension acts as interferers for the other dimensions. Some equalization structures for a 2×2 MIMO with V-Blast are discussed in

http://www.dsplog.com/tag/mimo

Hope this helps.

hi all

i have 3 for-loops here

for kk=users

for jj=receiver

for ii=Eb_N0_dB

nErr(jj,ii) = size(find([ip- ipHat]),2);

end

end

end

figure

semilogy(Eb_N0_dB,simBer(1,: ,’mo-’,'LineWidth’,2);

……..code……..

y-axis is bit_error, x-axis is dB

if i want to show what the error rate for each useris. What should i put in

nErr(??,??)????

@Martin: Maybe you should make a three dimensional matrix.

nErr(kk,jj,ii) = size(find([ip- ipHat]),2);

Both Matlab and Octave supports multi dimensional matrices.

Hello sir,

Sir, I am Student of Research area, working on MIMO Adaptive equalizer, How to start in this area?

pls, help me if you have any matlab based simulation plz, send me guidelines.

thanks

-pradip

@pradip: Sorry, I have not tried simulatiing MIMO with adaptive equalizers. Were you looking for a channel tracking scheme?

hi krishna we meet again ……..

I have problems,can you help me,…

1.what impact BPSK in high frequency?

2.in this script program have one error, I confused.can you fix the error?

this the program

% BER_sic=sim_sic_fn(SNR_dB, P, Nb, code_matrix, chan_type)

% return bit error rate of the sic receiver in AWGN or rayleigh fading

%

% PARAMETER:

% SNR_dB=signal-ti-nopise ratio in dB

% P=power control vector, P(i)=transmitted power of i-th user

% Nb=number of transmitted bit

% code_matrix=matrix of the spreading code used

% chan_type=channel type(1=AWGN, 2=rayleigh fading channel)

%

% output:

% BER_sic=bit error of the sic receiver

%

function BER_sic=sim_sic_fn(SNR_dB, P, Nb, code_matrix, chan_type)

BER_sic=zeros(size(SNR_dB));

K=length(P); %number of user

N=size(code_matrix,2); %spreading factor

G(:,1)=code_matrix(33,:)’;

G(:,2:K)=code_matrix(1:K-1,:)’;

%generate the crosscorrelation matrix R

R=G’*G;

for p=1:length(SNR_dB),

SNR=10^(SNR_dB(p)/10);

SNRchip=SNR/N;

disp(‘processing..’)

error_count=0;

%loop for Nb symbol bits

for n=1:Nb,

%generate a bernoulli symbol, i.e. b=1 or 0

b=input_symbols(K,1);

if chan_type==1

A=sqrt(P’);

elseif chan_type==2

A=rayleigh(1,K).*sqrt(P’);

else

disp(‘invalid channel type’);

dbquit;

end

%transmitted sinal

X=G*(A.*b);

%generate the observed data

r=awgn_chan(x,SNRchip);

%assume that the codes are known and that they are

%perfectly synchronized

%matched filter output

y=G’*r;

%short users according tho their received powers

[y_sorted,sort_order]=sort(y.^2);

user_index=find(sort_order==1);

y_sorted=y(sort_order);

A_sorted=A(sort_order);

R_sorted=R(sort_order,sort_order);

%detect users successively

b_hat=zeros(K,1);

b_hat(K)=sign(y_sorted(K));

for u=K-1:-1:1,

b_hat(u)=sign(y_sorted(u)-sum(A_sorted(u+1:K).*R_sorted(u+1:K,u).*b_hat(u+1:K)));

end

b_hat(1)=sign(y_sorted(1)-interference_estimate(1));

%if detected symbol is incerrect, increment the error_count

%if b_hat(user_index)~=b(1)

error_count=error_count+1;

end

end

%carculating Bit Error rate i.e. the percentage of erroneous

%symbol estimates

BER_sic(p)=error_count/Nb;

end

thank before….

@wap: My replies:

1. Hmm… I dont think there is impact. By high frequency, were you refering to carrier frequency.

2. Sorry, in general. I refuse to debug the code. If you have explicit q’s I can try to answer them.

Hope this helps.

thank khrishna for your answer……………..

how in the code_matrix so that the program is replaced with the gold code? gold code so that this function can be called by sic…..

thanks…

@ w@p: Sorry, am not familiar with the code which you provided. Thanks.

Hi~

Who can tell me how to establish

communicate with basestation 1*m error rate, it means 1 transmitter and m receivers

it is randomly pick one user

Please help, I am doing my final year project, i have 1 months only, i still have many tasks havent done yet

@Martin: You question is vague. Please state your requirement in bit more details to enable a reader to help.

i need to produce a program like this requirement

we have m receivers and n users

randonly pick one user and communicate with Basestation

we need to cal.the 1*m error rate

1 mean 1 transmitter, m is receivers

and then run it 10^6 times, find out the error rate

also in this simulation, we need to use MRC

it is a cellular system only.

please help me, because i need to produce distributed antenna system also.

Please, i have not enough time to produce

i try to produce it

am i right??

—————————–% Script for computing the BER for BPSK modulation in a

% Rayleigh fading channel

clear

N = 10^6 % number of bits or symbols

% Transmitter

ip = rand(1,N)>0.5; % generating 0,1 with equal probability

s = 2*ip-1; % BPSK modulation 0 -> -1; 1 -> 0

nRx = [1 2];

Eb_N0_dB = [-3:35]; % multiple Eb/N0 values, ratio of transmission power per bit to noise density

x1 = rand;

x2 = rand;

y1 = rand;

y2 = rand;

for jj = 1:length(nRx)

for ii = 1:length(Eb_N0_dB)

n = 1/sqrt(2)*[randn(nRx(jj),N) + j*randn(nRx(jj),N)]; % white gaussian noise, 0dB variance

h = 1/sqrt(2)*[randn(nRx(jj),N) + j*randn(nRx(jj),N)]; % Rayleigh channel

% Channel and noise Noise addition

d = sqrt((x1-x2).^2 + (y1-y2).^2);

sD = kron(ones(nRx(jj),1),s);

y = d^-2*h.*sD + 10^(-Eb_N0_dB(ii)/20)*n;

% equalization maximal ratio combining

yHat = sum(conj(h).*y,1)./sum(h.*conj(h),1);

% receiver – hard decision decoding

ipHat = real(yHat)>0;

% counting the errors

nErr(jj,ii) = size(find([ip- ipHat]),2);

end

end

simBer = nErr/N; % simulated ber

EbN0Lin = 10.^(Eb_N0_dB/10);

theoryBer_nRx1 = 0.5.*(1-1*(1+1./EbN0Lin).^(-0.5));

p = 1/2 – 1/2*(1+1./EbN0Lin).^(-1/2);

theoryBer_nRx2 = p.^2.*(1+2*(1-p));

% plot

close all

figure

semilogy(Eb_N0_dB,theoryBer_nRx1,’bp-’,'LineWidth’,2);

hold on

semilogy(Eb_N0_dB,simBer(1,:),’mo-’,'LineWidth’,2);

semilogy(Eb_N0_dB,theoryBer_nRx2,’rd-’,'LineWidth’,2);

semilogy(Eb_N0_dB,simBer(2,:),’ks-’,'LineWidth’,2);

axis([0 35 10^-5 0.5])

grid on

legend(‘nRx=1 (theory)’, ‘nRx=1 (sim)’, ‘nRx=2 (theory)’, ‘nRx=2 (sim)’);

xlabel(‘Eb/No, dB’);

ylabel(‘Bit Error Rate’);

title(‘BER for BPSK modulation with Maximal Ratio Combining in Rayleigh channel’);

—————————–

@Martin: Sorry, due to time constraints, its increasingly hard for me to debug the code. Good luck in your algorithm explorations.

@Martin: Sorry, I am not familiar with the system model which you proposing.

Hi~

Who can tell me how to establish

communicate with basestation 1*m error rate, m is receiver

it is randomly pick one user

Please help, I am doing my final year project, i have 1 months only, i still have many tasks havent done yet

@Martin: You question is vague. Typically, one would want to have the error rate as low as possible.

Do you still need the solution to simplify the Pb in Rayleigh Fading? I have the answer.

@Kenny: Thanks. However, Mr. Jose Antonio Urigüen kindly shared his findings. The URL:

http://www.dsplog.com/2009/01/22/derivation-ber-rayleigh-channel/

My program is fixed and reference to your program,I am trying to set a topology of M nodes which is 250,and then randomly pick 2 nodes and repeat 1000 times for each time,pick two nodes as tramitter and receiver and tranmsitt N=10^6 bits to the channel between them and calculate error performance.The error performance of the whole network is calculated and sum up all the d gives the total transmission power.

Now I would like to assume each node has a total pt= 100 so that for each pick up,each node used some of its energy to transmit and then therefore each node is different with the saved energy.i want to calculate the number of nodes that give zero in pt so that to plot a graph how many nodes alive against the time.my program is below: thank for you help

% Script for computing the BER for BPSK modulation in a Rayleigh fading channel

clear

N = 10^6 % number of bits or symbols

% Transmitter

ip = rand(1,N)>0.5; % generating 0,1 with equal probability

s = 2*ip-1; % BPSK modulation 0 -> -1; 1 -> 0

Eb_N0_dB = [-3:96]; % multiple Eb/N0 values

%Generate M nodes

M=250;

x1 = rand(1,M);

x2 = rand(1,M);

y1 = rand(1,M);

y2 = rand(1,M);

% initalize the total power

totalpwr = 0;

for ii = 1:length(Eb_N0_dB)

% Channel and noise Noise addition

n = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)]; % white gaussian noise, 0dB variance

h = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)]; % Rayleigh channel

% randomly pick two nodes and calculate the distance between them

i = randint(1,1,[1,100]);

u = randint(1,1,[1,100]);

d = sqrt((x1(i)-x2(u)).^2 + (y1(i)-y2(u)).^2);

d

totalpwr = totalpwr+d^4;

y = d.^-2.*h.*d.^2.*s + 10^(-Eb_N0_dB(ii)/20)*n;

% equalization

yHat = y./h;

% receiver – hard decision decoding

ipHat = real(yHat)>0;

% counting the errors

nErr(ii,1) = size(find([ip- ipHat]),2);

end

% error performance of the whole network

whole_error = sum(nErr)/N;

whole_error

totalpwr

simBer = nErr/N; % simulated ber

theoryBerAWGN = 0.5*erfc(sqrt(10.^(Eb_N0_dB/10))); % theoretical ber

EbN0Lin = 10.^(Eb_N0_dB/10);

theoryBer = 0.5.*(1-sqrt(EbN0Lin./(EbN0Lin+1)));

% plot

close all

figure

semilogy(Eb_N0_dB,theoryBer,’bp-’,'LineWidth’,2);

hold on

semilogy(Eb_N0_dB,simBer,’mx-’,'LineWidth’,2);

axis([-3 96 10^-5 0.5])

grid on

legend(‘Rayleigh-Theory’, ‘Rayleigh-Simulation’);

xlabel(‘Eb/No, dB’);

ylabel(‘Bit Error Rate’);

title(‘BER for BPSK modulation in Rayleigh channel’);

Hi Kirshna ;

could you give me some hints on how to solve the problem or give guideance on how to solve,thank you very much.your help is much considerable as i really frustrated about this question for long

@Vincent: Sorry, I am not familiar with the scenario which you are proposing (and hence unable to help).

Krishna,

I am glad to receive your reply,however i received your bad news about your reply,but i really don’t understand what is the relationship between your familiarity and my scenairo.Since i really don’t know how to continue write or fix my program,or may you give some any others useful link or reference so that i know more details about what happened?your reply is appreciate.

@vincent; I just realized a typo in my message. I wanted to say – I was not familiar with the scenario which you presented (and hence unable to help).

Khrisna,

Thank you for your explanation. But i think we should consider the Doppler shift in Rayleigh channel. What do you think? Thank you

@Filbert: I have not tried modeling Doppler. Let me try to do some reading and get back to you.

The fading considered in your program is Flat and Slow rayleigh fading where there is no doppler spread i.e. max doppler shift (fm=v/lamda) is zero. Am I right?

@Street hawk: Its flat fading, no doppler spread. However, the channel varies randomly from each symbol. So its fast fading.

I think that if channel is changing every symbol it is mean that it is slow fading (To>Ts, where To- channel coherence time; Ts -symbol duration). Fast fading is when channel is changing several times during symbol (To<Ts).

Khrisna,

Thank you for your response. There is a problem about your Matlab’s script i do not understand. For Rayleigh channel you write:

h = 1/sqrt(2)*[randn(1,N)+j*randn(1,N)]

I notice that the formula for h and n (AWGN) is the same. Why?

I think you just generate random complex number for Rayleigh channel. Is it suppose to have its own model, not just random complex number? Thank you

@Filbert: Well, the Rayleigh channel can be modeled as independent Gaussian random variables on the real and imaginary compnents. The factor 1/sqrt(2) to normalize the power of h to unity. You may look @

http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/

for bit more details.

Further, note that h is multiplied on the transmit signal, whereas n is added to the transmit signal. So the effects of h and n are different.

Hope this helps

Hi Krishna..

Good day. You told that sqrt(1/2) is used to make the variance unity. However, in Matlab. by default, when you use ‘randn’, it creates gaussian random variables with zero mean and variance 1. So does it really necessary to multiply with sqrt(1/2)?

Looking forward to your reply. Thanks.

-Jamil

@Jamil: Recall, we create two random variables, one for real and another one for imaginary. So to make the variance of the total variance to unity, the factor of 1/sqrt(2) is used.

Hi Krishna,

How are you? I hope you are fine. I tried to simulate BER for OFDM transmission via Rayleigh channel. I wonder if in OFDM receiver i need equalization? I did not use equalization and i got poor BER. Thank you

@Filbert: Yes, we need equalization. Without equalization, the random gain and phase error introduced by the channel results in high BER.

hi…….can u plz tell me how can i find the ber of a qpsk in rayleigh fading channel.

@ranjay: Did you want to try modify the Matlab code for BPSK and extend to QPSK case.

Hi Krishna,

How do you consider the effect of interference in this simulation ?. Suppose the interference is another BPSK modulated signal. Can you replace SNR with SINR to compute BER ?.

Your thoughts,

Arvind

@Arvind: Yes, the signal can be corrupted by another BPSK modulated signal. However, I would think that the BER performance will be poorer in that case.

IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. 56, NO. 2, FEBRUARY 2008 213

Soft Handover Overhead Reduction by

RAKE Reception with Finger Reassignment

Seyeong Choi, Member, IEEE, Mohamed-Slim Alouini, Senior Member, IEEE,

Khalid A. Qaraqe, Senior Member, IEEE, and Hong-Chuan Yang, Senior Member, IEEE

Abstract—We propose and analyze in this paper a new finger

assignment technique that is applicable for RAKE receivers when

they operate in the soft handover (SHO) region. This scheme

employs a new version of generalized selection combining (GSC).

More specifically, in the SHO region, the receiver uses by default

only the strongest paths from the serving base station (BS)

and only when the combined signal-to-noise ratio (SNR) falls

below a certain pre-determined threshold, the receiver uses more

resolvable paths from the target BS to improve the performance.

Hence, relying on some recent results on order statistics we

attack the statistics of two correlated GSC stages and provide

the approximate but accurate closed-form expressions for the

statistics of the output SNR. By investigating the tradeoff among

the error performance, the path estimation load, and the SHO

overhead, we show through numerical examples that the new

scheme offers commensurate performance in comparison with

more complicated GSC-based diversity systems while requiring

a smaller estimation load and SHO overhead.

Index Terms—Fading channels, diversity techniques, RAKE

receiver, generalized selection combining (GSC), performance

analysis.

Hi I am looking for generalized selection combining for QPSK (or for anything else). Any help please where can I find it?

@Leona: You can look at the post on Selection Diversity for BPSK case

http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/

Hi Kirshna Thanks for ur reply.

Please see this code , I am confused how i added the two version of signal received from two Taps with one sample delay.Initially I have done IFFT of the signal with Four subcarriers.Please guide me.

x=transdata1;

taps=2;

p1=1;

p2=0.8;

gain1=sqrt(p1/2)*[randn(1,length(x)) + j*randn(1,length(x))];

gain2=sqrt(p2/2)*[randn(1,length(x)) + j*randn(1,length(x))];

channel_ifft=ifft([gain1,gain2],length(x));

n=1:length(x);

delay1=1;

for n=delay1+1:length(x) % Delay of one sample in second tap

x1(n)=x(n-delay1);

end

transdata=channel_ifft.*(x+x1); % confusing here

%————————Addition of noise ——————————-

noise = 1/sqrt(2)*[randn(1,length(x)) + j*randn(1,length(x))]; % Noise

snr = [0:10]; % multiple Eb/N0 values

ps=mean(abs(x).^2); % Power of the transmitted signal

for i = 1:length(snr)

y = transdata + (sqrt(3*ps)*10^(-snr(i)/20))*noise; %Addition of Noise

%————————–Receiver —————————————

recdata=y;

rx0=recdata.*(1./(channel_ifft)); %Equilizatiom

Take care.

Hi!

For those who are looking for a more accurate model of the fading coefficients can read the paper by Christos Komninakis and download the program he wrote at:

http://www.ee.ucla.edu/~chkomn/RayleighFiles/files.html

It’s very accurate for fast fading channels!

@Fang: Thanks for the link

please help me to solve the following writing the matlab code about RCPC code family with memory M = 4 and puncturing period P = 8

and choose overall code rate R = 1/4 with source block size K = 128 bits.

it is very urgent for me. if u could do it for me i would be very thankful to you.

thanks

dean

@dean: Am sorry, I typically help people in their simulation model development rather than code myself. You may look at the post on convolutional coding, Viterbi decoding and hopefully modify them to suit your requirements.

Good luck.

Hi Kirshna

i have three taps and showing below the channel description.

I produced the delay of one and two samples in 2nd and third taps.

can u guide me how I equlize the channel now.

x=transdata1;

taps=3;

p1=1;

p2=0.8;

p3=0.5;

gain1=sqrt(p1/2)*[randn(1,length(x)) + j*randn(1,length(x))];

gain2=sqrt(p2/2)*[randn(1,length(x)) + j*randn(1,length(x))];

gain3=sqrt(p3/2)*[randn(1,length(x)) + j*randn(1,length(x))];

n=1:length(x);

delay1=1;

delay2=2;

for n=delay1+1:length(x)

x1(n)=x(n-delay1);

end

for n=delay2+1:length(x)

x2(n)=x(n-delay2);

end

transdata=gain1.*x+gain2.*x1+gain2.*x2;

%————————Addition of noise ——————————-

noise = 1/sqrt(2)*[randn(1,length(x)) + j*randn(1,length(x))]; % Noise

snr = [0:10]; % multiple Eb/N0 values

ps=mean(abs(x).^2); % Power of the transmitted signal

for i = 1:length(snr)

y = transdata + (sqrt(3*ps)*10^(-snr(i)/20))*noise; %

Thanks.

@mabrarbari: If I may, I find the code snippet which you have provided bit confusing. Some questions:

1. Why is each path gain1, gain2, gain3 having length of the input vector. I would have thought that each path is a single element.

2. Instead of using the for loop’s you could have used a convolutional operator.

In the code snippet, your channel is not a flat fading channel. Need to formulate the multipath channel as matrix and solve the equation.

Thanks Kirshna;

Yes I have created frequency selective fadded channel in which there is a one sample delay in taps.if i made power of tap 2 and 3 zero it becomes flat fadding and zeroforcing equlizer can be used. but if these two taps are used then how i equilize the channel.

1- Can I go for normalize the channel.

2- can u give me any othe alternate to represent frequency selective channel.

Take care.

@Abrar: May I try to provide you a simple Matlab code snippet for equalization in the case of a 3 tap channel.

clear all;close all

x = 2*(rand(1,7)>0.5)-1 % input sequence

h = [ 1 2 3]; % channel

op1 = conv(x,h); % output of channel

hM = toeplitz([h(1) zeros(1,length(x)-1) ], [h zeros(1,length(x)-1) ]); % writing the channel using toeplitz matrix

op2 = x*hM ; % finding the output using the matrix multiplication way

op1 – op2; % verifying that both the outputs are the same. This confirms that toeplitz way of representing h is correct.

xHat = op2*pinv(hM) % dividing received symbol with pseudo inverse of the channel to perform equalization

The code performs the simple zero forcing equalization. To represent mathematically, the convolution operation is represented as product a vector and a matrix. At the receiver, we divide by the matrix to get the equalized channel estimate. This is described in detail in Chapter 3.4.3 of Digital Communications: Fundamentals and Applications by Bernard Sklar

Hope this helps.

Dear sir

At first i want to give you thank for this interesting blog…..i go through ur posts of matlab code these are really helpfull for me. As i am a new user of matlab i faced some problem during coding. i will be very grate full to you if u help me to coding this using matlab:

Calculate path loss between 3 wireless node.By considering the following parameters: BPSk modulation, AWGN channel, transmitter power 15 dB, frequency 2.4 Ghz, patloss exponent 3, refercence distance 1 km. Assume that one node A is transmitting and 2nd node B is receiving signal,as the 3rd node C is located near to the node B it also recieve some signal . i have to calculate the recieved signal to noise ratio at node B and C.

It would be nice if u help me by giving some code or any suggestion. I can not understand how to start ……Please reply me…..i am waiting for your reply……

with regards

shoumi….

@Shoumi:

From your problem statement, it seems like you are looking for finding the path loss from between A to B and A to C.

Given that you know the transmit power, path loss exponent, and the distance, it should be reasonably easy to find the received signal power. Once the received power is known, based on the noise floor of each receiver (ideally depending on noise bandwidth), SNR can be computed.

The link in wiki on path loss might be helpful

http://en.wikipedia.org/wiki/Log-distance_path_loss_model

Hope this helps.

Thank u for ur reply. Yes i calculate the received power for each node. Now for calculating the SNR for each receiving node ..i have to know the noise at each each receiver. Please inform me what do you mean by noise floor…is it AWGN noise???..

Awgn noise is a complex value..and the recived signal is in dBm unit. for calculating the SNR these unit must be same…how can i conver the awgn noise into dBm..Please help me regarding this….

@shoumi: Yes, the noise is complex. But the noise power can be represented in dBm. For example, assume the receiver bandwidth is 20MHz. The noise floor in dBm = -174+10*log10(20e6) = -101dBm.

You may find bit more details on the equation @ http://en.wikipedia.org/wiki/Thermal_noise

Hope this helps.

Hi

Kirshna

I am away some time from country , now come back and start working.

As u made equilization by just dividing the received signal through channel, any other way to equilizing.

if channel become zero then wht happen if we used equlization by dividing.

I have defined the channel as

x=transdata1;

taps=3;

p1=1;

p2=0;

p3=0;

gain1=sqrt(p1/2)*[randn(1,length(x)) + j*randn(1,length(x))];

gain2=sqrt(p2/2)*[randn(1,length(x)) + j*randn(1,length(x))];

gain3=sqrt(p3/2)*[randn(1,length(x)) + j*randn(1,length(x))];

n=1:length(x);

delay1=1;

delay2=2;

for n=delay1+1:length(x)

x1(n)=x(n-delay1);

end

for n=delay2+1:length(x)

x2(n)=x(n-delay2);

end

transdata=gain1.*x+gain2.*x1+gain2.*x2;

%————————Addition of noise ——————————-

noise = 1/sqrt(2)*[randn(1,length(x)) + j*randn(1,length(x))]; % Noise

snr = [0:10]; % multiple Eb/N0 values

ps=mean(abs(x).^2); % Power of the transmitted signal

for i = 1:length(snr)

y = transdata + (sqrt(3*ps)*10^(-snr(i)/20))*noise; %

At receiver how I equlize it.

@Abrar: If the channel becomes close to zero, equalization by dividing by the channel results in noise amplification.

Sir,

Thank you for giving me the hint.

As you said, 2.TO NORMALIZE THE RECEIVED POWER,

% Calculation of normalization factor to

%make the average power of output signal to 1

norm_fact = sqrt(sum(tap_weights_ln.^2));

1.TO NORMALIZE EACH TAP WEIGHT:

I divided taps_weights by sqrt(length(tap_Weight_ln)).

tapWtdB = [ -5 -10 -15 -20 -25];

tap_weights_ln = 10.^(tapWtdB/20);

taps=1/sqrt(2)*(randn(10^4,length(tapWtdB)) + j*randn(10^4,length(tapWtdB)));

taps_weights = taps.*repmat(tap_weights_ln, 10^4,1);

norm_each_tap_weight=taps_weights/sqrt(length(tap_weights_ln));

norm_fact = sqrt(sum(tap_weights_ln.^2));

taps_norm = 1/norm_fact*norm_each_tap_weight;

z=mean(abs(taps_norm ).^2);

est_tap_weight_db=10*log10(z)

total_path_gain_of_channel=sum(mean(abs(taps_norm).^2,1))

o/p:

est_tap_weight_db = -8.6703 -13.5813 -18.6063 -23.6805 -28.6429

total_path_gain_of_channel = 0.1991

Here, the total path gain of the channel is not close to one.(very less comparatively)

(but the estimated tap weight is close to the original tapweight(dB) ,ofcourse,

its value is less because of normalisation)

since,total path gain of the channel is not close to 1.

Is it acceptable one Sir?

Waiting for your reply.

Thank you for your response sir.

In the last post You wrote:

Incase of flat fading,

Matlab/octave code snippet to illustrate:

tapWtdB = -15

h = 10^(tapWtdB/20)*1/sqrt(2)*[randn(1,10^4) + j*randn(1,10^4)];

estTapWt = 10*log10(h*h’/length(h))

estTapWt is obtained very close to -15dB.

I think length(h) in estTapWt is Normalisation factor.

In case of Multiple taps

tapWtdB = [0 -20 -21];

tapWtdB_linear=10.^(tapWtdB/20);

h = 1/sqrt(2)*[randn(10^4,length(tapWtdB)) + j*randn(10^4,length(tapWtdB))];

taps_weights = h.*repmat(tapWtdB_linear, 10^4,1);

z=mean(abs(taps_weights/sum(tapWtdB_linear)).^2);estTapWt=10*log10(z);estTapWt

estTapWt =

-1.5181 -21.4893 -22.5174

PowerdB level is someWhat close to the original tapWtdB. This happens becuase of the Normalisation (1/sqrt(2)) in h line in tha Matlab code.

But, if the tapWtdb is [-15 -20 -21];

tapWtdB = [-15 -20 -21];

tapWtdB_linear=10.^(tapWtdB/20);

h = 1/sqrt(2)*[randn(10^4,length(tapWtdB)) + j*randn(10^4,length(tapWtdB))];

taps_weights = h.*repmat(tapWtdB_linear, 10^4,1);

z=mean(abs(taps_weights/sum(tapWtdB_linear)).^2);estTapWt=10*log10(z);estTapWt

estTapWt =

-6.3169 -11.3482 -12.2874

Here the difference between the consecutive taps power level (db) is satisfied. estTapWt value is not closer to original original value.

My question is

Whether the first tap necessarily be 0dB or someother dB value.If we can have some other value of power dB for the first tap then Can U please correct the Matlab code inorder to get original tapWtdB.

@Ramya: In the case of multiple taps, there are two aspects:

1. Firstly, one would want to normalize each tap such that the relative power is as per the tap weight. This can be accomplished by scaling each tap as in the example which I have shown earlier.

2. Secondly, to normalize the received power, one may want to make the average power from many such channel realizations to unity.

Hope this helps.

Sir,

I understood & got the answer for the last two questions Which i have posted previously.

Please Check Whether the answers are right.

question :2

After modelling has absolute value of complex gaussian R.Vs, from the Power delay profile,

we have the Tap weights in dB. Since the variance of rayleigh R.V. h is 1.

Tap_weight_db here in our case is = 0 dB.

In linear scale its value is 1.

Tap_weight = h*(10^(Tap_weight_db)/10 )=h*1=h

As per theory,the received signal r is written as

r=(Tap_weight * s) + n;

question 3:

whether the tap (in flat fading)will occur at 0 microseconds 0r 0 ns i.e., without delay.or shall the tap occur after some particular delay?

The tap doesn’t actually occur at 0 microseconds.

These delays are relative to the delay of the tap to which demodulator synchronizes.

But my first question was

In BER_simulation_BPSK_Rayleigh_channel Matlab Code:

given two lines seems to be same

n = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)];% WGN,0dB variance

h= 1/sqrt(2)*[randn(1,N) + j*randn(1,N)];% Rayleigh channel

From theory, h=sqrt(X^2+Y^2)=abs(X+jY);

Why ‘abs’ have not been included ?

Waiting for your reply.

Sir,

I saw your Matlab code for BER_simulation_BPSK_Rayleigh_channel.

As Rayleigh channel is modelled as |h|

|h|=1/sqrt(2).sqrt(X^2+Y^2)=1/sqrt(2).abs(X+jY)

if X & Y are two gaussian R.Vs with zero mean and Variance of 1/2 each.

In the Matlab code,

% white gaussian noise, 0dB variance

n = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)];

% Rayleigh channel

h= 1/sqrt(2)*[randn(1,N) + j*randn(1,N)];

My question are:

1.Both the above lines are same ( even though complex guassian R.V.s were used)

Whether the line given below is correct.

Rayleigh R.V. h=1/sqrt(2).*abs([randn(1,N) + j*randn(1,N)]);

if it is not correct kindly explain me, why ‘abs’ should be used?

2.After modelling h, from the Power delay profile,

we get the Tap weights in dB .(In our case as theory suggest, since flat fading has been used no. of taps will be only one)

Shall we convert tap weights in dB to linear scale and multiply with h and txmitted signal s (since inourcase only one tap is being considered for flat fading, convolution becomes multiplication) and add gaussian noise(which is implemented by mutiplying with Guassian Noise with linear scale Eb_No)

3. For flat fading, since the no. of taps is one. my question is whether that particular will occur at 0 microseconds 0r 0 ns i.e., without delay.or shall that tap can occur after some particular delay.

Am i correct in this regard.

sorry for troubling U sir, as i am struggling to get hold of things right.This is my first sincere effort in learning digital comm conceps clearly.

waiting for your reply.

@Ramya:

1. With the Rayleigh channel, it is assumed that the phase distortion caused by the channel can be uniformly distributed between 0 to 2pi. Hence we model it as independent complex Gaussian random variable for real and imaginary component

Please look at http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/ for bit more details.

2. I agree that since we are assuming flat fading, the power delay profile has only one tap. We should scale the variance of channel h in accordance with the gain of the power power delay profile. Note that you should use 10^(-tapWeightdB/20) and not 10^(-tapWeightdB/10) . A quick Matlab/octave code snippet to illustrate:

tapWtdB = -15

h = 10^(tapWtdB/20)*1/sqrt(2)*[randn(1,10^4) + j*randn(1,10^4)];

estTapWt = 10*log10(h*h’/length(h))

3. Typically there will be small delay. However, typical receivers will be having a packet detection circuit which can estimate the start of the symbol. Hence for modeling purpose we can assume that there is zero delay.

hi there,

First, thanks for sharing the understanding of this subject. Regarding the missing part, one can find it on Page 823 in the book, CDMA System Engineering Handbook.

It is all about some mathematic manipulations, while the Q(.) function is used to replace the erfc function in the expression.

@Alex: Thanks. Let me try to get my hands on the book.

@Shefo666: I think the ZF equalizer is optimal in SISO systems.

With a simple Matlab/Octave script may I try to show that ZF and MMSE equalizer gives the same BER. Infact, am not sure whether we can call it as ZF/MMSE equalization – as there is no interference terms.

% Matlab/Octave code snippet for comparing zero forcing and MMSE equalization for SISI

clear

N = 10^4 % number of bits or symbols

ip = rand(1,N)>0.5; % generating 0,1 with equal probability

s = 2*ip-1; % BPSK modulation 0 -> -1;

Eb_N0_dB = 3;

n = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)]; % white gaussian noise, 0dB variance

h = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)]; % Rayleigh channel

% Channel and noise Noise addition

y = h.*s + 10^(-Eb_N0_dB/20)*n;

% equalization

yHat_zf = y./h; % zf equalization

yHat_mmse = 1./(h.*conj(h)+ 10^(-Eb_N0_dB/10)).*conj(h).*y; % mmse equalization

% receiver – hard decision decoding

ipHat_zf = real(yHat_zf)>0;

ipHat_mmse = real(yHat_mmse)>0;

% counting the errors

nErr_zf = size(find([ip- ipHat_zf]),2)

nErr_mmse = size(find([ip- ipHat_mmse]),2)

Can see that nErr_zf and nErr_mmse is the same. Hope this helps.

Hi Krishna

could you please refer :

http://books.google.com/books?id=3DY6OAIGu0kC&printsec=frontcover&source=gbs_v2_summary_r&cad=0#v=onepage&q=&f=false

(Introduction to Space Time Communications)

Which gives a comparison of ZF, MMSE,etc for SISO

It shows MMSE is better than ZF for SISO.

Regards

Page 141

@WirelessNewbie: Sorry, the page 141 is not available from the link you provided. Let me try to get the book from the library, and I will respond.

One query: Is the claim that “ZF is better than MMSE for SISO” for a flat fading channel?

Dear Krishna Pillai,

Thanks for the very interesting blog ..

I am wondering whether it makes any sense to apply MMSE and ML equalization on SISO communication systems in order to get over the ZF equalization major problem of noise enhancement (when channel is in deep fade).

In such a case, how is the theoretical BER performance of BPSK (for instance) over flat rayleigh (given by proakis, Digital communications) affected ? can u reccommend any papers that calculate the theoretical BER performance of ML and MMSE equalization over flat fading .. Most papers i got assume ZF equalization when it comes to SISO systems !

@Tahmid: Moving the simulation to 4-QAM should be reasonably simple. First thing is that we will be sending information on both real and imgainary arms. Correspondingly, use the information from both real and imaginary arms while demodulating.

You can use the Matlab model provided here

http://www.dsplog.com/2007/11/06/symbol-error-rate-for-4-qam/

as reference.

Hope this helps.

@Muhammad: Thanks. Glad that you are finding http://www.dsplog.com useful for your studies. Happy learning!

@Abrar: Thanks for the comment. Sorry for the delayed response. I was travelling.

Typically for channel estimation, we send a known sequence (called preamble). Using the knowledge of the known sequence and what we received, we can estimate the channel.

For eg, in a multi-carrier OFDM system, the system model can be written as

Y = HX + N where

Y – is the received symbol on subcarrier k

H – is the channel on subcarrier k

X – is the transmitted symbol on subcarrier k

N – is the noise on subcarrier k

Since X is known, the estimate of H is,

H^ = Y/X

Does this help?

Hi Krishna,

I am trying to follow this tutorial and extend it further to simulations for 4-QAM over Rayleigh fading channels. Could you please tell me how the code to apply the fading channel would change from the BPSK case. thanks

I am Muhammad from Pakistan,

I was in need to compute the Bit error rate of BPSK and QPSK in AWGN and Rayleigh fading channel.

After finding this document, i was very happy as it improve my knowledge and releive my burden.

Thanks for all

Hi Kirshna

thanks for ur reply and sorry for the t I am saying thanks very late bec i am busy in my exam.Hoping U r fine, may GOD bless u. I am doing Postgraduate diploma in Massey University New Zealand.

I am doing simulation on MC-CDMA. So I need ur help in modeling frequency selective channel. if u have any code about frequency selective chanel estimation, it will helpful for me. I am try to do so I will send u code if I have have succeded to some extent.

waiting for ur reply.

Thanks nad Take care.

abrar

@Sky: 1/sqrt(2) is for normalizing the channel gain to unity. I think it is required, as the derivations are for unity channel gain.

Hi,

Could you please explain on how you get

h = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)]; for Rayleigh?

Coz I think, it should be just randn(1,N) + j*randn(1,N)

Thanks

@Leon: Thanks. The derivation has been pending for a while from my side. Let me try to go over it by over the week.

Hi Krishna,

Jussi Poikonen is correct, you can determine the double integral in this fashion to obtain the desired result.

@davis: Sorry for the delayed response, I was down with viral fever.

I had written a post on BFSK some time in the past.

URI: http://www.dsplog.com/2007/08/30/bit-error-rate-for-frequency-shift-keying-with-coherent-demodulation/

The post describes BFSK in AWGN case. You can convolve the transmit signal x with the impulse response of the channel (lets say h) before adding the noise n.

@Abrar: Sorry for the delayed response, I got bogged down by viral fever.

Yes, this is a simple form equalization, which basically tries to make the effective channel as impulse. Since this a simple one tap channel, the equalization becomes just a division. This can be classified as a simple zero-forcing linear equalization. You can read more about equalization in Chapter 8 of

Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt.

Hope this helps.

Hello Kirshna

I doing project on bfsk with rayleigh multipath chanel

i wish to ask you if I use not an envelop but modulated fsk

0->cos(f1*t) and 1->cos(f2*t)

how can I simulate rayleigh multipath chanel for bfsk in this case?

thank you

Davis

Hi Kirshna:

hoping u r fine.

I just want to ask about equlization which u have done here just by dividing by channel. is it another method of equlization.

can u give some explation about equlization or refer some material about it, bec, I want to understand it.

Tahnks and take care.

abrar

yES SIR

IT IS USED IN OPTICAL COMMUNICATIONS.

@Sairam: Ok, thanks. Are there some practical systems taking advantage of non-orthogonality to improve spectral efficiency. Can you please point me to the same.

it can be obtained if the orthogonal condition is not satisfied. But the derivations become more complicated. the advantage is that it provide high spectral efficiency. ok sir thank u for ur fast response

@sairam: By non-coherent, hope you meant a scheme like DBPSK. Correct? I have not tried to discuss the deivation for BER with non-coherent modulation in any blog post yet. Something for me to look at.

Can you give an example for non-orthogonal modulation?

Sir

How to derive the probability error of Non-coherent Non-Orhogonal modulation. TRhanking you sir!!

@Jussi Poikonen: Thanks for the guidance. Let me try to go over the steps per your suggestion.

You will find that the missing steps in the derivation above can be completed by writing open the final (double) integral, and changing the order of integration. That is, integrate first with respect to gamma from 0 to x^2 (where x is the argument of the complementary error function), and then with respect to x from 0 to infinity. The integration problem is finally reduced to single integrals from 0 to infinity over zero-mean Gaussian distributions.

@Varun: Probability of error in Rayleigh fading channel.

@nano686: This post does the simulation with a Rayleigh channel. Further, there are quite a few posts on Rayleigh channel

URI: http://www.dsplog.com/tag/rayleigh/

Hope this helps.

Hi,

I don’t know simulate Rayleigh channel by Matlab. Can you help me?

Thanks and Regards.

@ Pillai

Which equation ? Is it the probability of error in a Rayleigh fading channel or the probability of error in a AWGN channel ? Let me know so that i can come up with a solution if i know

@Varun: Yes, here the random variable is instantaneous SNR. Btw, do you know how to complete the proof? I have not figured that out. Kindly let me know, if you can prove the equation.

Thanks.

@Krishna Pillai

Thanks a lot for giving me a reply. So this is similar to finding the mean of a random variable right !!! Here the random variable is instantaneous SNR right !!! Comment on this statement Sir !!!

@Varun: Because of Rayleigh fading, the instantaneous SNR for a particular channel realization is |h|^2*Eb/N0. So, to find the average bit error rate, we need to average over all channel realization. Hence we average over the conditional error probability. Does that help?

Hi,

In order to find probability of error in a rayleigh fading channel we will average the conditional probability of error of the particular modulation scheme in AWGN channel over the rayleigh probability density function. Why we do this ? Please give explanation.

Thanks and Regards,

V Varun.

thanks krishna

actually i have to show the performance comparison of different strategies of cooperative communication…..

different combining techniques at receiver may also be considered but not compulsory…..

it is quit similar to relaying…..

i have sent some matlab codes on krishna at dsplog com

it shows some error which i fail to removed.

please help me in removing these errors……

@khurshid: Sorry, I can help with your coding but not do coding for you.

As a first step, you may try to create a packet error rate simulation environment with the above mentioned parameters. Are you using OFDM?

hello!

i go through your posts of matlab code, these are very helpful. please help me to solve the following writing the matlab code.

Assume BPSK modulation with equal user transmit

powers. The channels between users (interuser) and to the

base station (uplink) are mutually independent with flat,

Rayleigh fading that remains constant over each source block

(quasi-static fading). Each receiver has perfect CSI and employs coherent detection, and we quantify the quality of each channel by its average received SNR. We use the RCPC code family with memory M = 4 and puncturing period P = 8

and choose overall code rate R = 1/4 with source block size K = 128 bits.

if more detail is needed please inform me…

it is very urgent for me. if u could do it for me i would be very thankful to you.

thanks

khurshid

Thank you! Kirishna, hope fully ur script is very very help full for my work but still i have doubt on how to write a matlab script to simulate a frequency selective fading channel.

for example a channel of 120 Hz Doppler frequency.

You may use the BEM channel model for it.It is great to simulate the Jake’s channel model.

@Ok, thanks

What is BEM channel ?

@Humayun: I have not tried to simulate BEM channel model

@Lealem: Yes, the above simulation is for a flat Rayleigh fading channel.

In the OFDM case, though the total channel is a frequency selective channel, the channel experienced by the individual subcarriers will still be flat fading channel. Am expecting that a rayleigh channel having a duration less than the cyclic prefix, the expected BER performance in OFDM will be the same as the results obtained in the above plot. I will run simulations to check and confirm.

您好 请问最后的误码率公式是怎么化简得到的？能给指点一下参考书吗 谢谢

@JERRY: You may refer to the post on Derivation of BER for BPSK in Rayleigh channel for the proof.

http://www.dsplog.com/2009/01/22/derivation-ber-rayleigh-channel/

Further, this is discussed in Chapter 11.3 of Digital Communications by Barry, Lee, Meesserschmit.

Hi mate… i went through your code and it was great. Im just learning about mimo equalization and all. Will you be able to help me generate a code for 2×2 MIMO equalised using svd?please mail me if you can.

thanks in advance.

cheers.

@balaji: I have written some posts on MIMO. You may look up them @ http://www.dsplog.com/tag/mimo

In those posts, you may use the MIMO with Zero Forcing Equalizer as the most simplest example. When you said, equalizing MIMO by SVD, did you mean pre-compensing the transmitted matrix by V matrix such that the channel experienced by the MIMO channel is diagonal?

[U S V] = svd(H)

It’s useful~~Thank you!

hi,

i need few matlab codes oF MIMO AND SISO SYSTEMS ON AWGN FADING CHANNEL AND RICIAN FADING CHANNEL.

MANY THANKS

@Ali: You may refer posts @ http://www.dsplog.com/tag/mimo

These are flat faiding Rayleigh channel

Hi Krishna ,

Thanks for the reply ! Well as we have detailed information about rayleigh fading channel on this link http://www.dsplog.com/2008/08/10/ber-bpsk-rayleigh-channel/

can i have a similar information about rician and awgn, as i have already gone through the link which you have posted me unfortunately i couldn’t find what exactly i am looking for .

Thanks

@Ali: The BER for BPSK in AWGN can be found @

http://www.dsplog.com/2007/08/05/bit-error-probability-for-bpsk-modulation/

Sorry, I have not written anything on Rician channel as on today.

@Ali: The BER for BPSK in AWGN can be found @

http://www.dsplog.com/2007/08/05/bit-error-probability-for-bpsk-modulation/

Sorry, I have not written anything on Rician channel as on today

I have seen the simulation code which u used to simulate a Rayleigh fading channel. But i would like u to ask a hit how to simulate a frequency selective Rayleigh fading channel for OFDM application. I think the above simulation out put is a flat Rayleigh fading.

Thank you!!!

One who is interested in including his/her name in 3 research conferences in international IEEE conference should reply on email address me_researcher@yahoo.com. Both papers are related to mobile communication. Total of 3 authors list will be included in each Paper. 2 author names have already been included. Interested candidates who want to include his/her name at 3rd position will be required to pay for the registration fee.

Matlab Code as well as the full paper will be sent to the individual after acceptance of paper from the conference.

Fee submission will be through Freelancer and elance. Candidate name will be included in Paper after milestone payment is released by him/her. Milestone payment is one which is in the custody of broker (Freelancer or Elance authority) and not in the custody of either party. After the paper is accepted and the client is conformed about the acheivement, only then he will be allowed to pay. Time is short so the policy of 1st come 1st serve will be entertained.

@giro82: In this simulation, the channel changes every symbol

Khishma, please, BER BPSK OFDM in raileigh channel should be better than BER OFDM in Rayleigh channel? Explain? Please. Thank you very much.

Khishma, please, BER BPSK OFDM in raileigh channel should be better than BER BPSK in Rayleigh channel? Explain? Please. Thank you very much.

Reply

@FERNANDO: Due to the cyclic prefix, OFDM can handle mutipath channels better

http://www.dsplog.com/2008/02/17/cylcic-prefix-in-orthogonal-frequency-division-multiplexing/

http://www.dsplog.com/2008/08/26/ofdm-rayleigh-channel-ber-bpsk/

{ 8 trackbacks }