Research question : how to compute chirality by DIV & CURL function

2 次查看（过去 30 天）
chi shing
chi shing 2014-1-21
评论： Roger Stafford 2014-1-23
----------------------------------
[X Y Z] = meshgrid(y,x,z);
[cx1,cy1,cz1] = curl(X,Y,Z,Ex,Ey,Ez);
div1 = divergence(X,Y,Z,cx1,cy1,cz1);
term1 = 8.854*10^-12*0.5*div1;
[cx2,cy2,cz2] = curl(X,Y,Z,Hx,Hy,Hz);
div2 = divergence(X,Y,Z,cx2,cy2,cz2);
term2 = (1/(4*pi*10^-7))*0.5*div2;
chirality = term1 + term2;
----------------------------------
The above are my script to calculate chirality However, all Ex,Ey,Ez,Hx,Hy,Hz are complex number, but I already to transfer to real number. Should I use complex or real? Any error in syntax?

请先登录，再回答此问题。

回答（1 个）

Roger Stafford
Roger Stafford 2014-1-21
As far as I can see, the equation you show does not contain any divergence operations. One term is the scalar (dot) product of the electric field, E, by its curl, and the other term the same with the magnetic field, B. I see no divergence here. The divergence of the curl would imply the presence of second partial derivatives and there are none in the equation. In any case, the divergence of the curl of any vector field would always be identically zero.
  2 个评论
chi shing
chi shing 2014-1-22
Sorry about the unclear question. In fact, I wanna to show that DIV of CURL of vector should be zero, so that it can ensure my syntax is correct. One more question is that if all Ex,Ey,Ez,Hx,Hy,Hz are complex number, the above result will still get zero?
Roger Stafford
Roger Stafford 2014-1-23
Yes, the divergence of the curl of a vector field with continuous second derivatives with respect to position coordinates will always be zero even if this vector field is complex-valued.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Code Performance 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by