How do I calculate the electromagnetic field radiated by lightning using MATLAB?

I want to calculate the electromagnetic field radiated by lightning using MATLAB.

• I calculated the channel base current using “Heidler function”

ib=(((i01/N1)*((t/t11).^2)./(1+(t/t11).^2)).*exp(-(t/t21)));
with N1=exp(-(t11/t21)*((n*t21/t11)^(1/n)));

• I calculated the return stroke current using MTLE model
• The decay constant which allows the current to reduce its amplitude with height is:

 i(z,t)=ib(0,t-z/v)*exp(-z/?)

Now I want to calculate the electromagnetic field radiated. The lightning channel is considered as a vertical antenna of height H divided into n dipoles of length dz,the expressions for electric and magnetic fields are the following

 dH (r, z, t) = (dz / 4 * pi) * [(r / R ^ 3) * i (z,t-R / c) + (r / c * R ^ 2) * di (z, t - R / c) / dt]
 dEr(r,z,t)=(dz / 4 * pi*? o) * [(3r(z'-z)/R^3)*?i(z,t-R/c)dt.......
 dEz(r,z,t)=...........

How can I implement this in MATLAB?

Answer 1

[Paulo Silva]

Take a look at Electromagnetic Field Demo

Answer 2

[kezon]

i have a same problem how to implemen a time depanding field equation in matlab and then integrate it...

i calculate the heidler function

%%heidlerova funkcija%%
clc
clear
%first stroke%
i01=10.7;  
t11=0.25;
t21=2.5;
n=2;
N1=exp(-(t11/t21)*((n*t21/t11)^(1/n)));
t=0:0.01:5;
ib=(((i01/N1)*((t/t11).^n)./(1+(t/t11).^n)).*(exp(-(t/t21))));
figure,plot(t,ib)

help nead