/* [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"$