正方形環形天線的長岡係數Ｋｓ與電感Ｌ

μ0 = 4*π*10^(-7);
k0 = p/w;
k1 = w/p;
k2 = k1^2;
k3 = k1³;
ksr1 = √(k2+1);
ksr2 = √(k2+2);
Ks = (2/π)*((k1-1/k1)*ln(ksr1/(ksr2+1)) + k1*ln((ksr1+1)/k1) + 2*atan(k1/ksr1) +
(1/k1)*(ksr1-ksr2-ln(sqrt(2)+1)+(1/3)*(k3+(k2+2)^(3/2)-2*(k2+1)^(3/2)+sqrt(2)-1)));
L = Ks*μ0*p^2*N^2/w;
LmicroH = L*10^6;
S = p^2;
print(k0);
print(Ks);
print(LmicroH);
print(S);
ω = 2*π*frequency;
XL = ω*L;
C = 1/(ω*XL);
CpF = C *10^12;
print(XL);
print(CpF);
LengthWire = 4*p*N;
print(LengthWire);
λ=(3*10^8)/frequency;
Rr = (8/3)*π^3*120*π*(N*S/λ^2)^2;
print(Rr);
ρ = 1.68*10⁻⁸;
δ = sqrt(2*ρ/(ω*μ₀));
radius = radius_mm/1000;
Rdc = ρ*LengthWire /(π*radius²);
Rac = ρ*導線長度 / (π*(半徑² - (半徑 - δ)²));
直徑_mm = 2*半徑_mm;
δ微米 = δ*10⁶;
print(Diameter_mm );
print(Rdc);
print(δmicron);
print(Rac);
η = Rr/(Rr+Rac)*100;
print(η);
heff = 2*pi*p^2/λ;
print(heff);
Q = XL/Rac;
print(Q);
print(ρ);
Vo = Q*heff*Eradio;
print(Vo);
Rp = (Rr+Rac)*Q^2;
print(Rp);