/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 16.04.1 ] */

/* [wxMaxima: title   start ]
Phase Controlled Rectifier, Fully Controlled
   [wxMaxima: title   end   ] */

/* [wxMaxima: comment start ]
RMS values
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
Vrms: Vm / sqrt(2);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
current, 1st harmonic
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
assume(Iout > 0);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Irms: Iout;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
I1s: 1 / %pi * (integrate(Iout * sin(x), x, alpha, %pi + alpha) + 
    integrate(- Iout * sin(x), x, %pi + alpha, 2 * %pi + alpha));
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
I1c: 1 / %pi * (integrate(Iout * cos(x), x, alpha, %pi + alpha) + 
    integrate(- Iout * cos(x), x, %pi + alpha, 2 * %pi + alpha));
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
I1: sqrt(I1c^2 + I1s^2);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
I1: trigsimp(I1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
I1rms: I1 / sqrt(2);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
DPF
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
DPF: I1s / I1;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
P
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
P: Vrms * I1rms * DPF;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
S
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
S: Vrms * Irms;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([P / S, S / S], [alpha, 0, %pi-1e-3], [box, false], [legend, false]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
PF
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
PF: P / S;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ev(PF, alpha = 0);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ev(PF, alpha = 0, numer);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d(PF, [alpha, 0, %pi-1e-3], [box, false], [legend, false]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
verify PF-DPF relation
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([PF, I1rms / Irms * DPF], [alpha, 0, %pi-1e-3], [box, false], [legend, false]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
THD
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
THD: sqrt((Irms / I1rms)^2 - 1) * 100;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ev(THD, numer);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d(THD, [alpha, 0, %pi-1e-3], [y, 0, 60], [box, false], [legend, false]);
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$