% calculo de R_ac de inductor debido a efecto de proximidad y efecto skin.
%datos del ejemplo de bobina usado en el teórico (Nl=32, a=1.2 pulgadas, seccion del alambre (2*r) = 1mm).

roCu = 1.69e-8;
sigmaCu=1/roCu;
u0=4*pi*1e-7;
anucleo = 1.2*25.4e-3 % diametro bobinado en m
r=0.5e-3 % radio del alambre en m
Nl=32; % nro de vueltas (per layer, 1 sola layer)
b=Nl*2*r; % ancho del bobinado 
L=0.2e-6*(anucleo/25.4e-3)^2*Nl^2/(3*(anucleo/25.4e-3)+9*(b/25.4e-3)) % L(Hy)
f=2e6;	% frecuencia
l=Nl*pi*anucleo % m, largo de la bobina

Rload = 5e3 % R carga oscilador colpitts
Qtanque = Rload/(2*pi*f*L)


Stotal = pi*r^2; % seccion total del conductor
RDC = roCu*l/Stotal % resistencia en DC

% resistencia considerando solo efecto skin (sin efecto de proximidad)

delta=sqrt(1/(u0*sigmaCu*f*pi)) % profundidad efecto skin

% seccion efectiva debido a efecto skin
if(delta < r) 
   Sskin= pi*(r^2 - (r-delta)^2);
else 
   Sskin= pi*(r^2);
end
  
% R considerando solo efecto skin
Rskinhf = roCu*l/Sskin


% R por efecto proximidad, método de Dowell
% http://en.wikipedia.org/wiki/Proximity_effect_(electromagnetism)
% Dowell, P.L. (August 1966). "Effects of Eddy Currents in Transformer Windings". Proceedings IEE 113 (8): 1387–1394.

a=sqrt(Stotal); % altura del conductor si fuera cuadrado tal que de igual area que el que tenemos 
h=a;% ancho del conductor si fuera cuadrado tal que de igual area que el que tenemos (= altura)

eta = Nl*(a/b);
u0 = 4e-7*pi;
alfa = sqrt(j*2*pi*f*u0*eta/roCu);
M=alfa*h*coth(alfa*h);
Rprox = RDC*real(M)





QL=2*pi*f*L/(Rprox+Rskinhf)

Rparalelo= QL^2*(Rprox+Rskinhf)