Clear Filters
Clear Filters

How code with for loop to get an array?

2 views (last 30 days)
Athira T Das
Athira T Das on 9 May 2024
Commented: Torsten on 9 May 2024
Ran in:
The code given below evaluates I which is 512*512 matrix. Using the for loop the I is evaluated for 757 values of z. Then I calculated SI and BER, but not an array. I need BER and SI for each I corresponding to each z value, or need an array of SI and BER of 1*757. Can anyone help me to with the code?
clc;close all;clear all;
lambda = 532;
k = 2*pi/lambda;
z=[1.03700000000000 1.39100000000000 1.59800000000000 1.84800000000000 1.88900000000000 2.01400000000000 2.34800000000000 2.45100000000000 2.57600000000000 2.72200000000000 2.78500000000000 2.88900000000000 2.97200000000000 3.16000000000000 3.20200000000000 3.36900000000000 3.41100000000000 3.51500000000000 3.55700000000000 3.57700000000000 3.66100000000000 3.66100000000000 3.68100000000000 3.70200000000000 3.74400000000000 3.82800000000000 3.93200000000000 4.05700000000000 4.18200000000000 4.55600000000000 4.70200000000000 4.91000000000000 5.01500000000000 5.16100000000000 5.24400000000000 5.26500000000000 5.28600000000000 5.28600000000000 5.28600000000000 5.24400000000000 5.26400000000000 5.22300000000000 5.39000000000000 5.57700000000000 5.72200000000000 5.93100000000000 6.16000000000000 6.41000000000000 6.63800000000000 6.88800000000000 7.03400000000000 7.24200000000000 7.34700000000000 7.45100000000000 7.53400000000000 7.63900000000000 7.68000000000000 7.80500000000000 7.82600000000000 7.95100000000000 8.09700000000000 8.20100000000000 8.30600000000000 8.47200000000000 8.66000000000000 8.82600000000000 8.97200000000000 9.07700000000000 9.18100000000000 9.24300000000000 9.28400000000000 9.32600000000000 9.34700000000000 9.34700000000000 9.36800000000000 9.36800000000000 9.40900000000000 9.47200000000000 9.51400000000000 9.57600000000000 9.65900000000000 9.74200000000000 9.84600000000000 9.99200000000000 10.1790000000000 10.5130000000000 10.5760000000000 10.9710000000000 11.1380000000000 11.3260000000000 11.4090000000000 11.4710000000000 11.5130000000000 11.5550000000000 11.6170000000000 11.7000000000000 11.8250000000000 11.8870000000000 12.2830000000000 12.4910000000000 12.6370000000000 12.8040000000000 12.9490000000000 13.1150000000000 13.2610000000000 13.5540000000000 13.6790000000000 13.8240000000000 13.9080000000000 14.3040000000000 14.3870000000000 14.4910000000000 14.5740000000000 14.6370000000000 14.7620000000000 14.8450000000000 14.9490000000000 15.0320000000000 15.1370000000000 15.2820000000000 15.3650000000000 15.5320000000000 15.6160000000000 15.7410000000000 15.8870000000000 16.0530000000000 16.1570000000000 16.3650000000000 16.5940000000000 16.7400000000000 16.8650000000000 16.9280000000000 17.0530000000000 17.0950000000000 17.1780000000000 17.2200000000000 17.3030000000000 17.3440000000000 17.5530000000000 17.6980000000000 17.8020000000000 17.9900000000000 18.1360000000000 18.3650000000000 18.5740000000000 18.8030000000000 18.9900000000000 19.1980000000000 19.5100000000000 19.6560000000000 19.8430000000000 19.9690000000000 20.0940000000000 20.1760000000000 20.3220000000000 20.4680000000000 20.5930000000000 20.7590000000000 20.9050000000000 21.0510000000000 21.1970000000000 21.4260000000000 21.5920000000000 21.7170000000000 21.9050000000000 22.0300000000000 22.2180000000000 22.3430000000000 22.5510000000000 22.8420000000000 23.0920000000000 23.1960000000000 23.4250000000000 23.5920000000000 23.7170000000000 23.8630000000000 24.0500000000000 24.1960000000000 24.3210000000000 24.5910000000000 24.7370000000000 24.8620000000000 25.0490000000000 25.2580000000000 25.5490000000000 25.7160000000000 25.8830000000000 26.0910000000000 26.2370000000000 26.3200000000000 26.6960000000000 27.4040000000000 27.6120000000000 27.7790000000000 27.9870000000000 28.1540000000000 28.3620000000000 28.5280000000000 28.6950000000000 28.8830000000000 29.2370000000000 29.4450000000000 29.5490000000000 29.7990000000000 29.9650000000000 30.1320000000000 30.3190000000000 30.5070000000000 30.6950000000000 30.8200000000000 30.9860000000000 31.1320000000000 31.3200000000000 31.3820000000000 31.4870000000000 31.5910000000000 31.7580000000000 31.8610000000000 31.9870000000000 32.0700000000000 32.3610000000000 32.4860000000000 32.6320000000000 32.8810000000000 33.0060000000000 33.0480000000000 33.0690000000000 33.0900000000000 33.0900000000000 33.0690000000000 33.0270000000000 32.9240000000000 32.9030000000000 32.8610000000000 32.8820000000000 32.9030000000000 32.9660000000000 33.0490000000000 33.0910000000000 33.1740000000000 33.2570000000000 33.3200000000000 33.3820000000000 33.4660000000000 33.4440000000000 33.4230000000000 33.4020000000000 33.4230000000000 33.3810000000000 33.4020000000000 33.4020000000000 33.4020000000000 33.4020000000000 33.4020000000000 33.4440000000000 33.4640000000000 33.5270000000000 33.5270000000000 33.6110000000000 33.6520000000000 33.6730000000000 33.6930000000000 33.6740000000000 33.6740000000000 33.6320000000000 33.5280000000000 33.4650000000000 33.3820000000000 33.4030000000000 33.3610000000000 33.3400000000000 33.3200000000000 33.2570000000000 33.3410000000000 33.3820000000000 33.4240000000000 33.4440000000000 33.4650000000000 33.4440000000000 33.5060000000000 33.5060000000000 33.5480000000000 33.5480000000000 33.5480000000000 33.5690000000000 33.6100000000000 33.7350000000000 33.8390000000000 33.9220000000000 34.0270000000000 34.1310000000000 34.3180000000000 34.5480000000000 34.7350000000000 34.8810000000000 35.0480000000000 35.2140000000000 35.4440000000000 35.5900000000000 35.7150000000000 36.0690000000000 36.1730000000000 36.3400000000000 36.5270000000000 36.5900000000000 36.7150000000000 36.9640000000000 37.0270000000000 37.1730000000000 37.3190000000000 37.4640000000000 37.6100000000000 37.7770000000000 37.9230000000000 38.0680000000000 38.1520000000000 38.2980000000000 38.4230000000000 38.5480000000000 38.5480000000000 38.5680000000000 38.6100000000000 38.6520000000000 38.7140000000000 38.7770000000000 38.9230000000000 39.0900000000000 39.2360000000000 39.3810000000000 39.7140000000000 39.8390000000000 40.0060000000000 40.1310000000000 40.2560000000000 40.3810000000000 40.4640000000000 40.5270000000000 40.6100000000000 40.6930000000000 40.7550000000000 40.8800000000000 41.0260000000000 41.1300000000000 41.2350000000000 41.3810000000000 41.5260000000000 41.6510000000000 41.7560000000000 41.9220000000000 42.0890000000000 42.2140000000000 42.4850000000000 42.6520000000000 42.8600000000000 42.9650000000000 43.1520000000000 43.2560000000000 43.4020000000000 43.5060000000000 43.7140000000000 43.7560000000000 43.8390000000000 43.9010000000000 43.9850000000000 44.1100000000000 44.2140000000000 44.3380000000000 44.5470000000000 44.6720000000000 44.8380000000000 44.9630000000000 45.1300000000000 45.3380000000000 45.4420000000000 45.6500000000000 45.7540000000000 46.0870000000000 46.1910000000000 46.3580000000000 46.4620000000000 46.5460000000000 46.6290000000000 46.6710000000000 46.7540000000000 46.8170000000000 46.8790000000000 47.0460000000000 47.1080000000000 47.2130000000000 47.3590000000000 47.6710000000000 47.9000000000000 48.0660000000000 48.2540000000000 48.4210000000000 48.5460000000000 48.7120000000000 48.7950000000000 48.8160000000000 48.9210000000000 48.9830000000000 49.0250000000000 49.0670000000000 49.0870000000000 49.2120000000000 49.2960000000000 49.4210000000000 49.4840000000000 49.7120000000000 49.9410000000000 50.0450000000000 50.6070000000000 50.7320000000000 50.8370000000000 50.9200000000000 51.0040000000000 51.1290000000000 51.1910000000000 51.2960000000000 51.4210000000000 51.4830000000000 51.5250000000000 51.6080000000000 51.7130000000000 51.7540000000000 51.8170000000000 51.9210000000000 52.0040000000000 52.0660000000000 52.1910000000000 52.2950000000000 52.4620000000000 52.5450000000000 52.6910000000000 52.8160000000000 53.0650000000000 53.2320000000000 53.2950000000000 53.3980000000000 53.4400000000000 53.5240000000000 53.6280000000000 53.6490000000000 53.7530000000000 53.8570000000000 53.9200000000000 54.0660000000000 54.1490000000000 54.2730000000000 54.4190000000000 54.5230000000000 54.6280000000000 54.7320000000000 54.8160000000000 54.8990000000000 55.0240000000000 55.0860000000000 55.1490000000000 55.2320000000000 55.2950000000000 55.4200000000000 55.5030000000000 55.7110000000000 55.8150000000000 56.0650000000000 56.1910000000000 56.4200000000000 56.5030000000000 56.6070000000000 56.6910000000000 56.7530000000000 56.8570000000000 56.9830000000000 57.3160000000000 57.4610000000000 57.6280000000000 57.7730000000000 57.9610000000000 58.1270000000000 58.5020000000000 58.6690000000000 58.8150000000000 58.9190000000000 59.0650000000000 59.1690000000000 59.2520000000000 59.4190000000000 59.4190000000000 59.5230000000000 59.5440000000000 59.6280000000000 59.7530000000000 59.8160000000000 59.9410000000000 60.1490000000000 60.2950000000000 60.4610000000000 61.1270000000000 61.2530000000000 61.5020000000000 61.9390000000000 62.0850000000000 62.2100000000000 62.3770000000000 62.7100000000000 62.8980000000000 63.0440000000000 63.2940000000000 63.4820000000000 63.7110000000000 63.8980000000000 64.1060000000000 64.3150000000000 64.5020000000000 64.6900000000000 64.8770000000000 65.0850000000000 65.2100000000000 65.3150000000000 65.4400000000000 65.6270000000000 65.7310000000000 65.8560000000000 66.0230000000000 66.1270000000000 66.2510000000000 66.3560000000000 66.6470000000000 66.7510000000000 66.8760000000000 66.9380000000000 67.1890000000000 67.3140000000000 67.3970000000000 67.5220000000000 67.5640000000000 67.7310000000000 67.8560000000000 67.9810000000000 68.0650000000000 68.1480000000000 68.2730000000000 68.3980000000000 68.4810000000000 68.6060000000000 68.6680000000000 68.8140000000000 68.8970000000000 69.0220000000000 69.2090000000000 69.2920000000000 69.6460000000000 69.6670000000000 69.8340000000000 70.0840000000000 70.1890000000000 70.3140000000000 70.4810000000000 70.6050000000000 70.7510000000000 70.8760000000000 71.0430000000000 71.1470000000000 71.2720000000000 71.4390000000000 71.7100000000000 71.8770000000000 71.9810000000000 72.1470000000000 72.2720000000000 72.3550000000000 72.5010000000000 72.5430000000000 72.6470000000000 72.7520000000000 72.8140000000000 72.9610000000000 73.0440000000000 73.1690000000000 73.3140000000000 73.4810000000000 74.0430000000000 74.1890000000000 74.3980000000000 74.5220000000000 74.6050000000000 74.7100000000000 74.7300000000000 74.7720000000000 74.8130000000000 74.8330000000000 74.8960000000000 74.9370000000000 75.0210000000000 75.0620000000000 75.2080000000000 75.3540000000000 75.7080000000000 75.8750000000000 76.0630000000000 76.2300000000000 76.4380000000000 76.5430000000000 76.7090000000000 76.8760000000000 77.0010000000000 77.1260000000000 77.2720000000000 77.3770000000000 77.4800000000000 77.6050000000000 77.8550000000000 77.9600000000000 78.0640000000000 78.1060000000000 78.2100000000000 78.3970000000000 78.5220000000000 78.6270000000000 78.8550000000000 79.0430000000000 79.2300000000000 79.6470000000000 79.8760000000000 80.5220000000000 80.6260000000000 80.7930000000000 80.8550000000000 80.8970000000000 80.9180000000000 80.9390000000000 80.9810000000000 81.0430000000000 81.1270000000000 81.2730000000000 81.4600000000000 81.6270000000000 81.8350000000000 82.1260000000000 82.6260000000000 82.8760000000000 83.1050000000000 83.2720000000000 83.4390000000000 83.5850000000000 83.6890000000000 83.7310000000000 83.7100000000000 83.6890000000000 83.6890000000000 83.7100000000000 83.8350000000000 83.9180000000000 84.0850000000000 84.4600000000000 84.6890000000000 84.9180000000000 85.3970000000000 85.5220000000000 85.6260000000000 85.8560000000000 85.8550000000000 85.8560000000000 85.8760000000000 85.9380000000000 85.9390000000000 86.0430000000000 86.1470000000000 86.2510000000000 86.4180000000000 86.6470000000000 86.8550000000000 87.0850000000000 87.3340000000000 87.5630000000000 87.7720000000000 88.0630000000000 88.1890000000000 88.4390000000000 88.5430000000000 88.6480000000000 88.6480000000000 88.7100000000000 88.7730000000000 88.8770000000000 88.8980000000000 89.0010000000000 89.1050000000000 89.2930000000000 89.4180000000000 89.6470000000000 89.7930000000000 90.0010000000000 90.2100000000000 90.4180000000000 90.6060000000000 90.7100000000000 90.8770000000000 91.0640000000000 91.1680000000000 91.2720000000000 91.3560000000000 91.4180000000000 91.5220000000000 91.6470000000000 91.7310000000000 91.8350000000000 91.9180000000000 92.1060000000000 92.2310000000000 92.3760000000000 92.6060000000000 92.7930000000000 93.0010000000000 93.4180000000000 93.6470000000000 93.8970000000000 94.4180000000000 94.5840000000000 94.7300000000000 95.0220000000000 95.1890000000000 95.3140000000000 95.5220000000000 95.8140000000000 95.9180000000000 96.0640000000000 96.2300000000000 96.3130000000000 96.5010000000000 96.6260000000000 96.7720000000000 96.8560000000000 97.0230000000000 97.1680000000000 97.2720000000000 97.3970000000000 97.5430000000000 97.6680000000000 97.7520000000000 97.8770000000000 97.9600000000000 98.0650000000000 98.3360000000000 98.4190000000000 98.5230000000000 98.6070000000000 98.7320000000000 98.8150000000000 98.9190000000000 99.1060000000000 99.2520000000000 99.4190000000000 99.5860000000000 99.7940000000000 99.9190000000000 100.044000000000];
z_diff = (diff(z)) ;
z_diff = [0 z_diff];
tilt_angle=10;
beta=(pi/180)*tilt_angle;
Lx = 0.1;
Ly = Lx;
Nx = 512;
Ny = Nx ;
dx = Lx/Nx;
dy = Ly/Ny;
x = (-Nx/2:Nx/2-1)*dx;
y = (-Ny/2:Ny/2-1)*dy;
[X,Y] = meshgrid(x,y);
r = sqrt(X.^2 + Y.^2);
kx = meshgrid(-Nx/2:Nx/2-1)/(Nx.*dx);
ky = meshgrid(-Ny/2:Ny/2-1)/(Ny.*dy);
kz = ones(Nx);
Kx=cos(beta).*kx + sin(beta).*kz;
Ky=ky;
Kz=-sin(beta).*kx + cos(beta).*kz;
K2 = Kx.^2 + Ky.^2 + Kz.^2;
delfx = 1/(Lx);
delfy = delfx;
R = randn(Nx,Ny) ;
R = R - mean(R(:,:));
R = R/std(R(:,:));
Phi = [2.50310494762771e-07 2.50946965543003e-07 4.74127533715154e-07 5.24210357233832e-06 2.42362427581638e-07 1.42653443399616e-07 5.16930235276850e-07 9.04590700562310e-07 2.29541191637123e-06 3.60105794008093e-07 9.17845719860523e-08 2.50390367004531e-07 1.09738615765455e-06 1.00311041468925e-06 6.50810613011221e-08 3.02640252523566e-07 5.31047019902531e-07 1.80313245178286e-06 3.04099822876624e-06 1.47008194185741e-07 4.08290582950848e-08 3.82118426499270e-06 5.29334397036188e-07 4.87810778598663e-06 1.18782521297988e-06 1.16666772914768e-06 1.69778706192608e-06 7.67426466736916e-06 5.27799451227444e-06 4.23216916272082e-07 2.66828750584071e-05 2.46205729358867e-06 2.80169965568206e-06 2.97378024567976e-07 2.21170607901861e-07 3.41285818450589e-07 3.85275104621375e-08 4.05141179263807e-07 1.17629321652161e-06 4.95233464403551e-06 1.80125151270180e-06 2.81483058944610e-07 1.24528365857859e-05 6.86781782456521e-06 2.95828059847895e-06 1.27035361174474e-06 6.34388599513081e-07 2.32012818083199e-06 2.81784487117724e-06 1.46096455213802e-06 1.60178082889894e-06 1.94791694377124e-07 6.35459501021705e-08 8.30940388112504e-08 2.53354933929590e-07 4.31247522543152e-06 9.63180625859285e-07 8.03726503455623e-07 2.98282116276067e-07 9.25402679823575e-08 1.61856374619784e-07 1.73242501417858e-07 2.31570160754109e-07 1.78600321810053e-07 1.25598618058698e-07 1.63418364388577e-06 1.02774042441542e-06 8.32018926678718e-08 2.53586554087675e-07 2.53508325030189e-07 4.07413642129177e-08 2.42096085408459e-07 5.91587468209778e-07 1.27761778023523e-07 1.65773731666773e-07 1.79565395840267e-07 3.78521468221561e-07 1.21471607158554e-07 3.25378687175317e-07 1.25569084575180e-07 5.49035746100893e-08 4.07733632799329e-08 4.79103771763424e-07 3.27135102365317e-07 2.10496092551324e-07 1.33641668807187e-06 1.68227872575447e-06 1.52265748820018e-06 2.53683706505161e-07 1.25630494723704e-07 1.12347756326375e-06 2.53692684133802e-07 2.10500428113819e-07 2.52934873659000e-07 6.77249701528946e-07 2.95386782896753e-07 3.16617852613830e-07 3.59879226900604e-07 2.81244834536819e-07 2.35438339516804e-07 3.40806542459708e-07 6.44683615098099e-08 5.05945401350734e-07 2.57387924440325e-07 7.82108421793457e-08 2.52513556939930e-07 7.41154472489689e-08 8.95664707278605e-08 3.44823854262242e-07 2.57498619329261e-07 2.07939465354519e-06 7.18570183290631e-07 3.25198024230522e-07 2.11933525065655e-07 8.30783371984955e-08 1.65710628194624e-07 5.48396993890243e-08 5.95342783012939e-08 2.28819020334242e-07 1.99067722101735e-07 2.40425372248607e-07 2.51916559208949e-07 2.98446361528618e-07 2.82958505602994e-07 2.10599141757522e-07 4.60997674322361e-08 2.59856693650578e-07 1.47642122461920e-07 3.93574284250877e-07 3.94729289526600e-08 6.37347851352551e-08 2.51536627057736e-07 2.28438571372065e-07 1.90144303446780e-07 4.73626330414131e-08 6.44142268235402e-08 1.72106272576708e-06 1.70721924857358e-06 2.94643811973548e-07 3.55943119776846e-08 9.36518317402136e-08 4.31201179335937e-08 1.56428758776995e-07 3.19975996745435e-07 3.34275554966361e-07 1.20064760586921e-06 4.30984033789707e-07 2.15616245567332e-07 3.03074660779611e-07 3.29934827134432e-07 -2.24864187842034e-08 2.33412514429474e-07 3.21092520408291e-07 3.38617740923693e-07 6.04579881670076e-07 1.91481555229955e-07 9.38472397526705e-08 6.13052684855612e-07 3.97260343618278e-06 2.45355935394897e-07 3.47585309986832e-07 3.38923893338132e-07 2.20926360286704e-07 5.25446559955282e-07 2.84158336921376e-07 4.04755442008443e-07 3.35541338245773e-07 3.76856346223440e-07 2.75568777774655e-07 2.09025928742013e-07 7.10170870233622e-08 2.95787195332815e-07 5.59078661720798e-07 2.51150081156080e-07 4.02707325242296e-08 8.24148001852386e-08 2.57363940284162e-07 3.55816103375846e-07 3.95777597622744e-07 3.27993323594138e-07 3.52690891300866e-07 3.52278318886910e-07 5.06712192595194e-08 1.23808524289852e-07 8.39347385990663e-07 2.02892942754780e-07 7.97641133582245e-08 3.64530615247902e-07 9.11185642215493e-07 2.44511041473120e-07 3.24963966176285e-07 1.11685979674259e-07 2.86261990645848e-07 3.75634147009314e-08 2.16925989302430e-07 3.84116828662608e-07 2.71791109589261e-07 1.06845659312697e-07 2.02225581949648e-07 3.02535099485185e-07 3.11957552477906e-08 9.15881915538412e-08 9.18891121954120e-07 3.70106500980225e-07 2.25779726927071e-07 2.86067288473873e-07 5.66382866246259e-07 4.26632108086086e-07 4.53045634671475e-07 1.17973378714372e-06 3.40969918417227e-07 2.48643473538260e-07 2.64505922988546e-07 2.38444406096982e-07 3.47101823144428e-07 3.37625374999722e-07 3.20790426924069e-07 5.02178252920690e-08 1.73576015527609e-07 8.57402376593969e-07 1.80131811857122e-07 2.76964595681315e-07 2.14838038341036e-07 2.21053263051776e-07 1.59406260405821e-07 4.59614264881235e-08 3.78789455759776e-07 2.58966590318316e-07 7.75696011578299e-08 3.39990416320565e-07 3.24572509489827e-07 3.52020828992088e-08 6.49912862168267e-08 2.22828739301052e-07 2.13473046057573e-07 1.70720547021041e-07 2.19582224668132e-07 2.05914598785446e-07 1.14899310738254e-06 1.03415047080757e-07 1.89177690522258e-06 2.24175047702725e-07 3.76150153666139e-07 9.18124634120261e-08 2.27678254839140e-07 8.61346207114181e-08 3.38947999645138e-06 1.71880093129416e-07 2.22461873169689e-07 3.72427593494344e-08 6.76102058093114e-08 1.13163515787091e-06 7.04052219453715e-07 3.28310827940334e-07 2.12060821359437e-07 2.68576701610855e-07 2.04744061114503e-07 1.14622616867334e-07 4.35675529886048e-08 5.29665131717990e-08 3.50574437154752e-08 5.32920831968263e-08 1.04801149974382e-07 2.08316009852988e-07 2.77349913208742e-07 1.48756450473382e-07 3.78517042358300e-07 1.19222278136494e-07 4.57902170141890e-08 2.03249187200450e-07 9.32437812713631e-08 6.61503340988668e-07 1.01572927594225e-06 3.26294032932256e-07 2.26313604682662e-07 7.01012285236915e-08 7.61666326540900e-08 4.16954886967803e-08 4.41719683698357e-08 4.02549730011751e-08 4.41603434482274e-08 3.50371730851861e-08 3.00718061774876e-07 3.50372423542925e-08 5.32741630316562e-08 2.09046711312522e-07 3.75611845531339e-06 2.26535469768828e-07 4.62710159517236e-07 5.55357457448765e-07 4.16702305525217e-07 2.05922476618442e-07 2.06490274401431e-07 2.33754328058605e-07 1.66615319203810e-07 2.13390792614979e-07 2.54792444311713e-07 4.20241426628409e-07 6.95047952441784e-07 6.02254152247896e-07 1.56218004764843e-07 1.99836446669293e-07 1.62974708502310e-07 2.93420921132526e-07 1.67909832369442e-07 1.53692619622575e-07 8.89709805825161e-08 1.64148574413349e-07 1.82871518257336e-07 1.99566615308843e-07 1.59784291140960e-07 1.87893662090374e-07 2.77364481028229e-07 1.10443572547260e-07 2.18556580353372e-07 1.83932131011210e-07 7.94652011731045e-08 2.43919458395486e-07 1.59781205439559e-07 1.94642054497664e-07 1.29989372686070e-07 1.90910393259075e-07 2.84452061601644e-07 3.30006526167654e-07 1.93485493531741e-07 2.40209860863182e-07 4.15333943239071e-07 2.91343651778087e-07 2.58377576473351e-07 1.81274762632482e-07 1.88893039486970e-07 2.19428417763272e-07 2.20994023873646e-07 2.27950728874121e-07 1.64400439994395e-07 1.00341668297253e-07 1.67577483763212e-07 5.61980802260245e-08 3.04539633123467e-08 3.19018898206730e-08 3.89504322796163e-08 3.19941498021468e-07 5.41874921040340e-06 1.01776790218898e-06 1.29079721421378e-06 5.77887058141474e-07 2.79939362285509e-07 5.55989804582382e-08 2.11144624551603e-07 5.73140726411827e-08 5.35992572193921e-08 2.80066251579433e-07 1.76964657359828e-07 3.25199602895421e-07 2.18271466195228e-07 1.80628024038326e-07 4.69550563668357e-08 1.48055689005769e-07 1.82317859358405e-07 1.33063533148151e-07 2.62672490298873e-08 2.52416541944235e-08 2.08244390466210e-07 3.06261128674104e-07 2.18183988930487e-07 1.32114003617060e-07 3.13731744862144e-07 2.78980205382911e-07 2.42110700965470e-08 1.28818767334643e-07 3.48416502648141e-07 1.43250842379859e-07 2.58790299254369e-08 2.36803673262810e-08 2.54679999490770e-08 1.53423624973069e-07 2.11044574009824e-07 2.45219073821534e-08 5.40329423903049e-08 4.07828910547019e-08 2.86866659397290e-08 5.74108869690081e-08 3.17780519921456e-08 2.89072270765441e-07 6.46668413187508e-08 3.89801393203052e-07 1.62322978572959e-07 2.69026419637955e-07 4.76120100475985e-07 5.32185045407108e-07 1.57575408820348e-07 2.08858600548085e-07 1.87565977660685e-07 2.01108310997773e-07 3.95620143987532e-08 5.24866345380986e-08 8.64930520826212e-08 9.74072997804584e-07 2.54439638933078e-08 1.00154777975012e-07 9.42429455464106e-07 2.23088555779751e-07 2.02097796370716e-07 1.99163165408265e-07 9.26536634692885e-08 1.31298738607173e-07 1.97726692783641e-07 2.28150683478771e-07 2.17170907865184e-07 2.69677497220439e-08 1.54765400027711e-07 6.26255760215395e-08 1.46771827369792e-07 1.81515822529873e-07 1.72918459779846e-07 8.15990146285526e-08 9.85514009220808e-08 7.03396855430348e-07 5.93657045438661e-07 1.51834395684521e-07 9.76927764035926e-08 1.18527867299870e-07 2.71435324116586e-08 5.68571003566081e-08 1.65983256102705e-07 2.29129793972797e-07 5.88646431511064e-08 3.08684402204026e-05 2.38265317337181e-06 4.30752744717656e-07 1.04689461805175e-07 2.32270648752108e-07 6.07211466407868e-07 1.95341105985221e-08 7.38231690963101e-08 1.95129857199575e-07 4.01671898332107e-08 4.75072551754564e-08 2.49046064060719e-06 1.78142428447807e-07 2.60935684174315e-07 4.29594269700864e-08 1.14070714922750e-07 5.26390027271695e-08 6.74132420131039e-08 2.11045265009953e-07 7.27021148761358e-07 2.83886081059763e-07 2.76343292567963e-08 1.61294742863238e-07 1.83304573271511e-07 2.45251911632262e-07 9.18095053702886e-07 9.40457085677553e-07 2.90451584499496e-07 1.04003175628026e-07 1.08284695915254e-07 5.44402482199441e-08 2.02578013696944e-08 5.88615942025308e-08 1.04406254724335e-07 1.20598226687063e-07 5.71264670628078e-07 3.15144191075623e-07 2.23052176438872e-08 2.98936409995515e-08 2.92428609360162e-07 2.47505388371339e-07 2.72866577624107e-07 2.69881858030521e-07 1.85519336468060e-07 2.56402576998181e-07 8.88653855726949e-08 1.23281438816645e-07 3.74059740200056e-08 3.63466577344003e-08 1.20887062612203e-06 5.31728280982040e-07 2.22245937870703e-08 1.24471844746278e-07 5.11735373064063e-07 3.13008213548308e-07 5.70247343660301e-08 7.23194564373494e-07 1.79616235944090e-07 6.78737172874271e-08 1.87999018332358e-07 1.56051038058590e-07 2.31029093824182e-08 3.65957708093627e-08 5.13581582166058e-07 2.55008648105981e-08 1.44166238137747e-07 3.81609462756446e-07 2.16242665285891e-07 2.34087134829968e-07 2.55408939395098e-07 5.31844204547311e-08 1.37230800917903e-07 2.24923680386313e-07 7.91242025449375e-08 1.56574476334133e-07 1.34926252415438e-07 1.70013819330808e-08 1.48130847036606e-07 1.10651211565701e-07 1.58123074045499e-07 1.47865663105776e-07 1.14316450201195e-07 2.02955526741989e-08 1.49518466142608e-07 2.99469635848321e-08 4.06175629714119e-07 2.32871553385763e-07 1.70948094081342e-07 1.86984334344187e-07 1.54774299548728e-07 1.27005955203776e-07 7.29605614875902e-08 4.26639317987837e-08 5.91930531596668e-08 1.50666658228087e-08 4.86748534643297e-08 1.64538615062307e-08 1.32252962938291e-07 1.06729078335828e-07 1.39785078669013e-07 1.39779800669610e-07 9.26422941859900e-08 1.13233915195479e-07 8.32946711657419e-08 1.67451405167575e-07 3.65621117891888e-08 2.21661117907828e-07 3.65652178907218e-08 1.41441861581090e-06 3.30075312403152e-07 2.05346680636497e-08 2.16432978289422e-07 1.92694668921994e-07 1.99379501668891e-07 2.60904461669433e-07 2.16709099204294e-07 1.33136064843243e-07 1.52413473502972e-07 7.53889226907776e-07 3.22036453180506e-07 1.68931528429636e-07 1.04094510285934e-07 1.60098934336117e-07 1.02657996503878e-07 6.83661821001145e-08 1.61931255758191e-07 2.01320024686194e-07 1.07534006069169e-07 1.07910634894973e-07 9.46522813533042e-08 1.66743437525985e-07 1.78365196465516e-07 1.90349497426267e-07 1.99051143944893e-07 1.10969246755729e-07 1.11264782143939e-07 1.02006294022338e-07 1.24767050443006e-07 1.04932154146443e-07 1.03659587531759e-07 4.90478970104965e-08 9.01107885624077e-08 9.30802687007756e-08 2.06842405795697e-07 1.53448550376276e-07 9.79061221486547e-08 8.26140384317179e-08 1.50146254483116e-07 1.40136146748372e-07 1.56837338087257e-07 1.76198808159403e-07 1.99122966368083e-07 1.23858764266714e-07 1.06727895602262e-07 8.15204687804434e-08 1.88227798887152e-08 2.31668781894091e-08 7.67325491571661e-07 1.97158231154847e-07 1.04232262587474e-07 8.81089549804027e-09 1.09695696650620e-08 3.67828021042941e-07 2.42873900856430e-07 9.30400519863338e-08 1.22815336391232e-08 1.04498084973760e-07 4.13726864269242e-07 1.28584375991528e-08 1.20376581165868e-07 1.07046398186240e-07 5.96207309583810e-08 3.66763563633301e-08 1.21422552113772e-07 1.39822102904246e-07 1.58475758239753e-07 8.54716186106630e-08 9.47107273229749e-08 1.57282457487923e-07 2.22341833330773e-07 1.47297286955243e-07 2.19435143761748e-07 1.39688617089499e-07 7.24259126481027e-08 7.09641458863346e-08 1.19226547380480e-07 7.87372467318109e-08 9.72314598010880e-08 7.74647075656411e-09 9.47412998435404e-08 1.16751285273194e-07 1.13487779361823e-07 1.13827758524977e-07 2.72522446715463e-07 8.16170870703966e-08 8.72321297381038e-08 9.95784392252335e-08 2.65989252708394e-09 1.06630932742284e-07 9.36921419305549e-08 8.41656329617571e-08 1.30269792124746e-08 4.86904164038956e-06 1.01469686899418e-07 1.47172552045969e-07 1.11949067863535e-07 9.34802107168088e-08 6.41098286680338e-08 6.68095594366830e-09 4.65659962399765e-08 1.74321722504626e-07 1.53177157686198e-07 1.72679939394816e-07 9.26791009109865e-08 7.38280843033631e-09 1.21431764801520e-08 1.83474589517062e-08 5.49758746890643e-08 5.49686356365556e-08 5.41396143667257e-07 1.02078716839878e-08 1.46280082306786e-07 9.20039209747377e-08 1.07464173701763e-08 4.37249887631680e-08 1.20696543908951e-07 7.37023067237918e-08 1.54370995646110e-07 8.75880053587604e-08 1.38972473988507e-07 1.50011697718075e-07 1.13900910702895e-07 7.30960487245964e-08 1.06352204500379e-07 7.15892552073097e-09 1.26382007766409e-08 2.34492999966371e-07 9.25652500580458e-08 8.91104525001921e-08 1.04346566674604e-08 7.96833204055166e-08 3.41626685510968e-07 1.46859973263572e-07 6.99097620123092e-08 5.60565018753087e-08 9.18716519344181e-08 5.82078430390537e-09 8.51757309259611e-08 1.41028716092150e-07 7.10210717526897e-09 3.02205722277794e-08 2.76205242367368e-07 9.08403872147274e-08 8.55641175181640e-08 8.07674139440669e-08 5.64950026938399e-09 1.47871928050032e-07 1.37623718221666e-07 9.60980796219824e-08 9.05395061229061e-09 3.58713661087815e-07 3.78246981154481e-07 1.97627162656081e-08 4.78463482905684e-08 5.47976120635378e-07 5.98012738422856e-07 1.38288635243591e-07 3.97209746106909e-08 1.04405402152745e-07 1.64230983323839e-08 2.14721179480984e-08 4.85083210883807e-08 1.22646011220379e-06 8.39456826640735e-09 8.27662116573820e-08 3.58891263221149e-07 9.16048888640066e-08 1.21278839928491e-07 2.31054125718236e-07 1.00889365297461e-07 1.19358834432565e-07 2.43427211256917e-07 9.11020950715102e-08 1.38162844699562e-07 8.00915682468241e-08 4.98258405386454e-09 7.30080219716715e-08 1.52701681888474e-07 2.93328698184570e-07 1.35353313498113e-07 4.59682925254202e-08 9.08782411223864e-08 2.20927979357592e-08 9.05503229704100e-08 6.68343694641667e-08 7.33025736057753e-08 7.49958147508063e-08 6.22824091009223e-09 8.02682015058001e-06 3.67877032868707e-07 7.19496678750717e-08 3.17758490462365e-09 5.75563157661276e-09 3.24272951280416e-09 1.00381201383279e-07 1.88260239049035e-07 7.07194739983604e-08 4.68947848174950e-09 4.16252458503278e-09 7.94576959843098e-08 2.74009776209153e-07 1.42430916804990e-07 1.26235913248699e-07 1.39626029880374e-07 1.16782761195100e-07 1.00256563902946e-07 -2.95101260456139e-09 6.83432403678117e-08 1.27348806727558e-07 5.96766339731236e-08 1.03604721020808e-07 1.21085249476580e-07 6.80561072346092e-08 2.39948064835149e-09 1.00688731186191e-08 2.43035685127838e-08 4.13254891501583e-08 9.63544005596039e-09 2.53883818143423e-07 7.44890597640385e-07];
w0=0.002;
F=2*(pi^2)*(k^2).*z_diff.*Phi;
x0 = linspace(-0.02,0.02,Nx);y0=x0;[X0,Y0] = meshgrid(x0,y0);
Uin=exp(-(X0.^2+Y0.^2)/w0^2);
x = x + z_diff(1)*tan(beta);
screen1 = exp(1i*((Kx.*x) + (Ky.*y))).*delfx.*delfy.*sqrt(F(1)).*R ;
Uout = Uin.*exp(1i.*screen1);
c=[0.740664357812470 0.715162315217582 0.689780546528124 0.664556687838984 0.639501229779350 0.614618312979130 0.589909677254809 0.565375929219662 0.541017064806704 0.516832720551119 0.492822307973243 0.468985091286926 0.445320235117717 0.421826835235231 0.398503939126698 0.375350560217132 0.352365687964806 0.329548295192595 0.306897343515746 0.284411787427262 0.262090577416753 0.239932662380474 0.217936991503006 0.196102515739251 0.174428188990141 0.152912969040777 0.131555818312352 0.110355704466664 0.0893116008928944 0.0684224870996102 0.0476873490298864 0.0271051793136574 0.00667497746850089 0.0136042499421828 0.0337334891863653 0.0537137192700918 0.0735459118688980 0.0932310312800026 0.112770034392025 0.132163870669524 0.151413482150099 0.170519803452171 0.189483761791847 0.208306277007520 0.226988261591059 0.245530620724636 0.263934252322323 0.282200047075789 0.300328888503445 0.318321653002534 0.336179209903699 0.353902421527639 0.371492143243500 0.388949223528713 0.406274504030009 0.423468819625392 0.440532998486862 0.457467862143721 0.474274225546305 0.490952897130006 0.507504678879476 0.523930366392895 0.540230748946220 0.556406609557339 0.572458725050048 0.588387866117797 0.604194797387153 0.619880277480917 0.635445059080876 0.650889888990134 0.666215508194993 0.681422651926368 0.696512049720699 0.711484425480339 0.726340497533411 0.741080978693105 0.755706576316408 0.770217992362264 0.784615923449129 0.798901060911951 0.813074090858528 0.827135694225279 0.841086546832386 0.854927319438336 0.868658677793846 0.882281282695167 0.895795790036781 0.909202850863477 0.922503111421817 0.935697213210995 0.948785793033077 0.961769483042646 0.974648910795838 0.987424699298783 1.00009746705545 1.01266782811491 1.02513639211797 1.03750376434332 1.04977054575300 1.06193733303735 1.07400471865938 1.08597329089862 1.09784363389434 1.10961632768826 1.12129194826677 1.13287106760249 1.14435425369543 1.15574207061353 1.16703507853271 1.17823383377643 1.18933888885472 1.20035079250267 1.21127008971854 1.22209732180122 1.23283302638735 1.24347773748787 1.25403198552414 1.26449629736359 1.27487119635488 1.28515720236266 1.29535483180181 1.30546459767133 1.31548700958769 1.32542257381781 1.33527179331163 1.34503516773420 1.35471319349739 1.36430636379121 1.37381516861470 1.38324009480640 1.39258162607448 1.40184024302645 1.41101642319850 1.42011064108442 1.42912336816425 1.43805507293243 1.44690622092571 1.45567727475060 1.46436869411057 1.47298093583279 1.48151445389465 1.48996969944983 1.49834712085410 1.50664716369079 1.51487027079592 1.52301688228302 1.53108743556759 1.53908236539135 1.54700210384610 1.55484708039729 1.56261772190731 1.57031445265850 1.57793769437581 1.58548786624928 1.59296538495611 1.60037066468257 1.60770411714556 1.61496615161393 1.62215717492952 1.62927759152798 1.63632780345927 1.64330821040796 1.65021920971321 1.65706119638862 1.66383456314170 1.67053970039319 1.67717699629614 1.68374683675467 1.69024960544263 1.69668568382193 1.70305545116068 1.70935928455115 1.71559755892740 1.72177064708287 1.72787891968756 1.73392274530517 1.73990249040993 1.74581851940328 1.75167119463032 1.75746087639608 1.76318792298160 1.76885269065981 1.77445553371120 1.77999680443936 1.78547685318630 1.79089602834756 1.79625467638721 1.80155314185264 1.80679176738913 1.81197089375438 1.81709085983270 1.82215200264920 1.82715465738367 1.83209915738443 1.83698583418192 1.84181501750219 1.84658703528021 1.85130221367304 1.85596087707287 1.86056334811985 1.86510994771484 1.86960099503199 1.87403680753115 1.87841770097023 1.88274398941727 1.88701598526258 1.89123399923051 1.89539834039132 1.89950931617273 1.90356723237146 1.90757239316461 1.91152510112087 1.91542565721170 1.91927436082231 1.92307150976254 1.92681740027765 1.93051232705897 1.93415658325440 1.93775046047891 1.94129424882478 1.94478823687186 1.94823271169764 1.95162795888723 1.95497426254328 1.95827190529575 1.96152116831153 1.96472233130413 1.96787567254306 1.97098146886325 1.97403999567435 1.97705152696986 1.98001633533631 1.98293469196218 1.98580686664685 1.98863312780943 1.99141374249746 1.99414897639556 1.99683909383401 1.99948435779718 2.00208502993198 2.00464137055607 2.00715363866619 2.00962209194621 2.01204698677524 2.01442857823560 2.01676712012071 2.01906286494297 2.02131606394144 2.02352696708957 2.02569582310279 2.02782287944605 2.02990838234126 2.03195257677470 2.03395570650435 2.03591801406711 2.03783974078601 2.03972112677735 2.04156241095768 2.04336383105087 2.04512562359498 2.04684802394913 2.04853126630031 2.05017558367013 2.05178120792145 2.05334836976503 2.05487729876606 2.05636822335072 2.05782137081251 2.05923696731876 2.06061523791685 2.06195640654055 2.06326069601621 2.06452832806891 2.06575952332861 2.06695450133616 2.06811348054936 2.06923667834885 2.07032431104407 2.07137659387907 2.07239374103836 2.07337596565262 2.07432347980443 2.07523649453394 2.07611521984447 2.07695986470806 2.07777063707106 2.07854774385950 2.07929139098462 2.08000178334822 2.08067912484798 2.08132361838281 2.08193546585808 2.08251486819086 2.08306202531507 2.08357713618662 2.08406039878854 2.08451201013601 2.08493216628138 2.08532106231914 2.08567889239087 2.08600584969017 2.08630212646748 2.08656791403493 2.08680340277114 2.08700878212594 2.08718424062516 2.08732996587524 2.08744614456794 2.08753296248490 2.08759060450230 2.08761925459532 2.08761909584274 2.08759031043135 2.08753307966046 2.08744758394630 2.08733400282638 2.08719251496389 2.08702329815202 2.08682652931821 2.08660238452849 2.08635103899165 2.08607266706347 2.08576744225093 2.08543553721630 2.08507712378131 2.08469237293121 2.08428145481885 2.08384453876873 2.08338179328099 2.08289338603538 2.08237948389527 2.08184025291152 2.08127585832643 2.08068646457760 2.08007223530177 2.07943333333869 2.07876992073487 2.07808215874743 2.07737020784777 2.07663422772538 2.07587437729149 2.07509081468279 2.07428369726507 2.07345318163688 2.07259942363312 2.07172257832864 2.07082280004183 2.06990024233815 2.06895505803365 2.06798739919850 2.06699741716044 2.06598526250827 2.06495108509530 2.06389503404273 2.06281725774310 2.06171790386362 2.06059711934957 2.05945505042762 2.05829184260918 2.05710764069365 2.05590258877174 2.05467683022873 2.05343050774771 2.05216376331278 2.05087673821232 2.04956957304210 2.04824240770851 2.04689538143169 2.04552863274867 2.04414229951646 2.04273651891521 2.04131142745123 2.03986716096009 2.03840385460964 2.03692164290309 2.03542065968197 2.03390103812916 2.03236291077188 2.03080640948460 2.02923166549207 2.02763880937220 2.02602797105897 2.02439927984538 2.02275286438630 2.02108885270135 2.01940737217775 2.01770854957318 2.01599251101859 2.01425938202100 2.01250928746633 2.01074235162217 2.00895869814052 2.00715845006058 2.00534172981146 2.00350865921493 2.00165935948813 1.99979395124624 1.99791255450517 1.99601528868428 1.99410227260896 1.99217362451332 1.99022946204283 1.98826990225688 1.98629506163147 1.98430505606170 1.98230000086442 1.98028001078079 1.97824519997880 1.97619568205582 1.97413157004115 1.97205297639850 1.96996001302854 1.96785279127132 1.96573142190885 1.96359601516744 1.96144668072028 1.95928352768979 1.95710666465011 1.95491619962948 1.95271224011267 1.95049489304337 1.94826426482658 1.94602046133099 1.94376358789132 1.94149374931069 1.93921104986296 1.93691559329504 1.93460748282926 1.93228682116561 1.92995371048410 1.92760825244698 1.92525054820109 1.92288069838009 1.92049880310668 1.91810496199492 1.91569927415240 1.91328183818251 1.91085275218661 1.90841211376627 1.90596002002548 1.90349656757275 1.90102185252339 1.89853597050162 1.89603901664270 1.89353108559514 1.89101227152279 1.88848266810698 1.88594236854864 1.88339146557039 1.88083005141866 1.87825821786576 1.87567605621199 1.87308365728765 1.87048111145519 1.86786850861116 1.86524593818833 1.86261348915771 1.85997125003053 1.85731930886031 1.85465775324484 1.85198667032819 1.84930614680269 1.84661626891091 1.84391712244766 1.84120879276190 1.83849136475875 1.83576492290140 1.83302955121309 1.83028533327899 1.82753235224814 1.82477069083540 1.82200043132329 1.81922165556395 1.81643444498099 1.81363888057138 1.81083504290734 1.80802301213817 1.80520286799215 1.80237468977836 1.79953855638852 1.79669454629885 1.79384273757188 1.79098320785825 1.78811603439854 1.78524129402510 1.78235906316378 1.77946941783578 1.77657243365938 1.77366818585177 1.77075674923075 1.76783819821655 1.76491260683352 1.76198004871193 1.75904059708966 1.75609432481398 1.75314130434320 1.75018160774846 1.74721530671537 1.74424247254577 1.74126317615938 1.73827748809551 1.73528547851471 1.73228721720050 1.72928277356095 1.72627221663041 1.72325561507115 1.72023303717496 1.71720455086485 1.71417022369661 1.71113012286052 1.70808431518289 1.70503286712770 1.70197584479823 1.69891331393861 1.69584533993546 1.69277198781941 1.68969332226677 1.68660940760100 1.68352030779434 1.68042608646938 1.67732680690054 1.67422253201569 1.67111332439765 1.66799924628574 1.66488035957729 1.66175672582918 1.65862840625934 1.65549546174827 1.65235795284053 1.64921593974624 1.64606948234258 1.64291864017527 1.63976347246002 1.63660403808406 1.63344039560754 1.63027260326505 1.62710071896701 1.62392480030119 1.62074490453409 1.61756108861240 1.61437340916446 1.61118192250163 1.60798668461974 1.60478775120050 1.60158517761291 1.59837901891464 1.59516932985347 1.59195616486864 1.58873957809224 1.58551962335059 1.58229635416566 1.57906982375634 1.57584008503991 1.57260719063331 1.56937119285456 1.56613214372404 1.56289009496590 1.55964509800932 1.55639720398990 1.55314646375097 1.54989292784488 1.54663664653434 1.54337766979374 1.54011604731041 1.53685182848594 1.53358506243751 1.53031579799910 1.52704408372284 1.52376996788025 1.52049349846351 1.51721472318676 1.51393368948733 1.51065044452699 1.50736503519323 1.50407750810048 1.50078790959138 1.49749628573799 1.49420268234301 1.49090714494105 1.48760971879982 1.48431044892133 1.48100938004314 1.47770655663956 1.47440202292281 1.47109582284426 1.46778800009561 1.46447859811008 1.46116766006355 1.45785522887583 1.45454134721172 1.45122605748225 1.44790940184585 1.44459142220945 1.44127216022968 1.43795165731401 1.43462995462187 1.43130709306582 1.42798311331269 1.42465805578467 1.42133196066046 1.41800486787641 1.41467681712758 1.41134784786894 1.40801799931637 1.40468731044786 1.40135582000454 1.39802356649183 1.39469058818047 1.39135692310765 1.38802260907807 1.38468768366504 1.38135218421152 1.37801614783119 1.37467961140955 1.37134261160493 1.36800518484960 1.36466736735075 1.36132919509161 1.35799070383243 1.35465192911156 1.35131290624645 1.34797367033472 1.34463425625514 1.34129469866868 1.33795503201952 1.33461529053606 1.33127550823193 1.32793571890697 1.32459595614830 1.32125625333121 1.31791664362027 1.31457715997020 1.31123783512696 1.30789870162865 1.30455979180654 1.30122113778599 1.29788277148748 1.29454472462753 1.29120702871964 1.28786971507532 1.28453281480497 1.28119635881886 1.27786037782809 1.27452490234548 1.27118996268656 1.26785558897046 1.26452181112088 1.26118865886698 1.25785616174431 1.25452434909574 1.25119325007237 1.24786289363441 1.24453330855215 1.24120452340678 1.23787656659136 1.23454946631169 1.23122325058718 1.22789794725176 1.22457358395477 1.22125018816181 1.21792778715567 1.21460640803712 1.21128607772588 1.20796682296139 1.20464867030376 1.20133164613454 1.19801577665766 1.19470108790021 1.19138760571334 1.18807535577308 1.18476436358118 1.18145465446596 1.17814625358312 1.17483918591663 1.17153347627945 1.16822914931448 1.16492622949528 1.16162474112694 1.15832470834688 1.15502615512566 1.15172910526776 1.14843358241246 1.14513961003454 1.14184721144515 1.13855640979258 1.13526722806305 1.13197968908148 1.12869381551230 1.12540962986023 1.12212715447104 1.11884641153231 1.11556742307426 1.11229021097045 1.10901479693858 1.10574120254124 1.10246944918668 1.09919955812955 1.09593155047166 1.09266544716271 1.08940126900108 1.08613903663452 1.08287877056091 1.07962049112901 1.07636421853917 1.07310997284407 1.06985777394945 1.06660764161482 1.06335959545421 1.06011365493683 1.05686983938786 1.05362816798911 1.05038865977972 1.04715133365691 1.04391620837666 1.04068330255439 1.03745263466568 1.03422422304698 1.03099808589625 1.02777424127369 1.02455270710241 1.02133350116911 1.01811664112476 1.01490214448529 1.01169002863224 1.00848031081345 1.00527300814371 1.00206813760545 0.998865716049361 0.995665760195108 0.992468286631939 0.989273311819364 0.986080852087796 0.982890923639205 0.979703542547760 0.976518724760467 0.973336486097822 0.970156842254431 0.966979808799657 0.963805401178249 0.960633634710965 0.957464524595207 0.954298085905637 0.951134333594803 0.947973282493756 0.944814947312659 0.941659342641416 0.938506482950264 0.935356382590394 0.932209055794548 0.929064516677626 0.925922779237284 0.922783857354535 0.919647764794513 0.916514553007848 0.913626501586297 0.957181910437377 1.18507611517074];
Ri=z_diff./cos(beta);
screens = zeros(Nx,Nx,length(F)-1);
for i = 2:length(F)
x = x + z_diff(i)*tan(beta);
screens(:,:,i)= screens(:, :, i-1) + exp(1i*((Kx.*x)+(Ky.*y))).*delfx.*delfy.*sqrt(F(i)).*R;
Hi = exp(1i*Kz*Ri(i) - c(i).*Ri(i));
Ui = fftshift(fft2(Uout));
Ui = Hi.*Ui;
Ui = ifft2(ifftshift(Ui));
Uio = Ui.*exp(1i.*screens(:,:,i));
I = Uio.*conj(Uio);
Uout = Uio;
end
% to calculate SigmaI
A1=0;A2=0;
for i1=1:Nx
for i2=1:Nx
I1 = I(i1,i2)^2 ;
A1 = A1 + I1 ;
end
end
A1
A1 = 3.1916e-257
for j1=1:Nx
for j2=1:Nx
I2 = I(j1,j2) ;
A2 = A2 + I2;
end
end
SI = abs(((A1/Nx)/((A2/Nx)^2)) - 1)
SI = 0.7497
%to calculate BER
syms u
PI = (1/(u*sqrt(2*pi*SI)))*exp(-((log(u) + SI )^2)/(2*SI)) ;
SNR0 = 1;
SNR = SNR0/sqrt(1+SI*SNR0) ;
fun = PI*erfc(real((SNR.*u)/(2*sqrt(2))));
ff = matlabFunction(fun,'Vars',u);
BER = (0.5*integral(ff,0,Inf))
BER = 0.4032

Sign in to answer this question.

Accepted Answer

Voss
Voss on 9 May 2024
Ran in:
You can make I a 3D array of size Nx-by-Nx-by-numel(F).
clc;close all;clear all;
lambda = 532;
k = 2*pi/lambda;
z=[1.03700000000000 1.39100000000000 1.59800000000000 1.84800000000000 1.88900000000000 2.01400000000000 2.34800000000000 2.45100000000000 2.57600000000000 2.72200000000000 2.78500000000000 2.88900000000000 2.97200000000000 3.16000000000000 3.20200000000000 3.36900000000000 3.41100000000000 3.51500000000000 3.55700000000000 3.57700000000000 3.66100000000000 3.66100000000000 3.68100000000000 3.70200000000000 3.74400000000000 3.82800000000000 3.93200000000000 4.05700000000000 4.18200000000000 4.55600000000000 4.70200000000000 4.91000000000000 5.01500000000000 5.16100000000000 5.24400000000000 5.26500000000000 5.28600000000000 5.28600000000000 5.28600000000000 5.24400000000000 5.26400000000000 5.22300000000000 5.39000000000000 5.57700000000000 5.72200000000000 5.93100000000000 6.16000000000000 6.41000000000000 6.63800000000000 6.88800000000000 7.03400000000000 7.24200000000000 7.34700000000000 7.45100000000000 7.53400000000000 7.63900000000000 7.68000000000000 7.80500000000000 7.82600000000000 7.95100000000000 8.09700000000000 8.20100000000000 8.30600000000000 8.47200000000000 8.66000000000000 8.82600000000000 8.97200000000000 9.07700000000000 9.18100000000000 9.24300000000000 9.28400000000000 9.32600000000000 9.34700000000000 9.34700000000000 9.36800000000000 9.36800000000000 9.40900000000000 9.47200000000000 9.51400000000000 9.57600000000000 9.65900000000000 9.74200000000000 9.84600000000000 9.99200000000000 10.1790000000000 10.5130000000000 10.5760000000000 10.9710000000000 11.1380000000000 11.3260000000000 11.4090000000000 11.4710000000000 11.5130000000000 11.5550000000000 11.6170000000000 11.7000000000000 11.8250000000000 11.8870000000000 12.2830000000000 12.4910000000000 12.6370000000000 12.8040000000000 12.9490000000000 13.1150000000000 13.2610000000000 13.5540000000000 13.6790000000000 13.8240000000000 13.9080000000000 14.3040000000000 14.3870000000000 14.4910000000000 14.5740000000000 14.6370000000000 14.7620000000000 14.8450000000000 14.9490000000000 15.0320000000000 15.1370000000000 15.2820000000000 15.3650000000000 15.5320000000000 15.6160000000000 15.7410000000000 15.8870000000000 16.0530000000000 16.1570000000000 16.3650000000000 16.5940000000000 16.7400000000000 16.8650000000000 16.9280000000000 17.0530000000000 17.0950000000000 17.1780000000000 17.2200000000000 17.3030000000000 17.3440000000000 17.5530000000000 17.6980000000000 17.8020000000000 17.9900000000000 18.1360000000000 18.3650000000000 18.5740000000000 18.8030000000000 18.9900000000000 19.1980000000000 19.5100000000000 19.6560000000000 19.8430000000000 19.9690000000000 20.0940000000000 20.1760000000000 20.3220000000000 20.4680000000000 20.5930000000000 20.7590000000000 20.9050000000000 21.0510000000000 21.1970000000000 21.4260000000000 21.5920000000000 21.7170000000000 21.9050000000000 22.0300000000000 22.2180000000000 22.3430000000000 22.5510000000000 22.8420000000000 23.0920000000000 23.1960000000000 23.4250000000000 23.5920000000000 23.7170000000000 23.8630000000000 24.0500000000000 24.1960000000000 24.3210000000000 24.5910000000000 24.7370000000000 24.8620000000000 25.0490000000000 25.2580000000000 25.5490000000000 25.7160000000000 25.8830000000000 26.0910000000000 26.2370000000000 26.3200000000000 26.6960000000000 27.4040000000000 27.6120000000000 27.7790000000000 27.9870000000000 28.1540000000000 28.3620000000000 28.5280000000000 28.6950000000000 28.8830000000000 29.2370000000000 29.4450000000000 29.5490000000000 29.7990000000000 29.9650000000000 30.1320000000000 30.3190000000000 30.5070000000000 30.6950000000000 30.8200000000000 30.9860000000000 31.1320000000000 31.3200000000000 31.3820000000000 31.4870000000000 31.5910000000000 31.7580000000000 31.8610000000000 31.9870000000000 32.0700000000000 32.3610000000000 32.4860000000000 32.6320000000000 32.8810000000000 33.0060000000000 33.0480000000000 33.0690000000000 33.0900000000000 33.0900000000000 33.0690000000000 33.0270000000000 32.9240000000000 32.9030000000000 32.8610000000000 32.8820000000000 32.9030000000000 32.9660000000000 33.0490000000000 33.0910000000000 33.1740000000000 33.2570000000000 33.3200000000000 33.3820000000000 33.4660000000000 33.4440000000000 33.4230000000000 33.4020000000000 33.4230000000000 33.3810000000000 33.4020000000000 33.4020000000000 33.4020000000000 33.4020000000000 33.4020000000000 33.4440000000000 33.4640000000000 33.5270000000000 33.5270000000000 33.6110000000000 33.6520000000000 33.6730000000000 33.6930000000000 33.6740000000000 33.6740000000000 33.6320000000000 33.5280000000000 33.4650000000000 33.3820000000000 33.4030000000000 33.3610000000000 33.3400000000000 33.3200000000000 33.2570000000000 33.3410000000000 33.3820000000000 33.4240000000000 33.4440000000000 33.4650000000000 33.4440000000000 33.5060000000000 33.5060000000000 33.5480000000000 33.5480000000000 33.5480000000000 33.5690000000000 33.6100000000000 33.7350000000000 33.8390000000000 33.9220000000000 34.0270000000000 34.1310000000000 34.3180000000000 34.5480000000000 34.7350000000000 34.8810000000000 35.0480000000000 35.2140000000000 35.4440000000000 35.5900000000000 35.7150000000000 36.0690000000000 36.1730000000000 36.3400000000000 36.5270000000000 36.5900000000000 36.7150000000000 36.9640000000000 37.0270000000000 37.1730000000000 37.3190000000000 37.4640000000000 37.6100000000000 37.7770000000000 37.9230000000000 38.0680000000000 38.1520000000000 38.2980000000000 38.4230000000000 38.5480000000000 38.5480000000000 38.5680000000000 38.6100000000000 38.6520000000000 38.7140000000000 38.7770000000000 38.9230000000000 39.0900000000000 39.2360000000000 39.3810000000000 39.7140000000000 39.8390000000000 40.0060000000000 40.1310000000000 40.2560000000000 40.3810000000000 40.4640000000000 40.5270000000000 40.6100000000000 40.6930000000000 40.7550000000000 40.8800000000000 41.0260000000000 41.1300000000000 41.2350000000000 41.3810000000000 41.5260000000000 41.6510000000000 41.7560000000000 41.9220000000000 42.0890000000000 42.2140000000000 42.4850000000000 42.6520000000000 42.8600000000000 42.9650000000000 43.1520000000000 43.2560000000000 43.4020000000000 43.5060000000000 43.7140000000000 43.7560000000000 43.8390000000000 43.9010000000000 43.9850000000000 44.1100000000000 44.2140000000000 44.3380000000000 44.5470000000000 44.6720000000000 44.8380000000000 44.9630000000000 45.1300000000000 45.3380000000000 45.4420000000000 45.6500000000000 45.7540000000000 46.0870000000000 46.1910000000000 46.3580000000000 46.4620000000000 46.5460000000000 46.6290000000000 46.6710000000000 46.7540000000000 46.8170000000000 46.8790000000000 47.0460000000000 47.1080000000000 47.2130000000000 47.3590000000000 47.6710000000000 47.9000000000000 48.0660000000000 48.2540000000000 48.4210000000000 48.5460000000000 48.7120000000000 48.7950000000000 48.8160000000000 48.9210000000000 48.9830000000000 49.0250000000000 49.0670000000000 49.0870000000000 49.2120000000000 49.2960000000000 49.4210000000000 49.4840000000000 49.7120000000000 49.9410000000000 50.0450000000000 50.6070000000000 50.7320000000000 50.8370000000000 50.9200000000000 51.0040000000000 51.1290000000000 51.1910000000000 51.2960000000000 51.4210000000000 51.4830000000000 51.5250000000000 51.6080000000000 51.7130000000000 51.7540000000000 51.8170000000000 51.9210000000000 52.0040000000000 52.0660000000000 52.1910000000000 52.2950000000000 52.4620000000000 52.5450000000000 52.6910000000000 52.8160000000000 53.0650000000000 53.2320000000000 53.2950000000000 53.3980000000000 53.4400000000000 53.5240000000000 53.6280000000000 53.6490000000000 53.7530000000000 53.8570000000000 53.9200000000000 54.0660000000000 54.1490000000000 54.2730000000000 54.4190000000000 54.5230000000000 54.6280000000000 54.7320000000000 54.8160000000000 54.8990000000000 55.0240000000000 55.0860000000000 55.1490000000000 55.2320000000000 55.2950000000000 55.4200000000000 55.5030000000000 55.7110000000000 55.8150000000000 56.0650000000000 56.1910000000000 56.4200000000000 56.5030000000000 56.6070000000000 56.6910000000000 56.7530000000000 56.8570000000000 56.9830000000000 57.3160000000000 57.4610000000000 57.6280000000000 57.7730000000000 57.9610000000000 58.1270000000000 58.5020000000000 58.6690000000000 58.8150000000000 58.9190000000000 59.0650000000000 59.1690000000000 59.2520000000000 59.4190000000000 59.4190000000000 59.5230000000000 59.5440000000000 59.6280000000000 59.7530000000000 59.8160000000000 59.9410000000000 60.1490000000000 60.2950000000000 60.4610000000000 61.1270000000000 61.2530000000000 61.5020000000000 61.9390000000000 62.0850000000000 62.2100000000000 62.3770000000000 62.7100000000000 62.8980000000000 63.0440000000000 63.2940000000000 63.4820000000000 63.7110000000000 63.8980000000000 64.1060000000000 64.3150000000000 64.5020000000000 64.6900000000000 64.8770000000000 65.0850000000000 65.2100000000000 65.3150000000000 65.4400000000000 65.6270000000000 65.7310000000000 65.8560000000000 66.0230000000000 66.1270000000000 66.2510000000000 66.3560000000000 66.6470000000000 66.7510000000000 66.8760000000000 66.9380000000000 67.1890000000000 67.3140000000000 67.3970000000000 67.5220000000000 67.5640000000000 67.7310000000000 67.8560000000000 67.9810000000000 68.0650000000000 68.1480000000000 68.2730000000000 68.3980000000000 68.4810000000000 68.6060000000000 68.6680000000000 68.8140000000000 68.8970000000000 69.0220000000000 69.2090000000000 69.2920000000000 69.6460000000000 69.6670000000000 69.8340000000000 70.0840000000000 70.1890000000000 70.3140000000000 70.4810000000000 70.6050000000000 70.7510000000000 70.8760000000000 71.0430000000000 71.1470000000000 71.2720000000000 71.4390000000000 71.7100000000000 71.8770000000000 71.9810000000000 72.1470000000000 72.2720000000000 72.3550000000000 72.5010000000000 72.5430000000000 72.6470000000000 72.7520000000000 72.8140000000000 72.9610000000000 73.0440000000000 73.1690000000000 73.3140000000000 73.4810000000000 74.0430000000000 74.1890000000000 74.3980000000000 74.5220000000000 74.6050000000000 74.7100000000000 74.7300000000000 74.7720000000000 74.8130000000000 74.8330000000000 74.8960000000000 74.9370000000000 75.0210000000000 75.0620000000000 75.2080000000000 75.3540000000000 75.7080000000000 75.8750000000000 76.0630000000000 76.2300000000000 76.4380000000000 76.5430000000000 76.7090000000000 76.8760000000000 77.0010000000000 77.1260000000000 77.2720000000000 77.3770000000000 77.4800000000000 77.6050000000000 77.8550000000000 77.9600000000000 78.0640000000000 78.1060000000000 78.2100000000000 78.3970000000000 78.5220000000000 78.6270000000000 78.8550000000000 79.0430000000000 79.2300000000000 79.6470000000000 79.8760000000000 80.5220000000000 80.6260000000000 80.7930000000000 80.8550000000000 80.8970000000000 80.9180000000000 80.9390000000000 80.9810000000000 81.0430000000000 81.1270000000000 81.2730000000000 81.4600000000000 81.6270000000000 81.8350000000000 82.1260000000000 82.6260000000000 82.8760000000000 83.1050000000000 83.2720000000000 83.4390000000000 83.5850000000000 83.6890000000000 83.7310000000000 83.7100000000000 83.6890000000000 83.6890000000000 83.7100000000000 83.8350000000000 83.9180000000000 84.0850000000000 84.4600000000000 84.6890000000000 84.9180000000000 85.3970000000000 85.5220000000000 85.6260000000000 85.8560000000000 85.8550000000000 85.8560000000000 85.8760000000000 85.9380000000000 85.9390000000000 86.0430000000000 86.1470000000000 86.2510000000000 86.4180000000000 86.6470000000000 86.8550000000000 87.0850000000000 87.3340000000000 87.5630000000000 87.7720000000000 88.0630000000000 88.1890000000000 88.4390000000000 88.5430000000000 88.6480000000000 88.6480000000000 88.7100000000000 88.7730000000000 88.8770000000000 88.8980000000000 89.0010000000000 89.1050000000000 89.2930000000000 89.4180000000000 89.6470000000000 89.7930000000000 90.0010000000000 90.2100000000000 90.4180000000000 90.6060000000000 90.7100000000000 90.8770000000000 91.0640000000000 91.1680000000000 91.2720000000000 91.3560000000000 91.4180000000000 91.5220000000000 91.6470000000000 91.7310000000000 91.8350000000000 91.9180000000000 92.1060000000000 92.2310000000000 92.3760000000000 92.6060000000000 92.7930000000000 93.0010000000000 93.4180000000000 93.6470000000000 93.8970000000000 94.4180000000000 94.5840000000000 94.7300000000000 95.0220000000000 95.1890000000000 95.3140000000000 95.5220000000000 95.8140000000000 95.9180000000000 96.0640000000000 96.2300000000000 96.3130000000000 96.5010000000000 96.6260000000000 96.7720000000000 96.8560000000000 97.0230000000000 97.1680000000000 97.2720000000000 97.3970000000000 97.5430000000000 97.6680000000000 97.7520000000000 97.8770000000000 97.9600000000000 98.0650000000000 98.3360000000000 98.4190000000000 98.5230000000000 98.6070000000000 98.7320000000000 98.8150000000000 98.9190000000000 99.1060000000000 99.2520000000000 99.4190000000000 99.5860000000000 99.7940000000000 99.9190000000000 100.044000000000];
z_diff = (diff(z)) ;
z_diff = [0 z_diff];
tilt_angle=10;
beta=(pi/180)*tilt_angle;
Lx = 0.1;
Ly = Lx;
Nx = 512;
Ny = Nx ;
dx = Lx/Nx;
dy = Ly/Ny;
x = (-Nx/2:Nx/2-1)*dx;
y = (-Ny/2:Ny/2-1)*dy;
[X,Y] = meshgrid(x,y);
r = sqrt(X.^2 + Y.^2);
kx = meshgrid(-Nx/2:Nx/2-1)/(Nx.*dx);
ky = meshgrid(-Ny/2:Ny/2-1)/(Ny.*dy);
kz = ones(Nx);
Kx=cos(beta).*kx + sin(beta).*kz;
Ky=ky;
Kz=-sin(beta).*kx + cos(beta).*kz;
K2 = Kx.^2 + Ky.^2 + Kz.^2;
delfx = 1/(Lx);
delfy = delfx;
R = randn(Nx,Ny) ;
R = R - mean(R(:,:));
R = R/std(R(:,:));
Phi = [2.50310494762771e-07 2.50946965543003e-07 4.74127533715154e-07 5.24210357233832e-06 2.42362427581638e-07 1.42653443399616e-07 5.16930235276850e-07 9.04590700562310e-07 2.29541191637123e-06 3.60105794008093e-07 9.17845719860523e-08 2.50390367004531e-07 1.09738615765455e-06 1.00311041468925e-06 6.50810613011221e-08 3.02640252523566e-07 5.31047019902531e-07 1.80313245178286e-06 3.04099822876624e-06 1.47008194185741e-07 4.08290582950848e-08 3.82118426499270e-06 5.29334397036188e-07 4.87810778598663e-06 1.18782521297988e-06 1.16666772914768e-06 1.69778706192608e-06 7.67426466736916e-06 5.27799451227444e-06 4.23216916272082e-07 2.66828750584071e-05 2.46205729358867e-06 2.80169965568206e-06 2.97378024567976e-07 2.21170607901861e-07 3.41285818450589e-07 3.85275104621375e-08 4.05141179263807e-07 1.17629321652161e-06 4.95233464403551e-06 1.80125151270180e-06 2.81483058944610e-07 1.24528365857859e-05 6.86781782456521e-06 2.95828059847895e-06 1.27035361174474e-06 6.34388599513081e-07 2.32012818083199e-06 2.81784487117724e-06 1.46096455213802e-06 1.60178082889894e-06 1.94791694377124e-07 6.35459501021705e-08 8.30940388112504e-08 2.53354933929590e-07 4.31247522543152e-06 9.63180625859285e-07 8.03726503455623e-07 2.98282116276067e-07 9.25402679823575e-08 1.61856374619784e-07 1.73242501417858e-07 2.31570160754109e-07 1.78600321810053e-07 1.25598618058698e-07 1.63418364388577e-06 1.02774042441542e-06 8.32018926678718e-08 2.53586554087675e-07 2.53508325030189e-07 4.07413642129177e-08 2.42096085408459e-07 5.91587468209778e-07 1.27761778023523e-07 1.65773731666773e-07 1.79565395840267e-07 3.78521468221561e-07 1.21471607158554e-07 3.25378687175317e-07 1.25569084575180e-07 5.49035746100893e-08 4.07733632799329e-08 4.79103771763424e-07 3.27135102365317e-07 2.10496092551324e-07 1.33641668807187e-06 1.68227872575447e-06 1.52265748820018e-06 2.53683706505161e-07 1.25630494723704e-07 1.12347756326375e-06 2.53692684133802e-07 2.10500428113819e-07 2.52934873659000e-07 6.77249701528946e-07 2.95386782896753e-07 3.16617852613830e-07 3.59879226900604e-07 2.81244834536819e-07 2.35438339516804e-07 3.40806542459708e-07 6.44683615098099e-08 5.05945401350734e-07 2.57387924440325e-07 7.82108421793457e-08 2.52513556939930e-07 7.41154472489689e-08 8.95664707278605e-08 3.44823854262242e-07 2.57498619329261e-07 2.07939465354519e-06 7.18570183290631e-07 3.25198024230522e-07 2.11933525065655e-07 8.30783371984955e-08 1.65710628194624e-07 5.48396993890243e-08 5.95342783012939e-08 2.28819020334242e-07 1.99067722101735e-07 2.40425372248607e-07 2.51916559208949e-07 2.98446361528618e-07 2.82958505602994e-07 2.10599141757522e-07 4.60997674322361e-08 2.59856693650578e-07 1.47642122461920e-07 3.93574284250877e-07 3.94729289526600e-08 6.37347851352551e-08 2.51536627057736e-07 2.28438571372065e-07 1.90144303446780e-07 4.73626330414131e-08 6.44142268235402e-08 1.72106272576708e-06 1.70721924857358e-06 2.94643811973548e-07 3.55943119776846e-08 9.36518317402136e-08 4.31201179335937e-08 1.56428758776995e-07 3.19975996745435e-07 3.34275554966361e-07 1.20064760586921e-06 4.30984033789707e-07 2.15616245567332e-07 3.03074660779611e-07 3.29934827134432e-07 -2.24864187842034e-08 2.33412514429474e-07 3.21092520408291e-07 3.38617740923693e-07 6.04579881670076e-07 1.91481555229955e-07 9.38472397526705e-08 6.13052684855612e-07 3.97260343618278e-06 2.45355935394897e-07 3.47585309986832e-07 3.38923893338132e-07 2.20926360286704e-07 5.25446559955282e-07 2.84158336921376e-07 4.04755442008443e-07 3.35541338245773e-07 3.76856346223440e-07 2.75568777774655e-07 2.09025928742013e-07 7.10170870233622e-08 2.95787195332815e-07 5.59078661720798e-07 2.51150081156080e-07 4.02707325242296e-08 8.24148001852386e-08 2.57363940284162e-07 3.55816103375846e-07 3.95777597622744e-07 3.27993323594138e-07 3.52690891300866e-07 3.52278318886910e-07 5.06712192595194e-08 1.23808524289852e-07 8.39347385990663e-07 2.02892942754780e-07 7.97641133582245e-08 3.64530615247902e-07 9.11185642215493e-07 2.44511041473120e-07 3.24963966176285e-07 1.11685979674259e-07 2.86261990645848e-07 3.75634147009314e-08 2.16925989302430e-07 3.84116828662608e-07 2.71791109589261e-07 1.06845659312697e-07 2.02225581949648e-07 3.02535099485185e-07 3.11957552477906e-08 9.15881915538412e-08 9.18891121954120e-07 3.70106500980225e-07 2.25779726927071e-07 2.86067288473873e-07 5.66382866246259e-07 4.26632108086086e-07 4.53045634671475e-07 1.17973378714372e-06 3.40969918417227e-07 2.48643473538260e-07 2.64505922988546e-07 2.38444406096982e-07 3.47101823144428e-07 3.37625374999722e-07 3.20790426924069e-07 5.02178252920690e-08 1.73576015527609e-07 8.57402376593969e-07 1.80131811857122e-07 2.76964595681315e-07 2.14838038341036e-07 2.21053263051776e-07 1.59406260405821e-07 4.59614264881235e-08 3.78789455759776e-07 2.58966590318316e-07 7.75696011578299e-08 3.39990416320565e-07 3.24572509489827e-07 3.52020828992088e-08 6.49912862168267e-08 2.22828739301052e-07 2.13473046057573e-07 1.70720547021041e-07 2.19582224668132e-07 2.05914598785446e-07 1.14899310738254e-06 1.03415047080757e-07 1.89177690522258e-06 2.24175047702725e-07 3.76150153666139e-07 9.18124634120261e-08 2.27678254839140e-07 8.61346207114181e-08 3.38947999645138e-06 1.71880093129416e-07 2.22461873169689e-07 3.72427593494344e-08 6.76102058093114e-08 1.13163515787091e-06 7.04052219453715e-07 3.28310827940334e-07 2.12060821359437e-07 2.68576701610855e-07 2.04744061114503e-07 1.14622616867334e-07 4.35675529886048e-08 5.29665131717990e-08 3.50574437154752e-08 5.32920831968263e-08 1.04801149974382e-07 2.08316009852988e-07 2.77349913208742e-07 1.48756450473382e-07 3.78517042358300e-07 1.19222278136494e-07 4.57902170141890e-08 2.03249187200450e-07 9.32437812713631e-08 6.61503340988668e-07 1.01572927594225e-06 3.26294032932256e-07 2.26313604682662e-07 7.01012285236915e-08 7.61666326540900e-08 4.16954886967803e-08 4.41719683698357e-08 4.02549730011751e-08 4.41603434482274e-08 3.50371730851861e-08 3.00718061774876e-07 3.50372423542925e-08 5.32741630316562e-08 2.09046711312522e-07 3.75611845531339e-06 2.26535469768828e-07 4.62710159517236e-07 5.55357457448765e-07 4.16702305525217e-07 2.05922476618442e-07 2.06490274401431e-07 2.33754328058605e-07 1.66615319203810e-07 2.13390792614979e-07 2.54792444311713e-07 4.20241426628409e-07 6.95047952441784e-07 6.02254152247896e-07 1.56218004764843e-07 1.99836446669293e-07 1.62974708502310e-07 2.93420921132526e-07 1.67909832369442e-07 1.53692619622575e-07 8.89709805825161e-08 1.64148574413349e-07 1.82871518257336e-07 1.99566615308843e-07 1.59784291140960e-07 1.87893662090374e-07 2.77364481028229e-07 1.10443572547260e-07 2.18556580353372e-07 1.83932131011210e-07 7.94652011731045e-08 2.43919458395486e-07 1.59781205439559e-07 1.94642054497664e-07 1.29989372686070e-07 1.90910393259075e-07 2.84452061601644e-07 3.30006526167654e-07 1.93485493531741e-07 2.40209860863182e-07 4.15333943239071e-07 2.91343651778087e-07 2.58377576473351e-07 1.81274762632482e-07 1.88893039486970e-07 2.19428417763272e-07 2.20994023873646e-07 2.27950728874121e-07 1.64400439994395e-07 1.00341668297253e-07 1.67577483763212e-07 5.61980802260245e-08 3.04539633123467e-08 3.19018898206730e-08 3.89504322796163e-08 3.19941498021468e-07 5.41874921040340e-06 1.01776790218898e-06 1.29079721421378e-06 5.77887058141474e-07 2.79939362285509e-07 5.55989804582382e-08 2.11144624551603e-07 5.73140726411827e-08 5.35992572193921e-08 2.80066251579433e-07 1.76964657359828e-07 3.25199602895421e-07 2.18271466195228e-07 1.80628024038326e-07 4.69550563668357e-08 1.48055689005769e-07 1.82317859358405e-07 1.33063533148151e-07 2.62672490298873e-08 2.52416541944235e-08 2.08244390466210e-07 3.06261128674104e-07 2.18183988930487e-07 1.32114003617060e-07 3.13731744862144e-07 2.78980205382911e-07 2.42110700965470e-08 1.28818767334643e-07 3.48416502648141e-07 1.43250842379859e-07 2.58790299254369e-08 2.36803673262810e-08 2.54679999490770e-08 1.53423624973069e-07 2.11044574009824e-07 2.45219073821534e-08 5.40329423903049e-08 4.07828910547019e-08 2.86866659397290e-08 5.74108869690081e-08 3.17780519921456e-08 2.89072270765441e-07 6.46668413187508e-08 3.89801393203052e-07 1.62322978572959e-07 2.69026419637955e-07 4.76120100475985e-07 5.32185045407108e-07 1.57575408820348e-07 2.08858600548085e-07 1.87565977660685e-07 2.01108310997773e-07 3.95620143987532e-08 5.24866345380986e-08 8.64930520826212e-08 9.74072997804584e-07 2.54439638933078e-08 1.00154777975012e-07 9.42429455464106e-07 2.23088555779751e-07 2.02097796370716e-07 1.99163165408265e-07 9.26536634692885e-08 1.31298738607173e-07 1.97726692783641e-07 2.28150683478771e-07 2.17170907865184e-07 2.69677497220439e-08 1.54765400027711e-07 6.26255760215395e-08 1.46771827369792e-07 1.81515822529873e-07 1.72918459779846e-07 8.15990146285526e-08 9.85514009220808e-08 7.03396855430348e-07 5.93657045438661e-07 1.51834395684521e-07 9.76927764035926e-08 1.18527867299870e-07 2.71435324116586e-08 5.68571003566081e-08 1.65983256102705e-07 2.29129793972797e-07 5.88646431511064e-08 3.08684402204026e-05 2.38265317337181e-06 4.30752744717656e-07 1.04689461805175e-07 2.32270648752108e-07 6.07211466407868e-07 1.95341105985221e-08 7.38231690963101e-08 1.95129857199575e-07 4.01671898332107e-08 4.75072551754564e-08 2.49046064060719e-06 1.78142428447807e-07 2.60935684174315e-07 4.29594269700864e-08 1.14070714922750e-07 5.26390027271695e-08 6.74132420131039e-08 2.11045265009953e-07 7.27021148761358e-07 2.83886081059763e-07 2.76343292567963e-08 1.61294742863238e-07 1.83304573271511e-07 2.45251911632262e-07 9.18095053702886e-07 9.40457085677553e-07 2.90451584499496e-07 1.04003175628026e-07 1.08284695915254e-07 5.44402482199441e-08 2.02578013696944e-08 5.88615942025308e-08 1.04406254724335e-07 1.20598226687063e-07 5.71264670628078e-07 3.15144191075623e-07 2.23052176438872e-08 2.98936409995515e-08 2.92428609360162e-07 2.47505388371339e-07 2.72866577624107e-07 2.69881858030521e-07 1.85519336468060e-07 2.56402576998181e-07 8.88653855726949e-08 1.23281438816645e-07 3.74059740200056e-08 3.63466577344003e-08 1.20887062612203e-06 5.31728280982040e-07 2.22245937870703e-08 1.24471844746278e-07 5.11735373064063e-07 3.13008213548308e-07 5.70247343660301e-08 7.23194564373494e-07 1.79616235944090e-07 6.78737172874271e-08 1.87999018332358e-07 1.56051038058590e-07 2.31029093824182e-08 3.65957708093627e-08 5.13581582166058e-07 2.55008648105981e-08 1.44166238137747e-07 3.81609462756446e-07 2.16242665285891e-07 2.34087134829968e-07 2.55408939395098e-07 5.31844204547311e-08 1.37230800917903e-07 2.24923680386313e-07 7.91242025449375e-08 1.56574476334133e-07 1.34926252415438e-07 1.70013819330808e-08 1.48130847036606e-07 1.10651211565701e-07 1.58123074045499e-07 1.47865663105776e-07 1.14316450201195e-07 2.02955526741989e-08 1.49518466142608e-07 2.99469635848321e-08 4.06175629714119e-07 2.32871553385763e-07 1.70948094081342e-07 1.86984334344187e-07 1.54774299548728e-07 1.27005955203776e-07 7.29605614875902e-08 4.26639317987837e-08 5.91930531596668e-08 1.50666658228087e-08 4.86748534643297e-08 1.64538615062307e-08 1.32252962938291e-07 1.06729078335828e-07 1.39785078669013e-07 1.39779800669610e-07 9.26422941859900e-08 1.13233915195479e-07 8.32946711657419e-08 1.67451405167575e-07 3.65621117891888e-08 2.21661117907828e-07 3.65652178907218e-08 1.41441861581090e-06 3.30075312403152e-07 2.05346680636497e-08 2.16432978289422e-07 1.92694668921994e-07 1.99379501668891e-07 2.60904461669433e-07 2.16709099204294e-07 1.33136064843243e-07 1.52413473502972e-07 7.53889226907776e-07 3.22036453180506e-07 1.68931528429636e-07 1.04094510285934e-07 1.60098934336117e-07 1.02657996503878e-07 6.83661821001145e-08 1.61931255758191e-07 2.01320024686194e-07 1.07534006069169e-07 1.07910634894973e-07 9.46522813533042e-08 1.66743437525985e-07 1.78365196465516e-07 1.90349497426267e-07 1.99051143944893e-07 1.10969246755729e-07 1.11264782143939e-07 1.02006294022338e-07 1.24767050443006e-07 1.04932154146443e-07 1.03659587531759e-07 4.90478970104965e-08 9.01107885624077e-08 9.30802687007756e-08 2.06842405795697e-07 1.53448550376276e-07 9.79061221486547e-08 8.26140384317179e-08 1.50146254483116e-07 1.40136146748372e-07 1.56837338087257e-07 1.76198808159403e-07 1.99122966368083e-07 1.23858764266714e-07 1.06727895602262e-07 8.15204687804434e-08 1.88227798887152e-08 2.31668781894091e-08 7.67325491571661e-07 1.97158231154847e-07 1.04232262587474e-07 8.81089549804027e-09 1.09695696650620e-08 3.67828021042941e-07 2.42873900856430e-07 9.30400519863338e-08 1.22815336391232e-08 1.04498084973760e-07 4.13726864269242e-07 1.28584375991528e-08 1.20376581165868e-07 1.07046398186240e-07 5.96207309583810e-08 3.66763563633301e-08 1.21422552113772e-07 1.39822102904246e-07 1.58475758239753e-07 8.54716186106630e-08 9.47107273229749e-08 1.57282457487923e-07 2.22341833330773e-07 1.47297286955243e-07 2.19435143761748e-07 1.39688617089499e-07 7.24259126481027e-08 7.09641458863346e-08 1.19226547380480e-07 7.87372467318109e-08 9.72314598010880e-08 7.74647075656411e-09 9.47412998435404e-08 1.16751285273194e-07 1.13487779361823e-07 1.13827758524977e-07 2.72522446715463e-07 8.16170870703966e-08 8.72321297381038e-08 9.95784392252335e-08 2.65989252708394e-09 1.06630932742284e-07 9.36921419305549e-08 8.41656329617571e-08 1.30269792124746e-08 4.86904164038956e-06 1.01469686899418e-07 1.47172552045969e-07 1.11949067863535e-07 9.34802107168088e-08 6.41098286680338e-08 6.68095594366830e-09 4.65659962399765e-08 1.74321722504626e-07 1.53177157686198e-07 1.72679939394816e-07 9.26791009109865e-08 7.38280843033631e-09 1.21431764801520e-08 1.83474589517062e-08 5.49758746890643e-08 5.49686356365556e-08 5.41396143667257e-07 1.02078716839878e-08 1.46280082306786e-07 9.20039209747377e-08 1.07464173701763e-08 4.37249887631680e-08 1.20696543908951e-07 7.37023067237918e-08 1.54370995646110e-07 8.75880053587604e-08 1.38972473988507e-07 1.50011697718075e-07 1.13900910702895e-07 7.30960487245964e-08 1.06352204500379e-07 7.15892552073097e-09 1.26382007766409e-08 2.34492999966371e-07 9.25652500580458e-08 8.91104525001921e-08 1.04346566674604e-08 7.96833204055166e-08 3.41626685510968e-07 1.46859973263572e-07 6.99097620123092e-08 5.60565018753087e-08 9.18716519344181e-08 5.82078430390537e-09 8.51757309259611e-08 1.41028716092150e-07 7.10210717526897e-09 3.02205722277794e-08 2.76205242367368e-07 9.08403872147274e-08 8.55641175181640e-08 8.07674139440669e-08 5.64950026938399e-09 1.47871928050032e-07 1.37623718221666e-07 9.60980796219824e-08 9.05395061229061e-09 3.58713661087815e-07 3.78246981154481e-07 1.97627162656081e-08 4.78463482905684e-08 5.47976120635378e-07 5.98012738422856e-07 1.38288635243591e-07 3.97209746106909e-08 1.04405402152745e-07 1.64230983323839e-08 2.14721179480984e-08 4.85083210883807e-08 1.22646011220379e-06 8.39456826640735e-09 8.27662116573820e-08 3.58891263221149e-07 9.16048888640066e-08 1.21278839928491e-07 2.31054125718236e-07 1.00889365297461e-07 1.19358834432565e-07 2.43427211256917e-07 9.11020950715102e-08 1.38162844699562e-07 8.00915682468241e-08 4.98258405386454e-09 7.30080219716715e-08 1.52701681888474e-07 2.93328698184570e-07 1.35353313498113e-07 4.59682925254202e-08 9.08782411223864e-08 2.20927979357592e-08 9.05503229704100e-08 6.68343694641667e-08 7.33025736057753e-08 7.49958147508063e-08 6.22824091009223e-09 8.02682015058001e-06 3.67877032868707e-07 7.19496678750717e-08 3.17758490462365e-09 5.75563157661276e-09 3.24272951280416e-09 1.00381201383279e-07 1.88260239049035e-07 7.07194739983604e-08 4.68947848174950e-09 4.16252458503278e-09 7.94576959843098e-08 2.74009776209153e-07 1.42430916804990e-07 1.26235913248699e-07 1.39626029880374e-07 1.16782761195100e-07 1.00256563902946e-07 -2.95101260456139e-09 6.83432403678117e-08 1.27348806727558e-07 5.96766339731236e-08 1.03604721020808e-07 1.21085249476580e-07 6.80561072346092e-08 2.39948064835149e-09 1.00688731186191e-08 2.43035685127838e-08 4.13254891501583e-08 9.63544005596039e-09 2.53883818143423e-07 7.44890597640385e-07];
w0=0.002;
F=2*(pi^2)*(k^2).*z_diff.*Phi;
x0 = linspace(-0.02,0.02,Nx);y0=x0;[X0,Y0] = meshgrid(x0,y0);
Uin=exp(-(X0.^2+Y0.^2)/w0^2);
x = x + z_diff(1)*tan(beta);
screen1 = exp(1i*((Kx.*x) + (Ky.*y))).*delfx.*delfy.*sqrt(F(1)).*R ;
Uout = Uin.*exp(1i.*screen1);
c=[0.740664357812470 0.715162315217582 0.689780546528124 0.664556687838984 0.639501229779350 0.614618312979130 0.589909677254809 0.565375929219662 0.541017064806704 0.516832720551119 0.492822307973243 0.468985091286926 0.445320235117717 0.421826835235231 0.398503939126698 0.375350560217132 0.352365687964806 0.329548295192595 0.306897343515746 0.284411787427262 0.262090577416753 0.239932662380474 0.217936991503006 0.196102515739251 0.174428188990141 0.152912969040777 0.131555818312352 0.110355704466664 0.0893116008928944 0.0684224870996102 0.0476873490298864 0.0271051793136574 0.00667497746850089 0.0136042499421828 0.0337334891863653 0.0537137192700918 0.0735459118688980 0.0932310312800026 0.112770034392025 0.132163870669524 0.151413482150099 0.170519803452171 0.189483761791847 0.208306277007520 0.226988261591059 0.245530620724636 0.263934252322323 0.282200047075789 0.300328888503445 0.318321653002534 0.336179209903699 0.353902421527639 0.371492143243500 0.388949223528713 0.406274504030009 0.423468819625392 0.440532998486862 0.457467862143721 0.474274225546305 0.490952897130006 0.507504678879476 0.523930366392895 0.540230748946220 0.556406609557339 0.572458725050048 0.588387866117797 0.604194797387153 0.619880277480917 0.635445059080876 0.650889888990134 0.666215508194993 0.681422651926368 0.696512049720699 0.711484425480339 0.726340497533411 0.741080978693105 0.755706576316408 0.770217992362264 0.784615923449129 0.798901060911951 0.813074090858528 0.827135694225279 0.841086546832386 0.854927319438336 0.868658677793846 0.882281282695167 0.895795790036781 0.909202850863477 0.922503111421817 0.935697213210995 0.948785793033077 0.961769483042646 0.974648910795838 0.987424699298783 1.00009746705545 1.01266782811491 1.02513639211797 1.03750376434332 1.04977054575300 1.06193733303735 1.07400471865938 1.08597329089862 1.09784363389434 1.10961632768826 1.12129194826677 1.13287106760249 1.14435425369543 1.15574207061353 1.16703507853271 1.17823383377643 1.18933888885472 1.20035079250267 1.21127008971854 1.22209732180122 1.23283302638735 1.24347773748787 1.25403198552414 1.26449629736359 1.27487119635488 1.28515720236266 1.29535483180181 1.30546459767133 1.31548700958769 1.32542257381781 1.33527179331163 1.34503516773420 1.35471319349739 1.36430636379121 1.37381516861470 1.38324009480640 1.39258162607448 1.40184024302645 1.41101642319850 1.42011064108442 1.42912336816425 1.43805507293243 1.44690622092571 1.45567727475060 1.46436869411057 1.47298093583279 1.48151445389465 1.48996969944983 1.49834712085410 1.50664716369079 1.51487027079592 1.52301688228302 1.53108743556759 1.53908236539135 1.54700210384610 1.55484708039729 1.56261772190731 1.57031445265850 1.57793769437581 1.58548786624928 1.59296538495611 1.60037066468257 1.60770411714556 1.61496615161393 1.62215717492952 1.62927759152798 1.63632780345927 1.64330821040796 1.65021920971321 1.65706119638862 1.66383456314170 1.67053970039319 1.67717699629614 1.68374683675467 1.69024960544263 1.69668568382193 1.70305545116068 1.70935928455115 1.71559755892740 1.72177064708287 1.72787891968756 1.73392274530517 1.73990249040993 1.74581851940328 1.75167119463032 1.75746087639608 1.76318792298160 1.76885269065981 1.77445553371120 1.77999680443936 1.78547685318630 1.79089602834756 1.79625467638721 1.80155314185264 1.80679176738913 1.81197089375438 1.81709085983270 1.82215200264920 1.82715465738367 1.83209915738443 1.83698583418192 1.84181501750219 1.84658703528021 1.85130221367304 1.85596087707287 1.86056334811985 1.86510994771484 1.86960099503199 1.87403680753115 1.87841770097023 1.88274398941727 1.88701598526258 1.89123399923051 1.89539834039132 1.89950931617273 1.90356723237146 1.90757239316461 1.91152510112087 1.91542565721170 1.91927436082231 1.92307150976254 1.92681740027765 1.93051232705897 1.93415658325440 1.93775046047891 1.94129424882478 1.94478823687186 1.94823271169764 1.95162795888723 1.95497426254328 1.95827190529575 1.96152116831153 1.96472233130413 1.96787567254306 1.97098146886325 1.97403999567435 1.97705152696986 1.98001633533631 1.98293469196218 1.98580686664685 1.98863312780943 1.99141374249746 1.99414897639556 1.99683909383401 1.99948435779718 2.00208502993198 2.00464137055607 2.00715363866619 2.00962209194621 2.01204698677524 2.01442857823560 2.01676712012071 2.01906286494297 2.02131606394144 2.02352696708957 2.02569582310279 2.02782287944605 2.02990838234126 2.03195257677470 2.03395570650435 2.03591801406711 2.03783974078601 2.03972112677735 2.04156241095768 2.04336383105087 2.04512562359498 2.04684802394913 2.04853126630031 2.05017558367013 2.05178120792145 2.05334836976503 2.05487729876606 2.05636822335072 2.05782137081251 2.05923696731876 2.06061523791685 2.06195640654055 2.06326069601621 2.06452832806891 2.06575952332861 2.06695450133616 2.06811348054936 2.06923667834885 2.07032431104407 2.07137659387907 2.07239374103836 2.07337596565262 2.07432347980443 2.07523649453394 2.07611521984447 2.07695986470806 2.07777063707106 2.07854774385950 2.07929139098462 2.08000178334822 2.08067912484798 2.08132361838281 2.08193546585808 2.08251486819086 2.08306202531507 2.08357713618662 2.08406039878854 2.08451201013601 2.08493216628138 2.08532106231914 2.08567889239087 2.08600584969017 2.08630212646748 2.08656791403493 2.08680340277114 2.08700878212594 2.08718424062516 2.08732996587524 2.08744614456794 2.08753296248490 2.08759060450230 2.08761925459532 2.08761909584274 2.08759031043135 2.08753307966046 2.08744758394630 2.08733400282638 2.08719251496389 2.08702329815202 2.08682652931821 2.08660238452849 2.08635103899165 2.08607266706347 2.08576744225093 2.08543553721630 2.08507712378131 2.08469237293121 2.08428145481885 2.08384453876873 2.08338179328099 2.08289338603538 2.08237948389527 2.08184025291152 2.08127585832643 2.08068646457760 2.08007223530177 2.07943333333869 2.07876992073487 2.07808215874743 2.07737020784777 2.07663422772538 2.07587437729149 2.07509081468279 2.07428369726507 2.07345318163688 2.07259942363312 2.07172257832864 2.07082280004183 2.06990024233815 2.06895505803365 2.06798739919850 2.06699741716044 2.06598526250827 2.06495108509530 2.06389503404273 2.06281725774310 2.06171790386362 2.06059711934957 2.05945505042762 2.05829184260918 2.05710764069365 2.05590258877174 2.05467683022873 2.05343050774771 2.05216376331278 2.05087673821232 2.04956957304210 2.04824240770851 2.04689538143169 2.04552863274867 2.04414229951646 2.04273651891521 2.04131142745123 2.03986716096009 2.03840385460964 2.03692164290309 2.03542065968197 2.03390103812916 2.03236291077188 2.03080640948460 2.02923166549207 2.02763880937220 2.02602797105897 2.02439927984538 2.02275286438630 2.02108885270135 2.01940737217775 2.01770854957318 2.01599251101859 2.01425938202100 2.01250928746633 2.01074235162217 2.00895869814052 2.00715845006058 2.00534172981146 2.00350865921493 2.00165935948813 1.99979395124624 1.99791255450517 1.99601528868428 1.99410227260896 1.99217362451332 1.99022946204283 1.98826990225688 1.98629506163147 1.98430505606170 1.98230000086442 1.98028001078079 1.97824519997880 1.97619568205582 1.97413157004115 1.97205297639850 1.96996001302854 1.96785279127132 1.96573142190885 1.96359601516744 1.96144668072028 1.95928352768979 1.95710666465011 1.95491619962948 1.95271224011267 1.95049489304337 1.94826426482658 1.94602046133099 1.94376358789132 1.94149374931069 1.93921104986296 1.93691559329504 1.93460748282926 1.93228682116561 1.92995371048410 1.92760825244698 1.92525054820109 1.92288069838009 1.92049880310668 1.91810496199492 1.91569927415240 1.91328183818251 1.91085275218661 1.90841211376627 1.90596002002548 1.90349656757275 1.90102185252339 1.89853597050162 1.89603901664270 1.89353108559514 1.89101227152279 1.88848266810698 1.88594236854864 1.88339146557039 1.88083005141866 1.87825821786576 1.87567605621199 1.87308365728765 1.87048111145519 1.86786850861116 1.86524593818833 1.86261348915771 1.85997125003053 1.85731930886031 1.85465775324484 1.85198667032819 1.84930614680269 1.84661626891091 1.84391712244766 1.84120879276190 1.83849136475875 1.83576492290140 1.83302955121309 1.83028533327899 1.82753235224814 1.82477069083540 1.82200043132329 1.81922165556395 1.81643444498099 1.81363888057138 1.81083504290734 1.80802301213817 1.80520286799215 1.80237468977836 1.79953855638852 1.79669454629885 1.79384273757188 1.79098320785825 1.78811603439854 1.78524129402510 1.78235906316378 1.77946941783578 1.77657243365938 1.77366818585177 1.77075674923075 1.76783819821655 1.76491260683352 1.76198004871193 1.75904059708966 1.75609432481398 1.75314130434320 1.75018160774846 1.74721530671537 1.74424247254577 1.74126317615938 1.73827748809551 1.73528547851471 1.73228721720050 1.72928277356095 1.72627221663041 1.72325561507115 1.72023303717496 1.71720455086485 1.71417022369661 1.71113012286052 1.70808431518289 1.70503286712770 1.70197584479823 1.69891331393861 1.69584533993546 1.69277198781941 1.68969332226677 1.68660940760100 1.68352030779434 1.68042608646938 1.67732680690054 1.67422253201569 1.67111332439765 1.66799924628574 1.66488035957729 1.66175672582918 1.65862840625934 1.65549546174827 1.65235795284053 1.64921593974624 1.64606948234258 1.64291864017527 1.63976347246002 1.63660403808406 1.63344039560754 1.63027260326505 1.62710071896701 1.62392480030119 1.62074490453409 1.61756108861240 1.61437340916446 1.61118192250163 1.60798668461974 1.60478775120050 1.60158517761291 1.59837901891464 1.59516932985347 1.59195616486864 1.58873957809224 1.58551962335059 1.58229635416566 1.57906982375634 1.57584008503991 1.57260719063331 1.56937119285456 1.56613214372404 1.56289009496590 1.55964509800932 1.55639720398990 1.55314646375097 1.54989292784488 1.54663664653434 1.54337766979374 1.54011604731041 1.53685182848594 1.53358506243751 1.53031579799910 1.52704408372284 1.52376996788025 1.52049349846351 1.51721472318676 1.51393368948733 1.51065044452699 1.50736503519323 1.50407750810048 1.50078790959138 1.49749628573799 1.49420268234301 1.49090714494105 1.48760971879982 1.48431044892133 1.48100938004314 1.47770655663956 1.47440202292281 1.47109582284426 1.46778800009561 1.46447859811008 1.46116766006355 1.45785522887583 1.45454134721172 1.45122605748225 1.44790940184585 1.44459142220945 1.44127216022968 1.43795165731401 1.43462995462187 1.43130709306582 1.42798311331269 1.42465805578467 1.42133196066046 1.41800486787641 1.41467681712758 1.41134784786894 1.40801799931637 1.40468731044786 1.40135582000454 1.39802356649183 1.39469058818047 1.39135692310765 1.38802260907807 1.38468768366504 1.38135218421152 1.37801614783119 1.37467961140955 1.37134261160493 1.36800518484960 1.36466736735075 1.36132919509161 1.35799070383243 1.35465192911156 1.35131290624645 1.34797367033472 1.34463425625514 1.34129469866868 1.33795503201952 1.33461529053606 1.33127550823193 1.32793571890697 1.32459595614830 1.32125625333121 1.31791664362027 1.31457715997020 1.31123783512696 1.30789870162865 1.30455979180654 1.30122113778599 1.29788277148748 1.29454472462753 1.29120702871964 1.28786971507532 1.28453281480497 1.28119635881886 1.27786037782809 1.27452490234548 1.27118996268656 1.26785558897046 1.26452181112088 1.26118865886698 1.25785616174431 1.25452434909574 1.25119325007237 1.24786289363441 1.24453330855215 1.24120452340678 1.23787656659136 1.23454946631169 1.23122325058718 1.22789794725176 1.22457358395477 1.22125018816181 1.21792778715567 1.21460640803712 1.21128607772588 1.20796682296139 1.20464867030376 1.20133164613454 1.19801577665766 1.19470108790021 1.19138760571334 1.18807535577308 1.18476436358118 1.18145465446596 1.17814625358312 1.17483918591663 1.17153347627945 1.16822914931448 1.16492622949528 1.16162474112694 1.15832470834688 1.15502615512566 1.15172910526776 1.14843358241246 1.14513961003454 1.14184721144515 1.13855640979258 1.13526722806305 1.13197968908148 1.12869381551230 1.12540962986023 1.12212715447104 1.11884641153231 1.11556742307426 1.11229021097045 1.10901479693858 1.10574120254124 1.10246944918668 1.09919955812955 1.09593155047166 1.09266544716271 1.08940126900108 1.08613903663452 1.08287877056091 1.07962049112901 1.07636421853917 1.07310997284407 1.06985777394945 1.06660764161482 1.06335959545421 1.06011365493683 1.05686983938786 1.05362816798911 1.05038865977972 1.04715133365691 1.04391620837666 1.04068330255439 1.03745263466568 1.03422422304698 1.03099808589625 1.02777424127369 1.02455270710241 1.02133350116911 1.01811664112476 1.01490214448529 1.01169002863224 1.00848031081345 1.00527300814371 1.00206813760545 0.998865716049361 0.995665760195108 0.992468286631939 0.989273311819364 0.986080852087796 0.982890923639205 0.979703542547760 0.976518724760467 0.973336486097822 0.970156842254431 0.966979808799657 0.963805401178249 0.960633634710965 0.957464524595207 0.954298085905637 0.951134333594803 0.947973282493756 0.944814947312659 0.941659342641416 0.938506482950264 0.935356382590394 0.932209055794548 0.929064516677626 0.925922779237284 0.922783857354535 0.919647764794513 0.916514553007848 0.913626501586297 0.957181910437377 1.18507611517074];
Ri=z_diff./cos(beta);
screens = zeros(Nx,Nx,length(F)-1);
% pre-allocate 3D array I
I = zeros(Nx,Nx,numel(F));
for i = 2:numel(F)
x = x + z_diff(i)*tan(beta);
screens(:,:,i)= screens(:, :, i-1) + exp(1i*((Kx.*x)+(Ky.*y))).*delfx.*delfy.*sqrt(F(i)).*R;
Hi = exp(1i*Kz*Ri(i) - c(i).*Ri(i));
Ui = fftshift(fft2(Uout));
Ui = Hi.*Ui;
Ui = ifft2(ifftshift(Ui));
Uio = Ui.*exp(1i.*screens(:,:,i));
% calculate page i of I
I(:,:,i) = Uio.*conj(Uio);
Uout = Uio;
end
whos I
Name Size Bytes Class Attributes I 512x512x757 1587544064 double
% to calculate SigmaI
A1 = reshape(sum(I.^2,[1 2]),1,[]); % sum over 1st two dimensions, and reshape to row vector
A2 = reshape(sum(I,[1 2]),1,[]); % same
SI = abs(((A1/Nx)./((A2/Nx).^2)) - 1) % use array (i.e., element-wise) operations
SI = 1x757
NaN 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503 0.7503
<mw-icon class=""></mw-icon>
<mw-icon class=""></mw-icon>
%to calculate BER
syms u
PI = (1./(u*sqrt(2*pi*SI))).*exp(-((log(u) + SI ).^2)./(2*SI)) ; % use array (i.e., element-wise) operations throughout
SNR0 = 1;
SNR = SNR0./sqrt(1+SI*SNR0) ;
fun = PI.*erfc(real((SNR.*u)/(2*sqrt(2))));
ff = matlabFunction(fun,'Vars',u);
BER = 0.5*integral(ff,0,Inf,'ArrayValued',true) % use ArrayValued
Warning: Inf or NaN value encountered.
BER = 1x757
NaN 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032 0.4032
<mw-icon class=""></mw-icon>
<mw-icon class=""></mw-icon>
  3 Comments
Show 1 older comment
Voss
Voss on 9 May 2024
I really don't know. Maybe ask a new question about it since the question of how to extend code to 3D and the question of why the results don't match the report are really very different questions. Be sure to include a link or other reference to the report whose results you're trying to replicate.
Torsten
Torsten on 9 May 2024
I'm confused that in the formula for SI, I seems to be a dependent variable and depend on x and y while in the expressions for BER and p_I, it turns into an independent variable of integration.

Sign in to comment.

More Answers (0)

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!