clear?all;
clc

lamda_FSO=1;

sigma_R2=10;

a1=（exp（0.49*sigma_R2/（（1+1.11*sigma_R2^（6/5））^（7/6）））-1）^（-1）
b1=（exp（0.51*sigma_R2/（（1+0.69*sigma_R2^（6/5））^（5/6）））-1）^（-1）
Gam_a1=gamma（a1）;
Gam_b1=gamma（b1）;

a=5;
b=1;

Gam_a=gamma（a）;
Gam_b=gamma（b）;

?x=1:1:100;
%%%%%%%%%%%%%%%%%%%????????負指數模型的pdf????
f1（x）=（lamda_FSO.*exp（-lamda_FSO.*（0.1.*x）））;

%%%%%%%%%%%%%%%%%%%????????K分布模型的pdf
f2（x）=（（2.*（（a.*b）^（（a+b）/2））.*（（0.1.*x）.^（（a+b）/2-1））.*（besselk（a-b2.*（（a.*b.*（0.1.*x））.^（1/2））））./（Gam_a.*Gam_b）））;

%%%%%%%%%%%%%%%%%%%????????Gamma-Gamma模型的pdf
f3（x）=（（2.*（（a1.*b1）^（（a1+b1）/2））.*（（0.1.*x）.^（（a1+b1）/2-1））.*（besselk（a1-b12.*（（a1.*b1.*（0.1.*x））.^（1/2））））./（Gam_a1.*Gam_b1）））;

plot（xf1（x）‘-b‘xf2（x）‘-k‘xf3（x）‘-m‘‘LineWidth‘1.5）;

%?plot（xf1（x）‘-b‘‘L