%BER-SNR
SNR?=?10:1:60;
n?=?size（SNR2）;
BERC?=?zeros（3n）;
BERL?=?zeros（3n）;
Lambda?=?1.55*10^（-6）;
k=2*pi/Lambda;
Cn?=?[0.51*10^（-13）?1.13*10^（-13）?2.10*10^（-13）];
L=[1000?3000?5000];

Rytov?=?1.23?*?Cn?*?k^（7/6）*?L（1）^（11/6）;
sigma?=?sqrt（Rytov）;
a=zeros（13）;b=zeros（13）;
for?i?=?1:3
????a（i）?=?（exp（0.49*Rytov（i）/（1+1.11*sigma（i）^（12/5））^（7/6））-1）^（-1）;
????b（i）?=?（exp（0.51*Rytov（i）/（1+0.69*sigma（i）^（12/5））^（5/6））-1）^（-1）;
????s=[‘[（1-‘num2str（a（i））‘）/2（2-‘num2str（a（i））‘）/2（1-‘num2str（b（i））‘）/2（2-‘num2str（b（i））‘）/21]‘];
????for?j?=?1:n
????????z?=（（2/（a（i）*b（i）））^2）*10^（SNR（j）/10）;
????????BERC（ij）?=?2^（a（i）+b（i）-3）/（pi^（3/2）*gamma（a（i））*gamma（b（i）））*double（feval（symengine‘meijerG‘‘2‘‘4‘s?‘[0?1/2]‘?num2str（z）））;?
????end
end
figure（1）;semilogy（SNRBERC（1:）‘-‘SNRBERC（2:）‘-x‘SNRBERC（3:）‘--‘）;
xlabel（‘Average?electrical?SNR（dB）‘）;ylabel（‘BER‘）;title（‘BER-SNR（Cn2）‘）;
legend（‘Cn2=0.51e-13‘‘Cn2=1.13e-13‘‘Cn2=2.10e-13‘‘Location‘‘NorthEast‘）;

Rytov?=?1.23?*?Cn（1）?*?k^（7/6）*?L.^（11/6）;
sigma?=?sqrt（Rytov）;
a=zeros（13）;b=zeros（13）;
for?i?=?1:3
????a（i）?=?（exp（0.49*Rytov（i）/（1+1.11*sigma（i）^（12/5））^（7/6））-1）^（-1）;
????b（i）?=?（exp（0.51*Rytov（i）/（1+0.69*sigma（i）^（12/5））^（5/6））-1）^（-1）;
????s=[‘[（1-‘num2str（a（i））‘）/2（2-‘num2str（a（i））‘）/2（1-‘num2str（b（i））‘）/2（2-‘num2str（b（i））‘）/21]‘];
????for?j?=?1:n
????????z?=（（2/（a（i）*b（i）））^2）*10^（SNR（j）/10）;
????????BERL（ij）?=?2^（a（i）+b（i）-3）/（pi^（3/2）*gamma（a（i））*gamma（b（i）））*double（feval（symengine‘meijerG‘‘2‘‘4‘s?‘[0?1/2]‘?num2str（z）））;?
????end
end
figure（2）;semilogy（SNRBERL（1:）‘-‘SNRBERL（2:）‘-x‘SNRBERL（3:）‘--‘）;
xlabel（‘Average?electrical?SNR（dB）‘）;ylabel（‘BER‘）;title（‘BER-SNR（L）‘）;
legend（‘L=1000‘‘L=3000‘‘L=5000‘‘Location‘‘NorthEast‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2730??2018-04-11?17:33??zernike.m

-----------?---------??----------?-----??----

?????????????????2730????????????????????1