Version pdf de ce texte : puissance_elec_scilab.pdf

lycee - Racine du site

Puissance électrique en sinusoïdal avec Scilab

On peut reconstituer l’allure de la puissance instantanée reçue par un dipôle en courant sinusoïdal avec quelques lignes de code.

Tout d'abord on génère un vecteur de points « temps ». Puisque nous utilisons fréquemment le réseau 50 Hz, prenons une durée de 2 périodes.
Si f = 50 hz alors T = 1/f donc T = 1/50 = 0,020 s = 20 ms Deux périodes font 0,040 s. Pour bien tracer les courbes on génère un vecteur de 1 000 valeurs également réparties entre 0 et 0,040.

--> t=linspace(0,0.040,1000);

Puis on donne les valeurs de la fréquence, du déphasage (ici le courant sera en retard : dipôle récepteur inductif)
et on en déduit la pulsation.

--> f=50;

--> phi=-%pi/6;

--> omega=2*%pi*f;

On fait calculer les matrices lignes de 1 000 valeurs représentant u(t) et i(t).

--> u=cos(omega*t);

--> i=cos(omega*t+phi);

Pour calculer le produit de u et i il faut faire un calcul membre à membre des éléments des deux matrices.

--> p=u.*i;

On trace u(t) en rouge.

--> plot2d(t,u,5)

On trace i(t) en bleu.

--> plot2d(t,i,2)

On trace p(t) en noir.

--> plot2d(t,p)

Et voilà. Ce n'est pas plus compliqué.

N.B. : On peut vérifier sur le graphique que si les valeurs moyennes de u(t) et i(t) sont bien nulles, ce n'est pas le cas de p(t).