
    fces=@(x) sin((pi*(x.^2))/2);
    fcec=@(x) cos((pi*(x.^2))/2);
    
    v1=-8:0.01:8;
    hodnoty=-8:0.01:8;
    hodnoty2=hodnoty;
    r=1;
    for x=v1
        
        E=((1+1i)/2)*(((1/2)-(integral(fcec,0,x)))-1i*((1/2)-(integral(fces,0,x))));  %(1/2)*((1/2)+(integral(fcec,0,x))^2-integral(fcec,0,x)+(integral(fces,0,x))^2-integral(fces,0,x));
        hodnoty(r)=20*log10(abs(E));
        hodnoty2(r)=20*log10(1/(pi*x*sqrt(2)));
        r=r+1;
        
    end
    x=1;
    E1=10*((1+1i)/2)*(((1/(2))-integral(fcec,0,x))-1i*((1/(2))-integral(fces,0,x)));
    abs(E1);
    20*log10(abs(E1));
    plot(v1,hodnoty)     
    grid on
    hold on
    %plot(v1,hodnoty2)
    xlabel('v [-]');
    ylabel('L [dB]');
    %%
f=1000;
h_l=1;
h_s=100;
h_hm=50;
vzd=880;
vzd_m=400;
E0=10;
R_v=-0.5;
if vzd<=vzd_m
    E=0;
else
vln=((3*10^(8))/(f*10^(6)));
k=((2*pi())/(vln));

[vzd_odr,d1,d2,d4,d3]=vzdalenosti_pro_odrazadifrakci(h_l,h_s,h_hm,vzd,vzd_m);
alfa=acos((((double(d2))^2)+((double(d4))^2)-((double(d3))^2))/(2*double(d2)*double(d4)));
h=(double(d2))*sin(alfa);
v2=h*sqrt((2/vln)*((1/double(d2))+(1/double(d3))));


%v=vypocek_koef_v(vzd_odr+vzd_m,0,h_s,f,h_hm);
Ep=pouze_primy(d3,E0,f); %intenzita v bode odrazu od primeho paprsku
Ed=intenzita_od_difrakce(vzd_odr+vzd_m,vzd_m,0,h_s,f,h_hm,v2);
Eo=Ep*Ed;  %intenzita difrakci v bode odrazu od vody
E=((Eo/double(d1))*exp(-1i*k*double(d1)))*(R_v*exp(-1i*pi()));
end
E
    
%%   
    
 vzdalenost = 5000; %[m]
    krok = 0.8;
    vyska_lodi = 1; %[m]
    vyska_satelitu = 100; %[m]
    vyska_mostu_1 = 50 %[m] vyska horni hrany mostu
    vyska_mostu_2 = 40 %[m] vyska spodni hrany mostu
    vzdalenost_mostu_1 = 500; %[m] vzdálenost mostu od kolmého průmětu satelitu na hladinu řeky
    vzdalenost_mostu_2 = 450;
    koef_odrazu_vody = -0.5;
    koef_odrazu_koryta = -0.8;
    frekvence = 1000; %[MHz]
    uhel_odchylky_koryta = (pi())/12; %[rad]
    vzdalenost_od_brehu = 10; %[m]
    E0 = 10; %[V/m]
   
v=-8:0.1:8;
hodnoty=-8:0.1:8;
r=1;
for x=v
E=intenzita_od_difrakce(vzdalenost,vzdalenost_mostu_2,vyska_lodi,vyska_satelitu,frekvence,vyska_mostu_2,x);
hodnoty(r)=abs(E);
r=r+1;
end

plot(v,hodnoty)
grid on



    
    
    
    
    