Matlab code for channel models for FSO

[Remo De Suza]

Hi All.
Does somebody have MATLAB code for any of the following channel models.
1. log normal
2. gamma-gamma
3. negative exponetial
4. diffuce( indoor)
5. double genralized gamma -gamma

I would appreciate your help. Looking forward for your responses. Thanking you.

Regards
Remo

[Remo De Suza]

Hi All.
Please can someone provide me with the matlab code.

Regards
Remo

[Remo De Suza]

For any of the above channel models. I would appreciate that. Thanking you!

[SAHIL SINGH]

Hi Remo,

I think we already had a discussion on this.. I would suggest you to find it on the internet… I did try to find it out for you but i am not sure whether it will work or not.. Try the following links
https://www.eecis.udel.edu/~bohacek/Classes/SensorNets2007/propagationhomework.pdf

Log Distance Path Loss or Log Normal Shadowing Model

Hope it helps

Regards
Sahil Singh

[SAHIL SINGH]

Hi Remo,

Do let me know if it was useful…
All the best

Regards
Sahil Singh

[Remo De Suza]

Hi Sahil.
I searched but did not get any substantial stuff about it. Anyone working with FSO will surely help and provide me the required help. Anyway thanks for your concern.
Regards
Remo

[SAHIL SINGH]

Hi Remo,

Well I hope you get your query resolved…
All the best

Regards
Sahil Singh

[Remo De Suza]

Hopefully.
Anyway Thanks Sahil for showing so much concern.
Regards
Remo

[Rajguru M. Mohan]

Hi Remo,
Here is the MATLAB code for LOg-Normal:
x = (10:1000:125010)’;
y = lognpdf(x,log(20000),1.0);

figure;
plot(x,y)
h = gca;
h.XTick = [0 30000 60000 90000 120000];
h.XTickLabel = {‘0′,’$30,000′,’$60,000’,…
‘$90,000′,’$120,000’};

Check this on Matlab software.

[Remo De Suza]

Hello Rajguru Mohan.
Thank you so much for the code. I was desperately looking for it.
You have posted two code. Are they different or it is the same code you have posted in two comments. Looking forward for your response. Cheers
Are you also working on FSO?

Thanks and regards
Remo

[Rajguru M. Mohan]

Hi ,
As we know that the exponential distribution is a special case of the gamma distribution.

mu = 1:5;

y = gampdf(1,1,mu)
y =
0.3679 0.3033 0.2388 0.1947 0.1637

y1 = exppdf(1,mu)
y1 =
0.3679 0.3033 0.2388 0.1947 0.1637
So, you can try this.

[nagesh kuknoor]

Hi Rajguru,

Can we directly make use of matlab code in optisystem?

Regards

[Domail Singh]

Hello Nagesh,

We can’t directly write code in optisystem. We have to co-simulate which means write code in either Matlab or C and then interface the code. Hope i answered your query.
Thank you
Regards
Domail

[nagesh kuknoor]

Hi Domail,
Thank you for the information. it is really helpful. I checked the binary switch example and i got some idea.

Regards

[Domail Singh]

Hello Nagesh
You are welcome. I am glad that it helped.

With Regards
Domail

[Remo De Suza]

Hello Rajguru.
Of course we can cosimulate with Matlab. Try finding out a video on youtube about it. That would help.

Regards
Remo De’Suza

[nagesh kuknoor]

Hi Remo,
Thanks Remo for the help. It was very helpful. Thanks again.

Regards

[Remo De Suza]

Hi Nagesh,
You are welcome
Regards
Remo

[Domail Singh]

Hello Rajguru,

Are these codes readily available ? I wonder.

[Remo De Suza]

Hello Domail.
It is not like that. You need to compile these codes and then interface. Hope you got my point.

Regards
Remo De’Suza

[Domail Singh]

Hi Remo
Thank you for that. I got the point. It was a bit confusing.

Regards
Domail

[Rajguru M. Mohan]

Hi remo,
clear all
close all
% subplot (nrows,ncols,plot_number)
x=0:.1:2*pi; % x vector from 0 to 2*pi, dx = 0.1

subplot(2,2,3) % plot negative exponential function
plot(x,exp(-x));

Plotting Negative exponential is one most simplest one in MATLAB.
Seeking your response.

[Remo De Suza]

Hello Rajguru.
I appreciate you for the Help. I shall try running thhis code. Hope it works.
Regards
Remo De’Suza

[Ankita Sharma]

Hi Rajguru,

First of all welcome to the forum…Thanks for providing the code…Surely this should work in the matlab software.. I hope it helps Remo… All the best..

Regards
Ankita

[Remo De Suza]

Thank you Ankita for the kind words. I appreciate your concern.

Regards
Remo De’Suza

[SAHIL SINGH]

Hi Rajguru,

Thanks for the code.. I hope this will be of help to remo…

Regards
Sahil Singh

[Remo De Suza]

Thank you all for your help. I appreciate your efforts especially Rajguru Mohan.

Regards
Remo

[prakash jat]

hello friends
i will post all of the matlab code some time later( 1 or 2 days)

[Remo De Suza]

Hello Prakash Jat.
Thank you so much for the help you have shown. I will wait it alright.
Do you work with FSO design too?

Thanks & regards
Remo

[prakash jat]

if someone do interfacing with optisytem then please post again so that everyone can get help

[prakash jat]

hello friends
if someone do interfacing with optisystem then post it again for everyone so that everyone can get help

[prakash jat]

hi remo
i am really happy you are working also on FSO.
1. we can mitigate power requirement by well-developed laser type( class of laser), linewidth, as well as channel estimation , coherence width type of modulation etc.
2. again second issue is that turbulence( which is more important than attenuation ) , which can be mitigated by aperture averaging , mimo,
divergence of beam and many more…………..

[Remo De Suza]

Hello Prakash.
Thank you so much buddy. I am really very very thankful to you for all the help you have provided to me.
What area are you exactly working on for FSO?
I am trying to implement some different channel model to improve the performance. Which model is Optisystem implementing when we design any FSO system?
Also what is the value of beam divergence you are taking for your design?
Looking forward for your response. Thanking you!

Thanks and regards
Remo

[prakash jat]

i am getting trouble in posting , why this is happning so i don’t know( net is fast), i replied you many time but my content is not posting well

[prakash jat]

matlab code for pdf calculation of by gamma-gamma channel model
%The Gamma-gamma pdf
clear
clc
alpha=4.2;
a=alpha;
beta=1.4;
b=beta;
k=(a+b)/2;
k1=a*b;
K =2*(k1ˆk)/(gamma(a)*gamma(b));
I=0.01:0.01:5;
K1=I.ˆ(k−1);
Z=2*sqrt(k1*I);
p=K.*K1.*besselk((a−b),Z);
plot(I,p)
xlabel(‘Irradiance, I’)
ylabel(‘Gamma gamma pdf, p(I)’)

[prakash jat]

gamma gamma channel is used by optisystem

[Remo De Suza]

Hi Prakash.
I am really thankful to you. What is the source from where you read optisystem uses gamma-gamma channel model?
Looking forward for your reply.

Thanks & regards
Remo

[Ankita Sharma]

Hi Prakash jat,

Thanks for the code…

Regards

[Remo De Suza]

Thank you All for the help and i appreciate it.

[nagesh kuknoor]

Hi All,
very helpful and constructive discussion indeed. All the best Remo!
Regards

[Remo De Suza]

Hi Nagesh,
Thank you for your concern.

Regards
Remo

[Ankita Sharma]

Hi remo,

Did the code work properly ?

Regards
Ankita

[Manoj Kumar]

Hi Remo,

Hope you got the desired information you were looking for…

Regards

[Domail Singh]

Hi Remo
Hope it all helped you. Thank you.
Regards

[Rahul Tiwari]

Hi all,

A very meaningful discussion… Hope remo got the desired help…

Regards
Rahul Tiwari

[GopiKrishna S]

Do any one has the diffuse channel model(Optisystem based matlab code)??

Thanks and Regards,
S.Gopikrishna.

[Author]

Posts