大气湍流相位屏仿真代码.zip

大气湍流相位屏仿真代码.zip

%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