%Corrigé de l'exercice "Tension induite de mouvement"
%Copyright C.K.

%paramètres
N1 = 1000
N2 = 1000
%On peut rendre la résistance dépendante du nombre de spires N1. Pour ce faire, on suppose que le volume de la bobine ne change pas
% et qu'on réduit le diamètre du fil lorsque N1 augmente => R1 = rho_Cu * lfil / Sfil avec Sfil = S_totale_bobine / N1 et l_fil = N1 * l_moyenne_1tour
% la résistance du fil est donc proportionelle à N1^2
% elle vaut 12 Ohm pour N1 = 1000
R1 = 12/(1000^2)*N1^2
%R2 n'apparait pas dans la réponse vu qu'il n'y a pas de courant dans la bobine 2
U1 = 12
mu0 = 4e-7*pi
Omega = 2*pi*50
Sf = 400e-6
lf=150e-3
murf = 500
delta_a = 0.2e-3
delta_b = 5e-3
%résultats selon corrigé
delta_c = lf/murf+delta_a+delta_b
f = @(t,i1)(U1-R1*i1)*(delta_c-delta_b*cos(Omega*t))/(N1^2*mu0*Sf)+ (i1*delta_b*Omega*sin(Omega*t) )/(delta_c-delta_b*cos(Omega*t))
%résolution de l'équation différentielle
[t,i1] = ode45(f,[0:0.001:0.08],0)

%Plots du courant dans la bobine 1 et de la tension induite dans la bobine 2
figure(1)
plot(t,i1)
xlabel('Temps [s]')
ylabel('Courant i1 [A]')
u2 = N2/N1*(U1-R1*i1)
grid()
figure(2)
plot(t,u2)
xlabel('Temps [s]')
ylabel('Tension U2 [V]')
grid()