Demo entry 6771634

aaa

   

Submitted by anonymous on Nov 19, 2018 at 21:58
Language: Matlab. Code size: 685 Bytes.

Li = [L(:,1), 2.*room_y-L(:,2), L(:,3)]; % Calculate the reflection point coordinate
d_reflect = sqrt((Li(:,1)-source_x).^2+(Li(:,2)-source_y).^2+(Li(:,3)-source_z).^2); %Total distance of the reflection path
t_reflect_1 = d_reflect./1000./c; %time delay
p_reflect = mod(t_reflect_1*2*pi/T, 2*pi); %phase
thet = abs(atan((Li(:,2)-source_y)./(Li(:,1)-source_x)));%Incident angle
reflect_coefficient = (sin(thet)-sqrt(epsilon_w-(cos(thet)).^2))./(sin(thet)+sqrt(epsilon_w-(cos(thet)).^2));% Cauculation of reflect coefficient
E_reflect = (lambda./(4.*pi.*d_reflect./1000)) .*  reflect_coefficient;
E1=E_reflect.* exp(1i.*(-p_reflect)); %Combine the delay phase to reflection phase

This snippet took 0.00 seconds to highlight.

Back to the Entry List or Home.

Delete this entry (admin only).