roldao
Newbie level 4
Help FDTD+CPML urgent
Hello,
I'm working in FDTD+CPML and I need urgent help.
I´m simulating the situation of the paper that I send in attach.
In the case of PEC buildings, I can´t get the results of power delay profiles shown in the paper. I use a different mode than they, I use TEz mode, with the same characteristics dt and dx that they use.
When I calculate the power using the expression:
10log(((Ex^2+Ey^2)*lambda^2*1)/(120*pi*4*pi))
and using the excitation like:
Jmax =1000;
alpha = (4/(32*dt))^2;
source=zeros(1,nSteps);
for n=1:88%(9 e 30)
source= Jmax*(exp(-alpha*(((dt*n)-(32*dt))^2)));
end
then I apply in the time loop the source like:
Hz(x,y)=source(x,y);
I´m using everything like they have in this document and I can´t reach the same results.
The plot aspect is similar but the power is bigger (50 dBw bigger like you can see in the file attached)
Please Help me.
Thanks a lot
Joao Roldao
Hello,
I'm working in FDTD+CPML and I need urgent help.
I´m simulating the situation of the paper that I send in attach.
In the case of PEC buildings, I can´t get the results of power delay profiles shown in the paper. I use a different mode than they, I use TEz mode, with the same characteristics dt and dx that they use.
When I calculate the power using the expression:
10log(((Ex^2+Ey^2)*lambda^2*1)/(120*pi*4*pi))
and using the excitation like:
Jmax =1000;
alpha = (4/(32*dt))^2;
source=zeros(1,nSteps);
for n=1:88%(9 e 30)
source= Jmax*(exp(-alpha*(((dt*n)-(32*dt))^2)));
end
then I apply in the time loop the source like:
Hz(x,y)=source(x,y);
I´m using everything like they have in this document and I can´t reach the same results.
The plot aspect is similar but the power is bigger (50 dBw bigger like you can see in the file attached)
Please Help me.
Thanks a lot
Joao Roldao