calculates self inductance and mutual inductance for wire loops in 3d. There are two functions:
L=self_inductance(x,y,r) calculates self indactance for plane wire loop, that is polygon. x,y is vertices of the polygon. r- wire radius. The L is calculated by numerical integration of magnetic field flux: F=L*I, then L=F/I. Considered high frequency case when no field inside wire. Note: small wire radius can not be neglected in self inductance.
calculates mutual inductance of two loops that is polygons in 3d. x1 y1 z1 -vertices of first loop. x2 y2 z2 -vertices of second loop. Note: wire radius considered much smaller then typical loops size, then no inputs for wire radius in the function. It calculated numerically by integration Newman formula:
test_self_inductance.m is test for self_inductance function it compare theoretical and calculated values for self inductance for ring.
calculated value: 7.2355e-006 H
real value: 7.9093e-006 H
i.e 10% error
guess for mutual inductance is more precise (not checked carefully).
Any way it is possible to make more small elements in meshes so increase precision (and computational time).
squares_in_polygon - auxiliary function.