f = 1e6:1e5:100e6; r=3.9904e-3; D=15.8e-3; delta = sqrt(1./pi./f./mu_c./sigma_c); R_solid = 1./pi./r./delta./sigma_c; R = (D./2./r)./sqrt((D./2./r).^2-1).*R_solid; Ls = R./2./pi./f; Lm = mu_c/pi*acosh(D/2/r); L = Ls+Lm; C = pi*epslon/acosh(D/2/r); G = 2.*pi.*f.*C.*tdelta; temp_a = complex(R, 2.*pi.*f.*L); temp_b = complex(G, 2.*pi.*f.*C); gama = sqrt(temp_a.*temp_b);将这代码全部转换为数学公式表达

$f = 1\times10^6:1\times10^5:100\times10^6$

$r=3.9904\times10^{-3}$

$D=15.8\times10^{-3}$

$\delta = \sqrt{\frac{1}{\pi f \mu_c \sigma_c}}$

$R_{solid} = \frac{1}{\pi r \delta \sigma_c}$

$R = \frac{D}{2r} \cdot \frac{R_{solid}}{\sqrt{\left(\frac{D}{2r}\right)^2-1}}$

$L_s = \frac{R}{2\pi f}$

$L_m = \frac{\mu_c}{\pi} \cdot \text{acosh}\left(\frac{D}{2r}\right)$

$L = L_s + L_m$

$C = \frac{\pi \epsilon}{\text{acosh}\left(\frac{D}{2r}\right)}$

$G = 2\pi f C \delta$

$\text{temp}_a = R + 2\pi f L i$

$\text{temp}_b = G + 2\pi f C i$

$\gamma = \sqrt{\text{temp}_a \cdot \text{temp}_b}$